作者zupo (外派越南工作者)
標題Re: [討論] 網格法求解問題
時間Tue May 31 22:00:42 2016
※ 引述《danny555361 ()》之銘言:
: 小弟的熱傳作業需要用網格法解每個node上的溫度,我現在手上有28條方程式&28個變數
: 方成組長的全部都想這樣 △x △y 是已知 T(m,n)是我想要解的
: (△y/△x)*T(1,2)+(△x/△y)*T(2,1)-((△y/△x)+(△x/△y))*T(1,1)=0
: 我把每個網格的溫度令成T(1,1) T(1,2).....T(6,5)
: 想說最後用solve來解28條方程式
: 最後想要解方程組的時候發現
: syms T(1,1) 無法編譯
: 所以有想到辦法是不是把T(m,n)都換成Tmn這樣似乎就行了
: 可是這樣我方程式全部都要重寫
: 想問大家有沒有其他好辦法
: 先謝謝各位的建議
用差分法的話是可以解熱傳方程式
不過建議還是把原始算式放上來比較好幫忙
另外邊界條件是? 沒有這些也沒有辦法解
你這是 dT/dx + dT/dy = f(T,x,y) 的形式吧?
總之還是先把原始方程式放上來看看
--
數值方法 : 求根值 . ODE系統 . DDE系統. DAE系統. PDE系統
非線性聯立系統. 化工程序控制系統 . 力學系統
RLC 系統 . MCK 系統 . Laplace 轉換 系統
其他 : 工程數學 .微積分.化工程序.輸送現象.化工程序線性控制
Simulator . MATLAB tool box . Simulink . FORTRAN
不會影像處理.....囧>
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 14.163.255.245
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1464703247.A.3FC.html
→ rex0707: 熱傳導是 dT/dt = k_x*(d^2T/dx^2) + k_y*(d^2T/d^2) 吧 06/01 00:50
→ rex0707: 後面可能再多一項source 06/01 00:51
→ zupo: 有一種是方向熱傳,解反矩陣即可求解 06/02 23:37
→ rex0707: 恩...我只聽過熱傳導、熱對流、熱輻射... 06/03 00:23
→ rex0707: 解矩陣是隱式方法 不用解矩陣是顯式方法 看原原PO用哪種 06/03 00:24