clcclear;
Im=imread（‘SVD去噪散斑10號.png‘）;
figureimshow（Im）;
title（‘原始模糊散斑‘）
Im=rgb2gray（Im）;
f=double（Im）;
g=fft2（f）;
lg=log（abs（g））;
g=fftshift（g）;
Og=abs（255*g/max（max（abs（g））））;
figureimshow（Og）;
title（‘模糊散斑頻譜‘）
cepg=ifft2（lg