24小时热门版块排行榜    

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

dingd

铁杆木虫 (职业作家)

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

upboxer

木虫 (正式写手)


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

change0618

铁杆木虫 (著名写手)

方丈大师

广告贴,绝对是广告贴
4楼2012-03-14 18:00:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 一志愿武理材料工程348求调剂 +3  ̄^ ̄゜汗 2026-03-19 4/200 2026-03-20 21:01 by zhukairuo
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
信息提示
请填处理意见