24小时热门版块排行榜    

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

akang9

铜虫 (初入文坛)

[求助] 求高手解答Matlab求解微分方程的初始值和边界值问题

为了一项工程模拟,有四个微分方程,有三个已知初始值,另一个只有终值,怎么编程求解?我也试过对未知初始值的方程假设一个初始值,用ode45求解验证计算的终值与实际值是否吻合的方法,换了很多假设值,都得不到合理的结果,甚至曲线趋势与实际完全相反。我该怎么办?

比如:
f(1)=-(y(1)-20)*3.8*pi+5*y(4);
f(2)=0.1*y(2)/(0.003*100);
f(3)=y(2)/(0.003*y(1));
f(4)=0.55*(1-y(4))/(y(3)/60)*exp(-4*log(2)*(y(1));

已知初值和终值为
y(1)0=220
y(2)0=14
y(3)0未知,y(3)终值为1000
y(4)0=0
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

akang9

铜虫 (初入文坛)

引用回帖:
3楼: Originally posted by dingd at 2013-02-03 09:20:28
两点边值问题,Matlab有个BVP4C专门求解这种问题,有点麻烦,不如1stOpt来的直观:

Variable  x=,y1=220,y2=12,y3=,y4=0;
Plot y1,y3;
ODEFunction y1'=-(y1-20)*3.8*pi+5*y4;
            y2'=0.1*y2/(0.003 ...

dingd您好!您给我的code直接能够在命令窗口运行吗?

直接运行的话,跳出这样的错误信息:
??? Undefined function or method 'Variable' for input arguments of type 'char'.

我的版本是R2010b。
4楼2013-02-04 21:49:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

dingd

铁杆木虫 (职业作家)

ODEFunction y1'=-(y1-20)*3.8*pi+5*y4;
            y2'=0.1*y2/(0.003*100);
            y3'=y2/(0.003*y1);
            y4'=0.55*(1-y4)/(y3/60)*exp(-4*ln(2)*y1);
[/code]

边值估算:
        y3(x=0): 11.2471932817486

t1.jpg

2楼2013-02-03 09:19:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
两点边值问题,Matlab有个BVP4C专门求解这种问题,有点麻烦,不如1stOpt来的直观:
CODE:
Variable  x=[0,3],y1=220,y2=12,y3=[,1000],y4=0;
Plot y1,y3[y2];
ODEFunction y1'=-(y1-20)*3.8*pi+5*y4;
            y2'=0.1*y2/(0.003*100);
            y3'=y2/(0.003*y1);
            y4'=0.55*(1-y4)/(y3/60)*exp(-4*ln(2)*y1);

3楼2013-02-03 09:20:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

akang9

铜虫 (初入文坛)

引用回帖:
4楼: Originally posted by akang9 at 2013-02-04 21:49:13
dingd您好!您给我的code直接能够在命令窗口运行吗?

直接运行的话,跳出这样的错误信息:
??? Undefined function or method 'Variable' for input arguments of type 'char'.

我的版本是R2010b。...

dingd您给我的代码是在1stOpt中的代码吧。
但是下载了1.5Pro运行了,报错信息为

x=[0中变量定义有误!
3]中变量定义有误!
y3=[中变量定义有误!
1000]中变量定义有误!
5楼2013-02-04 22:16:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 4/200 2026-02-08 08:47 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 4/200 2026-02-08 08:46 by vs90ilomwc
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 5/250 2026-02-08 08:32 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:26 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:12 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:07 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 08:06 by vs90ilomwc
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:52 by vs90ilomwc
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:46 by vs90ilomwc
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:32 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:27 by vs90ilomwc
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:26 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:07 by vs90ilomwc
[硕博家园] 博士延得我,科研能力直往上蹿 +8 偏振片 2026-02-02 8/400 2026-02-08 06:52 by liyeqik
[教师之家] 有院领导为了换新车,用横向课题经费买了俩车 +7 瞬息宇宙 2026-02-04 7/350 2026-02-07 21:47 by tfang
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[基金申请] 同年申请2项不同项目,第1个项目里不写第2个项目的信息,可以吗 +4 hitsdu 2026-02-06 4/200 2026-02-07 13:07 by jurkat.1640
[考博] 天津大学招2026.09的博士生,欢迎大家推荐交流(博导是本人) +4 a793625982 2026-02-05 5/250 2026-02-07 10:57 by a793625982
[公派出国] CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生|MSCA经费充足、排名优 +4 雨念 2026-02-01 6/300 2026-02-06 23:32 by MelissaPon
[教师之家] 遇见不省心的家人很难过 +18 otani 2026-02-03 22/1100 2026-02-04 11:06 by tangmnt
信息提示
请填处理意见