自己写的,如有不对敬请指正
以下的代码是用于某一参数改变时,三方主体演化博弈的仿真图,我不会使用for循环,所以是手动调整的参数的改变,但是比较容易理解和出图
简单来说,下列代码MATLAB运行的界面需要打开四个文件:一个主函数画图的文件+三个复制动态方程的文件。
以下代码出的是某一参数在高中低三种力度下三方主体的趋势图
代码如下
// 复制动态方程输入并保存为odefun1;odefun2;odefun3,因为主函数是要调用这三个文件所以保证文件名和主函数调用的名称一致
function dz=odefun1(t,z)
dz=zeros(3,1);
dz(1)=z(1)*(1-z(1))*(0.48*z(2)+0.39*z(3)-0.2);
dz(2)=z(2)*(1-z(2))*(-0.48*z(1)+3*z(3)-0.34);
dz(3)=z(3)*(1-z(3))*(0.4*z(1)-0.32);
end
function dy=odefun2(t,y)
dy=zeros(3,1);
dy(1)=y(1)*(1-y(1))*(0.48*y(2)+0.15