你能加入平滑因子写一个NLMS函数么

当然可以。以下是一个简单的NLMS函数，其中包含平滑因子参数mu和防止零除错误的epsilon参数。

import numpy as np

def nlms(x, d, M, mu, epsilon):
    """
    :param x: 输入信号
    :param d: 期望输出信号
    :param M: 滤波器长度
    :param mu: 步长因子
    :param epsilon: 防止零除错误的小常数
    :return: 滤波器系数向量
    """

    # 初始化滤波器系数向量
    w = np.zeros(M)

    # 初始化误差向量
    e = np.zeros_like(d)

    # 进行自适应滤波
    for i in range(M, len(x)):
        # 从输入信号中提取当前时刻的M个样本，作为当前子信号
        xi = x[i-M:i]

        # 计算当前子信号与滤波器系数向量的卷积，得到当前估计输出
        yi = np.dot(w, xi)

        # 计算当前时刻的误差
        e[i] = d[i] - yi

        # 更新滤波器系数向量
        w += mu * xi * e[i] / (np.dot(xi, xi) + epsilon)

    return w

在上面的代码中，输入信号x和期望输出信号d都是一维的numpy数组，滤波器长度M是一个整数，步长因子mu和防止零除错误的小常数epsilon都是浮点数。函数返回一个numpy数组，包含估计的滤波器系数。

需要注意的是，NLMS算法是一种自适应滤波算法，能够自动调整滤波器系数以适应不同的输入信号。但是，这种算法对于输入信号中存在的高频噪声比较敏感，可能会导致滤波器系数发散，因此在实际应用中需要进行一定的预处理和调试。