baipi 发表于 2009-1-16 16:04

请教:这个图如何绘制

传了几遍都没上来,图片太大了。
附件中是胡海岩《应用非线性动力学》中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);
    =ode45(@liuxing,,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 发表于 2009-1-17 14:28

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

tianyikillua 发表于 2009-1-30 17:06

帮你顶下。。
页: [1]
查看完整版本: 请教:这个图如何绘制