function??object?=?TV（blurimage）;

[N1N2]=size（blurimage）;

%%中心化算子
T=zeros（N1N2）;
x0=N1/2;y0=N2/2;
for?re=1:N1
?????for?im=1:N2
??????????ta=cos（??2*pi*（?（re-1）*x0+（im-1）*y0??）/N1?）?;
??????????tb=-sin（2*pi*（?（re-1）*x0+（im-1）*y0?）/N2?）?;
????????T（reim）=complex（tatb）;
?????end
end?

%一階差分算子
aa=0.66499;?????%?d0=[-aa?bb?0;
????????????????????????%?????????bb??0??0;
bb=0.17555;?????%?????????0???0??0];

%?********數據準備***********%
u=blurimage;
h=ones（N1N2）;
h=h/sum（sum（h））;

Z=fft2（blurimage）;

%*******************主要迭代過程*******************%
%?%規整化參數
lamda_u=2.0e-5;
lamda_h=1.0e-5;
%迭代控制參數
AM_Iteration=30;
CG=12;
%***AM迭代**
U=fft2（u）;

for?ll=1:AM_Iteration
????
????H=fft2（h）;
????HH=conj（H）.*H.*T;
????S=Z-H.*U.*T;
????B=conj（H）.*S.*T;
????b=re

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2012-11-18?19:41??TV\
?????文件????????3455??2011-07-11?16:42??TV\TV.m
?????文件????????1918??2011-11-02?09:40??TV\TV1.asv
?????文件????????2048??2011-11-08?15:52??TV\TV1.m