请教：这个图如何绘制

baipi

传了几遍都没上来，图片太大了。
附件中是胡海岩《应用非线性动力学》中54页的公式和相关的相轨迹图
我用下面的程序去模拟却怎么都得不到这样的图
请大家帮我看看是不是我的方法或者程序有问题
程序：
%建立函数liuxing
function dx=liuxing(t,x,mu,delta,beta,f)

dx=[-mu*x(1)+f./2*sin(x(2));delta+3/8*beta*x(1)^2+f./(2*x(1))*cos(x(2))];
%调用函数liuxing，用ODE45画相轨迹图
clear
tic
mu=0.1;
T=500;
f=0.1;delta=.3;A=rand(8,2);chuzhi=[];
for i=1:8
chuzhi(i,:)=A(i,:);
t0=0;                       %起始时间
t1=T*2*pi;
options=odeset('RelTol',1e-6,'AbsTol',1e-6);
    [t,x]=ode45(@liuxing,[0:0.01:t1],chuzhi(i,:),options,mu,delta,f);
figure(1);% subplot(2,1,1);
hold on
plot(x(20:end-10,1),x(20:end-10,2))  
xlabel('a');  ylabel('\phi');      %相图
end
toc

baipi

怎么没有人来啊,
只能自己顶了.

tianyikillua

帮你顶下。。