24小时热门版块排行榜    

查看: 3467  |  回复: 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的回帖

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

铜虫 (初入文坛)

引用回帖:
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的回帖

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的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

内容已删除
6楼2013-02-04 23:34:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 akang9 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 2rKfVQg2jQww 2026-11-06 15/750 2026-11-09 05:37 by l48hFd8Nwae2
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 5sJFDKp2kwuI 2026-11-06 13/650 2026-11-09 05:15 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 16:32 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
请填处理意见