作者obear7852 (歐巴耳)
看板MATLAB
標題[問題]矩陣內存入方程式再解聯立
時間Thu Nov 10 14:18:03 2011
要解的題目本身是N個方程式解N-2個未知數,
而這些未知數是有相關性的,
題目是二維熱傳。
我用矩陣的原因是每個矩陣內的元素跟相鄰的元素有相關,
相關就是上下左右項相加再除4,
一開始我將矩陣四周的邊界條件設定好,
用for迴圈列出了矩陣內每一項內含的方程式,
希望for迴圈列好後,然後解每一項的。
但是matlab輸出傾向數值解,
因此在每一項的方程式列好時他就幫我算出那項的解了,
拿4*4的矩陣舉例:
1 1 1 1
0 x1 x2 8
0 x3 x4 7
0 1 1 1
x1-x4就是該項上下左右相加在除四,
因此應該會有四條關於x1.x2.x3.x4的方程式可以解出x1.x2.x3.x4
但是跑出來卻是x1=1/4 我想是因為當x3.x2還沒設定時,算x2及x3都是0
他就先幫我算,但這不是我要的...
有可能說明的不太完整,希望有板友能幫幫忙,跨過這個障礙。
對題目有疑問我會再做說明,我的矩陣大約是50*50*50的三維矩陣。
謝謝!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.116.201.102
1F:→ Keiichi:重新將equation整理成 Ax=B 的形式 填入A跟B,解x 11/10 17:05