matlab 盲卷积算法,盲去卷积 图像复原 matlab代码解释

下面是一个人写的盲去卷积算法实现图像去模糊，其中有些代码不是很能理解，求解释！尤其每次都做这个循环：

for i=1:x_size         %进行乘(-1)^(i+j)处理

for j=1:y_size

flag=(-1)^(i+j);

tmp1(i,j)=tmp1(i,j)*flag;

end

end

到底是干嘛？最好能对整段代码加些注释再贴到回复里面。谢谢！

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%程序内容：盲迭代反卷积方法恢复图象

%%编程时间：2007年12月16日

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all;

gray_S=imread('cameraman.jpg');    %读入图像

figure;imshow(gray_S);title('original image')   %显示原始图

[x_size,y_size]=size(gray_S);

%求取PSF尺寸大小(圆形)

p_size=5;                       %定义psf尺寸大小

p_support=zeros(x_size,y_size); %psf掩模

p_support(1:p_size,1:p_size)=1;

beta=0.9;           %参数设置

g=im2double(gr