CN107505596A - 基于双扩展水声信道环境下的mimo主动探测信号设计与检测系统和方法 - Google Patents
基于双扩展水声信道环境下的mimo主动探测信号设计与检测系统和方法 Download PDFInfo
- Publication number
- CN107505596A CN107505596A CN201710605628.8A CN201710605628A CN107505596A CN 107505596 A CN107505596 A CN 107505596A CN 201710605628 A CN201710605628 A CN 201710605628A CN 107505596 A CN107505596 A CN 107505596A
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- msubsup
- msub
- msup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于双扩展水声信道环境下的MIMO主动探测信号设计与检测系统和方法，适用于快速时变的双扩展水声信道条件下的水声探测与定位领域。本发明提出的基于周期性自相关和互相关结构的ZCZ序列波形设计的优势在于不用过多的考虑和设计不同的正交结构，只利用波形序列的周期性循环移位相关特性即可实现干扰的完全消除。通过波束形成和匹配滤波处理，区分不同发射阵元回波信号的虚拟数据向量，扩展虚拟阵元孔径，提高目标分辨力，提升整个系统的空‑时分集增益。一次可以实现所有方位的检测，与现有传统波束形成技术相比，节省探测时间。
Description
技术领域
本发明涉及一种基于双扩展水声信道环境下的MIMO主动探测信号设计与检测系统和方法，适用于快速时变的双扩展水声信道条件下的水声探测与定位领域。
技术背景
声纳信号与通信信号的本质区别在于，通信信号包含通信的全部信息，而声纳信号毫无信息，只是信息的运载工具，当声纳的发射信号到达探测目标时，信号就会、发生反射，目标的全部信息都蕴藏在回波信号中。雷达信号的发射波形和他能运载信息有直接的关系，波形的选择会直接决定和影响声纳系统的性能参数，包括信噪比、距离分辨力、多普勒分辨力、时延和多普勒模糊度函数等。传输过程重点关注如何从声纳回波信号中尽可能多的提取目标的有用信息。而在波导水声环境中，信道具有快速时变特性和时延多普勒双扩展特性，给目标信号提取造成困难。因此要合理设计声纳信号发射波形和信号检测方法，克服时变双扩展水声信道的复杂环境，提高MIMO声纳系统空-时分集增益，实现多目标检测。
发明内容
本发明提出一种基于快时变、双扩展水声信道环境下的MIMO主动探测信号设计与检测系统和方法。
现有的MIMO主动探测信号设计，基本原理上都采用非周期的相关结构，这些发射信号虽然具有较好的非周期自相关与互相关特性，仍有一定的旁瓣干扰。本发明提出一种基于周期性自相关和互相关结构的ZCZ序列波形设计，这种波形的优势在于不用过多的考虑和设计不同的正交结构，只利用波形序列的周期性循环移位相关特性即可实现干扰的完全消除。通过波束形成和匹配滤波处理，区分不同发射阵元回波信号的虚拟数据向量，扩展虚拟阵元孔径，提高目标分辨力，提升整个系统的空-时分集增益，一次可以实现所有方位的检测。与现有传统波束形成技术相比，节省探测时间。
为克服波导环境水声信道快时变、大时延、严重多普勒扩展的缺点，实现多目标检测。本发明采取的技术方案为：周期性循环移位相关ZCZ序列设计和MIMO水声信号检测，系统主要装置包括：
ZCZ序列生成器，采用循环移位和循环前缀的方式，生成适用于MIMO主动探测的信号；
发射功率放大器组，与ZCZ序列生成器相连，用于发射信号功率放大；
接收功率放大器组，与接收水听器阵列相连，用于实现接收信号增益控制；
发射换能器阵列，与发射功率放大器组相连，用于将放大后的探测信号由电信号转换为声信号，并将声信号发射至探测水域；
接收水听器阵列，用于将接收到的回波声信号转换为电信号；
相干处理单元，将接收到的MIMO信号进行处理，得到MIMO波束模式，对MIMO波束模式进行估计，得到潜在目标方位，然后对匹配滤波输出的结果，利用广义似然比检测器检测估计的潜在目标方位上是否存在真实目标。
其中，umimo(θ)为匹配输出之后的MIMO波束模式，θ为可能存在目标的方位，H1表示估计出的潜在目标方位存在真实的目标，H0表示估计出的潜在目标方位不存在真实的目标，δ为广义似然比检测器的检测阈值。
与现有技术相比，本发明创新之处在于：
基于周期性循环移位相关的ZCZ发射信号设计，能够对抗波导时变水声信道，可以实现完全干扰消除；
通过MIMO阵处理，可以获得空间的分集增益，与传统的相控阵波束形成相比，扩展了虚拟孔径和空间自由度，能够获得更窄的主瓣和更低的旁瓣，提升了目标检测的分辨力和回波检测的信噪比；
与传统相控阵波束形成相比，本发明一次探测可实现对所有方位的检测，避免了相控扫描，节约了时间。
附图说明
图1是本发明中水声MIMO定位收发系统整体原理图。
图2是收发合置阵列示意图。
图3是本发明中的MIMO探测系统和相控阵探测系统波束模式图对比。
图4是ZCZ信号经过真实波导水声信道干扰抑制情况示意图。
具体实施方式
下面结合附图和具体实例对本发明做进一步的描述，但本发明的实施和保护范围不限于此。
本发明基于ZCZ探测信号设计的MIMO水声探测装置，包括：
ZCZ序列生成器，采用循环移位和循环前缀的方式，生成适用于MIMO主动探测的信号；
发射功率放大器组，与ZCZ序列生成器相连，用于发射信号功率放大；
接收功率放大器组，与接收水听器阵列相连，用于实现接收信号增益控制；
发射换能器阵列，与发射功率放大器组相连，用于将放大后的探测信号由电信号转换为声信号，并将声信号发射至探测水域；
接收水听器阵列，用于将接收到的回波声信号转换为电信号；
相干处理单元，将接收到的MIMO信号进行处理，得到MIMO波束模式，对MIMO波束模式进行估计，得到潜在目标方位，然后对匹配滤波输出的结果，利用广义似然比检测器检测估计的潜在目标方位上是否存在真实目标。
图1所示为本发明信号生成与探测方法的工作原理图，包括：ZCZ序列生成器、发射功率放大器组、接收功率放大器组、发射换能器阵列、接收水听器阵列、相干处理单元。
发射阵和接收阵布放如图2所示，发射阵元个数为Nt＝6，接收阵元个数为Nr＝7，发射阵元间距为10cm，接收阵元间距7cm。
需要说明的是，虽然图1为了更加简洁的表达，使用了Nt＝6，Nr＝7，但是本发明中，发射换能器阵列Nt与接收换能器阵列Nr不仅限于此特例，只要在收发条件允许的范围内，可以通过调整收发阵元个数，来得到合理的波束模式。
利用本发明的MIMO探测信号发生装置进行探测信号的生成，包括如下步骤：
步骤1：对于一个ZCZ序列结构，Fn代表具有M个ZCZ序列的簇，每个ZCZ序列的长度为L，可以进一步表示为F(L,M,ZCZ),ZCZ为零相关区间长度，当n为0时，用于生成ZCZ序列的基础矩阵可表示为：
步骤2：进一步推广，L0＝2为开始的序列长度，并且让起始n＝1，更长的ZCZ序列矩阵可以表示为:F(L,M,ZCZ)＝(22nL0,2n+1,2n+1)＝(8,4,3)。可以用式(4)来表示结构：
式(4)中矩阵表示对取反，F1的每一行都是一组ZCZ序列。
步骤3：推广ZCZ序列簇，对于Fn-1＝(22(n-1)L0,2n,2n-1+1)，一个更大的序列簇Fn＝(22nL0,2n+1,2n+1)可表示为：
式(5)中，和可以用式(6)处理方式来实现：
将式(5)生成的ZCZ信号进行以下处理：
步骤1：将生成的一路ZCZ序列进行循环移位处理，移位生成几路信号由发射换能器阵元个数来决定。
步骤2：为了保证生成信号序列具有良好的循环移位自相关特性，需要将每路生成的信号加循环前缀处理。
经过循环移位和循环前缀处理的ZCZ信号，经过发射驾驶向量整合，由发射功率放大器组将信号放大，最终由MIMO换能器将电信号转换成声信号传播出去。
整个信号收发和检测过程包括如下步骤：
步骤1：信号发射机发射一组M个经过循环移位+循环前缀处理的ZCZ循环移位序列，信号以向量形式表示为s[n]＝(s1[n],s2[n],…,sM[n])T，发射信号的协方差矩阵可以表示为：
式(7)中，M≥1的正整数，n表示正交信号第n时刻的采样点,1≤n≤L，L表示采样得到的码字序列长度，βij表示第i路发射信号si[n]和第j发射信号sj[n]之间的互相关系数，sH[n]是s[n]共轭转置运算，Rs可以通过SVD分解为：
Rs＝UΛUH (8)
其中U为酉矩阵，Λ为对角阵，UH为U的共轭转置。当发射信号完全正交时，Rs＝IM。
步骤2：由Nt个发射阵元和Nr个接收阵元组成的探测系统，发射阵列驾驶向量at(θ)和接收阵列驾驶向量ar(θ)分别表示为式(9)和式(10)：
其中at(θ)为发射换能器阵列的驾驶向量，ar(θ)为接收换能器阵列的驾驶向量，θ为潜在目标方位角度，Nt为发射换能器阵列的阵元个数，Nr为接收水听器阵列的阵元个数，dt为发射换能器阵元间距，dr为接收水听器阵元间距，f为发射机发射信号中心频率，c为水中声速，(·)T表示矩阵转置运算。
步骤3：经过驾驶向量波束形成的Nt个正交信号s[n]，经过功率放大器组发送至发射换能器阵列，发射换能器阵列将放大后的正交信号由电信号转化成声信号，并将声信号发送至待探测水域；信号经过待检测目标反射，由接收水听器阵列接收回波信号，并由接收水听器功率放大器组将回波声信号转化为电信号。接收到的信号可以表示为：
式(11)中，r[n]为接收水听器接收到的电信号，r[n]＝[r1[n],r2[n],…,rNr[n]]T，n＝1,2,…L为Nr元接收阵接收到的信号序列，α(θ)为信号传播衰减系数，w[n]为接收水听器阵列接收到的信号组与发射信号组不相关的加性噪声向量。服从的复高斯分布，其中是噪声功率，是秩为Nr的单位矩阵。
把长度为L的接收信号序列表示为矩阵的形式，采取阵列数据存储的形式：
其中R＝[r[1],r[2],…,r[L]],S＝[s[1],s[2],…s[L]]，则N＝[w[1],w[2],…，w[L]]，
步骤4：根据式(11)描述的MIMO接收信号模型r[n]来检测和定位目标，要通过与发射的ZCZ循环移位信号做匹配滤波处理来获得，通过匹配滤波处理得到的充分检验统计矩阵为：
发射信号具有完美的循环移位自相关特性，充分检验统计矩阵Ymimo进一步简化为：
将式(14)得到的充分统计量做列向量化处理，得到如式(15)所示的充分统计量的列向量化表示：
ymimo＝vec(Ymimo)＝α(θ0)d(θ0)+v (15)
式(15)中，θ0是目标潜在方位，α(θ0)为潜在目标方位角度的信号传播衰减系数，是长度为NtNr×1的接收水听器阵列匹配输出响应，是克罗内克积，是服从的复高斯噪声，其中是向量化后的噪声功率，是秩为MN的单位矩阵。
步骤5：根据式(15)得到的充分统计量ymimo，由式(16)得到系统收发波束模式图和含有潜在目标方位的最大似然估计
式(16)中，umimo(θ)为匹配输出之后的MIMO模式，θ为可能存在目标的方位，||表示求绝对值运算，||||表示向量求模运算，表示对ar(θ)做共轭转置运算，表示对at(θ)做共轭转置运算，表示潜在目标所在方位的最大似然估计。
使用式(16)表达式的波束模式，对于每一个-90°≤θ≤90°，计算其对应的波束模式数值，得到如图3实线所示的波束模式图，根据式(17)和波束模式图，搜索umimo(θ)的最大值所对应的角度θ，即为估计出的潜在目标的方位
步骤6：利用式(17)对MIMO波束模式进行估计，得到潜在的目标方位，然后利用式(18)的广义似然比检测器检潜在目标方位上是否真实存在目标：
式(18)中，是估计出的潜在目标方位上的波束模式，H1表示估计出的潜在目标方位存在真实的目标，H0表示估计出的潜在目标方位不存在真实的目标，δ为广义似然比检测器的检测阈值。由预先设定的虚警概率Pf来决定，当时认为潜在目标方位目标真实存在，否则认为潜在目标方位目标并不真实存在。
从图3可以看出，估计出来的目标方位与实际情况中待检目标方位相符。
与传统相控阵波束形成波束模式相比较，本发明的探测方法的波束模式具有更窄的主瓣和更窄的旁瓣，第一旁瓣比常规相控阵的第一旁瓣低13dB。
波束模式具有更窄的主瓣能够获得更高的目标分辨力，更低的旁瓣能够减小背景噪声和虚警的干扰，提高接收信号信噪比，提高目标检测概率。
本发明提出的基于循环移位的ZCZ序列的信号检测，能够对抗波导时变双扩展水声信道，能够实现干扰的完全消除，图4表示循环移位相关的ZCZ信号经过真实的水声信道。与传统的伪随机序列相比，循环移位相关的ZCZ序列能够实现干扰的完全消除，具有更好的检测性能。
Claims (5)
1.一种基于双扩展水声信道环境下的MIMO主动探测信号设计与检测系统，其特征在于包括:
ZCZ序列生成器，采用循环移位和循环前缀的方式，生成适用于MIMO主动探测的信号；
发射功率放大器组，与ZCZ序列生成器相连，用于发射信号功率放大；
发射换能器阵列，与发射功率放大器组相连，用于将放大后的探测信号由电信号转换为声信号，并将声信号发射至探测水域；
接收水听器阵列，用于将接收到的回波声信号转换为电信号；
接收功率放大器组，与接收水听器阵列相连，用于实现接收信号增益控制；
相干处理单元，将接收到的MIMO信号进行处理，得到MIMO波束模式，对MIMO波束模式进行估计，得到潜在目标方位，然后对匹配滤波输出的结果，利用广义似然比检测器检测估计的潜在目标方位上是否存在真实目标。
2.一种基于双扩展水声信道环境下的MIMO主动探测信号设计与检测方法，其特征在于包括如下步骤：
1)ZCZ序列生成器生成ZCZ序列；将生成的ZCZ序列进行循环移位处理，移位生成几路信号由发射换能器阵元个数来决定，将每路生成的信号加循环前缀处理；
2)经过循环移位和循环前缀处理的ZCZ信号，经过发射驾驶向量整合，由发射功率放大器组将信号放大，放大后的信号由发射功率放大器组发送至发射换能器阵列，由发射换能器阵列将电信号转换成声信号发送至待探测水域；
3)信号经过待检测目标反射，由接收水听器阵列接收回波信号，并由接收功率放大器组将回波声信号转化为电信号，相干处理单元对MIMO波束模式进行估计，得到潜在的目标方位，然后利用的广义似然比检测器检测潜在目标方位上是否真实存在目标。
3.根据权利要求2所述的基于双扩展水声信道环境下的MIMO主动探测信号设计与检测方法，其特征在于所述的步骤1)，具体为：
步骤1.1：对于一个ZCZ序列结构，Fn代表具有M个ZCZ序列的簇，每个ZCZ序列的长度为L，Fn进一步表示为F(L,M,ZCZ)，ZCZ为零相关区间长度，当n为0时，用于生成ZCZ序列的基础矩阵表示为：
F(L,M,ZCZ)＝(2,2,1) (1)
步骤1.2：以L0＝2为开始的序列长度，并且让起始n＝1，更长的ZCZ序列矩阵表示为:
F(L,M,ZCZ)＝(22nL0,2n+1,2n+1)＝(8,4,3)，用公式(2)来表示结构：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>F</mi>
<mn>1</mn>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>F</mi>
<mn>11</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mn>12</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mn>13</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mn>14</mn>
<mn>1</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>F</mi>
<mn>21</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mn>22</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mn>23</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mn>24</mn>
<mn>1</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>F</mi>
<mn>31</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mn>32</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mn>33</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mn>34</mn>
<mn>1</mn>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>F</mi>
<mn>41</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mn>42</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mn>43</mn>
<mn>1</mn>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mn>44</mn>
<mn>1</mn>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>F</mi>
<mn>11</mn>
<mn>0</mn>
</msubsup>
<msubsup>
<mi>F</mi>
<mn>11</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>F</mi>
<mn>12</mn>
<mn>0</mn>
</msubsup>
<msubsup>
<mi>F</mi>
<mn>12</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msubsup>
<mi>F</mi>
<mn>11</mn>
<mn>0</mn>
</msubsup>
<mo>)</mo>
<msubsup>
<mi>F</mi>
<mn>11</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msubsup>
<mi>F</mi>
<mn>12</mn>
<mn>0</mn>
</msubsup>
<mo>)</mo>
<msubsup>
<mi>F</mi>
<mn>12</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>F</mi>
<mn>21</mn>
<mn>0</mn>
</msubsup>
<msubsup>
<mi>F</mi>
<mn>21</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>F</mi>
<mn>22</mn>
<mn>0</mn>
</msubsup>
<msubsup>
<mi>F</mi>
<mn>22</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msubsup>
<mi>F</mi>
<mn>21</mn>
<mn>0</mn>
</msubsup>
<mo>)</mo>
<msubsup>
<mi>F</mi>
<mn>21</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msubsup>
<mi>F</mi>
<mn>22</mn>
<mn>0</mn>
</msubsup>
<mo>)</mo>
<msubsup>
<mi>F</mi>
<mn>22</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msubsup>
<mi>F</mi>
<mn>11</mn>
<mn>0</mn>
</msubsup>
<mo>)</mo>
<msubsup>
<mi>F</mi>
<mn>11</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msubsup>
<mi>F</mi>
<mn>12</mn>
<mn>0</mn>
</msubsup>
<mo>)</mo>
<msubsup>
<mi>F</mi>
<mn>12</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>F</mi>
<mn>11</mn>
<mn>0</mn>
</msubsup>
<msubsup>
<mi>F</mi>
<mn>11</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>F</mi>
<mn>12</mn>
<mn>0</mn>
</msubsup>
<msubsup>
<mi>F</mi>
<mn>12</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msubsup>
<mi>F</mi>
<mn>21</mn>
<mn>0</mn>
</msubsup>
<mo>)</mo>
<msubsup>
<mi>F</mi>
<mn>21</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msubsup>
<mi>F</mi>
<mn>22</mn>
<mn>0</mn>
</msubsup>
<mo>)</mo>
<msubsup>
<mi>F</mi>
<mn>22</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>F</mi>
<mn>21</mn>
<mn>0</mn>
</msubsup>
<msubsup>
<mi>F</mi>
<mn>21</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>F</mi>
<mn>22</mn>
<mn>0</mn>
</msubsup>
<msubsup>
<mi>F</mi>
<mn>22</mn>
<mn>0</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mo>-</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>+</mo>
<mo>+</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>+</mo>
<mo>-</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mo>-</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mo>-</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>+</mo>
<mo>-</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>+</mo>
<mo>-</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mo>-</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mo>+</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mo>-</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>+</mo>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mo>-</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>+</mo>
<mo>-</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mo>-</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mo>-</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mn>4</mn>
<mo>&times;</mo>
<mn>8</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式(2)中，矩阵表示对取反，F1的每一行都是一组ZCZ序列；
步骤1.3：推广ZCZ序列簇，对于Fn-1＝(22(n-1)L0,2n,2n-1+1),Fn＝(22nL0,2n+1,2n+1)，一个更大的序列簇Fn＝(22nL0,2n+1,2n+1)表示为：
<mrow>
<msup>
<mi>F</mi>
<mi>n</mi>
</msup>
<mo>=</mo>
<msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>F</mi>
<mn>11</mn>
<mi>n</mi>
</msubsup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mn>1</mn>
<mi>M</mi>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>M</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>F</mi>
<mn>21</mn>
<mi>n</mi>
</msubsup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mn>2</mn>
<mi>M</mi>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>M</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mi>M</mi>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<mi>M</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mi>M</mi>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>M</mi>
<mo>)</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>M</mi>
<mo>)</mo>
<mi>M</mi>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>M</mi>
<mo>)</mo>
<mo>(</mo>
<mi>M</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>M</mi>
<mo>)</mo>
<mi>M</mi>
</mrow>
<mi>n</mi>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<msup>
<mi>M</mi>
<mo>&prime;</mo>
</msup>
<mo>&times;</mo>
<msup>
<mi>L</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤1.4：式(3)中，和用式(4)的处理方式来实现
<mrow>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>M</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>M</mi>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>M</mi>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>M</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>n</mi>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>M</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>n</mi>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>M</mi>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>M</mi>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>+</mo>
<mi>M</mi>
<mo>)</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>M</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>n</mi>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>+</mo>
<mi>M</mi>
<mo>)</mo>
<mi>M</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>M</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>n</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>+</mo>
<mi>M</mi>
<mo>)</mo>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>M</mi>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>+</mo>
<mi>M</mi>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>M</mi>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>F</mi>
<mrow>
<mi>i</mi>
<mi>M</mi>
</mrow>
<mi>n</mi>
</msubsup>
<mo>.</mo>
</mrow>
4.根据权利要求2所述的基于双扩展水声信道环境下的MIMO主动探测信号设计与检测方法，其特征在于所述的步骤2)，具体为：
步骤2.1：信号发射机发射一组M个经过循环移位+循环前缀处理的ZCZ循环移位序列，信号以向量形式表示为s[n]＝(s1[n],s2[n],…,sM[n])T，发射信号的协方差矩阵表示为：
式(5)中，M≥1的正整数，n表示正交信号第n时刻的采样点,1≤n≤L，βij表示第i路发射信号si[n]和第j发射信号sj[n]之间的互相关系数，sH[n]是s[n]共轭转置运算，Rs通过SVD分解为：
Rs＝UΛUH (6)
式(6)中U为酉矩阵，Λ为对角阵，UH为U的共轭转置，当发射信号完全正交时，Rs＝IM；
步骤2.2：发射阵列驾驶向量at(θ)和接收阵列驾驶向量ar(θ)分别表示为(7)和(8)：
<mrow>
<msub>
<mi>a</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>,</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;fd</mi>
<mi>t</mi>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>/</mo>
<mi>c</mi>
</mrow>
</msup>
<mo>,</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mn>2</mn>
<msub>
<mi>d</mi>
<mi>t</mi>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>/</mo>
<mi>c</mi>
</mrow>
</msup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>d</mi>
<mi>t</mi>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>/</mo>
<mi>c</mi>
</mrow>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>,</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;fd</mi>
<mi>t</mi>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>/</mo>
<mi>c</mi>
</mrow>
</msup>
<mo>,</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mn>2</mn>
<msub>
<mi>d</mi>
<mi>t</mi>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>/</mo>
<mi>c</mi>
</mrow>
</msup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mi>r</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>d</mi>
<mi>r</mi>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>/</mo>
<mi>c</mi>
</mrow>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中at(θ)为发射换能器阵列的驾驶向量，ar(θ)为接收换能器阵列的驾驶向量，θ为潜在目标方位角度，Nt为发射换能器阵列的阵元个数，Nr为接收水听器阵列的阵元个数，dt为发射换能器阵元间距，dr为接收水听器阵元间距，f为发射机发射信号中心频率，c为水中声速，(·)T表示矩阵转置运算；
步骤2.3：经过驾驶向量波束形成的Nt个正交信号s[n]，经过功率放大器组发送至发射换能器阵列，发射换能器阵列将放大后的正交信号由电信号转化成声信号，并将声信号发送至待探测水域；信号经过待检测目标反射，由接收水听器阵列接收回波信号，并由接收水听器功率放大器组将回波声信号转化为电信号。
5.根据权利要求2所述的基于双扩展水声信道环境下的MIMO主动探测信号设计与检测方法，其特征在于所述的步骤3)，具体为：
步骤3.1：接收水听器接收到的信号表示为：
<mrow>
<mi>r</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>a</mi>
<mi>t</mi>
<mi>T</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<mo>+</mo>
<mi>w</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>n</mi>
<mo>&le;</mo>
<mi>L</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式(9)中，r[n]为接收水听器接收到的电信号，n＝1,2,…L为Nr元接收阵接收到的信号序列，α(θ)为信号传播衰减系数，w[n]为接收水听器阵列接收到的信号组与发射信号组不相关的加性噪声向量，服从的复高斯分布，其中是噪声功率，是秩为Nr的单位矩阵；
步骤3.2：利用式(9)描述的MIMO接收信号模型r[n]来检测和定位目标，要通过与ZCZ循环移位信号做匹配滤波处理来获得，通过匹配滤波处理得到的充分检验统计矩阵为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>m</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<mi>r</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>&Lambda;</mi>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<msup>
<mi>U</mi>
<mi>H</mi>
</msup>
<mi>s</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<mo>|</mo>
<msup>
<mi>s</mi>
<mi>H</mi>
</msup>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<mi>s</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>a</mi>
<mi>t</mi>
<mi>T</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<msup>
<mi>U&Lambda;</mi>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mo>+</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<mi>w</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>&Lambda;</mi>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<msup>
<mi>U</mi>
<mi>H</mi>
</msup>
<mi>s</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<mo>|</mo>
<msup>
<mi>s</mi>
<mi>H</mi>
</msup>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<mi>s</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
发射信号具有完美的循环移位自相关特性，充分检验统计矩阵Ymimo进一步简化为：
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>m</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>a</mi>
<mi>t</mi>
<mi>T</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<mi>w</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<msup>
<mi>s</mi>
<mi>H</mi>
</msup>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<mo>|</mo>
<msup>
<mi>s</mi>
<mi>H</mi>
</msup>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<mi>s</mi>
<mo>&lsqb;</mo>
<mi>n</mi>
<mo>&rsqb;</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
将式(11)得到的充分统计量做列向量化处理，得到如式(12)所示的充分统计量的列向量化表示：
ymimo＝vec(Ymimo)＝α(θ0)d(θ0)+v (12)
式(12)中，θ0是目标潜在方位，α(θ0)为潜在目标方位角度的信号传播衰减系数，是长度为NtNr×1的接收水听器阵列匹配输出响应，是克罗内克积，是服从的复高斯噪声，其中是噪声功率，是秩为NtNr的单位矩阵；
步骤3.3：根据步骤3.2得到的充分统计向量ymimo，由式(13)得到系统收发波束模式图和含有潜在目标方位的最大似然估计
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>m</mi>
<mi>o</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msup>
<mi>d</mi>
<mi>H</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>|</mo>
<mo>|</mo>
<mi>d</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msubsup>
<mi>a</mi>
<mi>t</mi>
<mi>H</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>|</mo>
<msubsup>
<mi>a</mi>
<mi>t</mi>
<mi>H</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mrow>
<msup>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>N</mi>
<mi>r</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mrow>
<mi>&theta;</mi>
<mo>,</mo>
<mi>&alpha;</mi>
</mrow>
</munder>
<msub>
<mi>u</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>m</mi>
<mi>o</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式(13)中，umimo(θ)为匹配输出之后的MIMO波束模式，||表示求绝对值运算，|| ||表示向量求模运算，表示对ar(θ)做共轭转置运算，表示对at(θ)做共轭转置运算，表示潜在目标所在方位的最大似然估计；
步骤3.4：利用式(13)对MIMO波束模式进行估计，得到潜在的目标方位，然后利用式(14)的广义似然比检测器检测潜在目标方位上是否真实存在目标：
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>m</mi>
<mi>o</mi>
</mrow>
</msub>
<msubsup>
<mrow>
<mo>(</mo>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
<mtable>
<mtr>
<mtd>
<mo><</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>H</mi>
<mn>0</mn>
</msub>
</mtd>
</mtr>
</mtable>
<mtable>
<mtr>
<mtd>
<msub>
<mi>H</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>&GreaterEqual;</mo>
</mtd>
</mtr>
</mtable>
</msubsup>
<mi>&delta;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
式(14)中，为估计出的潜在目标方位上的波束模式，H1表示估计出的潜在目标方位存在真实的目标，H0表示估计出的潜在目标方位不存在真实的目标，δ为广义似然比检测器的检测阈值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710605628.8A CN107505596B (zh) | 2017-07-24 | 2017-07-24 | 基于双扩展水声信道环境下的mimo主动探测信号设计与检测系统和方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710605628.8A CN107505596B (zh) | 2017-07-24 | 2017-07-24 | 基于双扩展水声信道环境下的mimo主动探测信号设计与检测系统和方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107505596A true CN107505596A (zh) | 2017-12-22 |
| CN107505596B CN107505596B (zh) | 2020-10-16 |
Family
ID=60688783
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710605628.8A Active CN107505596B (zh) | 2017-07-24 | 2017-07-24 | 基于双扩展水声信道环境下的mimo主动探测信号设计与检测系统和方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107505596B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111147156A (zh) * | 2019-12-27 | 2020-05-12 | 浙江大学 | 一种基于水声通信数据再利用的多址接入交叉估计方法 |
| CN114500199A (zh) * | 2022-04-13 | 2022-05-13 | 鹏城实验室 | 一种基于空时簇处理的水声通信方法及相关设备 |
| CN115085838A (zh) * | 2022-04-22 | 2022-09-20 | 哈尔滨工程大学 | 基于虚拟发射阵列的水声探测通信一体化波形验证方法 |
| WO2023160714A1 (zh) * | 2022-02-28 | 2023-08-31 | 华为技术有限公司 | 多流感知信号生成方法、装置及系统 |
| CN117081894A (zh) * | 2023-07-17 | 2023-11-17 | 中国科学院声学研究所 | 一种利用信道稀疏特性的水声信号检测方法 |
| WO2024032562A1 (zh) * | 2022-08-09 | 2024-02-15 | 华为技术有限公司 | 一种通信方法及装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100767821B1 (ko) * | 2006-10-25 | 2007-10-17 | 한국해양연구원 | 음향 탐지 신호 은닉 방법 및 그 장치 |
| CN102608609A (zh) * | 2011-12-19 | 2012-07-25 | 浙江大学 | 一种基于mimo的高分辨水下目标探测装置及方法 |
| WO2016096186A1 (de) * | 2014-12-18 | 2016-06-23 | Thyssenkrupp Marine Systems Gmbh | Verfahren zur ortung eines objekts unter wasser |
| CN106330251A (zh) * | 2016-08-29 | 2017-01-11 | 东南大学 | 基于零相关带序列的水声通信系统多普勒扩展估计方法 |
-
2017
- 2017-07-24 CN CN201710605628.8A patent/CN107505596B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100767821B1 (ko) * | 2006-10-25 | 2007-10-17 | 한국해양연구원 | 음향 탐지 신호 은닉 방법 및 그 장치 |
| CN102608609A (zh) * | 2011-12-19 | 2012-07-25 | 浙江大学 | 一种基于mimo的高分辨水下目标探测装置及方法 |
| WO2016096186A1 (de) * | 2014-12-18 | 2016-06-23 | Thyssenkrupp Marine Systems Gmbh | Verfahren zur ortung eines objekts unter wasser |
| CN106330251A (zh) * | 2016-08-29 | 2017-01-11 | 东南大学 | 基于零相关带序列的水声通信系统多普勒扩展估计方法 |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111147156A (zh) * | 2019-12-27 | 2020-05-12 | 浙江大学 | 一种基于水声通信数据再利用的多址接入交叉估计方法 |
| CN111147156B (zh) * | 2019-12-27 | 2021-03-26 | 浙江大学 | 一种基于水声通信数据再利用的多址接入交叉估计方法 |
| WO2023160714A1 (zh) * | 2022-02-28 | 2023-08-31 | 华为技术有限公司 | 多流感知信号生成方法、装置及系统 |
| CN114500199A (zh) * | 2022-04-13 | 2022-05-13 | 鹏城实验室 | 一种基于空时簇处理的水声通信方法及相关设备 |
| CN115085838A (zh) * | 2022-04-22 | 2022-09-20 | 哈尔滨工程大学 | 基于虚拟发射阵列的水声探测通信一体化波形验证方法 |
| CN115085838B (zh) * | 2022-04-22 | 2023-09-29 | 哈尔滨工程大学 | 基于虚拟发射阵列的水声探测通信一体化波形验证方法 |
| WO2024032562A1 (zh) * | 2022-08-09 | 2024-02-15 | 华为技术有限公司 | 一种通信方法及装置 |
| CN117081894A (zh) * | 2023-07-17 | 2023-11-17 | 中国科学院声学研究所 | 一种利用信道稀疏特性的水声信号检测方法 |
| CN117081894B (zh) * | 2023-07-17 | 2024-05-24 | 中国科学院声学研究所 | 一种利用信道稀疏特性的水声信号检测方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107505596B (zh) | 2020-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107505596A (zh) | 基于双扩展水声信道环境下的mimo主动探测信号设计与检测系统和方法 | |
| CN104297734B (zh) | 基于频率分集阵列的mimo雷达的欺骗式干扰抑制方法 | |
| CN102608609B (zh) | 一种基于mimo的高分辨水下目标探测方法 | |
| Hassanien et al. | Transmit/receive beamforming for MIMO radar with colocated antennas | |
| JP5631763B2 (ja) | Ｍｉｍｏレーダシステム、送信装置、受信装置及びｍｉｍｏレーダ信号処理方法 | |
| CN106338723A (zh) | 一种基于互质脉冲重复间隔的空时自适应处理方法及装置 | |
| Zhu et al. | Target detection performance analysis of FDA-MIMO radar | |
| JP2007507691A (ja) | ソナー・システムおよびプロセス | |
| CN109521430B (zh) | 一种抑制窄带干扰的距离扩展目标检测方法 | |
| CN102183435A (zh) | 一种基于多路径反射理论的海底密度和声速测量方法 | |
| CN107272005A (zh) | 基于可靠声路径下目标回波到达时延和到达角度的主动定位方法 | |
| Li et al. | Super-resolution time delay estimation for narrowband signal | |
| CN101387701A (zh) | 一种基于前向预测的被动时间反转混响抑制方法 | |
| CN103713289B (zh) | 基于分布式Phased-MIMO联合处理的目标探测方法 | |
| CN105181800A (zh) | 基于多通道空时逆滤波技术的声学覆盖层回声降低测量方法 | |
| CN105301580A (zh) | 一种基于分裂阵互谱相位差方差加权的被动探测方法 | |
| Gong et al. | A machine learning-based approach for auto-detection and localization of targets in underwater acoustic array networks | |
| CN102967852B (zh) | 基于dps序列的mimo-oth雷达波形生成方法 | |
| US8400875B2 (en) | Active sonar system and active sonar method using a pulse sorting transform | |
| CN109283497A (zh) | 双基地fda-mimo雷达距离欺骗干扰识别方法 | |
| CN106324602A (zh) | 一种mimo声呐系统 | |
| US8116169B2 (en) | Active sonar system and active sonar method using noise reduction techniques and advanced signal processing techniques | |
| CN108896974A (zh) | 一种改进的mimo阵列高分辨空间谱估计方法 | |
| EP2317335A1 (en) | Improved beamforming method for analysing signals received by a transducer arrray, and relative detection system | |
| CN101470201B (zh) | 一种多发多收网络干扰抑制系统及干扰抑制方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |