直流分量去除方法
- 直流分量即信号的均值,因此将原始数据直接减去其均值可消去。
y_data=y1-mean(y1); % 原始信号减去均值得到去直流分量后的信号
subplot(121);plot(y_data)
y=fft(y_data,N);
mag=abs(y)*2/N;
subplot(1,2,2);plot(f(1:N/2),mag(1:N/2));
- 利用频谱消除:频谱中第一个数据对应的是0HZ,也即直流分量,将0HZ对应的数据置为0即可消去直流分量。(令f(1)=0)
y_new=y;
y_new(1)=0;
x_new=real(ifft(y_new)); % 求得去除直流分量后的原始信号
subplot(121);plot(x_new)
y=fft(x_new,N);
mag=abs(y)*2/N;
subplot(122);plot(f(1:N/2),mag(1:N/2));