RLS最小二乘法自适应滤波器m文件

Fs = 40000;                                                     %设置采样频率
t = 0:1/Fs:3.5;  
t = t‘;
Size_t = size（t1）;
F1 = 20000;
F2 = 0;
F3 = 0;
F4 = 1000;
Signal = sin（2*pi*F1*t） + 0.5*sin（2*pi*F2*t） + 0.25*sin（2*pi*F3*t）; %生成信号
noise_amp = 1;                                           %定义噪声的标准差
noise1 = noise_amp*randn（Size_t1）;                      %生成高斯白噪声
noise2 = noise_amp*randn（Size_t1）;
noise3 = 5*sin（2*pi*F4*t+pi/2）;

noise = noise2;
Signal_noise = Signal + 0.2*noise;                           %加入高斯白噪声
%Signal_noise（2:end） = Signal_noise（2:end） + 0.15*noise（1:end-1）;
%Signal_noise（3:end） = Signal_noise（3:end） + 0.1*noise（1:end-2）;

subplot（211）;
plot（tSignal）;
title（‘原始信号‘）;
subplot（212）;
plot（tSignal_noise）;
title（‘加入干扰噪声的信号‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%计算rls算法
M = 3;                       %定义FIR滤波器阶数
lamda = 1;                %定义遗忘因子
Signal_Len = Size_t - M -1;   %定义信号数据的个数
I = eye（M）;                   %生成对应的单位矩阵
c = 1;