LMS自适应滤波器(共5页)

1、精选优质文档-倾情为你奉上LMS自适应滤波器是使滤波器的输出信号与期望响应之间的误差的均方值为最小，因此称为最小均方（LMS）自适应滤波器。function yn,W,en=LMS(xn,dn,M,mu,itr)% LMS(Least Mean Squre)算法% 输入参数:%     xn   输入的信号序列      (列向量)%     dn 

2、;  所期望的响应序列    (列向量)%     M    滤波器的阶数        (标量)%     mu   收敛因子(步长)      (标量)     要求大于0,小

3、于xn的相关矩阵最大特征值的倒数    %     itr  迭代次数            (标量)     默认为xn的长度,M<itr<length(xn)% 输出参数:%     W    滤波器的权

4、值矩阵     (矩阵)%          大小为M x itr,%     en   误差序列(itr x 1)    (列向量)  %     yn   实际输出序列 

5、60;           (列向量)% 参数个数必须为4个或5个if nargin = 4                 % 4个时递归迭代的次数为xn的长度     itr = length(x

6、n);elseif nargin = 5             % 5个时满足M<itr<length(xn)    if itr>length(xn) | itr<M        error('迭代次数过大或过小!'); &

7、#160;  endelse    error('请检查输入参数的个数!');end% 初始化参数en = zeros(itr,1);             % 误差序列,en(k)表示第k次迭代时预期输出与实际输入的误差W  = zeros(M,itr);     &#

8、160;       % 每一行代表一个加权参量,每一列代表-次迭代,初始为0% 迭代计算for k = M:itr                  % 第k次迭代    x = xn(k:-1:k-M+1); 

9、60;      % 滤波器M个抽头的输入    y = W(:,k-1).' * x;        % 滤波器的输出    en(k) = dn(k) - y         %

10、 第k次迭代的误差        % 滤波器权值计算的迭代式    W(:,k) = W(:,k-1) + 2*mu*en(k)*x;end% 求最优时滤波器的输出序列yn = inf * ones(size(xn);for k = M:length(xn)    x = xn(

11、k:-1:k-M+1);    yn(k) = W(:,end).'* x;end            LMS函数的一个实例：%function main()close  all% 周期信号的产生 t=0:99;xs=10*sin(0.5*t);figure;subplot(2,1,1);plot(t,xs);grid;ylabel('幅值'

12、);title('it输入周期性信号');% 噪声信号的产生randn('state',sum(100*clock);xn=randn(1,100);subplot(2,1,2);plot(t,xn);grid;ylabel('幅值');xlabel('时间');title('it随机噪声信号');% 信号滤波xn = xs+xn;xn = xn.'    % 输入信号序列dn =

13、0;xs.'    % 预期结果序列M  = 20      % 滤波器的阶数rho_max = max(eig(xn*xn.');   % 输入信号相关矩阵的最大特征值mu = rand()*(1/rho_max)       % 收敛因子 0 <

14、; mu < 1/rhoyn,W,en = LMS(xn,dn,M,mu);% 绘制滤波器输入信号figure;subplot(2,1,1);plot(t,xn);grid;ylabel('幅值');xlabel('时间');title('it滤波器输入信号');% 绘制自适应滤波器输出信号subplot(2,1,2);plot(t,yn);grid;ylabel('幅值');xlabel('时间');title('it自适应滤波器输出信号');% 绘制自适应滤波器输出信号,预期输出信号和两者的误差figure plot(t,yn,'b&#