数字图像处理DCT矩阵分块变换matlab程序

function re=hw122（）
%说明：hw122为对图像进行系数20%33%50%随机保留的DCT以及IDCT变换
I=imread（‘aaa2.jpg‘）;
temp=I;

s=size（I）;

A=randint（s（1）s（2）[04]）;		
%生成和s大小一致的随机整数矩阵A，其中0至4的五个整数值各占所有元素的20%
for k=1:1:s（1）
   for t=1:1:s（2）
     if A（kt）~=1 
         A（kt）=0;	%对于A中非1元素全部置零
     end
   end
end
i=dct2（I）;
I=i.*A;					%按20%的比例随机保留dct矩阵中的系数
I2=idct2（I）;I3=idct2（i）;
I2=uint8（I2）;I3=uint8（I3）;	
%转换数据类型
figure（1）;
imshow（temp）;title（‘原始图像‘）;
figure（2）;
imshow（I3）;title（‘DCT系数全部保留的重构图像‘）;
  imwrite（I3‘DCT系数全部保留的重构图像.jpg‘）;
figure（3）;
imshow（I2）;title（‘DCT系数20%随机保留的重构图像‘）;%显示图像
  imwrite（I2‘DCT系数20%随机保留的重构图像.jpg‘）;

A=randint（s（1）s（2）[02]）;		
%生成和s大小一致的随机整数矩