跪求一个H无穷滤波matlab程序代码！万分感谢！

用一个简单的代码，为了AR模型Y（K）=写LMS自适应滤波 - AY（K-1）+ X（k）的

对于k = 1：M

X = 0.2 * randn（1， N）

Y（1）= 1％初始

对于i = 2：N

Y（I）= - A * Y（I-1）+ X（I）％按AR模型生成信号

结束

对于i = 2：N

E（K，I）= Y（I）-H（K，I）* Y（I-1）％误差

H（K ，I + 1）= H（K，I）+增量* Y（I-1）* E（K，I）％一估计值

端

端

% H 为已知

% Img 为已知

Img_FFT = fft2(Img)

Img_FFT_Shift = fftshft(Img_FFT)

Img_Fre = Img_FFT_Shift .* H % 频域的结果

Img_C = ifft2( ifftshift(Img_Fre) )% Complex

Img_Low = ImgC.abs()

---