%?wien.m
%Image?restoration?by?Wiener?filter
%?call?from:?x=wien（yh）
%y--observer?image;h--blur;x--restored?image
%參考文獻(xiàn)：反卷積和信號(hào)復(fù)原[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