| 5 | 1/1 | 返回列表 |
| 查看: 1787 | 回復(fù): 5 | ||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
dingd鐵桿木蟲 (職業(yè)作家)
|
[交流]
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); 積分區(qū)間: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求解代碼: 結(jié)果: 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 |
北大科研分享 |
鐵桿木蟲 (職業(yè)作家)
木蟲 (正式寫手)

鐵桿木蟲 (著名寫手)
方丈大師
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程(專)一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-17 | 8/400 |
|
|
[考研] 317求調(diào)劑 +8 | 申子申申 2026-03-19 | 13/650 |
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 一志愿中國海洋大學(xué),生物學(xué),301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
[考研] 生物學(xué)調(diào)劑招人。! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|