改进的维纳滤波的matlab代码

% wien.m
%Image restoration by Wiener filter
% call from: x=wien（yh）
%y--observer image;h--blur;x--restored image
%参考文献：反卷积和信号复原[M]6.2，6.4
function x=wien（yh）
[mn]=size（y）;
y=double（y）;
[mh nh]=size（h）;
m2=2*m;
n2=2*n;
yy=zeros（m2n2）;
yy（1:m1:n）=y;
yy（1:mn+1:n2）=y（1:mn:-1:1）;
yy（m+1:m21:n2）=yy（m:-1:11:n2）;
Y=fft2（yy）;
he=zeros（m2n2）;
he（1:mh1:nh）=h/sum（sum（h