矩阵A的特征值分解:
a=[1 2;3 4];
[u,v]=eig(a);
%% eig test
test1_1=u*v*inv(u);
test1_2=u*v*u';
test1_3=sqrt(sum(u(:,1).^2));
test1_4=u'*u;
test1_5=u(:,1)'*u(:,2);
test1_1 =
1.0000 2.0000
3.0000 4.0000
test1_2 =
0.6765 2.2059
2.2059 4.3235
test1_3 =1
test1_4 =
1.0000 -0.1715
-0.1715 1.0000
test1_5 =
-0.1715
其中u不是标准正交的。
%% svd test
[u,s,v]=svd(a);
test2_1=u*s*v';
test2_2=v*v';
test2_3=u*u';
test2_1 =
1.0000 2.0000
3.0000 4.0000
test2_2 =
1.0000 -0.0000
-0.0000 1.0000
test2_3 =
1.0000 -0.0000
-0.0000 1.0000