| 查看: 426 | 回复: 3 | ||
[求助]
编写控制程序
|
|
M文件代码如下： function dx=guidance(t,x,u) if 40 elseif 20 elseif 0 end dx=[-10.0211*u*x(1)^2/2200-3.716*3395^2*sin(x(2))/(3395+x(3))^2; -0.837*u*x(1)/2200-3.716*3395^2*cos(x(2))/x(1)/(3395+x(3))^2+x(1)*cos(x(2))/(3395+x(3)); x(1)*sin(x(2)); 3395*x(1)*cos(x(2))/(3395+x(3))]; end 命令行如下： x=ones(1,4);t=1;u=2;dx=guidance(t,x,u); >> x0=[6000;-15/180;125;0]; >> [t,y]=ode45(@guidance,[0,100],x0); 结果报错： Warning: Failure at t=6.283114e-220. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.258037e-234) at time t. > In ode45 at 371 请问该如何解决？ 还有一部分是如何在终点时间 tf 时，使x3（tf）=7，x1（tf）=400，请问该如何编写代码？ |
回复此楼» 猜你喜欢
274求调剂
已经有20人回复
-
265求调剂
已经有16人回复
-
材料与化工304求B区调剂
已经有6人回复
-
321求调剂
已经有5人回复
-
085600材料与化工306
已经有7人回复
-
085602化学工程求调剂。
已经有4人回复
-
312求调剂
已经有10人回复
-
材料考研求调剂
已经有6人回复
-
调剂推荐
已经有5人回复
-
中国科学院深圳先进技术研究院-光纤传感课题组招生-中国科学院大学、深圳理工大学联培
已经有4人回复
ifly21cn
银虫 (初入文坛)
- 应助: 7 (幼儿园)
- 金币: 764.7
- 散金: 10
- 帖子: 35
- 在线: 324.1小时
- 虫号: 1634933
- 注册: 2012-02-22
- 性别: GG
- 专业: 控制理论与方法
2楼2012-11-06 09:14:47
3楼2012-11-08 08:28:51
ifly21cn
银虫 (初入文坛)
- 应助: 7 (幼儿园)
- 金币: 764.7
- 散金: 10
- 帖子: 35
- 在线: 324.1小时
- 虫号: 1634933
- 注册: 2012-02-22
- 性别: GG
- 专业: 控制理论与方法
4楼2012-11-08 10:38:43
