matlab软件怎么做演化博弈模型的仿真分析

应该是微分方程的不同initial condition下的解，如果是微分方程的数值解法google上培猛桐matlab的code到处都配坦是，而且matlab自带知锋多个built-in微分方程解法，ode113,ode 45等等

1、在MATLAB plots中，没有找到与老型胡复制动态相位图相关的绘图工具。

2、复制动态相位图是根据复制动态方程和雅克比矩阵的计算结果画出来的，若再用Matlab做出相应的复制相位图，也仅侍拦仅是用另一种方式将之前的结果重新表述一遍，意义不大。

3、我觉得应该从MATLAB的常微分与偏微分方程入手，找到表达方程的plot，能够表达复制动态相位图的含义却又不与之前的复制动态相位图完全相同租凳。

这个代码有瞎凳蠢两个错误：

在定义function方程代码中：

楼主原：dxdt=[x(1)*(1-x(1))*(2*x(3)+1x(2)*(1-x(2))*(3*x(3)+1x(3)*(1-x(3))*(-1*x(1)-x(2)-x(1)*x(2)+4)]

改正后：dxdt=[x(1)*(1-x(1))*(2*x(3)+1)x(2)*(1-x(2))*(3*x(3)+1)x(3)*(1-x(3))*(-1*x(1)-x(2)-x(1)*x(2)+4)]

错误原因：在前两个分号‘；’前，缺少右括号‘）’

function的x需要传递三个参数想x(1) x(2) x(3)，而楼主在主磨陪程序的调用ode45函数时，只传递了两个参数i和j，可引入变量k，通过x-y-z为例，粗裤改成如下（修改部分已经加粗），其它部分修改相同：

%z-y-x

for i=0:0.1:1

for j=0:0.1:1

for k=0:0.1:1

[T,Y]=ode45('differential',[0 5],[i j k])

figure (1)

grid on

plot3(Y(:,1),Y(:,2),Y(:,3))

hold on

end

end

end

---