利用matlab中的ode45解一个线性方程问题,在第四阶频率是出现了这种图,不明白为什么会这样?
程序如下。
function dxdt=guan1(t,x,flag,omega,f1)%k=20000
dxdt=[x(2);
-((1.32)*x(2)+(0.19+4766.59)*x(1)--0.80*f1*cos(omega*t))/1.49;
x(4);
-((1.28)*x(4)+(-3.04+174123.42)*x(3)--0.43*f1*cos(omega*t))/1.45;
x(6);
-((1.32)*x(6)+(-10.43+1368005.25)*x(5)--0.26*f1*cos(omega*t))/1.45;
x(8);
-((0.32)*x(8)+(-23.85+4993849.88)*x(7)--0.27*f1*cos(omega*t))/1.38];
end
clear
clc
options=odeset;
ts=300;
omega=1900.9;
f1=0.01;
Xx=[0,0,0,0,0,0,0,0];
tspan=0:0.1:ts;
x=1;
[t1,x1]=ode45('guan1',tspan,Xx,options,omega,f1);
u1=(x1(:,1)*(cos(1.8953*x)-0.7342*sin(1.8953*x)-1.0000*cosh(1.8953*x)+0.7342*sinh(1.8953*x))+x1(:,3)*(cos(4.6954*x)-1.0185*sin(4.6954*x)-1.0000*cosh(4.6954*x)+1.0185*sinh(4.6954*x))+x1(:,5)*(cos(7.8550*x)-0.9993*sin(7.8550*x)-1.0001*cosh(7.8550*x)+0.9993*sinh(7.8550*x))+x1(:,7)*(cos(10.9955*x)-1.0003*sin(10.9955*x)-1.0003*cosh(10.9955*x)+1.0003*sinh(10.9955*x)));
plot(t1,u1,'k')
hold on