wh.m

维纳滤波器 通过求解Wiwner-Hopf方程设计滤波器

function [h,e]=WH(Rss,Rww,M)
                                       %编制最优化维纳霍夫滤波器，其中M为信号的点长,h为冲激响应，e为误差，Rss,Rww分别是真实信 
                                       %号和噪声的自相关函数
e1=10;e0=0;N=0;                        %给e1,e0,N赋初值
while abs(e0-e1)>1e-6                  %e0,e1不很接近则执行循环体
    N=N+1;
    e0=e1;
    Rxs=Rss(M:M+N-1);                  %s(n)与v(n)不相关，则Rxs=Rss
    Rxx=Rww(M:M+N-1)+Rss(M:M+N-1);
    R_xx=zeros(N);                     %R_xx*H=Rxs
    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                                  %循环结束并显示N的最终结果 
e=e1;