LMS自适应滤波器MATLAB实现

%filter tap:k
%data length:N

g=100;  %统计仿真次数为g
N=1024;  %输入信号抽样点数N
k=128;   %时域抽头LMS算法滤波器阶数
pp=zeros（gN-k）;  %将每次独立循环的误差结果存于矩阵pp中，以便后面对其平均
u=0.0001;

for q=1:g
    t=1:N;
    a=1;
    s=a*sin（0.05*pi*t）; %input signal
    figure（1）;
    subplot（311）
    plot（treal（s））;  %信号s的时域波形
    title（‘信号s时域波形‘）;
    xlabel（‘n‘）;
    ylabel（‘s‘）;
    axis（[0N-a-1a+1]）;
    xn=awgn（s5）     %加入零均值高斯白噪声，信噪比为3dB
    
    %设置初值
    y=zeros（1N）;   %output signal
    y（1:k）=xn（1:k）;  %将输入信号xn的前k个值作为输出y的前k个值