| 查看: 1776 | 回复: 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求解代码： 结果： 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 北大科研分享 |
» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
-
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
物理学I论文润色/翻译怎么收费?
已经有119人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有4人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求微分方程 最大值问题
已经有4人回复
- 一个1stOpt的初级问题：复数方程组计算问题（有代码）
已经有5人回复
- 急求，常微分方程的参数计算问题，有1stopt的帮忙解答一下
已经有5人回复
- 1stopt 计算问题
已经有7人回复
- 谁用1stopt这个软件啊，离散数据拟合方程，方程未知
已经有5人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- 一个二阶常微分方程的求解方法？
已经有5人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- 谁有1stopt，帮忙跑一下啊
已经有6人回复
- 求高版本1stopt，解决一个问题，微分方程组参数拟合
已经有3人回复
- 机械系统动力学微分方程有了以后，如何计算固有频率呢？
已经有4人回复
- 用Wofram Mathematica解方程的问题
已经有9人回复
- 1stopt4阶非线性常微分方程，帮小弟看一下，跪求。
已经有7人回复
- 有关磁流体方程计算的询问
已经有8人回复
- 最优化计算求助！（1stopt或者Matlab）
已经有16人回复
- 【求助】求助一个一阶常微分方程的初值问题
已经有16人回复
- 【求助】微分方程问题求解
已经有5人回复
upboxer
木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 2880.9
- 散金: 37
- 红花: 4
- 帖子: 498
- 在线: 135.5小时
- 虫号: 1169321
- 注册: 2010-12-13
- 性别: GG
- 专业: 水文地质学(含地热地质学)
2楼2012-03-14 14:10:04
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
3楼2012-03-14 14:54:51
change0618
铁杆木虫 (著名写手)
方丈大师
- 计算强帖: 1
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
4楼2012-03-14 18:00:35
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
5楼2012-03-16 09:04:55
zhdgzhdg
至尊木虫 (著名写手)
- 应助: 32 (小学生)
- 金币: 22701.4
- 红花: 6
- 帖子: 1750
- 在线: 973.9小时
- 虫号: 95566
- 注册: 2005-10-05
- 专业: 化学反应工程
6楼2013-03-05 06:50:26
