服务商与监管部门决策演化复制动态方程matlab仿真——matlab2016a版本

1.输入常微分方程组

function dxdt=differential3(t,x)
dxdt=[x(1)*(1-x(1))*(-1.6-2.8*x(2));x(2)*(1-x(2))*(-1.4-2*x(1))];
end

保存为differential3.m文件

2.主函数

clear
%y-x
for i=0.1:0.2:0.9
    for j=0.1:0.2:0.9
        [T,Y]=ode45('differential3',[0 20],[i,j]);
        figure(1)
        grid on
        plot(T,Y(:,1),'r*-','lineWidth',1);    %T：时序，Y(:,1)：1指的是x。总表示dxdt
        hold on   
        plot(T,Y(:,2),'b--','lineWidth',1);    %T：时序，Y(:,2)：2指的是y。总表示dydt
        hold on
    end
end
axis([0 20 -0.1 1.1]);    %为了标出坐标轴最大值最小值，横坐标[0 20] ,纵坐标[-0.1 1.1]
set(gca,'XTick',[0:2:20],'YTick',[-0.1:0.1:1.1]);    %gca（ca即current axes）返回当前axes对象的句柄值;把当前的图形x轴标刻度，区间[0 20],间隔为2；Y同理
xlabel('evolutionary steps');    %标x轴名
ylabel('polution');     %标y轴名
legend('服务商','监管部门')    %标曲线名
title('服务商与监管部门决策演化路径')    %加标题
axis([0,20,0,1])

保存为for3.m文件

3.run结果

run出来的图形跟原视频作者的图不一致，各位大神有什么发现吗，如何修改能得到原视频结果图