matlab实现nlms算法

以下是一个简单的MATLAB代码实现NLMS算法：

function [y, h] = nlms(x, d, mu, M)
% x: 输入信号
% d: 目标输出信号
% mu: 步长
% M: 滤波器阶数

% 初始化权值向量
h = zeros(1, M);

% 初始化输出信号向量
y = zeros(size(d));

% 对输入信号进行滤波
for i = 1:length(d)
% 提取当前输入信号的子序列
x_sub = x(max(1, i-M+1):i);

% 计算当前输出信号
y(i) = h * x_sub';

% 计算误差信号
e = d(i) - y(i);

% 更新权值向量
h = h + mu * e * x_sub / (x_sub * x_sub');
end

end

在这个代码中，我们首先初始化了权值向量h和输出信号向量y。然后，我们对输入信号进行滤波，对于每个输入信号样本，我们提取当前输入信号的子序列，计算当前输出信号，计算误差信号，然后使用NLMS算法更新权值向量。最终，我们返回输出信号向量和更新后的权值向量。