应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stOpt求解常微分方程边值问题
61/1返回列表
查看: 1778  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

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

北大科研分享

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-02-21 14:40:44
已阅   回复此楼   关注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万多美元买的！
想做项目，发文章（尤其是国外），购买正版是唯一途径，即帮别人更是帮自己。
赞一下(1人) 回复此楼
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 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 296求调剂 +5 大口吃饭 身体健 2026-03-13 5/250 2026-03-17 21:05 by 不惑可乐
[考研] 0703化学调剂 ，六级已过，有科研经历 +8 曦熙兮 2026-03-15 8/400 2026-03-17 20:31 by xilongliang
[考研] 301求调剂 +5 yy要上岸呀 2026-03-17 5/250 2026-03-17 20:20 by peike
[考研] 281求调剂（0805） +3 烟汐忆海 2026-03-16 7/350 2026-03-17 20:16 by peike
[考研] 材料专硕306英一数二 +7 z1z2z3879 2026-03-16 9/450 2026-03-17 17:31 by ccjequ
[考研] 机械专硕325，寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 一志愿中科院，化学方向，295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
信息提示
关闭
请填处理意见
关闭