wienerdec.m

斯坦福大学的小波的工具箱 比matlab自带的工具箱更完善

function [H,out] = Wienerdec(g,in,sig,sigma)
% wienerdec -- Wiener filter for deconvolution.
% Usage
%   [H,out] = wienerdec(g,in,sigma)
% Inputs
%   g		Convolution filter.
%   in		Degraded signal.
%   sig		Original signal.
%   sigma		S.d. of the additive GWN.
% Outputs
%   H		Wiener filter
%   out		Restored signal.
%

Y=fft(sig);
FN=fft(in);
n=length(in);
PY=abs(Y.*Y);
G=fft(g);
Gsq=G.*conj(G);
H=(PY.*conj(G))./(PY.*Gsq + n*sigma^2);
FN=FN.*H;
out=real(ifft(FN));


% Written by Maureen Clerc and Jerome Kalifa, 1997
% [e-mail], [e-mail]






    
    
%   
% Part of WaveLab Version 802
% Built Sunday, October 3, 1999 8:52:27 AM
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%