matlab求解振动微分方程组中出现的问题

huzhaolin · 发表于 2006-3-30 17:49

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

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

x

振动方程组简化如下：x''+x'+y'-f(x,y)(x+y)=g; y''+x'+y'+f(x,y)(x-y)=k; g,k为已知量，if (x+y)>=0 f(x,y)=1 else f(x,y)=0; 如何用ode（45）进行求解！！！！！！！！

happy · 发表于 2006-3-30 17:50

把方程转化到状态空间，然后写成function，其中包含对f(x,y)的判断 然后直接调用ode45就行了

huzhaolin · 发表于 2006-3-30 18:34

function F=fun(t,f) if(f(1)+f(2)>=0); f(e)=1 else f(e)=0; F=[f(4); f(3); g-f(3)-f(4)+f(e)*(f(2)+f(1)); k-f(3)-f(4)-f(e)*(f(2)-f(1))] [t,y]=ode(@fun,0,10,0,0,0,0) 请问这样可以吗？

happy · 发表于 2006-3-30 18:37

<DIV class=quote>以下是引用huzhaolin在2006-3-30 18:34:21的发言： 
function F=fun(t,f) if(f(1)+f(2)>=0); f(e)=1 else f(e)=0; F=[f(4); f(3); g-f(3)-f(4)+f(e)*(f(2)+f(1)); k-f(3)-f(4)-f(e)*(f(2)-f(1))] [t,y]=ode(@fun,0,10,0,0,0,0) 请问这样可以吗？ </DIV>
 你这里的f是什么东西？

huzhaolin · 发表于 2006-3-30 18:43

令f(1)=y,f(2)=x,f(3)=x',f(4)=y'; function F=fun(t,f) if(f(1)+f(2)>=0); fe=1 else fe=0; F=[f(4); f(3); g-f(3)-f(4)+fe*(f(2)+f(1)); k-f(3)-f(4)-fe*(f(2)-f(1))] [t,y]=ode(@fun,0,10,0,0,0,0)

happy · 发表于 2006-3-30 18:49

<DIV class=quote>以下是引用huzhaolin在2006-3-30 18:43:24的发言： 
令f(1)=y,f(2)=x,f(3)=x',f(4)=y'; function F=fun(t,f) if(f(1)+f(2)>=0); fe=1 else fe=0; F=[f(4); f(3); g-f(3)-f(4)+fe*(f(2)+f(1)); k-f(3)-f(4)-fe*(f(2)-f(1))] [t,y]=ode(@fun,0,10,0,0,0,0)</DIV>
？？？？？？

huzhaolin · 发表于 2006-3-30 18:55

我这是大致步骤，是个大致草稿，请问是不是这样进行求解的。谢谢教授！

happy · 发表于 2006-3-30 18:57

是的

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

matlab求解振动微分方程组中出现的问题

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

回复：（huzhaolin）matlab求解振动微分方程组中出现的...

回复：（huzhaolin）function F=fun(t,f)if(f(1)+f(2)...

回复：（huzhaolin）令f(1)=y,f(2)=x,f(3)=x',f(4)=y'...

回复：（huzhaolin）matlab求解振动微分方程组中出现的...