xin_shou 发表于 2007-4-22 21:07

为什么增大激励频率相图会成这样？

function aa()

= ode23s(@vdp,,');
plot(Y(:,1),Y(:,2));
function dy=vdp(t,y)
R0=25e-3;
R1=33.0e-3;
Ig1=2.17e-3;
I0=400;
I1=2;
k1=4.0711e+004;
kg1=4.7599e+005;
   
    c1=1.1445;
    cg1=5.573;
dy=zeros(4,1);
dy(1)=y(2);
dy(2)=sin(0.5*t)-((1/I0+1/I1)*k1*y(1)+(1/I0+1/I1)*c1*y(2)...
    -(kg1*R0)/(R1*I1)*y(3)-(cg1*R0)/(R1*I1)*y(4));
dy(3)=y(4);
dy(4)=(k1*R1)/(R0*I1)*y(1)+(c1*R1)/(R0*I1)*y(2)...
    -((1/I1+1/Ig1)*kg1)*y(3)-((1/I1+1/Ig1)*cg1)*y(4);
运行后图：

xin_shou 发表于 2007-4-22 21:10

增大为dy(2)=sin(10*t).......
运行后：

gghhjj 发表于 2007-4-23 07:55

很显然是你画图的问题，你画出来的结果中的比较乱的线都是没有收敛前的结果

画图的时候要略去没有收敛前的结果

gghhjj 发表于 2007-4-23 07:57

第一种情况

= ode23s(@vdp,,');
plot(Y(500:end,1),Y(500:end,2));

gghhjj 发表于 2007-4-23 07:59

第二种情况

= ode23s(@vdp,,');
plot(Y(5000:end,1),Y(5000:end,2));

xin_shou 发表于 2007-4-23 09:58

非常谢谢！:handshake
页: [1]
查看完整版本: 为什么增大激励频率相图会成这样？