[发明专利]噪声功率谱线谱检测方法

摘要

一种噪声功率谱线谱检测方法：设噪声信号序列为s(n)，对其做功率谱估计，获取功率谱p1(n)及对数功率谱p2(n)，n为非负整数；利用滑动窗正交多项式相关拟合获取平滑谱，对功率谱p1(n)及对数功率谱p2(n)进行滑动窗正交多项式拟合，得到噪声信号的功率平滑谱ps1(n)和噪声信号的对数功率平滑谱ps2(n)，计算功率谱和对数功率谱的差值谱并规范化，pd1(n)＝[p1(n)-ps1(n)]/std{[p1(n)-ps1(n)]}pd2(n)＝[p2(n)-ps2(n)]/std{[p2(n)-ps2(n)]}，pd1(n)，pd2(n)分别为功率谱规范差值谱和对数功率谱规范差值谱，设定幅度门限G1和对数门限G2，提取功率谱和对数功率谱线谱pl1(n)，pl2(n)：，，综合功率谱和对数功率谱获取目标线谱pl(n)，谱值以对数形式给出：。

主权项

1.一种噪声功率谱线谱检测方法，其特征在于， A.设噪声信号序列为s(n)，对其做功率谱估计，获取功率谱p1(n)及对数功率谱p2(n)，n为非负整数， B.利用滑动窗正交多项式相关拟合获取平滑谱 设功率谱数据p(n)长度N点，取一长为M点的窗，M＜N，步长m点，m＜M，设窗的初始起点w＝0与功率谱的起点重合，N、M及m为正整数，w非负整数，w为窗的起点位置， a)预先计算5组规范止交序列，计算方法： 其中i＝0，1，2，3，4，n＝0，1，...，M-1，对x_i(n)进行施密特Schimidt止交化，获得规范止交向量y_i(n)，i＝0，1，2，3，4， b)在功率谱数据p(n)中取w，w+1，...，w+M-1点作为待处理的窗口功率谱数据s_M(n)，即： s_M(n)＝p(n+w)，n＝0，1，...，M-1 c)计算窗口功率谱数据s_M(n)和规范止交向量y_i(n)的相关系数r_i， d)计算窗口拟合谱取窗口拟合谱s_M(n)前m点数据填入平滑谱数据ss(n)中，即ss(n+w)＝s_M(n)，n＝0，1，...，m-1 e)当w+m+M＜N-1时，令w＝w+m，返回b)，否则，进入f)， f)当w+m+M＞＝N-1时，记w0＝w，再令w＝N-1-M，在依次执行b)、c)后，计算 取ss(w0+m+n)＝s_M(w0+m-w+n)，其中n＝0，1，...，M-1+w-w0-m， C.根据B方法，对功率谱p1(n)及对数功率谱p2(n)进行滑动窗正交多项式拟合，得到噪声信号的功率平滑谱ps1(n)和噪声信号的对数功率平滑谱ps2(n)， D.计算功率谱和对数功率谱的差值谱并规范化， pd1(n)＝[p1(n)-ps1(n)]/std{[p1(n)-ps1(n)]} pd2(n)＝[p2(n)-ps2(n)]/std{[p2(n)-ps2(n)]} pd1(n)，pd2(n)分别为功率谱规范差值谱和对数功率谱规范差值谱， E.设定幅度门限G1和对数门限G2，提取功率谱和对数功率谱线谱pl1(n)，pl2(n)： F.综合功率谱和对数功率谱获取目标线谱pl(n)，谱值以对数形式给出： 。