声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1052|回复: 1

[求助]请大家帮我看一下这个错误是怎么回事啊

[复制链接]
发表于 2006-3-29 15:20 | 显示全部楼层 |阅读模式

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

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

x
<P><STRONG>M文件:<BR></STRONG>function fun=xuanjian(t,x,xd)<BR>m1=2;<BR>m2=5;<BR>L=5;<BR>p=100;<BR>k=1000;<BR>u0=0.1;<BR>fun=[xd(1)-x(5);<BR>    xd(2)-x(6);<BR>    xd(3)-x(7);<BR>    xd(4)-x(8);<BR>    -2*x(5)-3*x(1)-2*xd(5);<BR>    -2*xd(6)-2*x(6)-3*x(2);<BR>    3/2*(3*L*m1*cos(x(3))*x(2)-3*L*m1*sin(x(3))*x(1)+2*L*m1*cos(x(3))*x(6)-2*L*m1*sin(x(3))*x(5)-L*m1*sin(x(3))*xd(5)+L*m1*cos(x(3))*xd(6)-4*m2*x(8)*x(7)*x(4)+2*p+6*cos(x(3))*m2*u0*x(2)+2*cos(x(3))*m2*x(4)*xd(6)+4*cos(x(3))*m2*x(4)*x(6)+6*cos(x(3))*m2*x(4)*x(2)-4*m2*x(8)*x(7)*u0-4*sin(x(3))*x(5)*m2*x(4)-6*sin(x(3))*x(1)*m2*x(4)-6*sin(x(3))*x(1)*m2*u0-2*sin(x(3))*xd(5)*m2*x(4)-4*sin(x(3))*x(5)*m2*u0-2*sin(x(3))*xd(5)*m2*u0+2*cos(x(3))*m2*u0*xd(6)+4*cos(x(3))*m2*u0*x(6))/(m1*L^2+3*m2*x(4)^2+6*m2*x(4)*u0+3*m2*u0^2)-xd(7);<BR>    -(k*x(4)-2*sin(x(3))*m2*x(6)-3*sin(x(3))*m2*x(2)-3*cos(x(3))*m2*x(1)-sin(x(3))*m2*xd(6)-m2*x(7)^2*x(4)-m2*x(7)^2*u0-cos(x(3))*m2*xd(5)-2*cos(x(3))*m2*x(5))/m2-xd(8)];<BR><STRONG>执行</STRONG><BR>x0=[1.0;0;0;0;10;15;12;20];xd0=[10;15;12;20;8;13;36;25];  %ding yi x0 he xd0<BR>res=ode15i('xuanjian',[0,10],x0,xd0)<BR><STRONG>结果<BR></STRONG>res = </P>
<P>     solver: 'ode15i'<BR>    extdata: [1x1 struct]<BR>          x: [1x1223 double]<BR>          y: [8x1223 double]<BR>      stats: [1x1 struct]<BR>      idata: [1x1 struct]</P>[em03]
回复
分享到:

使用道具 举报

 楼主| 发表于 2006-3-29 15:55 | 显示全部楼层
呵呵,不好意思,误会了,这不是错误,而是一个报告结果.我第一次见到不知道.还以为错了呢.其实是没问题的
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-9-26 01:18 , Processed in 0.057176 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表