一种基于OFDM的水声通信系统脉冲噪声抑制方法技术方案

本发明专利技术公开了一种基于OFDM的水声通信系统脉冲噪声抑制方法，其构建空子载波矩阵，并利用OFDM符号中的各个子载波之间的正交性，从接收信号中提取出仅包含脉冲噪声和背景噪声的信号，通过这种方式不仅能准确分离出接收信号中的噪声部分，而且计算复杂度低；基于压缩感知理论，通过引入L1/2范数迭代求解，首先将脉冲噪声估计的NP难问题转化为对脉冲噪声L1/2范数最小值的正则化问题，然后将问题转化为含两步迭代结果的形式，并利用迭代求解方法进行求解，使得本发明专利技术方法在不同空子载波个数和信噪比条件下具有较小的误码率；通过采用阈值函数对迭代结果进行处理，使得本发明专利技术方法在未受到脉冲噪声干扰的采样点处具有较小的均方误差。

【技术实现步骤摘要】
一种基于OFDM的水声通信系统脉冲噪声抑制方法
本专利技术涉及一种脉冲噪声抑制技术，尤其是涉及一种基于OFDM(OrthogonalFrequencyDivisionMultiplexing，多载波正交频分复用)的水声通信系统脉冲噪声抑制方法。
技术介绍
水声通信(UnderwaterAcousticCommunication，UAC)是一种利用声波将信源信息经过水声信道高速、可靠地传输到信宿的通信方式。水声通信是海洋水下通信的主要手段，其与传统的无线通信相比有很大不同，水声通信的带宽只有几十kHz，与无线通信的几十GHz带宽相比，可用带宽资源十分有限。在远海中，水声通信中的噪声一般可以近似为高斯白噪声；但是在近海区域，除了来自水声通信系统的高斯白噪声外，还有来自周围环境的脉冲噪声，脉冲噪声的自然来源有潮汐、洋流、海面波浪、热噪声、海洋生物活动以及地震等，脉冲噪声的人为来源由航运、航道爆破、石油和天然气勘探与生产以及声纳设备运行等产生。由于脉冲噪声的功率谱密度要远远高于背景噪声的功率谱密度，因此脉冲噪声是影响水声通信系统性能的主要因素，抑制脉冲噪声对于水声通信中信号的可靠传输具有重要意义。在水声通信系统中，由于受水底水面反射、折射和海洋环境的漫射效应等影响，水声信道的多径效应明显。多径效应会引起频率选择性衰落，造成码间串扰。与陆上的无线通信相比，声速比光速低很多个数量级，因此水声通信中的水声信道的多径效应引起的码间串扰对水声通信系统性能的影响更加严重。为了对抗水声信道的多径效应，多载波正交频分复用(OrthogonalFrequencyDivisionMultiplexing，OFDM)技术在水声通信中已被广泛采用。OFDM技术能够有效地应对频率选择性衰落信道，因此对脉冲噪声的敏感度要比单载波小很多，但是，实测表明，当脉冲噪声的能量超过某个门限值时，例如比背景噪声高10～20dB左右时，传统的OFDM接收机由于脉冲噪声的出现性能将急剧恶化，信号衰落严重，因此必须采用相应的技术来抑制脉冲噪声。目前，常见的水声通信系统脉冲噪声抑制方法可以分为两类，一类是参数辅助法，另一类是非参数辅助法。参数辅助法主要有限幅法和消隐法，此类方法针对随机到达、持续时间短、功率高的脉冲噪声，处理时可以采用非线性操作法，通过设定一个门限值来判断接收到的采样信号是否被脉冲噪声干扰，如果被干扰则将该采样信号置为0(blanking)或将该采样信号的幅值限制为门限值(clipping)；此类方法虽然实现简单，但是通常难以选择适当的非线性操作门限值来使其足以抑制脉冲噪声而不会使OFDM信号失真，导致误码率和均方误差较大。非参数辅助法利用脉冲噪声随机到达、持续时间短、在时域上可视为稀疏信号这一特性，将压缩感知(CompressedSensing，CS)技术用于脉冲噪声抑制，已被证明优于传统的限幅法和消隐法，其中比较常见的非参数辅助法有正交匹配追踪(OrthogonalMatchingPursuit，OMP)法和最小二乘法(LeastSquares，LS)，这两种方法可以较为准确地估计出脉冲噪声信号，从而在接收信号中减去脉冲噪声的估计值，完成对脉冲噪声的抑制，但是这两种方法所适用的脉冲噪声模型较为单一，在更加接近实际情况的混合高斯模型下，估计性能较差。
技术实现思路
本专利技术所要解决的技术问题是提供一种基于OFDM的水声通信系统脉冲噪声抑制方法，其计算复杂度低，且在更加接近实际情况的混合高斯模型下，具有较小的误码率和均方误差。本专利技术解决上述技术问题所采用的技术方案为：一种基于OFDM的水声通信系统脉冲噪声抑制方法，其特征在于包括以下步骤：步骤一：在基于OFDM的水声通信系统的发送端，将发送端发出的初始二进制数据序列记为d，将d经过编码器编码和正交相移键控映射得到的频域OFDM信号记为D，将D通过离散傅里叶反变换转换得到的对应的离散时域信号记为X，X＝FHD＝[x1,x2,…,xN]T；然后在X中的OFDM符号的前面插入长度大于水声信道的最大时延扩展的循环前缀，得到加有循环前缀的离散时域信号；再将加有循环前缀的离散时域信号通过水声信道传输给基于OFDM的水声通信系统的接收端；其中，X的维数为N×1，F表示维数为N×N的离散傅里叶变换范德蒙德矩阵，FH为F的厄米特变换，D的维数为N×1，符号“[]”为向量表示符号，[x1,x2,…,xN]T为[x1,x2,…,xN]的转置，x1,x2,…,xN对应表示X中的OFDM符号中的第1个子载波携带的信号元素、第2个子载波携带的信号元素、…、第N个子载波携带的信号元素，N表示X中的OFDM符号中的子载波的总个数；步骤二：在基于OFDM的水声通信系统的接收端，将接收端接收到的带有脉冲噪声信号的加有循环前缀的离散时域信号中的OFDM符号的前面的循环前缀去掉，得到带有脉冲噪声干扰的离散时域信号，记为r，r＝E(ε)HX+i+n；其中，r的维数为N×1，ε表示基于OFDM的水声通信系统的发送端与接收端相对移动产生的多普勒频移，E(ε)表示维数为N×N的多普勒频移矩阵，E(ε)＝diag(1,ej2πε,…,ej2πε(N-1))，diag()为对角矩阵表示，e表示自然基数，j为虚数单位，H表示水声信道，H为一个维数为N×N的信道循环卷积矩阵，i表示服从混合高斯分布的脉冲噪声信号，i的维数为N×1，i的概率密度函数为1≤s≤S，S表示i的稀疏度，也即表示i中包含的分量的总个数，ps表示i中的第s个分量出现的概率，且表示具有均值为0、方差为的高斯分布，表示i中的第s个分量的方差，n表示均值为0、方差为σn2的高斯白噪声信号，n的维数为N×1；步骤三：构造一个维数为K×N的空子载波矩阵，记为Φ，Φ由F中的第N-K行至第N行构成；然后根据OFDM符号中的各个子载波之间的正交性，利用Φ将r＝E(ε)HX+i+n转换成Φr＝ΦE(ε)HX+Φi+Φn＝Φi+Φn；接着令nΦ＝Φn，将Φr＝Φi+Φn改写成Φr＝Φi+nΦ；再令y＝Φr＝Φi+nΦ，根据y＝Φr＝Φi+nΦ构造得到求解i的原问题，描述为：约束条件为||y-Φi||2≤ζ；其中，K表示r中的OFDM符号中的空子载波的总个数，1<K<N，nΦ表示水声信道的背景噪声，y为引入的中间变量，y的维数为K×1，表示i的估计值，的维数为N×1，符号“||||0”为求矩阵的L0范数符号，表示使||i||0取最小值时i的值，符号“||||2”为求矩阵的L2范数符号，ζ为与n的方差σn2相关的一个正数；步骤四：引入i的Lq范数，将求解i的原问题转化为正则化问题，描述为：然后对使用一阶最优性条件，得到其中，q＝1/2，符号“||||1/2”为求矩阵的L1/2范数符号，λ为正则化参数，表示取使得的值最小时的i的值，I是维数为N×N的单位矩阵，表示惩罚因子的梯度算子，为的逆，ΦT为Φ的转置；步骤五：利用迭代法对进行求解，得到的最优值，具体过程为：1)、令p表示迭代的次数，p的初始值为1；令θp表示第p次迭代的判决残差；令表示设定的判决阈值；2)、在第p次迭代时，将改写成其中，表示第p次迭代后得到的的值，的维数为N×1，当p≠1时表示第p-1次迭代后得到的本文档来自技高网...

