function?image=bldconv（g）


a1=0.1;
a2=0.01;

PQ=paddedsize（size（g））;
G=fft2（gPQ（1）PQ（2））;

[y?x]=size（G）;
htemp=ones（3）;
h0=freqz2（htempPQ（1）PQ（2））;
R=Rcreat（yx）;

H=h0;
for?k=1:10?%計算IMG和psf
iMG=（conj（H）.*G）./（（conj（H）.*H）+a1*R）;
H=conj（iMG）.*G./（（conj（iMG）.*iMG）+a2*R）;
end

IMG=mat2gray（real（ifft2（iMG）））;
image=IMG（1:size（g1）1:size（g2））;
imwrite（image‘testfile.tif‘）;
imshow（image[]）;
hold?on;

end

%===========================
%計算R矩陣的函數Rcreat
%===========================

function?R=Rcreat（yx）
%R矩陣生成子函數
%by?Realasking
%為bdeconv.m編制

i=1:y;
j=1:x;

RI=zeros（yx）;
RJ=zeros（yx）;
R=zeros（yx）;

for?k=1:y?%向量化代碼生成R的矩陣
RI（ki）=-2*cos（2*pi*i./y）;
end

for?k=1:x
RJ（jk）=-2*cos（2*pi*j‘./x）;
end

R=4+RI+RJ;

end

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

?????文件????????816??2010-05-13?11:28??新建文件夾?（2）\H1模型.m

?????文件???????3144??2010-05-27?11:05??新建文件夾?（2）\PSNR.txt

?????文件???????1129??2010-05-13?16:08??新建文件夾?（2）\TV去噪.m

?????文件???????6252??2010-05-19?15:09??新建文件夾?（2）\weixing.jpg

?????文件???????1564??2010-05-11?15:20??新建文件夾?（2）\改進的TV模型.m

?????文件???????1185??2010-05-05?15:06??新建文件夾?（2）\自適應TV去噪模型.m

?????目錄??????????0??2010-05-31?12:00??新建文件夾?（2）

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

????????????????14090????????????????????7