[問題] 矩陣裡面放矩陣已回收

看板MATLAB作者 (GOGOGO)時間16年前 (2009/03/28 19:10), 編輯推噓2(204)
留言6則, 4人參與, 最新討論串1/1
現在想解一個2D的poisson equation 需要一個矩陣A大小是 (n^2,n^2) 一個矩陣B大小是(n,n) 然後B包含在A裡面: A=[B -I 0 0 ... 0 -I B -I 0 ... 0 0 -I B -I ... 0 . . . . . . . . . . -I B -I 0 -I B ] 數字的話我會寫個迴圈,但是現在裡面都是n*n的矩陣 苦惱了很久 想請大大們幫個忙 謝謝:) -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 118.160.67.37

03/28 20:34, , 1F
help spdiags ?
03/28 20:34, 1F

03/28 22:22, , 2F
先用blkdiag把B塞進去,I的部份用diag去處理
03/28 22:22, 2F

03/29 15:38, , 3F
A=kron(eye(n),B)-kron(diag(ones(9,1),1)+...
03/29 15:38, 3F

03/29 15:39, , 4F
diag(ones(9,1),-1),I);
03/29 15:39, 4F

03/29 15:39, , 5F
抱歉,把9換成n-1即可
03/29 15:39, 5F

03/29 22:24, , 6F
謝謝各位:)
03/29 22:24, 6F
文章代碼(AID): #19pWMBtq (MATLAB)
文章代碼(AID): #19pWMBtq (MATLAB)