自适应噪声抵消

weixin_43468843

已于 2022-06-01 15:49:49 修改

阅读量492

点赞数

文章标签： matlab 算法开发语言

于 2022-06-01 15:29:22 首次发布

本文链接：https://blog.csdn.net/weixin_43468843/article/details/125084145

版权

自适应噪声抵消

声源：CW信号 15kHz
干扰：宽带、连续信号 3k~30kHz

信号和干扰设置

clear all
close all
clc
%% 信号
f0=15e3;
T=1.0000e-02;%时宽
B=1/T;
fs=1e5;%采样频率
Ts=1/fs;%采样时间
N=round(T/Ts);%采样点个数
n=0:N-1;
t=linspace(0,T,N);
xs=exp(1i*2*pi*f0*t);%cw信号
noise=wgn(1,N,0);% 噪声信号
wp1=[3000,30000]*2/fs; %通带边界频率 ，单位为rad/s
ws1=[1000,35000]*2/fs; %阻带边界频率 ，单位为rad/s
Rp1=3;           %通带最大波纹度 ，单位dB (不要太小)
Rs1=10;          %表示阻带最小衰减，单位dB
[Nn1,Wn1]=buttord(wp1,ws1,Rp1,Rs1);%巴特沃斯数字滤波器的阶数n和-3dB归一化截止频率Wn
[B1,A1]=butter(Nn1,Wn1);%得到n阶巴特沃斯滤波的分子分母
vn=filter(B1,A1,noise);                    
ff=n*fs/N;
y=fft(vn,N);    %对信号进行快速Fourier变换
mag=abs(y);     %求得Fourier变换后的振幅
% % %-----------------------画图------------------------
figure()
subplot(221)
plot(t,real(xs));
% xlim([0,0.005]);
grid on;hold on;
title('CW脉冲信号时域波形');
xlabel('时间/s');ylabel('幅度');
subplot(222)
f=linspace(0,2*f0,N);
Y=T*(sin(pi*(f-f0)*T))./(pi*(f-f0)*T);
plot(f/1000,Y);grid on;hold on;
title('CW脉冲信号频谱');
xlabel('频率/kHz');ylabel('幅度');
subplot(223)
plot(t,real(vn));
grid on;hold on;
title('干扰信号时域波形');
xlabel('时间/s');ylabel('幅度');
subplot(224)
plot(ff(1:N/2)/1000,mag(1:N/2));%绘出Nyquist频率之前随频率变化的振幅
grid on;hold on;
title('干扰信号频谱');
xlabel('频率/kHz');ylabel('幅度');