Question

matlab画图，程序没有报错，但是画出的图形是空白的

clc
clear all
%画系统的庞加莱图主程序；
T=2*pi;%周期；
tspan=0:T/100:200*T;%求解区间；
y0=[0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005];%初始值；
options=odeset('RelTol',10^-3,'AbsTol',10^-3);%选项优化；
[t,y]=ode45('fun',tspan,y0,options);
Y1=y(15000:end,1);
Y2=y(15000:end,2);
Y3=y(15000:end,3);
Y4=y(15000:end,4);
Y5=y(15000:end,5);
Y6=y(15000:end,6);
Y7=y(15000:end,7);
Y8=y(15000:end,8);
Y9=y(15000:end,9);
Y10=y(15000:end,10);
figure(1)
plot(y(15000:end,2),y(15000:end,3),'.');
%调用子程序；
function d=fun(t,y)
m1=220;
m2=20;
m3=440;
R1=0.092;
R2=0.46;
cb=200;
kb=7.85e9;
ks=1.43e8;
T1=8000;
T2=15000;
J1=0.5*m2*R1.^2;
J2=0.5*m3*R2.^2;
b=0.0125;
xig=0.1;
z1=23;
z2=120;
xn1=0.3622;
xn2=0.5129;
k0=0.04723+0.15551/z1+0.25791/z2-0.00635*xn1-0.11654*xn1/z1+0.00193*xn2...
+0.24188*xn2/z2+0.00529*xn1^2+0.00182*xn2^2;
epsilona=0.12;
km=(0.75*epsilona+0.25)*k0;
cm=2*xig*sqrt(km*m2*m3/(m2+m3));
e0=0;
ea=1.08e-3;
fai=0;
n1=2000;
w1=n1/30;
wh=z1*w1;
et=e0+ea*sin(wh*t+fai);
zz=y(3)+R1*y(5)-y(7)+R2*y(9)-et;
if (abs(zz)>=b)
z=zz-b;
else z=zz+b;
end
dz=y(4)+R1*y(6)-y(8)+R2*y(10)+wh*ea*cos(wh*t);
df=cm*dz+km*z;
d= [y(2);
-(cb/m1)*y(2)-( (kb+ks)/m1 )*y(1)+(ks/m1)*y(3);
y(4);
-df/m2+(ks/m2)*y(3)-(ks/m2)*y(1);
y(6);
-(df/J1)*R1-T1/J1;
y(8);
-df/m3;
y(10);
-(df/J2)*R2-T2/J2;];

Answer 1

输入y回车 看看是多少

追问

是NAN

NAN