PDE 偏微分方程图像复原代码

偏微分方程图像复原代码，其中有频域迭代盲解卷积H1模型，TV去噪模型，自己编的自适应去噪模型，以及改进的自适应去噪模型，令附卫星图片，PSNR代码

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