请教高手：这个非线性方程组该如何求解？

vicky_boy · 发表于 2006-4-18 10:47

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

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

x

方程如下： mx''+cx'+∑k(x*cos(θi)+y*sin(θi))1.5*cos((θi)=W+Fucos(wt) my''+cy'+∑k(x*cos(θi)+y*sin(θi))1.5*sin((θi)=Fucos(wt) 式中m，c，k，W，Fu，w都是已知，θi＝2*π*（i－1)/z＋Wc*t，z，Wc也是已知。 在matlab下编写的程序如下： xdot=[x(2);(Fx+Fu*cos(w*t)-c*x2-sum((ko*(x(1)*cos(ai)+x(3)*sin(az)).^1.5).cos(az)))/m;... x(4);(Fy+Fu*sin(w*t)-c*x2-sum((ko*(x(1)*cos(ai)+x(3)*sin(az)).^1.5).sin(az)))/m] 然后用ode45求解，这样编写是否正确？

xj2070 · 发表于 2006-4-19 21:59

it is suggested that the simultanous equations 'd better be storaged in one dependent file; and usually it is expressed in the form of the state equations; then empoly the ode45 or other functions to solve the equations

flybaly · 发表于 2006-4-20 19:12

<DIV class=quote>以下是引用vicky_boy在2006-4-18 10:47:34的发言： 
方程如下： mx''+cx'+∑k(x*cos(θi)+y*sin(θi))1.5*cos((θi)=W+Fucos(wt) my''+cy'+∑k(x*cos(θi)+y*sin(θi))1.5*sin((θi)=Fucos(wt) 式中m，c，k，W，Fu，w都是已知，θi＝2*π*（i－1)/z＋Wc*t，z，Wc也是已知。 在matlab下编写的程序如下： xdot=[x(2);(Fx+Fu*cos(w*t)-c*x(2)-sum((ko*(x(1)*cos(ai)+x(3)*sin(az)).^1.5).cos(az)))/m;... x(4);(Fy+Fu*sin(w*t)-c*X(4)-sum((ko*(x(1)*cos(ai)+x(3)*sin(az)).^1.5).sin(az)))/m] 然后用ode45求解，这样编写是否正确？ </DIV>
X=X(1);X'=X(2) Y=X(3);Y'=X(4) 然后编写xdot=[dx;dx';dy;dy']就可以了

[此贴子已经被作者于2006-4-20 19:13:12编辑过]

vicky_boy · 发表于 2006-4-23 15:07

这个错误我已经改过来了，但还是不行，错误还是一样的 是不是sum用的不对？

vicky_boy · 发表于 2006-4-23 15:09

错误为：")" expected, "." found，不得其解

flybaly · 发表于 2006-4-23 19:40

<DIV class=quote>以下是引用vicky_boy在2006-4-23 15:09:08的发言： 错误为：")" expected, "." found，不得其解</DIV> 呵呵,cos(az)和sin(az)前少了乘号了 认真检查检查吧

[此贴子已经被作者于2006-4-23 19:41:23编辑过]

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

[分形与混沌] 请教高手：这个非线性方程组该如何求解？

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

回复：（flybaly）以下是引用vicky_boy在2006-4-18 10...

回复：（flybaly）以下是引用vicky_boy在2006-4-18 10...

回复：（vicky_boy）回复：（flybaly）以下是引用vick...