运动模糊图像倒谱程序

 %% 读入并显示图像
 filename = ‘lena1.jpg‘;
 I = imread（filename）;
figure
  imshow（uint8（I））;
 title（‘原图‘）; 
 %% 生成运动模糊图像
 PSF = fspecial（‘motion‘80 150）;
 g = imfilter（I PSF ‘circular‘）;
 figure 
 imshow（uint8（g））;
 title（‘运动模糊图‘）; 
 %% 对运动模糊图像进行灰度化，并进行二维快速傅里叶变换，生成其频谱图
 gb = rgb2gray（g）;
 figure
 imshow（uint8（gb））;
 PQ = paddedsize（size（gb））;
 F = fft2（gb PQ（1） PQ（2））;
 figure
 imshow（uint8（F））; 
 %% 作出倒频谱
 F1 = log（1+abs（F））;
 F2 = abs（F1）.^2;
 F3 = real（ifft2（F2））;
 figure
 imshow（uint8（F3））; 
 %% 将倒频谱压缩，居中
 H = log（1+abs（