怎样用MATLAB计算my''+cy'+ky+A*sgn(y')=B*sin(t)运动微分方程？

sdpxzcz · 发表于 2007-8-20 10:20

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

怎样用MATLAB计算运动微分方程？

运动微分方程为:my''+cy'+ky+A*sgn(y')=B*sin(t)
其中y是t的函数,m为质量,,c阻尼系数,,k结构刚度,A,B是常数,他们都是已知量.
y''表示加速度,y'表示速度,y表示位移,sgn()表示符号函数;
怎样用MATLAB计算这个运动微分方程呀?
请各位高手帮帮忙呀!!!

[ 本帖最后由 sdpxzcz 于 2007-8-20 10:42 编辑 ]

咕噜噜 · 发表于 2007-8-20 10:31

用ode45计算，具体你可以参考matlab帮助文件

sdpxzcz · 发表于 2007-8-20 10:36

关键是里面有关于y'(速度)的符号函数呀!
不知道怎么用,
咕噜噜你可以帮帮我吗?
我的QQ:119452259
电子邮箱:sdpxzcz@163.com

无水1324 · 发表于 2007-8-20 10:43

直接加一个判断语句。

你可以参考分段线性的例子，原来的都是判断位移，你把它改成速度就可以了

sdpxzcz · 发表于 2007-8-20 10:49

我是Matlab新手,可不可以把程序编给我看看呀!非常感谢!

咕噜噜 · 发表于 2007-8-20 11:18

速度也可以啊，在ode45里面
function doty=myff(t,y)
c=1;k=20;A=13;B=2;
doty=[y(2);-(c*y(2)+k*y(1)+A*sign(y(1))+B*sin(t))];
[t,y]=ode45('myff',[0,1],[0 0])
相应参数自己修改

sdpxzcz · 发表于 2007-10-7 22:48

谢谢各位呀!!!

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 怎样用MATLAB计算my''+cy'+ky+Asgn(y')=Bsin(t)运动微分方程？

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

评分

回复 #3 sdpxzcz 的帖子

评分

回复 #3 sdpxzcz 的帖子

评分

[编程技巧] 怎样用MATLAB计算my''+cy'+ky+A*sgn(y')=B*sin(t)运动微分方程？

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

评分

回复 #3 sdpxzcz 的帖子

评分

回复 #3 sdpxzcz 的帖子

评分

[编程技巧] 怎样用MATLAB计算my''+cy'+ky+Asgn(y')=Bsin(t)运动微分方程？