24小时热门版块排行榜    

查看: 2109  |  回复: 5

dingd

铁杆木虫 (职业作家)

[交流] 1stOpt求解常微分方程边值问题 已有3人参与

1stOpt 5.0已出了,在解决常微分方程边值问题功能方面有了很大提高,远比Matlab自带的BVP4C等求解器简单好用的多,下面是一实例:
微分方程组:
           dy1/dx=y3;
            dy2/dx=-y4;
            dy3/dx=-exp(-x)+y1*y2/(y2+y1);
            dy4/dx=exp(-x)-y1*y2/(y2+y1);
积分区间:x = [0,1]
已知边界条件:
1:  x = 0 时,y1 = 1,y4 = 0 或者:y1(x=0) = 1, y4(x=0) = 0;
2:  x = 1 时,y3 = 0 或者: y3(x=1) = 0
3:  y3(x=0) = y4(x=1)
4:  dy2/dx(x=0.5) + y1(x=0.4) = 0.6

1stOpt求解代码:
CODE:
ODEStep = 0.1;
Variable x,y1,y3,y4;
SubjectTo y3[0]=y4[1], y2'[0.5]+y1[0.4]=0.6;
ODEFunction y1'=y3;
            y2'=-y4;
            y3'=-exp(-x)+y1*y2/(y2+y1);
            y4'=exp(-x)-y1*y2/(y2+y1);

结果:
y3初值        6.12634551613304
y2初值        -0.988421760677377

微分方程拟合约束(SubjectTo):
y3[初值]-(y4[1]): -8.88178419700125E-16
y2'[0.5]+y1[0.4]-(0.6): 5.32885367776625E-16


x        计算y1        计算y3        计算y4        计算y2'
0.1        1.55036161834342        5.22444856226826        0.901896953864782        -0.901896953864782
0.2        2.05364638237305        4.85196325368909        1.27438226244395        -1.27438226244395
0.3        2.52162389318013        4.50818981022655        1.61815570590649        -1.61815570590649
0.4        2.95496040832548        4.15507786835396        1.97126764777908        -1.97126764777908
0.5        3.35162040499198        3.77138510780756        2.35496040832548        -2.35496040832548
0.6        3.70752818818257        3.33619776453884        2.7901477515942        -2.7901477515942
0.7        4.01623535985374        2.82175078568221        3.30459473045083        -3.30459473045083
0.8        4.267774566885        2.18291509351099        3.94343042262205        -3.94343042262205
0.9        4.44578141831917        1.32889337172753        4.79745214440551        -4.79745214440551
1        4.51847089112852        4.16333634234434E-17        6.12634551613304        -6.12634551613304
回复此楼

» 收录本帖的淘帖专辑推荐

北大科研分享

» 猜你喜欢

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

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

upboxer

木虫 (正式写手)


小木虫(金币+0.5): 给个红包,谢谢回帖
LZ的1stOpt 5.0在哪找的啊?
轻轻的我来了正如我轻轻的走,像一阵风,希望给朋友带来的幸福,带走的是烦恼...
2楼2012-03-14 14:10:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

只有买啊!
3楼2012-03-14 14:54:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

广告贴,绝对是广告贴
4楼2012-03-14 18:00:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

手头还有正版的GAMS,4万多美元买的!
想做项目,发文章(尤其是国外),购买正版是唯一途径,即帮别人更是帮自己。
5楼2012-03-16 09:04:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhdgzhdg

至尊木虫 (著名写手)

5.0要上万啊
6楼2013-03-05 06:50:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dingd 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +9 majunge000 2026-08-11 11/550 2026-08-15 23:20 by botar
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见