njuhhf 发表于 2009-10-21 15:00

matlab 中ODE45 问题

各位好，
       对于一个带阻尼的单振子系统，在外简谐力作用下的运动方程为   M*加速度+R*速度+K*位移=Fsin（wt）（F,R,K,M,W,t已知）

小弟用MATLAB 的ode45函数 求解方程组
A       dy(1)=速度；
B       dy2=（Fsin（wt）-R*速度-K*位移）/M

这样得出的dy（1）就是位移，dy（2）就是速度。但是发现这时候速度与位移之间不满足微分关系.这是为何？
另外如果我要求加速度怎么办？是对速度微分还是要把速度，位移值代入方程B呢？

yufeng 发表于 2009-10-22 08:22

y（1)是位移 dy（1）是速度

octopussheng 发表于 2009-10-22 09:51

求加速度的话，把位移、速度求解结果再代入原方程即可
页: [1]
查看完整版本: matlab 中ODE45 问题