l4_1.m

DigitalImageProcessing_base_on_Matlab 基于Matlab的数字图像处理

%程序L4_1.m: SVD变换%
A=[ 1 2 1;2 3 2;1 2 1 ];      %宣告一个矩阵A
[u d v]=svd(A);               %对A矩阵进行SVD变换 
[nx,ny]=size(A);              %算出A 矩阵的大小
I=eye(nx,nx);                 %取出一个单位矩阵，其大小为nx×nx
for i=nx:-1:2                
  I(i,i)=0;                   %将对应到的主对角元素设为0
end      
B=u*I*d*v';                   %忽略第二和第三个特征值所还原的矩阵 
B                             %将还原后的矩阵显示出来 

%对实际图像进行SVD变换%
[X,map]=bmpread('L4_1.bmp');  %读取图像
A=X(120:220,120:220);         %将图像的一部分取出，并命名为A 
[nx,ny]=size(A);              %算出A矩阵的大小