24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1610  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

驿边桥CQU

新虫 (小有名气)

[求助] matlab求解边界值问题 已有1人参与

请问大家  我用bvp4c求解边界值问题时发现估测值对微分方程求解有很大影响,应该怎样确定合适的估测值呢

下面是程序:
function dy=wenxian(x,y)
Pr=7;Da=0.02;Ra=10^4;Ja=0.15;Sw=0.0;
dy=[y(2);(Ja/Da/Ra*(Sw*y(1)-1)-y(1)*y(2)^2)/y(1)^2];
function bc=wenxianbcr(ya,yb)
Pr=7;Da=0.02;Ra=10^4;
bc=[ya(2);yb(2)^2/yb(1)-Pr/Da^2/Ra];

clear all   %估测值为[10 65]时
solinit=bvpinit(linspace(0,1,10000),[10 65]);
sol=bvp4c(@wenxian,@wenxianbc,solinit);
x=linspace(0,1,10000);
y=deval(sol,x);figure(12);
plot(x,y(1,: ),'-o');title('10000,10,65')

clear all   %估测值为[10 40]时
solinit=bvpinit(linspace(0,1,10000),[10 40]);
sol=bvp4c(@wenxian,@wenxianbc,solinit);
x=linspace(0,1,10000);
y=deval(sol,x);figure(12);
plot(x,y(1,: ),'-o');title('10000,10,40')

matlab求解边界值问题
无标题.png


matlab求解边界值问题-1
无标题2.png

下面为原方程的描述:
matlab求解边界值问题-2
原方程描述.png
变形后:
matlab求解边界值问题-3
变形后.png
回复此楼

» 猜你喜欢

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

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

驿边桥CQU

新虫 (小有名气)

我确实当时没弄明白文章的意思,,谢谢了
5楼2015-04-08 21:11:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
确认方程和系数都无误吗?似乎没有精确解。好好再检查下吧。
2楼2015-04-08 11:26:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

驿边桥CQU

新虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2015-04-08 11:26:43
确认方程和系数都无误吗?似乎没有精确解。好好再检查下吧。

我是从一篇英文文献里找到的这个式子,他是算了出来的,但我总是重复不出来,,,它上面只是说用的是龙格库塔法跟射击法,,但具体就没说了  所以过来请大家帮下忙
3楼2015-04-08 14:54:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
驿边桥CQU: 金币+5, ★★★★★最佳答案 2015-04-08 21:13:01
xzhdty: 谢谢应助 2015-04-08 22:51:39
xzhdty: 金币+2 2015-04-08 22:51:55
第二个边值条件很难满足:Ja/Da/Ra*(Sw*y(1)-1)-y(1)*y(2)^2)/y(1)^2
如果变成:
Ja/Da/Ra*(Sw*y(1)-1)+y(1)*y(2)^2)/y(1)^2
就有比较好的解。
4楼2015-04-08 17:01:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂,一志愿郑州大学材料与化工专硕,英二数二342分,求老师收留 +8 v12abo 2026-04-02 8/400 2026-04-02 20:53 by 1104338198
[考研] 265求调剂 +12 林深温澜 2026-04-01 15/750 2026-04-02 19:14 by 一只好果子?
[考研] 环境工程297分求调剂一志愿杭高院 +15 GENJIOW 2026-03-31 16/800 2026-04-02 17:56 by cyh—315
[考研] 求调剂 302分初试 0854 +5 伶可乐 2026-04-02 5/250 2026-04-02 17:53 by 笔落锦州
[考研] 292分,材料与化工,申请调剂 +19 程晴之 2026-04-01 23/1150 2026-04-02 11:59 by 程晴之
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
[考研] 385分 生物学(071000)求调剂 +6 qf626 2026-04-01 6/300 2026-04-02 10:39 by guoweigw
[考研] 11408 321分求调剂 +3 huchun12138 2026-03-30 4/200 2026-04-01 22:48 by guanxin1001
[考研] 339求调剂,想调回江苏 +7 烤麦芽 2026-03-27 10/500 2026-04-01 21:35 by 495374996
[考研] 【求调剂】085601材料工程专硕 | 总分272 | +10 脚滑的守法公民 2026-03-27 10/500 2026-04-01 17:23 by pies112
[考研] 求调剂 +4 DADA怪 2026-03-31 4/200 2026-04-01 14:30 by ZXlzxl0425
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 11:35 by wangjy2002
[考研] 求调剂,一志愿北林食品与营养095500,301分,已过六级,有科研经历 +4 快乐储蓄罐 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
[考研] 329求调剂,一志愿西北工业大学,材料工程(085601) +6 小小机灵虫 2026-03-29 12/600 2026-03-31 16:58 by 记事本2026
[考研] 315求调剂 +6 akie... 2026-03-28 7/350 2026-03-31 16:48 by asdfzly
[考研] 085601 329分调剂 +6 yzsa12 2026-03-31 6/300 2026-03-31 15:23 by yanflower7133
[考研] 356求调剂 +4 gysy?s?a 2026-03-28 4/200 2026-03-29 10:32 by 唐沐儿
[考研] 求调剂 +7 争取九点睡 2026-03-28 8/400 2026-03-28 21:07 by 争取九点睡
[考研] 复试调剂 +3 raojunqi0129 2026-03-28 3/150 2026-03-28 15:27 by 落睿可思
[考研] 295求调剂 +5 1428151015 2026-03-27 6/300 2026-03-28 04:04 by fmesaito
信息提示
请填处理意见