www.pudn.com > speechenhancement.rar > WH.m

function [h,e]=WH(Rss,Rww,M) 
%Compliment Wiener_Hoff Optimal filter 
%Rss denotes auto-correlation of signal and Rww denotes auto-correlation 
%of White noise generated by computer 
%M denotes the point-number of signal 
e1=10;e0=0;N=0;%Value e1,e0 and convice to execute 
while abs(e0-e1)>1e-6 
    N=N+1; 
    e0=e1; 
    Rxs=Rss(M:M+N-1); 
    Rxx=Rww(M:M+N-1)+Rss(M:M+N-1); 
    R_xx=zeros(N); 
    for j=1:N 
        for n=1:N 
            R_xx(j,n)=Rxx(abs(j-n)+1); 
        end 
    end 
    h=inv(R_xx)*Rxs'; 
    e1=Rss(M)-h'*Rxs'; 
end,N 
e=e1;