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

vicky_boy

方程如下：<BR>mx<SUP>''</SUP>+cx<SUP>'</SUP>+∑k(x*cos(θ<SUB>i</SUB>)+y*sin(θ<SUB>i</SUB>))<SUP>1.5</SUP>*cos((θ<SUB>i</SUB>)=W+F<SUB>u</SUB>cos(wt)<BR>my<SUP>''</SUP>+cy<SUP>'</SUP>+∑k(x*cos(θ<SUB>i</SUB>)+y*sin(θ<SUB>i</SUB>))<SUP>1.5</SUP>*sin((θ<SUB>i</SUB>)=F<SUB>u</SUB>cos(wt)<BR>式中m，c，k，W，Fu，w都是已知，θ<SUB>i</SUB>＝2*π*（i－1)/z＋Wc*t，z，Wc也是已知。<BR>在matlab下编写的程序如下：<BR>ai=2*π*（i－1)/z＋Wc*t<BR>xdot=[x(2);(Fx+Fu*cos(w*t)-c*x2-sum((ko*(x(1)*cos(ai)+x(3)*sin(ai)).^1.5).cos(ai)))/m;...<BR>        x(4);(Fy+Fu*sin(w*t)-c*x2-sum((ko*(x(1)*cos(ai)+x(3)*sin(ai)).^1.5).sin(ai)))/m]<BR>然后用ode45求解，运行程序后，总是提示：Error：")" expected, "." found

happy

<DIV class=quote><B>以下是引用<I>vicky_boy</I>在2006-4-19 9:18:02的发言：</B><BR>方程如下：<BR>mx<SUP>''</SUP>+cx<SUP>'</SUP>+∑k(x*cos(θ<SUB>i</SUB>)+y*sin(θ<SUB>i</SUB>))<SUP>1.5</SUP>*cos((θ<SUB>i</SUB>)=W+F<SUB>u</SUB>cos(wt)<BR>my<SUP>''</SUP>+cy<SUP>'</SUP>+∑k(x*cos(θ<SUB>i</SUB>)+y*sin(θ<SUB>i</SUB>))<SUP>1.5</SUP>*sin((θ<SUB>i</SUB>)=F<SUB>u</SUB>cos(wt)<BR>式中m，c，k，W，Fu，w都是已知，θ<SUB>i</SUB>＝2*π*（i－1)/z＋Wc*t，z，Wc也是已知。<BR>在matlab下编写的程序如下：<BR>ai=2*π*（i－1)/z＋Wc*t<BR><FONT color=#ff0000>xdot=[x(2);(Fx+Fu*cos(w*t)-c*x2-sum((ko*(x(1)*cos(ai)+x(3)*sin(ai)).^1.5).cos(ai)))/m;...<BR>        x(4);(Fy+Fu*sin(w*t)-c*x2-sum((ko*(x(1)*cos(ai)+x(3)*sin(ai)).^1.5).sin(ai)))/m]<BR></FONT>然后用ode45求解，运行程序后，总是提示：Error：")" expected, "." found</DIV>
<P>这句有语法错误，自己先查括号之类的是否对应</P>