自适应滤波器matlab仿真

是关于自适应滤波器的列子，matlab实现。可以运行。适合初学者

%双频陷波器
%%% 改变第一级滤波器参考输入的幅值由1→2 %%%
clear all
clc
 
N=2000;                     %总采样长度
t=linspace（02N）; 
s=10*sin（4*pi*t）;           %原始正弦信号                   
n1=3*sin（100*pi*t）;         %50Hz干扰信号
n2=5*sin（50*pi*t）;          %25Hz干扰信号
x=s+n2+n1;                  %引入正弦双频干扰的原始输入
x1=2*cos（100*pi*t）;         %参考输入1
x2=2*sin（100*pi*t）;         %参考输入2
w1=0.1;
w2=0.1;                     %权矢量初值
e1=zeros（1N）;
y1=0;
u=0.01;                  
for i=1:N                   %LMS算法
    y1=w1*x1（i）+w2*x2（i）;
    e1（i）=x（i）-y1; 
    w1=w1+u*e1（i）*x1（i）; 
    w2=w2+u*e1（i）*x2（i）; 
end

x1=cos（50*pi*t）;
x2=sin（50*pi*t）;
w1=0.1;
w2=0.1;  
e2=zeros（1N）;
y2=0;
for i=1:N                   %LMS算法
    y2=w1*x1（i）+w2*x2（i）;
    e2（i）=e1（i）-y2; 
    w1=w1+u*e2（i）*x1（i）; 
    w2=w2+u*e2（i）*x2（i）; 
end

figure（1） 
subplot（211）;
plot（ts）;
title（‘原始正弦信号‘）;
axis（[0 2 -20 20]）;
subplot（212）; 
plot（tx）;
title（‘加入双频干扰的原始信号‘）;

figure（2）
subplot（311）;
plot（ts+n2）;
title（‘1级自适应陷波器输出理论值‘）;
subplot（312）;
plot（te1）;
title（‘1级自适应陷波器输出实验值‘）;
subplot（313）;
plot（ts+n2-e1）;
title（‘1级信号误差‘）;


figure（3）
subplot（311）;
plot（ts）;
title（‘2级自适应陷波器输出理论值‘）;
axis（[0 2 -20 20]）;
subplot（312）;
plot（te2）;
title（‘2级自适应陷波器输出实验值‘）;
subplot（313）;
axis（[0 2 -1 1]）;
plot（ts-e2）;
title（‘2级信号误差‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1022  2010-01-23 17:30  single frequency--different u.m

     文件       1499  2010-01-23 17:32  double frequency--change first noise amplitude.m

     文件       1067  2010-01-23 17:31  single frequency--different △ f.m

     文件       1013  2010-01-25 10:23  single frequency--different amplitude.m

     文件       1013  2010-01-25 10:55  single frequency--different Fudu.m

     文件       1022  2010-01-23 17:31  single frequency--different theta.m

----------- ---------  ---------- -----  ----

                 6636                    6