[发明专利]一种辐射噪声线谱源深度辨识方法

权利要求书

1.一种辐射噪声线谱源深度辨识方法，其特征在于包括如下步骤：

(1)初始化处理参数；

(2)读入辐射噪声线谱幅度值，更新线谱幅度序列A(n)；

(3)判断辐射噪声线谱幅度序列对应的时间长度T是否满足要求，若满足，进行步骤(4)，否则返回步骤(2)；

(4)将辐射噪声线谱幅度序列A(n)划分为Q个子序列，计算每个子序列中辐射噪声线谱幅度的起伏指数，并构成线谱幅度起伏指数序列LSAFI(q)；

(5)统计步骤(4)得到的线谱幅度起伏指数序列LSAFI(q)中各指数超过门限λ的概率p，若p0.5，则判断对应的线谱源为水面目标；若p0.5，则判断对应的线谱源为水下目标；p＝0.5时，则对应的线谱源属性无法判断；

其中，步骤(4)具体包括如下步骤：

(4.1)将辐射噪声线谱幅度序列A(n)划分为Q个子序列，Q≥300，每个子序列的长度为L，L≥50，相邻两个子序列数据重叠长度为K，L/2≤KL，第q个子序列为A_q(k_q)，1+(L-K)(q-1)≤k_q≤2L+(L-K)(q-1)，q＝1,2,…,Q；

(4.2)根据式(1)、(2)，分别计算每个子序列的一阶幂平均函数和三阶幂平均函数，构成序列PMF₁(q)、PMF₃(q)；

(4.3)根据得到的序列PMF₁(q)、PMF₃(q)，计算辐射噪声线谱幅度序列的每个子序列的线谱幅度起伏指数，并构成辐射噪声线谱幅度起伏指数序列LSAFI(q)，如式(3)所示：

2.根据权利要求1所述的一种辐射噪声线谱源深度辨识方法，其特征在于，步骤(1)初始化如下参数：

(1.1)设定用于保存辐射噪声线谱幅度值的序列为A(n)，n＝1,2,…,N，并令N＝1；

(1.2)设定辐射噪声线谱幅度值更新时间间隔为Δt；

(1.3)设定辐射噪声线谱幅度序列更新累积时间长度Γ；

(1.4)设定线谱幅度起伏指数门限值为λ。

3.根据权利要求1所述的一种辐射噪声线谱源深度辨识方法，其特征在于，步骤(2)具体包括如下步骤：读入一个辐射噪声线谱幅度值并赋给变量a，并令A(N)＝a。

4.根据权利要求2所述的一种辐射噪声线谱源深度辨识方法，其特征在于，步骤(3)具体包括如下步骤：

(3.1)根据步骤(2)得到的辐射噪声线谱幅度序列A(n)，n＝1,2,…,N，计算A(n)对应的时间长度T＝N×Δt；

(3.2)若T≥Γ，进行步骤(4)；否则，令N＝N+1，并返回步骤(2)。

5.根据权利要求1所述的一种辐射噪声线谱源深度辨识方法，其特征在于，步骤(5)具体包括如下步骤：

(5.1)根据步骤(4)得到的LSAFI(q)，统计其中各起伏指数超过门限λ的概率p，如式(4)所示：

式中M{·}表示辐射噪声线谱幅度起伏指数序列中大于等于门限λ的起伏指数值个数；

(5.2)若p0.5，判断该辐射噪声线谱源为水面目标；若p0.5，判断该辐射噪声线谱源为水下目标；若p＝0.5，则无法判断该辐射噪声线谱源为于水面目标还是水下目标。