clear?all;close?all;
Image_Input=imread（‘1.bmp‘）;
dims=ndims（Image_Input）;?????????????????????%獲取矩陣的維數
if?dims==3
???Image_Gray=rgb2gray（Image_Input）;???????????
elseif?dims~=2
????return;
end
[hw]=size（Image_Gray）;??????????????????
xr=8;
xl=8;
yt=8;
yb=8;
width=w-xr-xl;
height=h-yt-yb;
Image_Cut=imcrop（Image_Gray[xrytwidth-1height-1]）;???%??切除黑邊
figure（1）;
imshow（Image_Cut[]）;
[RowColumn]=size（Image_Cut）;
Image_FFT=fftshift（fft2（Image_Cut））;
if?rem（Row2）==0
????for?u=-Row/2:Row/2-1
????????if?rem（Column2）==0
????????????for?v=-Column/2:Column/2-1
????????????????H（u+Row/2+1v+Column/2+1）=exp（-0.001*（u^2+v^2）^（5/6））;
????????????end
????????else?
????????????for?v=-（Column-1）/2:（Column-1）/2
????????????????H（u+Row/2+1v+（Column-1）/2+1）=exp（-0.001*（u^2+v^2）^（5/6））;
????????????end
????????end
????end
else
????for?u=-（Row-1）/2:（Row-1）/2
????????if?rem（Column2）==0
????????????for?v=-Column/2:Column/2-1
????????????????H（u+（Row-1）/2+1v+Column/2+1）=exp（-0.001*（u^2+v^2）^（5/6））;
????????????end
????????else
????????????for?v=-（Column-1）/2:（Column-1）/2
????????????????H（u+（Row-1）/2+1v+（Column-1）/2+1）=exp（-0.001*（u^2+v^2）^（5/6））;
????????????end
????????end
????end
end
Degeneration=Image_FFT.*H;
I=ifft2（ifftshift（Degeneration））;
Image_Deg=0.5*（abs（real（I））+real（I））;
Image_Deg=uint8（Image_Deg）;
figure（2）;
imshow（Image_Deg[]）;??
Image_fft=fftshift（fft2（Image_Deg））;????%傅里葉變換平移
Image_abs=abs（Image_fft）;????%求頻譜
Image_norm=Image_abs/max（max（Image_abs））;?????%頻譜歸一化
Image_ln=log10（Image_norm）;?????%對數縮放處理
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?Image_array=Image_ln（fix（Row/2）1:end）;??????%獲取水平中心處平頻譜
%?if?rem（Column2）==0????????????????????????????%使用中心零點對稱坐標顯示
%?????X=-fix（Column/2）+1:1:fix（Column/2）;
%?else
%?????X=-fix（Column/2）:1:fix（Column/2）;
%?end
%?figure（3）;
%?plot（XImage_array）;
%?xlabel（‘離散頻率‘）
%?ylabel（‘傅里葉頻率歸一化后取自然對數‘）;
%?Datablur=Image_array（fix（Column/2）:end）;
%?figure（4）;
%?semilogx（Datablur‘b:*‘）;
%?xlabel（‘離散頻率的自然對數‘）;
%?ylabel（‘傅里葉頻率歸一化后取自然對數‘）;
%?hold?on;
%?X_half=fix（Column/2）:Column;
%?count=size（X_half）;
%?k=1:count（2）;
%?p=polyfit（log10（k）Datablur1）;
%?line=polyval（plog10（k））;
%?semilogx（line‘r:*‘）;
%?hold?off;
if?rem（Row2）~=0?????????%獲取水平中心處平頻譜
????Image_array=Image_ln（（Row+1）/21:end）;??????
else
????Image_array=（Image_ln（Row/21:end）+Image_ln（Row/2+11:end））/2;
end
if?rem（Column2）==0????????????????????????????%使用中心零點對稱坐標顯示
????Data_X=-Column/2:Column/2-1;
else
????Data_X=-（Column-1）/2:（Column-1）/2;
end
figure（3）;
plot（Data_XImage_array）;
xlabel（‘離散頻率‘）
ylabel（‘傅里葉頻率歸一化后取自然對數‘）;
if?rem（Column2）==0
????Data_Y=Image_array（Column/2+1:end）;
else
????Data_Y=Image_array（（Column+1）/2:end）;
end
figure（4）;
semilogx（Data_Y‘b:*‘）;
xlabel（‘離散頻率的自然對數‘）;
ylabel（‘傅里葉頻率歸一化后取自然對數‘）;
hold?on;
if?rem（Column2）==0
????k=1:Column/2;
else
????k=1:（Column+1）/2;
end
p=polyfit（log10（k）Da

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????3253854??2012-04-19?23:14??仿真代碼\1.bmp

?????文件?????265231??2012-04-24?16:05??仿真代碼\10.bmp

?????文件?????814356??2012-05-02?22:52??仿真代碼\11.bmp

?????文件?????786488??2012-05-02?22:53??仿真代碼\12.bmp

?????文件?????559160??2012-05-02?22:53??仿真代碼\13.bmp

?????文件?????536856??2012-05-02?22:54??仿真代碼\14.bmp

?????文件????3253856??2012-05-02?22:56??仿真代碼\15.bmp

?????文件????2233400??2012-05-02?22:56??仿真代碼\16.bmp

?????文件????2143352??2012-05-02?22:56??仿真代碼\17.bmp

?????文件????3145784??2012-05-02?22:56??仿真代碼\18.bmp

?????文件????1189246??2011-06-30?23:41??仿真代碼\19.bmp

?????文件????2307126??2012-04-19?23:17??仿真代碼\2.bmp

?????文件????1244214??2011-06-30?20:05??仿真代碼\20.bmp

?????文件????1118626??2012-04-19?23:24??仿真代碼\3.bmp

?????文件????2098230??2012-04-24?15:26??仿真代碼\4.bmp

?????文件????2098230??2012-04-24?15:58??仿真代碼\5.bmp

?????文件????3145782??2012-04-24?15:59??仿真代碼\6.bmp

?????文件????2233398??2012-04-24?15:59??仿真代碼\7.bmp

?????文件????2359350??2012-04-24?16:03??仿真代碼\8.bmp

?????文件????2143350??2012-04-24?16:03??仿真代碼\9.bmp

?????文件???????8368??2012-05-31?17:36??仿真代碼\blur_deblur.m

?????文件???????1519??2012-05-03?01:45??仿真代碼\clear_image.m

?????文件???????7577??2012-05-31?17:27??仿真代碼\deblur.asv

?????文件???????7401??2012-05-31?17:34??仿真代碼\deblur.m

?????文件????3478377??2012-03-31?16:37??仿真代碼\論文.pdf

?????目錄??????????0??2013-04-01?16:43??仿真代碼

-----------?---------??----------?-----??----

?????????????40433131????????????????????26