PTT數位生活區 / MATLAB

[討論] 網格法求解問題

看板MATLAB作者時間9年前 (2016/05/29 23:47), 編輯推噓0(0推 0噓 12→)
留言12則, 2人參與, 最新討論串1/2 (看更多)
小弟的熱傳作業需要用網格法解每個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這樣似乎就行了 可是這樣我方程式全部都要重寫 想問大家有沒有其他好辦法 先謝謝各位的建議 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.113.250.44 ※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1464536827.A.13A.html
→
05/29 23:48, , 1F
寫成矩陣的形式 x=A\b
05/29 23:48, 1F
→
05/29 23:49, , 2F
你是用有限差分嗎? 怎麼覺得離散式寫的怪怪的
05/29 23:49, 2F
→
05/29 23:52, , 3F
沒錯 是有限差分 可能是第一個點兩邊絕熱 看起來怪怪
05/29 23:52, 3F
→
05/29 23:53, , 4F
感謝大大的建議想問有辦法可以用方程式直接求解嗎
05/29 23:53, 4F
→
05/29 23:54, , 5F
因為每個變數都可能會出現在好幾條方程式中
05/29 23:54, 5F
→
05/29 23:54, , 6F
我懂你的意思 正常的作法就是寫成矩陣的形式再解線性系統
05/29 23:54, 6F
→
05/29 23:54, , 7F
用矩陣的話我天生眼殘 蠻可能漏看的
05/29 23:54, 7F
→
05/29 23:55, , 8F
可以用Gauss-Seidel方法求解
05/29 23:55, 8F
→
05/29 23:57, , 9F
不過28條式子寫成矩陣 直接法求解就很好了
05/29 23:57, 9F
→
05/29 23:58, , 10F
ok 感謝 原來一般就是用換成矩陣來解
05/29 23:58, 10F
→
05/29 23:59, , 11F
直接法是x=A\b嗎
05/29 23:59, 11F
→
05/30 00:01, , 12F
我不知道在matlab裡x=A\b是用甚麼方法解
05/30 00:01, 12F
‣ 更多 danny555361 的文章
文章代碼(AID): #1NImxx4w (MATLAB)
討論串 (同標題文章)
以下文章回應了本文：
完整討論串 (本文為第 1 之 2 篇)：
排序： 最新先 | 最舊先 | 留言數
在新視窗開啟完整討論串 (共2篇)
MATLAB 近期熱門文章
更多 近期熱門文章 >>
PTT數位生活區 即時熱門文章
更多 即時熱門文章 >>
‣ 更多 danny555361 的文章
文章代碼(AID): #1NImxx4w (MATLAB)