作者diva (Tremendo!)
看板MATLAB
標題Re: [討論]矩陣分解 用for loop
時間Mon May 21 22:15:58 2012
% 建立你原文說的矩陣
M=1:16 ;
M=reshape(M,4,4) ;
% 把每個2*2矩陣的第1項提出
A=M(1:2:end,1:2:end) ;
A=reshape(A,1,numel(A)) ;
% 把每個2*2矩陣的第2項提出
B=M(2:2:end,1:2:end) ;
B=reshape(B,1,numel(B)) ;
% 把每個2*2矩陣的第3項提出
C=M(1:2:end,2:2:end) ;
C=reshape(C,1,numel(C)) ;
% 把每個2*2矩陣的第4項提出
D=M(2:2:end,2:2:end) ;
D=reshape(D,1,numel(D)) ;
% 重新整理為一個比較好求平均值的矩陣
MM=[A;B;C;D] ;
% 求平均
mean(MM)
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 59.115.98.36
1F：→ diva:8*8的矩陣應該也可以用 05/21 22:20
2F：→ diva:如果你一定要用loop請再回報 05/21 22:21
3F：推 simon0216117:感謝diva的解惑,不過是一個100*100以上的矩陣,我想要 05/22 00:29
4F：→ simon0216117:需用loop吧!我已經改成此形式,A=magic(4), 05/22 00:31
5F：→ simon0216117:B=reshape(A,2,[]),可分成奇數列行和偶數行,在整理合 05/22 00:32
6F：→ simon0216117:併,不過還沒寫出來,還是要寫出loop才行,謝謝你的幫! 05/22 00:33