[发明专利]基于希尔伯特黄变换和自适应陷波的非平稳干扰抑制方法

权利要求书

1.基于希尔伯特黄变换和自适应陷波的非平稳干扰抑制方法，其特征在于具体步骤如下：

步骤一、对基带接收数据的采样结果进行归一化；

步骤二、对归一化后的数据，按照希尔伯特黄变换中的经验模式分解算法进行分解和希尔伯特变换，得到希尔伯特谱，具体实现步骤为：

A、根据希尔伯特黄变换的经验模式分解，将经步骤一得到的归一化后的数据 分解为M个内蕴模式函数(IMF)c_i(n)分量之和，即

其中，r_M(n)为分解后剩余分量，M为分解得到的内蕴模式函数的个数，M的值根据对r_M(n)的大小要求来确定；

B、对每一个内蕴模式函数分量进行希尔伯特变换，得到相应的解析信号，将所得到的M个解析信号带入 中， 表示为

其中，H[c_i(n)]为c_i(n)的布尔伯特变换，a_i(n)为c_i(n)解析信号的瞬时幅度，f_i(n)为c_i(n)的瞬时频率，θ_i(n)为c_i(n)的瞬时相位，T_s为采样间隔，i为内蕴模式函数分量的序号，即i＝1～M；

b为采样点序号，b＝1～n；n为第n个采样点；

C、将 表示成时间n、瞬时频率ω、幅值a_i(n)的三维谱图，即希尔伯特谱H(ω，n)：

其中，b_i为开关因子，当ω＝ω_i时，b_i＝1，否则b_i＝0；M代表内蕴模式函数分量的个数，i为内蕴模式函数分量的序号，即i＝1～M；

步骤三、根据得到的希尔伯特谱确定抑制噪声和直扩信号谱的门限电平，将高于门限电平的谱保留，将低于门限电平的谱清零，得到非平稳干扰的希尔伯特谱，具体实现步骤为：

(一)通过连续平均抑制算法，确定抑制噪声和直扩信号谱的门限电平，即

①将所有希尔伯特谱线作为初始谱线，并把它们放入集合I，得到所有希尔伯特谱线值的累加和A_m，将A_m除以希尔伯特谱线总数，求出希尔伯特谱线的均值E[|H(ω，n)|]：

E[|H(ω，n)|]＝A_m/(N×M)

其中，M为内蕴模式函数分量的个数，N为采样点数；

②根据希尔伯特谱线的均值E[|H(ω，n)|]，采用下式估计希尔伯特谱线的标准差σ[|H(ω，n)|]：

根据得到的标准差σ[|H(ω，n)|]，确定出门限η

η＝E[|H(ω，n)|]+3σ[|H(ω，n)|]

③将集合I中的希尔伯特谱线的模与门限η进行比较，大于该门限的希尔伯特谱线被认为是被干扰的希尔伯特谱线，从集合I中去除，集合I被更新；

④更新集合I后，更新A_m和E[|H(ω，n)|]，之后重复②和③，直到没有希尔伯特谱线的模大于门限η为止，最终得到没有希尔伯特谱线的模大于门限η的集合I′；

⑤以集合I′中的希尔伯特谱线，估计最终的门限η＝E[|H(ω，n)|]+3σ[|H(ω，n)|]；

(二)将高于最终门限η电平的希尔伯特谱保留，将低于门限η电平的希尔伯特谱清零，从而得到非平稳干扰的希尔伯特谱；

步骤四、依照非平稳干扰的希尔伯特谱，求取出非平稳干扰的瞬时频率，由该瞬时频率决定自适应无限冲激响应格型陷波器的参数，由此得到自适应陷波器冲激响应；

步骤五、将步骤一中的采样结果与自适应陷波器冲激响应进行卷积运算，得到抑制非平稳干扰后的数据，将该数据和同步的PN码序列进行相关运算，完成接收数据的解扩。