【水声信号分析】水声信道中通信信号的循环平稳特性【附MATLAB代码】

文章来源：微信公众号 EW Frontier

摘要：对水声传播对通信信号循环平稳特性的影响进行了建模和分析。两种信道被认为是：多尺度多滞后信道，在其中移动的和宽带声学系统通常进行通信，和色散信道导致低频模态传播在浅水中。它示出了多尺度多滞后信道变换循环平稳信号的速度和加速度依赖的时间弯曲的循环平稳过程的总和。这种时间扭曲被仔细考虑，以有效地恢复循环平稳特征。另一方面，它被发现，低频色散信道保持原来的周期性，但衰减较短的周期和传播的相关性。为了说明理论结果，模拟和真实的数据的应用也被提出。具体地说，估计时变多普勒尺度的问题是解决多尺度多滞后信道以及检测独特的循环平稳特性的信号。最后给出了盲码速率估计在色散信道隐蔽通信中的应用实例。特别注意相移键控（PSK），正交幅度调制（QAM），正交频分复用（OFDM），和直接序列扩频（DSSS）信号。附带的补充材料提供了用于估计和检测示例的MATLAB代码。

声学通信，循环平稳性（CS），模态色散，多尺度多滞后信道，水下作战。

为了促进在UAC环境中使用基于CS的方法并说明理论结果，还提供了模拟和真实的数据的应用。在多尺度多滞后MSML信道中进行了两个实验，一个在色散信道中进行。

1)在第一个实验中，考虑了强加速度引起的时变多普勒尺度的估计。本文提出了一种在多径信道中估计多个多普勒尺度的方法，并在模拟的相移键控（PSK）/正交幅度调制（QAM）信号上进行了评估。该方法利用了由于多普勒缩放而变得时间扭曲的CS特征。更具体地，时间的二次函数被用来对运动引起的多普勒进行建模，并且基于自相关特性导出成本函数。然后，该代价被迭代优化以估计与每条路径的速度和加速度相关的多项式系数。它不影响数据速率，适用于非合作场景。

2)第二应用处理检测具有特定CS签名的UAC信号的问题。这样的签名可能是唯一的，因此与检测水下战争类应用中的给定通信系统相关。一个统计测试，明确考虑随时间变化的多普勒尺度的推导。该测试提供了合理的复杂性，其性能与真实的正交频分复用（OFDM）的浅水MSML信道中收集的数据说明。它被证明提供更好的性能比现成的RF方法。

3)最后一个实验是专门的问题，盲符号率估计适用于直接序列扩频（DSSS）传输在色散信道。这样一个问题可能是关键的重要性窃听低频隐蔽communications.As相对于MSML信道，分散的信道不扭曲CS功能，使估计已经在文献中需要很少的适应。提出了一种估计器，并利用真实的低频单输入多输出（SIMO）浅水脉冲响应进行了重放模拟。结果表明，沿阵列的空间相干性沿着的结果在正确的估计率显着提高。

大多数通信信号是循环平稳的，其周期等于符号周期Ts。这在图1中用水下声学通信中常用的四种类型的信号的循环自相关函数来说明：PSK/QAM [14]，[15]（左上），短码DSSS [16]，[17]，[18]（右上），循环前缀OFDM（CP-OFDM）[19]，[20]，[21]（左下）和具有叠加导频的无保护OFDM（SIP-OFDM）[22]（右下）。每个信号都有自己的CS模式，这对盲调制识别很有用，每个图案都揭示了这个符号周期Ts的价值。如附录A中进一步讨论和说明的，由于重复使用相同的脉冲整形滤波器、相同的扩频或导频序列，或者由于循环前缀引起的冗余，这样的信号表现出周期性的二阶统计。其他类型的UAC波形，如MFSK，跳频FSK（如JANUS标准[23]中使用的）或啁啾扩频也是二阶CS [24]，[25]。虽然本文没有讨论，但下一节中介绍的应用示例可以很容易地转换为这些非相干调制方案。最后，除了数据符号之外，UAC信号通常包括信令模式，例如可生成附加CS特征的伪码、后同步码或时频导频[26]、[27]。在这种情况下，循环频率可能是不相称的，并且信号被称为几乎是CS [4]。为了不依赖于特定的帧结构，本文其余部分介绍的方法仅利用了符号率倍数的循环频率提供的信息。

图 1.循环自相关函数的例子。(a)带根升余弦滤波器的PSK/QAM。(b)具有7码片长扩频序列和根升余弦滤波器的短码DSSS。(c)具有根升余弦窗的CP-OFDM。(d)具有矩形窗口的SIP-OFDM。

图7。在用Bellhop模拟的MSML信道中获得的成本函数Jr（j）的演变的示例。(a)迭代j =0，（B）迭代j = 1。成本表示为基于模型（8）的速度和加速度的函数，其中c = 1500 m.s-1，即，v = 1500（1 − μ1），a = −3000μ2。Tobs =1 s，Ts =0.2 ms，fc =20 kHz，Eb/N0 =5 dB。

见图9。接收到的OFDM数据包的频谱图示例。(a)带有音调噪声的高SNR数据包。(b)具有脉冲噪声的中等SNR数据包。(c)低SNR数据包。

图19.具有均方根升余弦滤波器g的PSK/QAM信号的循环平稳特征的图示，其中g =0.5。(a)自相关函数，（B）循环自相关函数，（c）谱相关函数，和（d）循环谱。

图20.具有Ts = Ts/4、N =32和矩形窗口的ZP-OFDM信号的循环平稳特征的图示。(a)自相关函数，（B）循环自相关函数，（c）谱相关函数，和（d）循环谱。

图21。具有最大长度扩展序列的DSSS信号的循环平稳特征的图示，Nc =7，根升余弦滤波器g，其中k =0.5。(a)自相关函数，（B）循环自相关函数，和（c）谱相关函数，（d）循环谱。