是这样的
我要做二维图形的运算
先把它做成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点跑到现在还没跑完.....
请问有什么比较简略的方法吗 谢谢大家