发明名称 一种适用于浅海低频条件的水声信号盲解卷方法
摘要 本发明涉及一种适用于浅海低频条件的水声信号盲解卷方法,技术特征在于:利用波数估计方法从水平阵接收信号中估计出波导中传播的简正波模态的水平波数;然后,根据估计的简正波模态水平波数信息构造一组权向量,并利用该权向量对水平阵接收信号进行模态分解处理,从中恢复出源信号波形。本发明给出的这种约束最优权向量也可以用于进一步改善Sabra等人提出的水平阵人工时反盲解卷方法的解卷性能。
申请公布号 CN102333052B 申请公布日期 2014.04.16
申请号 CN201110139594.0 申请日期 2011.05.26
申请人 西北工业大学 发明人 孙超;郭国强;杨益新;卓颉
分类号 H04L25/02(2006.01)I;H04L27/26(2006.01)I 主分类号 H04L25/02(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种适用于浅海低频条件的水声信号盲解卷方法,其特征在于步骤如下:步骤1:利用水平线列阵接收水下远程目标辐射的低频信号,并对每个阵元接收的信号进行快速傅立叶变换FFT,将各阵元接收的时域信号转换到频域上,得到接收信号的频域向量形式:P(r,z<sub>r</sub>,z<sub>s</sub>,ω<sub>i</sub>),其中:i=1,…,N<sub>FFT</sub>,N<sub>FFT</sub>为FFT变换点数,z<sub>r</sub>为水平阵位于水下深度,z<sub>s</sub>为目标深度,ω<sub>i</sub>为第i点傅立叶变换的角频率,r=[r<sub>1</sub>,r<sub>2</sub>,…,r<sub>N</sub>]为声源与每个阵元间的水平距离,N为水平阵阵元个数;步骤2:采用高分辨的MUSIC参数估计方法从水平线列阵接收的每个频率对应的目标辐射信号中估计出简正波水平波数<img file="FDA0000415212900000011.GIF" wi="179" he="85" />M为高分辨的MUSIC参数估计方法中涉及的波导有效传播的模态阶数,m=1,2,…,M;步骤3:根据估计出的简正波水平波数<img file="FDA0000415212900000012.GIF" wi="170" he="85" />以一组最优权值<img file="FDA0000415212900000013.GIF" wi="398" he="157" />对水平阵接收信号进行模态分解处理,其中:<img file="FDA0000415212900000014.GIF" wi="128" he="79" />为解相关处理后的满秩的协方差矩阵,<img file="FDA0000415212900000015.GIF" wi="458" he="141" />R<sub>l</sub>(ω<sub>i</sub>)为第l个子阵的协方差矩阵,L为子阵的个数;<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>e</mi><mo>^</mo></mover><mi>m</mi></msub><mo>=</mo><msup><mrow><mo>[</mo><mn>1</mn><mo>,</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><msub><mover><mi>k</mi><mo>^</mo></mover><mi>m</mi></msub><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>d</mi></mrow></msup><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mover><mi>k</mi><mo>^</mo></mover><mi>m</mi></msub><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>d</mi></mrow></msup><mo>]</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths>d为水平线阵阵元间距;步骤4:分别对每个频点ω<sub>i</sub>,其中i=1,…,N<sub>FFT</sub>,上的水平阵接收信号P(r,z<sub>r</sub>,z<sub>s</sub>,ω<sub>i</sub>)重复步骤2和步骤3,得到模态分解后的<img file="FDA0000415212900000017.GIF" wi="550" he="84" />将所有频点的<img file="FDA0000415212900000018.GIF" wi="574" he="97" />进行逆傅立叶变换IFFT,得到源信号波形y(t)=IFFT(Y(ω)),实现水平阵接收信号的盲解卷处理。
地址 710072 陕西省西安市友谊西路127号