从信号中去除 60Hz 杂声
美国和其他几个国家/地区的交流电以 60 Hz 的频率振荡。这些振荡通常会破坏测量结果,必须将其减去。
% 在存在 60 Hz 电力线噪声的情况下,研究模拟仪器的输入的开环电压。
% 电压采样频率为 1 kHz。
load openloop60hertz, openLoop = openLoopVoltage;
Fs = 1000; % 采样频率
t = (0:length(openLoop)-1)/Fs;
plot(t,openLoop)
ylabel('Voltage(V)')
xlabel('Time(s)')
title('Open-Loop Voltage with 60 Hz Noise')
grid
使用 Butterworth 陷波滤波器消除 60 Hz 噪声。使用 designfilt
设计该滤波器。陷波的宽度定义为 59 至 61 Hz 的频率区间。滤波器至少去除该范围内频率分量的一半功率。
% % % 设计一个2阶带阻IIR滤波器,低3db频率59hz,高3db频率61hz,
% 采样率为Fs=1000hz。可视化滤波器的幅值响应。
% 用它来过滤openloop60hertz数据样本。
d = designfilt('bandstopiir','FilterOrder',2, ...
'HalfPowerFrequency1',59,'HalfPowerFrequency2',61, ...
'DesignMethod','butter','SampleRate',Fs);%设计Butterworth 陷波滤波器
% 绘制滤波器的频率响应。请注意,此陷波滤波器提供高达 45 dB 的衰减。
fvtool(d,'Fs',Fs)
% 用 filtfilt 对信号进行滤波,以补偿滤波器延迟。注意振荡是如何显著减少的。
figure
buttLoop = filtfilt(d, openLoop);%使用数字滤波器d对输入数据x进行零相位滤波。
plot(t,openLoop,t,buttLoop)
ylabel('Voltage (V)')
xlabel('Time (s)')
title('Open-Loop Voltage')
legend('Unfiltered','Filtered')
grid
% 使用周期图可以看到 60 Hz 的“峰值”已去除。
figure
[popen,fopen] = periodogram(openLoop,[],[],Fs);
[pbutt,fbutt] = periodogram(buttLoop,[],[],Fs);
plot(fopen,20*log10(abs(popen)),fbutt,20*log10(abs(pbutt)),'--')
ylabel('Power/frequency (dB/Hz)')
xlabel('Frequency (Hz)')
title('Power Spectrum')
legend('Unfiltered','Filtered')
grid
【我是小蜜蜂,知识的搬运工!】