Matlab - 演化博弈论实现

小白入门演化博弈论，演化博弈论的复制动态方程一些推导问题都已经搞明白。下面在学习如何使用matlab进行防止实验。
记录一下，第一次跟着B站视频敲的代码：视频来源
回头写一个如何推导复制动态方程的博客

动态微分方程函数

function dxdt = differential(t, x) 
    dxdt = [x(1)*(1-x(1))*(-2.3 - 2.8*x(2));x(2)*(1 - x(2))*(1.3 - 2 * x(1))];
end

绘图代码

clear
for i = 0.1 : 0.2 : 0.9
    for j = 0.1 : 0.2 : 0.9
        [T, Y] = ode45('differential', [0 20], [i j]);
        figure(1)
        grid on
        plot(T, Y(:, 1), 'r*--', 'lineWidth', 1);
        hold on;
        plot(T, Y(:, 2), 'b--', 'lineWidth', 1);
        hold on;
    end
end

axis([0 20 -0.1 1.1]);
set(gca, 'XTick', [0:2:20], 'YTick', [-0.1:0.1:1.1])
xlabel('evolutionary steps');
ylabel('populations');
legend('数据拥有者', '数据消费者');

结果图