|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
clear all;close all;clc;
N=512;%采样点数
fs=50; %系统采样频率 Ts=125e-6;
t=(0:N-1)/fs;
f=(0:N-1)*fs/N;
xs=exp(j*pi*(0.5*t.^2));
e=randn(size(t));
sn=xs+e;
% %----------LMS-----------------------------------------
xn=sn';%%自适应滤波起的参考输入,即自适应滤波器的输入
dn=xs; %期望信号(对的不??)
M=5;%%滤波器介数
mu=0.00001;
N=512;%%信号长度
en=zeros(N,1);%%自适应滤波系统输出
W=zeros(M,1); %LMS滤波器系数
%------------迭代算法-------------------
for k=M:N
x=xn(k:-1:k-M+1);
yn(k)=W'*x;
en(k)=dn(k)-yn(k);
out(k)=en(k);
W=W+mu*en(k)*x;
k=k+1;
end
figure(1);
subplot(211);
plot(t,sn);title('原始带噪信号');grid;
subplot(212);
% hold on;
% plot(t,xs','r');title('qiwang');grid;
plot(t,out);title('去噪后信号');grid; |
|