博弈仿真算法

function dxdt=differential(t,x)

dxdt=[x(1)*(1-x(1))*(5-6*x(2));x(2)*(1-x(2))*(1-6*x(1))]; end

clear

%y-x

for i=0:0.1:1

for j=0:0.1:1

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

figure(1)

grid on

plot(Y(:,1),Y(:,2));

hold on

end

end

%dx/dt-t

for i=0:0.1:1

for j=0:0.1:1

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

figure(2)

grid on

plot(T,Y(:,1));

hold on

end

end

%dy/dt-t

for i=0:0.1:1

for j=0:0.1:1

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

figure(3)

grid on

plot(T,Y(:,2));

hold on

end

end

[t,x]=ode45(@vdp,[0,5],[0.2,0.1],1)

figure;

plot(t,x,'r')

hold on;

[t,x]=ode45(@vdp,[0,5],[0.5,0.1],1)

plot(t,x,'b')

[t,x]=ode45(@vdp,[0,5],[0.8,0.1],1)

plot(t,x,'k') hold off