| 查看: 495 | 回复: 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，请问该如何编写代码？ |
回复此楼» 猜你喜欢
综述论文作为代表作会不会影响评审专家的印象分？
已经有12人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有29人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有9人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
3楼2012-11-08 08:28:51
ifly21cn
银虫 (初入文坛)
- 应助: 7 (幼儿园)
- 金币: 764.7
- 散金: 10
- 帖子: 35
- 在线: 324.1小时
- 虫号: 1634933
- 注册: 2012-02-22
- 性别: GG
- 专业: 控制理论与方法
2楼2012-11-06 09:14:47
ifly21cn
银虫 (初入文坛)
- 应助: 7 (幼儿园)
- 金币: 764.7
- 散金: 10
- 帖子: 35
- 在线: 324.1小时
- 虫号: 1634933
- 注册: 2012-02-22
- 性别: GG
- 专业: 控制理论与方法
4楼2012-11-08 10:38:43
