24小时热门版块排行榜    

查看: 429  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

657801288

铁虫 (初入文坛)

[求助] 编写控制程序

M文件代码如下:
function dx=guidance(t,x,u)
if 40     u=3.7*10e-4*exp(-1.371*10e-4*(x(3)-40*10e+3));
elseif 20     u=2.8*10e-3*exp(-1.92*10e-4*(x(3)-20*10e+3));
elseif 0     u=0.01417*exp(-0.8109*10e-4*x(3));
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,请问该如何编写代码?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

657801288

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by ifly21cn at 2012-11-06 09:14:47
=ode45(@guidance,,x0);
改成
=ode45(@guidance,,x0);
试试看。
积分的步长太短了

这样还是出现了问题。
??? Out of memory. Type HELP MEMORY for your options.
3楼2012-11-08 08:28:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

ifly21cn

银虫 (初入文坛)

【答案】应助回帖

[t,y]=ode45(@guidance,[0,100],x0);
改成
[t,y]=ode45(@guidance,[0:0.000001:100],x0);
试试看。
积分的步长太短了
2楼2012-11-06 09:14:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ifly21cn

银虫 (初入文坛)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 657801288 at 2012-11-07 18:28:51
这样还是出现了问题。
??? Out of memory. Type HELP MEMORY for your options....

步长小,内存不够用了。增加matlab的虚拟内存
4楼2012-11-08 10:38:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286求调剂 +12 Faune 2026-03-21 12/600 2026-03-26 19:49 by 不吃魚的貓
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +7 丹青奶盖 2026-03-26 7/350 2026-03-26 19:18 by macy2011
[考研] 266分求材料化工冶金矿业等专业的调剂 +3 哇呼哼呼哼 2026-03-26 3/150 2026-03-26 19:16 by JourneyLucky
[考研] 085602 289分求调剂 +8 WWW西西弗斯 2026-03-24 8/400 2026-03-26 16:33 by 不吃魚的貓
[考研] 考研调剂 +8 小蜡新笔 2026-03-26 8/400 2026-03-26 16:18 by dick_runner
[考研] 求调剂 +3 芦lty 2026-03-25 4/200 2026-03-25 23:25 by 芦lty
[考研] 调剂310 +3 温柔的晚安 2026-03-25 4/200 2026-03-25 23:16 by peike
[考研] 考研调剂 +6 来好运来来来 2026-03-21 7/350 2026-03-25 22:43 by 418490947
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研经历 +6 cccchenso 2026-03-23 6/300 2026-03-25 22:25 by 544594351
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +7 niko- 2026-03-22 7/350 2026-03-25 20:14 by qingfeng258
[考研] 333求调剂 +6 87639 2026-03-21 11/550 2026-03-25 16:17 by 87639
[考研] 293求调剂 +7 加一一九 2026-03-24 7/350 2026-03-25 12:02 by userper
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 生物学一志愿985,分数349求调剂 +6 zxts12 2026-03-21 9/450 2026-03-23 18:37 by macy2011
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 一志愿华中科技大学071000,求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
信息提示
请填处理意见