将目标算法用MATLAB实现（本人新手，请多多添加标注，以免我看不懂）

参考文献：万方数据 《车内噪声FXFU-LMS声振混合主动控制算法》

（FXLMS）

u ( k ) = x ( k )*W ( k ) “*”代表线性卷积。

 y ( k ) = u ( k )*s( k )

e( k ) = d ( k ) - y ( k )

滤 波 器 的 权 系 数 W ( k ) = [ W0 ( k ),W1 ( k ),⋯,WL - 1 ( k ) ]自适应更新为 W ( k + 1 ) = W ( k ) + μ∙xs ( k )*e( k )

xs( k ) = x ( k )*ŝ( k )

•其中：ŝ( k )是Ŝ ( z )的脉冲响应，μ为迭代步长，使得 滤波器系数W ( k )向目标函数的梯度反方向更新迭 代一步。μ的选择将会影响算法的收敛速度和最小 均方误差（SSE）

（FULMS） 

B( k ) = SS ( k )*pp ( k )

Y ( k ) = X͂( k ) + Y͂( k )

X͂( k ) = X ( k )*WW ( k )

 Y͂( k ) = Y ( k )*D ( k )

 WW ( k + 1 ) = WW ( k ) + 2δE ( k )*R( k )

 D( k + 1 ) = D( k ) + 2αE ( k )*G ( k )

R( k ) = X ( k )*ĥ( k )

 G ( k ) = Y ( k )*ĥ( k )

Z ( k ) = ( WW ( k )*X ( k ) + D( k )Y ( k ) )*h ( k ) = WW ( k )*R( k ) + D( k )*G ( k )

  E ( k ) = B( k ) - Z ( k ) = B( k ) - WW ( k )*R( k ) - D( k )*G ( k )

X ( k ) = [ x ( k ), x ( k - 1 ),⋯, x ( k - p ),⋯, x ( k - P + 1 ) ]

 WW ( k ) = [ ww0 ( k ),ww1 ( k ),⋯,wwp ( k ),⋯, wwP - 1 ( k ) ]

D( k ) = [ d0 ( k ), d1 ( k ),⋯, dq ( k ),⋯, dQ - 1 ( k ) ]

E_fc( k ) = IFFT ( FFT ( E ( k ) )*ZSFRF )

X_fc( k ) = IFFT ( FFT ( X ( k ) )*ZSFRF )  

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，目前超出我们的服务范围，暂时无法为您解答。

首次提问人员可免费体验一次有问必答服务。目前首次提问的问题服务范围为：编程语言、Java开发、python、数据库、前端开发 领域专业技术问题，为您提供问题的解决思路和指导。不提供源码代写、项目文档代写、论文代写、安装包资源发送或安装、软件使用指导等服务。

我们后续会持续优化，扩大我们的服务范围，为您带来更好地服务。