PTT數位生活區 / Programming

Fw: [求救 懸賞5000p] 一個7x7的矩陣 =.=

看板Programming作者 (綠扁帽)時間6年前 (2018/06/22 16:53), 5年前編輯推噓2(2推 0噓 8→)
留言10則, 3人參與, 6年前最新討論串1/1
各位請看，這是有限差分法的隱式解法 https://i.imgur.com/9leImX6.png
我們可以把它 拆解成下面的矩陣 https://i.imgur.com/WwrrD7A.png
我設定 係數A=1 係數B=0.8 係數C=1 然後最右邊那一列的C1*、C2*、C3*.....到C7* 都是1 所以最後矩陣相乘的結果Ｄ，每一個的計算結果都是(1*1)+(1*0.8)+(1*1)＝2.8 =================================== 這是一個迭代的矩陣，計算出每一組Ｄ以後， 就把Ｄ當成C1*、C2*～～C7*代入，再計算出下一組Ｄ； 然後再把Ｄ當成C1*、C2*～～C7*代入，再計算出下一組Ｄ；一直迭代下去 重點來了， 如果我把中間的那一行(第四行)的A、B、C 改成 0、1、0 這樣子，會導致Ｄ的數據，跟原來的不一樣 然後Ｄ接下來又被當作C1*、C2*～～C7*代入，又影響到下一組的Ｄ 然後下一組的Ｄ又當作C1*、C2*～～C7*代入，又影響到下下一組的Ｄ 一直影響下去。 ==================================== 總共要迭代幾次、或者要帶什麼數字，您們都可以自行決定。 我只是想要瞭解這個程式的寫法而已，數字不是重點，可自行決定數字大小。 感謝 如果看不懂我文中的"迭代"的意思，可以參閱下面兩張圖。 https://imgur.com/a/hlehNcU (求出最左邊的那一列 0,0,1,8,1,0,0) 把那一列代入到下面這張圖的最右邊 https://imgur.com/a/2enLKV3 求出最左邊的那一列(0 , 0.1 , 1.6 , 6.6 , 1.6 , 1 , 0.1 , 0) 以此類推 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 60.91.73.136 ※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1529653193.A.53B.html
→
06/22 16:44, 6年前 , 1F
while ess>1e-3;x = W*x0;ess = norm(x-x0);x0 = x;end
06/22 16:44, 1F
→
06/22 16:47, 6年前 , 2F
起始猜值x0跟係數矩陣W 在while外要先給好
06/22 16:47, 2F
→
06/22 16:48, 6年前 , 3F
W可以用diag創建
06/22 16:48, 3F
→
06/22 16:50, 6年前 , 4F
誤差也要先給一個能跑的起始值 大概就這樣
06/22 16:50, 4F
※ 編輯: GreenBeret (114.38.108.68), 06/22/2018 16:54:10
推
06/28 08:12, 6年前 , 5F
不好意思，請問你到問題是？
06/28 08:12, 5F
HI~~ e大大 後來在matlab版有問到解法了XD https://www.ptt.cc/bbs/MATLAB/M.1529667860.A.921.html 我敘事能力有點差QQ 您看他的解法 會比我用說的清楚XD
推
06/29 00:41, 6年前 , 6F
一樓有matlab的解法了
06/29 00:41, 6F
→
06/29 00:42, 6年前 , 7F
06/29 00:42, 7F
→
06/29 00:42, 6年前 , 8F
python with numpy
06/29 00:42, 8F
→
06/29 00:44, 6年前 , 9F
矩陣需存在要如下的檔案
06/29 00:44, 9F
→
06/29 00:44, 6年前 , 10F
06/29 00:44, 10F
! 其實1樓的方法 我現在才驚覺看懂 之前沒注意囧 ※ 編輯: GreenBeret (111.246.27.144), 06/29/2018 02:04:20 ※ 編輯: GreenBeret (60.91.73.136), 05/17/2019 15:58:58
‣ 更多 GreenBeret 的文章
文章代碼(AID): #1RBBZ_yt (Programming)
Programming 近期熱門文章
更多 近期熱門文章 >>
PTT數位生活區 即時熱門文章
更多 即時熱門文章 >>
‣ 更多 GreenBeret 的文章
文章代碼(AID): #1RBBZ_yt (Programming)