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

xin_shou

function aa()

[t,Y] = ode23s(@vdp,[0 100],[0 0 0 0]');
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

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

gghhjj

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

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

gghhjj

第一种情况

1. [t,Y] = ode23s(@vdp,[0 100],[0 0 0 0]');
   
2. plot(Y(500:end,1),Y(500:end,2));

复制代码

gghhjj

第二种情况

1. [t,Y] = ode23s(@vdp,[0 100],[0 0 0 0]');
   
2. plot(Y(5000:end,1),Y(5000:end,2));

复制代码

xin_shou

非常谢谢！:handshake