??shadow? 2021-09-24 11:09 采纳率: 0%
浏览 181
已结题

利用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的解得到的时程曲线而已

    评论

报告相同问题?

问题事件

  • 系统已结题 10月2日
  • 赞助了问题酬金 9月24日
  • 创建了问题 9月24日

悬赏问题

  • ¥15 CNVcaller安装后无法找到文件
  • ¥15 visual studio2022中文乱码无法解决
  • ¥15 关于华为5g模块mh5000-31接线问题
  • ¥15 keil L6007U报错
  • ¥15 webapi 发布到iis后无法访问
  • ¥15 初学者如何快速上手学习stm32?
  • ¥15 如何自动更换布娃娃图片上的衣服
  • ¥15 心理学eprime编程
  • ¥15 arduino esp8266开发
  • ¥15 stm32单片机通过485发送命令给驱动器控制电机转动,同样的代码f103可以控制电机转动,换到f407不能动了,但是用串口助手调试407显示发送的命令都是正确的,卡了好久了这是发送规则