Re: [討論]矩陣分解 用for loop

[diva (Tremendo!)]

% 建立你原文說的矩陣
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
→ diva:8*8的矩陣應該也可以用
→ diva:如果你一定要用loop請再回報
推 simon0216117:感謝diva的解惑,不過是一個100*100以上的矩陣,我想要
→ simon0216117:需用loop吧!我已經改成此形式,A=magic(4),
→ simon0216117:B=reshape(A,2,[]),可分成奇數列行和偶數行,在整理合
→ simon0216117:併,不過還沒寫出來,還是要寫出loop才行,謝謝你的幫!