[問題] 迴圈值輸出成矩陣

[feather0911 (我親愛的偏執狂)]

是這樣的
我要做二維圖形的運算
先把它做成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點跑到現在還沒跑完.....

請問有什麼比較簡略的方法嗎  謝謝大家

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 111.255.46.190
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1422372392.A.AD9.html
推 sunev: 都diagonalize了還用logm?