请问自适应陷波程序
请问这个自适应滤波程序是lms算法还是rls算法？我需要原理部分的说明 谢谢！N=2000; %总采样长度
t=0:N-1; %时间的变化范围
x=load('ecg.txt'); %输入信号
subplot(211); %作第1子图
plot(t,x(1:N));% x is the original signal mixed with 50Hz noise.
title('original'); grid;
xx=x(t+1)';
for k=1 :9
x1=cos(2*pi*t*50*k/1000);% 50Hz noise 采样频率1000hz
x2=sin(2*pi*t*50*k/1000);% 50Hz noise
% 50Hz notch filt starting
%初始化
w1=0;
w2=1;
e=zeros(1,N);
y=0.1;
u=0.1;%迭代步长
for i=1:N
y=w1*x1(i)+w2*x2(i);
e(i)=xx(i)-y;%误差信号 output
w1=w1+u*e(i)*x1(i);%迭代方程
w2=w2+u*e(i)*x2(i);%迭代方程
end
xx=e;
end
subplot(212); %作第2子图
plot(t,e); grid;
title('output')
页:
[1]