task2.m

来自「自己编的LMS自适应算法的MATLAB源程序 本程序内容： 产生一个信号+白」· M 代码 · 共 28 行

% Adaptive filter 
close all; clear all; clc
f=[200;300];A=[1,2];fs=1000;t=(0:fs)/fs;
s=A*sin(2*pi*f*t);
u=wgn(1,length(t),0);
m=length(t);
subplot(311);
plot(t,s);
title('窄带周期信号2sin(2*pi*200*t)+sin(2*pi*300*)');
subplot(312);
plot(t,u);
title('宽带高斯白噪声');
x=s+u;
d=zeros(1,m);
d=x;

delay=3;
xshift=zeros(1,m);
for i=delay+1:m
    xshift(i)=x(i-delay);
end
    
[h,y,e] = lmsx(xshift,d,0.001,50);
subplot(313);
plot(t,y);
title('经过滤波的窄带周期信号');