使用svd分解将图片压缩

奇异值分解
原理：其中U是一个m * m的矩阵，Σ是一个m * n的矩阵，除了对角线上的元素以外全为0，主对角线上的每个元素都称为奇异值，V是一个n * n的矩阵。U和V都是酉矩阵。

具体的代码如下

k=144;
I=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part1\image.jpg’);
h=size(I,1);
w=size(I,2);
R=I(:,:,1);
G=I(:,:,2);
B=I(:,:,3);

Rd=im2double®;
Gd=im2double(G);
Bb=im2double(B);%转换为双精度

[Ur,Sr,Vr]=svd(Rd);%svd分解
[Ug,Sg,Vg]=svd(Gd);
[Ub,Sb,Vb]=svd(Bb);

Rt=zeros(h,w);%生成svd重构后的矩阵
Gt=zeros(h,w);
Bt=zeros(h,w);

for i=1:k
Rt=Rt+Sr(i,i)*Ur(:,i)*Vr(:,i)’;
Gt=Gt+Sg(i,i)*Ug(:,i)*Vg(:,i)’;
Bt=Bt+Sb(i,i)*Ub(:,i)*Vb(:,i)’;%分量求和
end
I2(:,:,1)=im2uint8(Rt);
I2(:,:,2)=im2uint8(Gt);
I2(:,:,3)=im2uint8(Bt);%转换为uint8
imshow(I2);
结果如下：