是這樣的
我要做二維圖形的運算
先把它做成4*4的矩陣 矩陣內的每個數值都是401*401的矩陣
然後因為有用到特徵值 所以我是寫迴圈讓他每個點逐項運算
再把它算結果出來的mSL也變成一個4*4的矩陣 矩陣內的每個數值都是401*401的矩陣
參考了一下前面的留言 這是我目前的式子
for i=1:401
for j=1:401
MSL=[m11(i,j) m12(i,j) m13(i,j) m14(i,j);
m21(i,j) m22(i,j) m23(i,j) m24(i,j);
m31(i,j) m32(i,j) m33(i,j) m34(i,j);
m41(i,j) m42(i,j) m43(i,j) m44(i,j)]/m11(i,j);
[VMSL,DMSL]=eig(MSL);
MSL_lamda=logm(DMSL);
mSL=VMSL*MSL_lamda*inv(VMSL);
Ans(i,j)={mSL}
end
end
for i= 1:401
for j=1:401
mSL11(i,j)=Ans{i,j}(1,1);
mSL12(i,j)=Ans{i,j}(1,2);
...從11~44
end
end
結果他就從10點跑到現在還沒跑完.....
請問有什麼比較簡略的方法嗎 謝謝大家