?%%?讀入并顯示圖像
?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（