发明名称 一种对原始数据低失真的谐波噪声干扰和白噪声干扰的去除方法
摘要 本发明公开一种对原始信号低失真的谐波噪声干扰和白噪声干扰的去除方法;通过对观测数据信号进行分析,提取出其所包含的谐波噪声干扰的频率值;采用与谐波噪声干扰的基波频率成正整数倍的采样频率对观测数据信号进行重采样,利用基波频率与采样频率的比值对重采样后的观测数据信号的高频率分量在时域进行下采样处理,使各个下采样子带中所包含的谐波噪声干扰分量为常数值;应用基线偏移的方法将谐波噪声干扰去除;对去除谐波噪声干扰的信号进行去除白噪声干扰的处理;本发明可以低失真的将观测数据信号中包含的谐波噪声干扰和白噪声干扰去除,提高数据处理的准确度。
申请公布号 CN102624349A 申请公布日期 2012.08.01
申请号 CN201210068379.0 申请日期 2012.03.15
申请人 北京航空航天大学 发明人 徐立军;孔德明;李小路
分类号 H03H7/01(2006.01)I 主分类号 H03H7/01(2006.01)I
代理机构 代理人
主权项 1.一种对原始信号低失真的谐波噪声干扰和白噪声干扰去除方法,其特征在于包括以下四个步骤:步骤一、根据谐波噪声干扰信号的性质,对观测数据信号进行分析,得出其中所包含的谐波噪声干扰信号的谐波频率值;步骤二、对观测数据信号y进行重采样;当观测数据信号y的采样频率与谐波噪声干扰信号r的基波频率f<sub>1</sub>为正整数倍关系时,直接进行步骤三,当观测数据信号y的采样频率与谐波噪声干扰信号r的基波频率f<sub>1</sub>之间不为正整数倍关系时,采用与谐波噪声干扰信号r的基波频率f<sub>1</sub>成正整数倍且满足采样定理的最小值作为采样频率f<sub>s</sub>对观测数据信号y进行重采样;步骤三、通过低通滤波器将观测数据信号y中不包含谐波噪声干扰的低频率分量y<sub>LP</sub>和包含谐波噪声干扰的高频率分量y<sub>HP</sub>分别提取出来,并对高频率分量y<sub>HP</sub>进行参数为f<sub>1</sub>/f<sub>s</sub>的下采样处理,通过基线偏移的方法去除其各个下采样子带中的谐波噪声干扰信号分量r<sub>j</sub>[k],再采用上采样处理的方式对去除谐波噪声干扰后的信号进行重构;谐波噪声干扰信号r的离散数学模型为:<img file="FDA0000143577400000011.GIF" wi="811" he="120" />其中a<sub>i</sub>,f<sub>i</sub>和<img file="FDA0000143577400000012.GIF" wi="36" he="37" />分别为第i阶谐波噪声干扰的幅值、频率和初始相位变量,观测数据信号y的离散数学模型表示为:<img file="FDA0000143577400000013.GIF" wi="882" he="120" />其中x为原始数据信号,z为白噪声干扰信号,当观测数据信号的采样频率f<sub>s</sub>与谐波噪声干扰的基波频率f<sub>1</sub>的比值f<sub>s</sub>/f<sub>1</sub>为一个正整数值时,在时域内对观测数据信号y进行参数为f<sub>1</sub>/f<sub>s</sub>的下采样处理,从y的第j个采样点开始进行,该下采样子带y<sub>j</sub>中采样点之间的离散时间间隔为t<sub>j</sub>[k]=k/f<sub>1</sub>+j/f<sub>s</sub>,其中k=0,1,2,...,K-1,j=0,1,2,...,f<sub>s</sub>/f<sub>1</sub>-1,K为该子带的长度,j为该子带的阶数,将所得到的离散时间间隔t<sub>j</sub>[k]代入谐波噪声干扰信号r的离散数学模型中得:<img file="FDA0000143577400000014.GIF" wi="812" he="120" />即观测数据信号y的第j阶下采样子带中所包含的谐波噪声干扰信号分量r<sub>j</sub>[k]为一个常数值,通过基线偏移的方法将其去除,由于观测数据信号y的直流分量<img file="FDA0000143577400000015.GIF" wi="29" he="60" />会影响对各个下采样子带中的谐波噪声分量r<sub>j</sub>[k]进行准确的估计,因此采用截止频率为f<sub>c</sub>(f<sub>c</sub><f<sub>1</sub>)的低通滤波器F<sub>LP</sub>对观测数据信号进行低通滤波,将观测数据信号的低频率分量y<sub>LP</sub>提取出来,对应的高频率分量为y<sub>HP</sub>=y-y<sub>LP</sub>,y的直流分量<img file="FDA0000143577400000021.GIF" wi="31" he="62" />仅包含于其低频率分量y<sub>LP</sub>中,对高频率分量y<sub>HP</sub>进行参数为f<sub>1</sub>/f<sub>s</sub>的下采样处理,各个下采样子带的基线为其算术平均值<img file="FDA0000143577400000022.GIF" wi="458" he="71" />通过基线偏移归零的方法:<img file="FDA0000143577400000023.GIF" wi="504" he="82" />将各个下采样子带中所包含的谐波噪声干扰信号分量去除,得到去除谐波噪声干扰后的下采样子带信号<img file="FDA0000143577400000024.GIF" wi="183" he="83" />对去除谐波噪声干扰后的下采样子带信号<img file="FDA0000143577400000025.GIF" wi="160" he="82" />进行参数为f<sub>s</sub>/f<sub>1</sub>的上采样重构处理,得到高频率分量<img file="FDA0000143577400000026.GIF" wi="107" he="74" />将其与观测数据信号y的低频率分量y<sub>LP</sub>叠加,得到去除谐波噪声干扰后的信号<img file="FDA0000143577400000027.GIF" wi="316" he="74" />将<img file="FDA0000143577400000028.GIF" wi="40" he="74" />的采样频率f<sub>s</sub>还原为观测数据信号y的采样频率初始值,完成谐波噪声干扰的去除过程;步骤四、对去除谐波噪声干扰后的信号<img file="FDA0000143577400000029.GIF" wi="40" he="74" />进行去除白噪声干扰的处理。
地址 100191 北京市海淀区学院路37号