tusvd.asv

图像的奇异值分解

A=[1 2 1;2 3 2;1 2 1 ];
[u d v] = svd(A);
[nx,ny] = size(A);
I=eye(nx,nx);
for i=nx:-1:2
    I(i,i) = 0;
end
B=u*I*d*v';
B

[X,map] = imread('G:\cameraman.bmp');
 X= rgb2gray(X);
X=double(X);
A=X(100:250 ,100:250);
[nx,ny] = size(A);
I=eye(nx,nx);
for i=nx:-1:21
    I(i,i) = 0;
end
[u d v]=svd(A);
B=u*I*d*v';
subplot(211)
imshow(A,map);
subplot(212)
imshow(B,map);