|
滤除工频信号
本帖最后由 wdhd 于 2016-9-14 10:37 编辑
肌电数据中始终有一个50赫的工频信号
图a中兰色为N=2048阶fft振幅谱a1,红色为apfft振幅谱a2,a1和a2的频谱相同的信号为稳态信号,即数据从头到尾都含有这个频率成份的信号.在图a中,k=21处a1和a2重合(相当21*5000/2048=51.2赫).只有这个工频信号是稳态信号,其它频率成分是随机的. 取原数据二通道中任一段作fft都有k=21这条谱线
用校正法时将k=21处的频率校正值(49.8918Hz),相位校正值(8.1753度)和振幅校正值(18.9734)算出,重构时域工频干扰信号,从原信号(图b兰色)中减去工频信号,得滤除工频的肌电信号(图b红色),
滤除工频的肌电信号再作谱分析,图c中兰色为FFT振幅谱a1’,红色为apFFT振幅谱a2’,a1’和a2’频谱中 k=21处峰值消除.
时域中原肌电信号(图b兰色)和滤除工频的肌电信号(图b红色)相差不大,频域中a1’和a2’不重合,表明信号非稳态信号,从时域图和频谱图中要看出毛病来还要下大功夫.有一文说用频谱肌电图中可得肌电中间频率,平均能量频率,不知求时有否改善
不过这个例子说明校正法可以从实际信号中滤除某一频率,它好像把某一频率挑了出来,做到单频陷波器,若用FIR滤波器难以做到.
[ 本帖最后由 zhwang554 于 2009-3-12 18:55 编辑 ] |
|