推 nemies : 那我想請教一下，那種跑好幾天才跑出結果的是甚麼 10/22 23:13
→ nemies : 樣的數學式呢？ 10/22 23:13
推 rex0707 : 我寫程式解Navier-Stokes 網格多就要跑數天到數月 10/23 00:42
※ 引述《nemies (...)》之銘言：
: 宣告 diff(x,t,2) 為 x(t) 被t微分兩次的寫法
: diff(x,t,1) 一次
: diff(y,t,2) y(t)
: 想請問有兩個微分方程式
: diff(x,t,2) = diff(x,t,1) + diff(y,t,1) + x + y + 1
: diff(y,t,2) = diff(y,t,1) + x*y*diff(x,t,1) + x*sin(y) + 5
: 要怎麼用數值的方法去解呢？
: 感謝......
照你的描述
原式
x" = x' +y' +x +y +1
y" = xy x' +y' +sin(y) x +5
為二階方程組
每個函數需要兩個"初始條件" 或"邊界條件"
你要解的問題是 IVP, Initial Value Problem
所以你要給定
x(0)= , x'(0)= , y(0)= , y'(0)=
一共四個初始條件 才有辦法求解
數值方法
令
x' = u
y' = v
u' = x" = x' +y' +x +y +1
v' = y" = xy x' +y' +sin(y) x +5
一階方程組及 初始條件 x(0)= , u(0)= , y(0)= , v(0)=
MATLAB有很好用的ode函式
應該不會花太多時間
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 70.95.164.2
※ 文章網址: https://www.ptt.cc/bbs/Math/M.1477126991.A.23E.html