这篇文章其实也是矩阵分析课程的一个作业,当时也是觉得矩阵奇异值分解挺有意思的,所以就在网上查找奇异值分解的相关应用。另外自己对图像处理也很感兴趣,因此就选择了图像压缩这个应用。写博客是个很好的学习方法,坚持写博客也是一个好习惯。下面就直接把论文贴上来。
附录1:奇异值分解几何意义程序
t=[0:0.01:2*pi];
x=cos(t);
y=sin(t);
%变换矩阵
A = [1 -2;1 2];
A = double(A);
%对A进行奇异值分解
[U S V] = svd(A)
s = svd(A);
%第一个分力
S1 = [s(1) 0;0 0];
A1 = U*S1*V'
%第二个分力
S2 = [0 0; 0 s(2)];
A2 =U*S2*V'
%画单位圆
subplot(2,2,1);
plot(x,y)
title('单位圆')
axis equal;
xlabe