方程如下:<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