【技术保护点】
1.一种基于OFDM的水声通信系统脉冲噪声抑制方法，其特征在于包括以下步骤：步骤一：在基于OFDM的水声通信系统的发送端，将发送端发出的初始二进制数据序列记为d，将d经过编码器编码和正交相移键控映射得到的频域OFDM信号记为D，将D通过离散傅里叶反变换转换得到的对应的离散时域信号记为X，X＝FHD＝[x1,x2,…,xN]T；然后在X中的OFDM符号的前面插入长度大于水声信道的最大时延扩展的循环前缀，得到加有循环前缀的离散时域信号；再将加有循环前缀的离散时域信号通过水声信道传输给基于OFDM的水声通信系统的接收端；其中，X的维数为N×1，F表示维数为N×N的离散傅里叶变换范德蒙德矩阵，FH为F的厄米特变换，D的维数为N×1，符号“[]”为向量表示符号，[x1,x2,…,xN]T为[x1,x2,…,xN]的转置，x1,x2,…,xN对应表示X中的OFDM符号中的第1个子载波携带的信号元素、第2个子载波携带的信号元素、…、第N个子载波携带的信号元素，N表示X中的OFDM符号中的子载波的总个数；步骤二：在基于OFDM的水声通信系统的接收端，将接收端接收到的带有脉冲噪声信号的加有循环前缀的离散时域信号中的OFDM符号的前面的循环前缀去掉，得到带有脉冲噪声干扰的离散时域信号，记为r，r＝E(ε)HX+i+n；其中，r的维数为N×1，ε表示基于OFDM的水声通信系统的发送端与接收端相对移动产生的多普勒频移，E(ε)表示维数为N×N的多普勒频移矩阵，E(ε)＝diag(1,e...

【技术特征摘要】
1.一种基于OFDM的水声通信系统脉冲噪声抑制方法，其特征在于包括以下步骤：步骤一：在基于OFDM的水声通信系统的发送端，将发送端发出的初始二进制数据序列记为d，将d经过编码器编码和正交相移键控映射得到的频域OFDM信号记为D，将D通过离散傅里叶反变换转换得到的对应的离散时域信号记为X，X＝FHD＝[x1,x2,…,xN]T；然后在X中的OFDM符号的前面插入长度大于水声信道的最大时延扩展的循环前缀，得到加有循环前缀的离散时域信号；再将加有循环前缀的离散时域信号通过水声信道传输给基于OFDM的水声通信系统的接收端；其中，X的维数为N×1，F表示维数为N×N的离散傅里叶变换范德蒙德矩阵，FH为F的厄米特变换，D的维数为N×1，符号“[]”为向量表示符号，[x1,x2,…,xN]T为[x1,x2,…,xN]的转置，x1,x2,…,xN对应表示X中的OFDM符号中的第1个子载波携带的信号元素、第2个子载波携带的信号元素、…、第N个子载波携带的信号元素，N表示X中的OFDM符号中的子载波的总个数；步骤二：在基于OFDM的水声通信系统的接收端，将接收端接收到的带有脉冲噪声信号的加有循环前缀的离散时域信号中的OFDM符号的前面的循环前缀去掉，得到带有脉冲噪声干扰的离散时域信号，记为r，r＝E(ε)HX+i+n；其中，r的维数为N×1，ε表示基于OFDM的水声通信系统的发送端与接收端相对移动产生的多普勒频移，E(ε)表示维数为N×N的多普勒频移矩阵，E(ε)＝diag(1,ej2πε,…,ej2πε(N-1))，diag()为对角矩阵表示，e表示自然基数，j为虚数单位，H表示水声信道，H为一个维数为N×N的信道循环卷积矩阵，i表示服从混合高斯分布的脉冲噪声信号，i的维数为N×1，i的概率密度函数为S表示i的稀疏度，也即表示i中包含的分量的总个数，ps表示i中的第s个分量出现的概率，且表示具有均值为0、方差为的高斯分布，表示i中的第s个分量的方差，n表示均值为0、方差为σn2的高斯白噪声信号，n的维数为N×1；步骤三：构造一个维数为K×N的空子载波矩阵，记为Φ，Φ由F中的第N-K行至第N行构成；然后根据OFDM符号中的各个子载...

【专利技术属性】
技术研发人员：李娜娜，李有明，余明宸，卢倩倩，
申请(专利权)人：宁波大学，
类型：发明
国别省市：浙江,33