[chenzheyi]

我自己失败的程序：
N=13;
[fx fy]=gradient(f);
[fxx fxy]=gradient(fx);
[fyx fyy]=gradient(fy);
left=fft2(fxx+fyy);%计算等式左边的值。

for u=0:N-1
for v=0:N-1
h(u+1,v+1)=-1*(u^2+v^2);
end
end %计算求得滤波方程。

fz2=left./h;
fz2(1,1)=0;
f2=ifft2(fz2); %求得重建的f'。
不知道为什么反正总是不对……