LMS算法 matlab编程

用matlab 编写的LMS算法。其中包含了滤波器的建立及仿真，以及滤波前后的比较

function [ynWen]=LMS（xndnMmuitr）
% LMS（Least Mean Squre）算法
% 输入参数:
%     xn   输入的信号序列      （列向量）
%     dn   所期望的响应序列    （列向量）
%     M    滤波器的阶数        （标量）
%     mu   收敛因子（步长）      （标量）     要求大于0小于xn的相关矩阵最大特征值的倒数    
%     itr  迭代次数            （标量）     默认为xn的长度M% 输出参数:
%     W    滤波器的权值矩阵     （矩阵）
%          大小为M x itr
%     en   误差序列（itr x 1）    （列向量）  
%     yn   实际输出序列             （列向量）

% 参数个数必须为4个或5个
if nargin == 4                 % 4个时递归迭代的次数为xn的长度 
    itr = length（xn）;
elseif nargin == 5             % 5个时满足M    if itr>length（xn） | itr        error（‘迭代次数过大或过小!‘）;
    end
else
    error（‘请检查输入参数的个数!‘）;
end


% 初始化参数
en = zeros（itr1）;             % 误差序列en（k）表示第k次迭代时预期输出与实际输入的误差
W  = zeros（Mitr）;             % 每一行代表一个加权参量每一列代表-次迭代初始为0

% 迭代计算
for k = M:itr                  % 第k次迭代
    x = xn（k:-1:k-M+1）;        % 滤波器M个抽头的输入
    y = W（:k-1）.‘ * x;        % 滤波器的输出
    en（k） = dn（k） - y ;        % 第k次迭代的误差
    
    % 滤波器权值计算的迭代式
    W（:k） = W（:k-1） + 2*mu*en（k）*x;
end

% 求最优时滤波器的输出序列
yn = inf * ones（size（xn））;
for k = M:length（xn）
    x = xn（k:-1:k-M+1）;
    yn（k） = W（:end）.‘* x;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1457  2008-06-02 16:58  lms算法\LMS.m

     文件       1091  2008-06-03 10:03  lms算法\main.m

     文件        766  2008-06-03 09:28  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\1.jpg

     文件       1372  2008-06-03 09:28  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\2.jpg

     文件       1386  2008-06-03 09:28  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\3.jpg

     文件       1420  2008-06-03 09:28  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\4.jpg

     文件       1417  2008-06-03 09:28  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\5.jpg

     文件      42849  2008-06-03 09:38  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\common.js

     文件       1840  2008-06-03 09:38  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\fgetccode_btn.gif

     文件       1453  2008-06-03 09:38  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\fronthelper.pack.js

     文件      21070  2008-06-03 09:38  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\inf.js

     文件       1890  2008-06-03 09:39  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\lightbox.css

     文件       6385  2008-06-03 09:38  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\lightbox.pack.js

     文件      27051  2008-06-03 09:39  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\mozine.css

     文件        415  2008-06-03 09:28  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\rss.png

     文件      22213  2008-05-28 15:59  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\show_ads.js

     文件       6299  2008-06-03 09:38  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\T1z_haXoVbSIyCtQg6_062423.jpg_sum.jpg

     文件       1975  2008-06-03 09:38  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\widget.js

     文件       3528  2008-06-03 09:38  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\wz.js

     文件          0  2008-06-03 09:39  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files\ystat.js

     文件      30107  2008-06-03 09:39  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.htm

     目录          0  2009-03-02 16:14  lms算法\rand函数\Matlab rand 随机函数 - 开发经营管理 - 陌路阑珊.files

     目录          0  2009-03-02 16:14  lms算法\rand函数

     目录          0  2009-05-13 22:41  lms算法

----------- ---------  ---------- -----  ----

               175984                    24