CN107769862A - 一种仿生低截获通信方法 - Google Patents
一种仿生低截获通信方法 Download PDFInfo
- Publication number
- CN107769862A CN107769862A CN201710979305.5A CN201710979305A CN107769862A CN 107769862 A CN107769862 A CN 107769862A CN 201710979305 A CN201710979305 A CN 201710979305A CN 107769862 A CN107769862 A CN 107769862A
- Authority
- CN
- China
- Prior art keywords
- low
- signal
- interception
- sonar
- communication interception
- 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.)
- Pending
Links
- 238000004891 communication Methods 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 28
- 241001481833 Coryphaena hippurus Species 0.000 claims abstract description 34
- 230000000739 chaotic effect Effects 0.000 claims abstract description 18
- 238000013461 design Methods 0.000 claims abstract description 18
- 150000001875 compounds Chemical class 0.000 claims abstract description 11
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 10
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 10
- 230000005540 biological transmission Effects 0.000 claims description 15
- 230000000694 effects Effects 0.000 claims description 15
- 238000013507 mapping Methods 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 5
- 238000005094 computer simulation Methods 0.000 claims description 3
- 230000001427 coherent effect Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 11
- 230000014509 gene expression Effects 0.000 description 9
- 238000001514 detection method Methods 0.000 description 8
- 239000002131 composite material Substances 0.000 description 5
- 238000005314 correlation function Methods 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 5
- 238000005311 autocorrelation function Methods 0.000 description 3
- 108010003272 Hyaluronate lyase Proteins 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241000283153 Cetacea Species 0.000 description 1
- 206010019133 Hangover Diseases 0.000 description 1
- 206010044565 Tremor Diseases 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000000556 factor analysis Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012421 spiking Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04K—SECRET COMMUNICATION; JAMMING OF COMMUNICATION
- H04K3/00—Jamming of communication; Counter-measures
- H04K3/80—Jamming or countermeasure characterized by its function
- H04K3/82—Jamming or countermeasure characterized by its function related to preventing surveillance, interception or detection
- H04K3/827—Jamming or countermeasure characterized by its function related to preventing surveillance, interception or detection using characteristics of target signal or of transmission, e.g. using direct sequence spread spectrum or fast frequency hopping
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04K—SECRET COMMUNICATION; JAMMING OF COMMUNICATION
- H04K3/00—Jamming of communication; Counter-measures
- H04K3/40—Jamming having variable characteristics
- H04K3/44—Jamming having variable characteristics characterized by the control of the jamming waveform or modulation type
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B13/00—Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
- H04B13/02—Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
- H04B17/3912—Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04K—SECRET COMMUNICATION; JAMMING OF COMMUNICATION
- H04K3/00—Jamming of communication; Counter-measures
- H04K3/40—Jamming having variable characteristics
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04K—SECRET COMMUNICATION; JAMMING OF COMMUNICATION
- H04K2203/00—Jamming of communication; Countermeasures
- H04K2203/10—Jamming or countermeasure used for a particular application
- H04K2203/12—Jamming or countermeasure used for a particular application for acoustic communication
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Electromagnetism (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
一种仿生低截获通信方法，涉及一种水声对抗。主动声呐波形设计；Alpha稳态对水下环境噪声建模；混沌序列的低截获通信；海豚信号合成的低截获通信。依据低截获声呐的低截获因子和模糊函数，采用低频信号、大时宽带宽积、复合频(码)制、随机或非线性体制、时频捷变以及低功率高能量效率六个原则，设计低截获声呐波形。利用混沌映射m序列控制法产生随机相位和混沌映射的基函数的随机性及相关性，实现混沌序列的低截获通信。计算出whistles脉冲串的包络，高斯函数作为基函数，拟合海豚脉冲，重构海豚信号。
Description
技术领域
本发明涉及一种水声对抗，尤其是涉及一种仿生低截获通信方法。
背景技术
随着水声对抗技术的需求，潜艇隐身技术得到进一步的发展。隐身雷达中的低截获(Low Probability of Intercept，LPI)技术被引入到声呐界。1983年，伦敦大学的J.R.Forest首次提出了LPI雷达方程。目前，LPI技术在雷达领域中已用于实际装备，而声呐中的LPI技术还在起步阶段。在声呐界，广义的低截获性能不仅是指不容易被截获，还表示即使信号被截获也不容易被识别。
从上世纪90年代起，LPI技术开始被逐步引入主动声呐研究领域。然而水声环境比空中环境要更加恶劣。水下环境中存在很多无规则散射体，例如：海洋生物、泥沙、气团、水泡等，这些散射体使水下信道存在严重的多经效应；海面和海底等对入射信号的散射造成主动声呐的接收端信号重叠形成混响；这些原因使常规信号在水下环境中衰减严重，造成LPI声呐发展较为缓慢。
然而，现在较为常用的主动声呐发射信号主要有CF信号、伪随机码、LFM信号等，在不同的水声条件下，不同的信号有各自的优点，但是这些信号形式较为简单，抗截获和解调能力不强，这造成声呐通信系统收发并不理想。
发明内容
本发明的目的在于提供一种仿生低截获通信方法。
本发明包括以下步骤：
1)主动声呐波形设计；
在步骤1)中，所述主动声呐波形设计的方法可为：依据主动声呐实现低截获通信的原理，推导低截获概率因子的数学表达式，根据低截获因子和模糊函数提出低截获主动声呐波形设计原则。
2)Alpha稳态对水下环境噪声建模；
在步骤2)中，所述Alpha稳态对水下环境噪声建模的具体方法可为：在进行波形设计时，Alpha稳态对海洋环境噪声进行建模，确立计算机仿真的信道模型；分析海洋环境的多普勒效应、多径效应、传播损耗以及海洋环境噪声，产生Alpha稳定分布及概率密度。
3)混沌序列的低截获通信；
在步骤3)中，所述混沌序列的低截获通信的具体方法可为：用混沌信号作为低截获声呐通信的载波信号，实现混沌键控和混沌调相的低截获通信；低截获声呐波形设计，分析低截获声呐的低截获因子和模糊函数，并提出低截获声呐波形设计的6个原则：低频信号、大时宽带宽积、复合频(码)制、随机或非线性体制、时频捷变以及低功率高能量效率；根据三种典型的低截获声呐信号，提出一种形式较为简单，在特定的环境中低截获性能较好的信号作为低截获声呐的发射信号；混沌键控和混沌调相的低截获通信系统，采用DCSK调制，将要发送的信息比特隐藏在两段信号的相位差内，克服CSK在非相干解调时存在门限漂移的问题，采用改进型Logistic映射和等间隔的方式来生成初值序列，产生个复值随机相位，可以将产生的混沌序列依次移位进入一个n级移位寄存器，选取n级移位寄存器中的ri个抽头构成随机相位映射器，混沌序列中的相位长度为基函数采样点数N。
4)海豚信号合成的低截获通信。
在步骤4)中，所述海豚信号合成的低截获通信的具体方法可为：实现高斯函数的导函数来对海豚信号进行建模，将实际通信信息伪装成海豚信号，以达到信息隐藏的目的；海豚信号合成的低截获通信，计算出whistles脉冲串的包络，就可以将海豚脉冲进行拟合，采用高斯信号及其各阶导函数模拟海豚信号，通过高斯函数及各阶导数来进行海豚信号的重构，可以将信息隐藏在高斯函数的各阶导数中，这样可以大大提高了信息的保密能力。
本发明依据低截获声呐的低截获因子和模糊函数，采用低频信号、大时宽带宽积、复合频(码)制、随机或非线性体制、时频捷变以及低功率高能量效率六个原则，设计低截获声呐波形。利用混沌映射m序列控制法产生随机相位和混沌映射的基函数的随机性及相关性，实现混沌序列的低截获通信。计算出whistles脉冲串的包络，高斯函数作为基函数，拟合海豚脉冲，重构海豚信号。
本发明设计了基于主动声呐的低截获通信系统发射波形；根据水下环境复杂多变的情况，实现了信道仿真，用Alpha稳态对水下环境噪声建模；实现混沌序列的低截获通信。实现海豚信号合成的低截获通信。
本发明具有以下有益效果：
本发明以海洋中鲸鱼信号作为信息载体，以合成信号作为传输运载信号，提高了仿生低截获通信的性能，改善仿生隐蔽通信。将信息隐藏在高斯函数的各阶导数中，并进行有效合成，可以大大提高了信息的保密能力，并实现发射的仿生信号可以包含多种载波、调制样式以及不同的带宽范围的基带信号，很大程度上提高信息传输量。
附图说明
图1为本发明的仿生低截获通信方法框图；
图2为本发明的信号的多普勒效应图；
图3为本发明的信号的多径效应图；
图4为本发明的低截获主动声呐；
图5为本发明中的DCSK调制器；
图6为本发明中的DCSK解调器；
图7为本发明中基于混沌映射m序列控制法的随机相位产生器图；
图8为本发明中高斯拟合信号合成流程图。
具体实施方式
以下实施例将结合附图对本发明作进一步的说明。
1、仿生低截获通信方法
请参考图1，本发明主要包括主动声呐波形的设计、海洋环境信道仿真和Alpha稳态对海洋环境噪声建模、混沌序列的低截获通信系统和海豚信号合成的低截获通信四大模块。
考虑水下多普勒效应、多径效应、传播损耗及海洋环境噪声因素，产生Alpha稳定分布的产生及概率密度，用Alpha稳态建立海洋环境噪声模型，计算机仿真的信道模型。
设计低频信号技术、大时宽带宽积技术、复合频(码)制、随机或非线性体制、时频捷变和低功率高能量效率六大原则，依据LFM+CF信号、LFM+Barker码信号和HFM信号，计算低截获因子和模糊函数，设计并实现了低截获主动声呐波形。
依据混沌通信理论，通过混沌Logistic映射产生混沌序列。实现混沌键控(DCSK)的低截获通信，DCSK调制中，每个比特时间被分为两个时隙，第一个时隙来发送一段混沌序列作为参考信号，第二个时隙发送经信息信号调制的混沌信号，若发送的信息比特为“+1”，则两个时隙的信号是相同的混沌信号；若发送的信息比特为“-1”，则信息段和参考信号段相位相反。换言之，DCSK将要发送的信息比特隐藏在两段信号的相位差内。实现混沌调相的低截获通信，利用改进型改进型Logistic映射，产生个复值随机相位，可以将产生的混沌序列依次移位进入一个n级移位寄存器，选取n级移位寄存器中的ri个抽头构成随机相位映射器，混沌序列中的相位长度为基函数采样点数N。混沌映射m序列控制法中混沌序列与m序列的移位速率不同。混沌序列每移位N次，m序列移位一次。所以每一个TDCS数据符号由m序列寄存器的取值和混沌序列的N次相位映射所决定。
采用高斯信号及其各阶导函数模拟海豚信号，实现海豚信号合成的低截获通信。使用N个高斯导函数的线性组合，高斯每个导函数认为是N维空间中的基函数，提取高斯信号包络，求海豚信号极值点，提取两个极值点之间的波形，产生N个高斯函数的导函数，与极值之间波形进行相关并取最值，得到N个导数的线性组合给海豚信号波形。
本发明将信息隐藏在海豚信号中，解决了直接将信息调制在海豚信号中的低速率问题，而且信息隐藏更为隐蔽，提高了低截获通信的保密能力，达到了低截获通信的目的。
2、海洋环境信道仿真和噪声建模
请参考图2和3，本发明中海洋信道仿真，考虑了海洋环境的多普勒效应、多径效应、传播损耗和海洋环境噪声，并采用Alpha稳定分布对海洋环境噪声建模，保证系统通信质量。
(1)多普勒效应
在频率域，信道的多普勒效应是由于信道的时变特性所引起的频谱扩展。当发送端和接收端固定时，唯一的时变特性就是由于海洋表面的风浪与对应引起的海洋表面波浪运动。对于水下通信系统而言，带宽的起伏度表示为：
其中，w是风速，单位为m/s，fw＝2/w是波浪频率，单位为Hz，hw＝0.005×w/2是波浪长度，单位为m，f0是载波频率，单位为Hz，θ是随机相位角度，单位为度，c是声速，单位为m/s。
(2)多径效应
多径特性使得原始发送信号的延迟混响在不同的时间时刻达到接收端。延迟主要是由于海面与海底的反射造成的。多径效应造成了接收信号强度与响应的浮动。多径信道可以表示为：
其中，wi是每一径所对应的幅度，对应的第i条作用于δ的kronecker函数。
(3)传播损耗
传播损耗主要分为两部分，衰减损耗和扩散损耗：
扩散损耗：
TLcylind＝10·log10(r)
衰减损耗：
TLatten＝αr×10-3
其中，α为衰减系数，
r为传播距离，单位为m，f为频率，单位为Hz。
水声信道划分为三个频段：低频段、中频段和高频段。低频段在0～15KHz之间，较适用于远距离通信；中频段在15～150KHz之间，较适用于中距离通信；高频段在150～1.5MHz之间，主要用于告诉的近距离通信。由于传播损耗的存在，水声通信实验表明，水声通信的距离与传输速率的乘积上限可表示为：
通信距离(s)·通信速率(v)＝40Km·Kbps
(4)海洋环境噪声建模
自α稳定分布产生的80年来，α稳定分布理论得到了极大的重视和发展，并形成了较为完善的理论体系。通过α稳定分布进行海洋环境建模，Alpha稳定分布相关数学定义如下：
定义1：若对于任意正数A和B，存在着正数C和一个实数D，满足
则称随机变量X是稳定分布的。式中的X1和X2为X的独立样本，符号表示分布相同。若当D＝0时上式成立，则称X为严格稳定的。若X和-X具有的分布相同，则称X为对称稳定的，可以看出，对称稳定分布总是满足严格稳定分布的。
若X1,X2,…Xn是相互独立的随机变量，均服从稳定分布，且具有相同的参数α和β，则它们的线性组合也一定服从稳定分布，且具有相同的参数α和β。
定理1：对于任何随机变量X，存在α∈(0,2]，使C满足
Cα＝Aα+Bα
定义2：设X1,X2,…Xn为随机变量X的独立样本，对于任意n＞1，若存在一个整数Cn和一个实数Dn，满足
则认为变量X具有稳定分布，并且定义1与定义2等价。
本发明主要依据标准参数系S，获得Alpha稳定分布特征函数的表征。
若随机变量X存在四个参数α,β,σ,μ，其中，0＜α≤2,-1≤β≤1,σ,μ≥0为实数，其特征函数如下式所示：
式中
则称随机变量t服从Alpha稳定分布，记为X～S(α,β,σ,μ)。式中，称参数α为特征指数，它决定该分布脉冲特性的程度。α值越小，拖尾越厚，脉冲特性越显著。相反，α值越大，拖尾越薄，脉冲特性减弱。称参数β为对称参数，主要用于确定分布的斜度。β＝0对应对称分布，简称SαS。0＜β≤1对应于右偏斜分布，-1≤β＜0对应于左偏斜分布。称参数α为分散系数，又称为尺度参数，它用来度量样本相对均值的分散程度。称参数μ为位置参数。若μ,σ满足μ＝0,σ＝1，则称Alpha稳定分布为标准Alpha稳定分布。
Alpha稳定分布为具有尖峰脉冲特性的噪声过程提供了行之有效的模型。利用Alpha稳定分布来代替传统的高斯分布，这样在系统海洋环境建模仿真时，可以使仿真结果更加具有实用价值，设计出的信号处理算法有较好的鲁棒性。
3、低截获声呐波形
请参考图4，计算低截获声呐的低截获因子和模糊函数，设计低截获声呐波形的6个原则：低频信号、大时宽带宽积、复合频(码)制、随机或非线性体制、时频捷变以及低功率高能量效率。在特特定的环境中低截获性能较好的信号来作为低截获声呐的发射信号
(1)声呐截获因子
根据雷达截获概率因子的定义，将低截获标准引入声呐，计算声呐截获概率因子，这里的低截获主要针对信号被截获：
α＝RI/RD
其中，RI是声呐截获接收机能发现发射信号的最大距离，RD是主动声呐能够探测到的最大距离。
根据声呐方程，声呐截获接收机的优质因素写为：
FOMI＝TL(RI)＝SLI-NLI+GSI+GTI-DTI
其中，TL为传播损失，SL为声源级，NL为环境噪声级，GS为空间增益，GT为时间增益，DT为检测阈。
主动声呐接收机的优质因素可以写为：
FOMD＝TL(RD)＝(SLD+TSD-NLD+GSD+GTD-DTD)/2
若截获接收机和主动声呐接收机的检测阈根据低截获的定义满足下列条件：DTI/DTD＝ρ，其中ρ是根据战术需求确定的检测阈比例，可以得到：
TL(RI)-TL(RD)＝10*β*lg(RI/RD)+σ(RI-RD)
＝(SL-NL)/2+(1/2-ρ)DTD-TSD/2+
(2*GSI-GSD+2*GTI-GTD)/2
这里不考虑RI-RD＞＞0或＜＜0的情况(这时低截获性能不明显)，上式可以写为：
αdB＝(A+P-T+G)/β
式中，
αdB＝10*lgα,A＝(SL-NL)/2,P＝(1/2-ρ)DTD,T＝TSD/2,
G＝(2*GSI-GSD+2*GTI-GTD)/2
从声呐的截获因子分析可以得到：当αdB＞0时，表明截获声呐可以探测到主动声呐的存在，而主动声呐无法探测到截获声呐；αdB＜0时，表明截获声呐不能探测到主动声呐，而主动声呐可以探测到截获声呐，具有这样性质的主动声呐成为低截获(LPI)声呐。
(2)模糊函数
声呐模糊函数表示匹配滤波器的输出，描述目标的距离和多普勒频移对回波信号的影响，信号μ(t)的声呐模糊函数通常被定义为二维互相关函数的模的平方|x(τ,ε)|2。具体表达式为：
(3)设计低截获声呐六大原则
①低频信号技术：采用低频信号进行发射波形设计有利于声呐的隐蔽以及反截获；
②大时宽带宽积技术：为了弥补单个脉冲BT积的不足，雷达和国外一些声呐中常采用脉冲串的形式来等效增大时宽积。这类形式的波形的模糊函数是板钉型，频谱呈梳状。它们不仅克服了单个长时宽带脉冲的许多缺点、保持了较高的时间处理增益，并且具有较好的多普勒分辨能力和距离分辨能力。其组合形式较为多样、处理方式比较灵活，是一种比较好的低截获波形形式。
③复合频(码)制：采用复合频制可以避免由于信道对不同频段的影响所造成的性能损失，这里的复合频段指的是频段与频段之间的频率跳变较大的情况。一些简单的信号，跳频信号(FHOP)，BPSK等，都具有复合频(码)制的特点。而正交码、Costas码所进行得复合调频脉冲串(OC-PTFM、Costas-PTFM)，都具有较好的低截获性能。
④随机或非线性体制：在水声低截获通信中，随机特性和混沌特性的编码调制信号具有较优的低截获性能，在使用随机和混沌信号进行信号调制时，必须严格控制发射功率，同时必须考虑水声信道起伏影响。
⑤时频捷变：在波形设计时，可以在一个“短”的时间周期中对信号的各个有效特征按一定的统计规律或者随机的方式进行调制。这对声呐信号的低截获性能的提高是非常有利的。
⑥低功率高能量效率：发射信号能量效率可以定义为待测脉冲信号辐射能量与以待测信号频带为中心的单频信号的比值。不同的信号具有不同的发射效率，发射效率高的波形有利于信号的检测，更容易被截获接收机截获。
(4)低截获声呐波形
组合发射两种信号，并采用波形融合检测的方法检测出声呐信号的特征。组合信号和仿生信号不仅可以更加适应海洋环境，而且由于信号形式较为复杂，使信号的截获概率较低。
①LFM+Barker码组合信号
线性调频信号有较好的抗多普勒性能，将二相码信号与线性调频信号进行组合处理，用于低截获主动声呐系统中。
线性调频信号的表达式为：
二相编码信号的表达式为：
其中，δ(t)为冲激函数，T为子脉冲宽度，P为码长，K为线性调频信号斜率，cm为巴克码，混合脉冲信号为线性调频信号与二相编码脉冲函数的卷积，即：
组合信号的带宽近似为线性调频信号带宽，不再受到巴克码子码宽度的限制，因此可以得到较大的时宽带宽积。
线性调频信号模糊函数表达式为：
LFM+Barker组合信号的模糊函数为：
令t'＝t-mT
②LFM+CF组合信号
LFM+CF组合波形具有有较好的抗频选特性，在一定程度上可以克服浅海声信道中的频率选择性衰落所引起的检测不稳问题，实现较大速度范围内的目标检测。
假设在一个宽度为T的脉内，CF的持续时间为0～kT,(0＜k＜1)，FM的持续时间为kT～T，则恒频信号为：
其模糊函数：
线性调频信号为：
其模糊函数为：
③双曲调频信号
双曲调频(Hyperbolic Modulation Frequency)信号是多普勒不变信号，作为导频信号来计算恶劣水声通信环境中的多普勒偏移，从而实现系统同步。
HFM信号表达式为：
式中，rect(t)＝1，f0为中心频率。
调频率为：
根据HFM的函数表达式，用包络α(t)代替信号的矩形窗，得到模糊函数为：
4、混沌序列的低截获通信
请参考图5～7，根据混沌映射，产生混沌序列。利用混沌信号，实现混沌序列的低截获通信，包括混沌键控和混沌调相的低截获通信。
(1)混沌键控的低截获通信
请参考图5和6，DCSK调制中，DCSK将要发送的信息比特隐藏在两段信号的相位差内。
以双极性二进制比特信息为例，DCSK系统发送端的信号可以表示为：
其中β表示每个时隙内传输信号的采样点数，根据扩频因子(Spreading Factor，SF)的定义，DCSK系统中SF＝2β。
为恢复bk，接收端将接受的信号ri与自身延时β后的信号ri-β进行相关运算。
相关器输出为：
假定传输信号在信道中只受到加性高斯白噪声，则上式中ri可表示为：
ri＝si+ξi
式中ξi服从高斯分布。
接收端得到输出，判决器需要判决的变量Zk：
其中，第一项的符号由当前传输的信息比特的符号决定。若传输的信息比特为“+1”，则第一项为正；若传输的信息比特为“-1”，则第一项为负。其余项数学期望均为零，将判决器的判决门限设置为零，则恢复信号为：
其中，sign[·]表示符号函数。
(2)混沌调相的低截获通信
改进的Logistic映射定义为：
xn+1＝1-2xn 2；-1＜xn＜1
其具有以下性质：
①不变测度
②均值
由于其均值为零，因此改进的Logistic映射又称零均值Logistic映射，序列二值化判决门限ξ＝0。
③自相关函数
④互相关函数
考虑两个混沌序列{x1,k|k＝0,1,2,…}和{x2,k|k＝0,1,2,…}，其互相关函数为：
利用不同的初值产生的混沌序列分配给不同的用户来作为多址接入的地址码。对此，可以采用等间隔的方式来生成初值序列。同时，由于混沌序列无限长的特点，要产生2ri个复值随机相位，可以将产生的混沌序列依次移位进入一个n级移位寄存器，选取n级移位寄存器中的ri个抽头构成随机相位映射器，混沌序列中的相位长度为基函数采样点数N，通过上述混沌映射法来提高基函数相位谱的相关性。
基函数相位的随机性主要和相位映射器的输入r有关，且r是固定的。一旦r给定，在单位圆上一个周期内的随机相位总数就可以固定，而相位的随机性则直接决定了基函数的随机性。这里利用m序列所确定的r抽头来控制混沌序列，通过改变相位映射器的输入r，基函数相位的随机性得以提高。
请参考图7，m为产生的移位寄存器的级数；映射级数ri由m序列的m级移位寄存器的值决定。若m序列前的m个移位寄存器的值为ri＝4，则混沌序列有4级映射抽头，抽头所对应的混沌序列前ri个移位寄存器的值决定相位映射，可选相位有点。
对于给定的r，则可以映射成种基数。若随机相位映射器的输入为r，则产生的基函数个数为：
若混沌序列的可选数量为Nc，m序列的移位寄存器级数为m，则相位映射总数为：
Ns＝(2m-1)Nc＝NmNc
检测概率为1/(2m-1)，即随着m的增加，检测概率呈指数下降，通过分析可以得到，引入混沌序列后，基函数的随机性得到大大提高，这样就增加了信号的保密性，即使信号被截获，也不容易被破译。
两个长度相同的时域基函数bl(n)、其相关函数为：
其中，N为基函数采样点函数，m为基函数之间的时延，Au、是两个基函数频谱的幅值，分别为两个基函数第u和v个子载波相位，Sl和Sq分别表示两个基函数上的调制数据。
当l≠q时，且m＝0时，可以得到基函数的互相关函数最大值为：
当l＝q，m＝0时，可以得到基函数的自相关函数为：
当m＝0时，可以得到基函数的自相关函数最大值为：
5、海豚信号合成的低截获通信
请参考图8，海豚声呐对浅海和浑浊海域的适应能力较强，借鉴海豚通过发射一组滴答脉冲串的方式进行通信或探测，设计仿生低截获通信发射信号。
(1)高斯函数及导函数
通过高斯函数及各阶导数进行海豚信号的拟合，本发明对对高斯函数及其导函数限定：
一个高斯脉冲p(t)可以用下列表达式进行描述：
其中，σ为标准差。
高斯函数的导函数表达式为：
其中n表示求导的阶数。
高斯脉冲的n阶导数的幅度谱为：
设当f＝fM时|xn(f)|取最大值，则fM称为最大峰值发射频率，fM可通过对上式求导并令其为零求出：
(2)高斯拟合过程
使用N个高斯导函数的线性组合，高斯每个导函数认为是N维空间中的基函数。对线性组合的权重则与所拟合信号本身有关。
用数学式表达：
其中，Si为权重系数。gi(t)表示第i阶导数。
F(f)为f(t)的傅立叶变换。
N个导数的线性组合给海豚信号波形的产生带来了很大的灵活性，但是为了使信号得到更好的拟合。
Claims (5)
1.一种仿生低截获通信方法，其特征在于包括以下步骤：
1)主动声呐波形设计；
2)Alpha稳态对水下环境噪声建模；
3)混沌序列的低截获通信；
4)海豚信号合成的低截获通信。
2.如权利要求1所述一种仿生低截获通信方法，其特征在于在步骤1)中，所述主动声呐波形设计的方法为：依据主动声呐实现低截获通信的原理，推导低截获概率因子的数学表达式，根据低截获因子和模糊函数提出低截获主动声呐波形设计原则。
3.如权利要求1所述一种仿生低截获通信方法，其特征在于在步骤2)中，所述Alpha稳态对水下环境噪声建模的具体方法为：在进行波形设计时，Alpha稳态对海洋环境噪声进行建模，确立计算机仿真的信道模型；分析海洋环境的多普勒效应、多径效应、传播损耗以及海洋环境噪声，产生Alpha稳定分布及概率密度。
4.如权利要求1所述一种仿生低截获通信方法，其特征在于在步骤3)中，所述混沌序列的低截获通信的具体方法为：用混沌信号作为低截获声呐通信的载波信号，实现混沌键控和混沌调相的低截获通信；低截获声呐波形设计，分析低截获声呐的低截获因子和模糊函数，并提出低截获声呐波形设计的6个原则：低频信号、大时宽带宽积、复合频(码)制、随机或非线性体制、时频捷变以及低功率高能量效率；根据三种典型的低截获声呐信号，提出一种形式较为简单，在特定的环境中低截获性能较好的信号作为低截获声呐的发射信号；混沌键控和混沌调相的低截获通信系统，采用DCSK调制，将要发送的信息比特隐藏在两段信号的相位差内，克服CSK在非相干解调时存在门限漂移的问题，采用改进型Logistic映射和等间隔的方式生成初值序列，产生个复值随机相位，将产生的混沌序列依次移位进入一个n级移位寄存器，选取n级移位寄存器中的ri个抽头构成随机相位映射器，混沌序列中的相位长度为基函数采样点数N。
5.如权利要求1所述一种仿生低截获通信方法，其特征在于在步骤4)中，所述海豚信号合成的低截获通信的具体方法为：实现高斯函数的导函数来对海豚信号进行建模，将实际通信信息伪装成海豚信号，以达到信息隐藏的目的；海豚信号合成的低截获通信，计算出whistles脉冲串的包络，将海豚脉冲进行拟合，采用高斯信号及其各阶导函数模拟海豚信号，通过高斯函数及各阶导数来进行海豚信号的重构，将信息隐藏在高斯函数的各阶导数中。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710979305.5A CN107769862A (zh) | 2017-10-19 | 2017-10-19 | 一种仿生低截获通信方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710979305.5A CN107769862A (zh) | 2017-10-19 | 2017-10-19 | 一种仿生低截获通信方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107769862A true CN107769862A (zh) | 2018-03-06 |
Family
ID=61268924
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710979305.5A Pending CN107769862A (zh) | 2017-10-19 | 2017-10-19 | 一种仿生低截获通信方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107769862A (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110109087A (zh) * | 2019-05-07 | 2019-08-09 | 中国科学院声学研究所 | 一种声纳不规则探测范围显示方法及系统 |
| CN110429954A (zh) * | 2019-07-12 | 2019-11-08 | 厦门大学 | 一种基于混沌扩频的变换域水声通信方法 |
| CN111246463A (zh) * | 2020-02-10 | 2020-06-05 | 浙江机电职业技术学院 | 基于alpha稳定分布随机过程偏斜度参数的隐蔽通信方法 |
| CN112104439A (zh) * | 2020-11-23 | 2020-12-18 | 北京中超伟业信息安全技术股份有限公司 | 一种自适应录音干扰方法及系统 |
| CN112182868A (zh) * | 2020-09-21 | 2021-01-05 | 天津大学 | 一种基于鲸/海豚click叫声的浅海主动声呐探测方法 |
| CN112630739A (zh) * | 2020-11-30 | 2021-04-09 | 海鹰企业集团有限责任公司 | 一种梳状谱信号的宽带条件任意速度副本产生方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106503336A (zh) * | 2016-10-21 | 2017-03-15 | 哈尔滨工程大学 | 一种海豚嘀嗒声信号建模与合成的方法 |
| CN107231228A (zh) * | 2017-07-11 | 2017-10-03 | 厦门声戎科技有限公司 | 混合混沌的低截获通信系统设计方法 |
-
2017
- 2017-10-19 CN CN201710979305.5A patent/CN107769862A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106503336A (zh) * | 2016-10-21 | 2017-03-15 | 哈尔滨工程大学 | 一种海豚嘀嗒声信号建模与合成的方法 |
| CN107231228A (zh) * | 2017-07-11 | 2017-10-03 | 厦门声戎科技有限公司 | 混合混沌的低截获通信系统设计方法 |
Non-Patent Citations (4)
| Title |
|---|
| 张路蔚: "隐蔽水声通信技术研究", 《中国优秀硕士学位论文全文数据库》 * |
| 李玉强等: "几种新型主动声呐发射信号性能分析研究", 《舰船科学技术》 * |
| 葛轶洲等: "Alpha稳定分布噪声下的水下通信调制识别", 《通信对抗》 * |
| 魏玉华等: "低截获主动声呐信号特征研究", 《兵工自动化》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110109087A (zh) * | 2019-05-07 | 2019-08-09 | 中国科学院声学研究所 | 一种声纳不规则探测范围显示方法及系统 |
| CN110429954A (zh) * | 2019-07-12 | 2019-11-08 | 厦门大学 | 一种基于混沌扩频的变换域水声通信方法 |
| CN111246463A (zh) * | 2020-02-10 | 2020-06-05 | 浙江机电职业技术学院 | 基于alpha稳定分布随机过程偏斜度参数的隐蔽通信方法 |
| CN111246463B (zh) * | 2020-02-10 | 2023-04-18 | 浙江机电职业技术学院 | 基于alpha稳定分布随机过程偏斜度参数的隐蔽通信方法 |
| CN112182868A (zh) * | 2020-09-21 | 2021-01-05 | 天津大学 | 一种基于鲸/海豚click叫声的浅海主动声呐探测方法 |
| CN112182868B (zh) * | 2020-09-21 | 2023-03-24 | 天津大学 | 一种基于鲸/海豚click叫声的浅海主动声呐探测方法 |
| CN112104439A (zh) * | 2020-11-23 | 2020-12-18 | 北京中超伟业信息安全技术股份有限公司 | 一种自适应录音干扰方法及系统 |
| CN112630739A (zh) * | 2020-11-30 | 2021-04-09 | 海鹰企业集团有限责任公司 | 一种梳状谱信号的宽带条件任意速度副本产生方法 |
| CN112630739B (zh) * | 2020-11-30 | 2024-02-20 | 海鹰企业集团有限责任公司 | 一种梳状谱信号的宽带条件任意速度副本产生方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107769862A (zh) | 一种仿生低截获通信方法 | |
| Jia-jia et al. | Bio-inspired steganography for secure underwater acoustic communications | |
| Liu et al. | Biologically inspired covert underwater acoustic communication by mimicking dolphin whistles | |
| CN110535537A (zh) | 一种水下通信探测一体化方法 | |
| CN104753638B (zh) | 一种混沌扩频水声通信方法 | |
| CN105227246A (zh) | 一种利用分段lfm信号仿海豚哨声信号的水声通信方法 | |
| CN105391500A (zh) | 一种基于超宽带信号的仿海豚嘀嗒声水声通信方法 | |
| CN107231228B (zh) | 混合混沌的低截获通信系统设计方法 | |
| Liu et al. | Covert underwater acoustic communication using whale noise masking on DSSS signal | |
| CN109039506A (zh) | 一种水下移动信道仿真方法 | |
| Han et al. | Experimental demonstration of underwater acoustic communication using bionic signals | |
| CN111580110B (zh) | 一种基于浅海多径时延的复合码水声测距方法 | |
| Jiang et al. | Covert underwater communication based on combined encoding of diverse time-frequency characteristics of sperm whale clicks | |
| CN109347568A (zh) | 一种仿海豚哨声连续相位多元频率调制水声通信方法 | |
| Qiao et al. | A frequency hopping pattern inspired bionic underwater acoustic communication | |
| CN110736977A (zh) | 一种适用于声呐通探的一体化信号产生方法 | |
| Demirors et al. | Chirp-based LPD/LPI underwater acoustic communications with code-time-frequency multidimensional spreading | |
| Jiang et al. | A sonar-embedded disguised communication strategy by combining sonar waveforms and whale call pulses for underwater sensor platforms | |
| CN104062641B (zh) | 一种ofdm雷达距离像处理的新方法 | |
| CN110224958A (zh) | 一种基于混沌序列的正交宽带调制解调方法 | |
| CN106910508A (zh) | 一种仿海洋打桩声源的隐蔽水声通信方法 | |
| CN103944845A (zh) | 全频谱载波调制的相关检测方法 | |
| CN108279418B (zh) | 利用鲸声的组合探测与通信装置和方法 | |
| En et al. | Multiuser underwater acoustic communication based on multicarrier-multiple chirp rate shift keying | |
| Jia et al. | Bionic Covert Underwater Acoustic Communication Based on Dolphin Whistle with Time-delay |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180306 |
|
| RJ01 | Rejection of invention patent application after publication |