??shadow?
2021-09-24 11:09
采纳率: 0%
浏览 175

利用matlab中的ode45解一个线性方程问题,在第四阶频率是出现了这种图,不明白为什么会这样?

利用matlab中的ode45解一个线性方程问题,在第四阶频率是出现了这种图,不明白为什么会这样?
程序如下。

img

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 

  • 收藏

1条回答 默认 最新

  • 技术专家团-Joel 2021-09-24 11:21

    同学你好,你这个图跟频率没有什么关系啊,这是一个幅值时程曲线。你求解了一个ode,然后根据ode的解得到的时程曲线而已

    打赏 评论

相关推荐 更多相似问题