%程序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矩陣的大小