| 查看: 424 | 回复: 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，请问该如何编写代码？ |
回复此楼» 猜你喜欢
寻找调剂
已经有8人回复
-
289求调剂
已经有16人回复
-
312求调剂
已经有8人回复
-
263求调剂
已经有10人回复
-
一志愿C9材料与化工专业总分300求调剂
已经有6人回复
-
296求调剂
已经有5人回复
-
0703化学/290求调剂/本科经历丰富/工科也可
已经有5人回复
-
284求调剂
已经有4人回复
-
化学工程085602 305分求调剂
已经有12人回复
-
考研调剂
已经有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
