作者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