自适应陷波器
链接:
1、自适应陷波滤波器
自适应算法分两种,一种是对本身信号中的某个信息感兴趣,一个是对信号的预期感兴趣(即降噪)
一、信号波形
二、频谱分析
三、代码实现
clear all;close all;
N = 128;
T = 1;
t = linspace(0,T,N);
x = 12*sin(2*pi*10*t+pi/4)+5*cos(2*pi*40*t);
figure,plot(x)
dt = t(2)-t(1);
f = 1/dt;//采样频率
X= fft(x);
F = X(1:N/2+1);
f = f*(0:N/2)/N;
figure,plot(f,abs(F),'-');
xlabel('Frequency');
ylabel('|F(k)|');
窄带信号产生wav文件
fs = 55000; %采样频率,决定fft的横坐,采样点数增加信号的保真率更高
fc = 1100; %载波1.1kHz
T = 1/fc; %载波周期
t = 0:1/fs:11000*T;%11000/1100=10s钟wav信号
YY = cos(2*pi*fc.*t);
audiowrite(sprintf('cos1100.wav'),YY,fs);