CN116707661A - 一种仿生隐蔽水声通信探测一体化方法 - Google Patents
一种仿生隐蔽水声通信探测一体化方法 Download PDFInfo
- Publication number
- CN116707661A CN116707661A CN202310490069.6A CN202310490069A CN116707661A CN 116707661 A CN116707661 A CN 116707661A CN 202310490069 A CN202310490069 A CN 202310490069A CN 116707661 A CN116707661 A CN 116707661A
- Authority
- CN
- China
- Prior art keywords
- signal
- whistle
- bionic
- detection
- original
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 110
- 239000011664 nicotinic acid Substances 0.000 title claims abstract description 96
- 238000004891 communication Methods 0.000 title claims abstract description 79
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000005236 sound signal Effects 0.000 claims abstract description 33
- 238000005070 sampling Methods 0.000 claims abstract description 7
- 238000012216 screening Methods 0.000 claims abstract description 4
- 238000013507 mapping Methods 0.000 claims description 27
- 238000005314 correlation function Methods 0.000 claims description 12
- 238000001228 spectrum Methods 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 8
- 230000021615 conjugation Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000003592 biomimetic effect Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000000827 velocimetry Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000003672 processing method Methods 0.000 abstract 1
- 238000013461 design Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 241000283153 Cetacea Species 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Classifications
-
- 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
- H04B14/00—Transmission systems not characterised by the medium used for transmission
- H04B14/02—Transmission systems not characterised by the medium used for transmission characterised by the use of pulse modulation
- H04B14/026—Transmission systems not characterised by the medium used for transmission characterised by the use of pulse modulation using pulse time characteristics modulation, e.g. width, position, interval
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/32—Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种仿生隐蔽水声通信探测一体化方法。现有方法没有联合考虑通信和探测的隐蔽性问题。本发明方法首先提取出嘀嗒声信号和哨声信号，进行功率和相关性筛选；然后对原始嘀嗒声信号进行波形编码和脉冲位置调制，对原始哨声信号进行时频形状编码和波形参数调制。利用仿生已调信号的良好探测性能，生成仿生通信探测一体化信号。在回波信号接收端，采用不同的回波信号处理方法。在通信接收端，采用组合解调方法。本发明生成的仿生哨声信号时域上与原始信号具有高相似度，且时频轮廓连续。本发明采样组合调制的方式，可以提高通信速率，具有较好的隐蔽性，不仅能够实现对目标距离、速度和方位角的准确估计，且具有良好的通信性能。
Description
技术领域
本发明涉及水声通信和目标探测一体化技术领域，特别涉及一种仿生隐蔽水声通信探测一体化方法。
背景技术
现代军事和民用平台需要同时实现通信和探测功能。按照传统方式进行分立式设计会增加系统的体积、功耗和成本，同时也会降低系统的兼容性能。采用通信探测一体化的设计方法可以让通信和探测模块共享硬件资源，从而解决上述问题。
为了防止用于军事目的水声通信信号在传输过程中被识别和破译，通信必须秘密进行。与传统的降低传输功率或增加调制带宽的方法不同，仿生隐蔽水声通信通过模仿鲸目动物的声音，诱导非通信合作方将鲸目动物的声音与仿生通信信号混淆，从而排除接收到的信号，实现隐蔽效果。由于仿生通信信号不需要降低传输能量，因此，仿生隐蔽水声通信可以获得比现有的隐蔽水下通信方法更好的传输距离和误码率性能。而主动声呐为了探测水下噪声环境中的远程目标，需要发射高功率探测信号以获得所需的信噪比。因此，主动声呐也具有安全性和隐蔽性需求。综合考虑水声通信隐蔽性和水声探测隐蔽性，对于水下作战平台等应用场景具有十分重要的意义。
目前已有的仿生隐蔽通信方法，如：申请号为202110111520.X的中国发明专利申请，公开了一种基于分段时频轮廓时延调制的仿鲸目哨声通信方法，该方法对哨声信号进行分段调制，局限在于调制后生成的仿生哨声信号时频轮廓产生不连续性，影响了仿生伪装性能。目前已有的水下通信探测一体化方法，如：申请号为201910909461.3的中国发明专利申请，公开了一种水下通信探测一体化方法，该方法将通信信息调制在相邻的pattern码之间的时延差值中，选用广义正弦调频信号作为pattern码，实现通信的同时，可以估计目标距离和速度，但是没有考虑对目标方位的估计。
综上可知，目前公开的专利中还没有关于仿生隐蔽水声通信探测一体化的方法。在仿生水声隐蔽通信方法中，需要考虑通信速率和隐蔽性的矛盾问题，使之既能提高通信速率，又不影响隐蔽性。在水下通信探测一体化方法中，需要提高对目标位置的估计准确度。
发明内容
本发明针对水声通信与目标主动探测皆存在隐蔽性问题，设计了一种仿生隐蔽水声通信与目标探测一体化方法，实现较高速率的水声通信，以及目标探测。
本发明方法具体是：
步骤(1)提取原始叫声信号：采用频带方差端点检测方法，提取原始嘀嗒声信号；采用双门限端点检测方法，提取原始哨声信号。
步骤(2)根据功率与相关性，筛选原始嘀嗒声信号和哨声信号：保留大于设定功率阈值，且大于自相关阈值、小于互相关阈值的原始嘀嗒声信号和原始哨声信号。
步骤(3)对筛选出的原始嘀嗒声信号和哨声信号进行组合信息编码和调制，具体如下：
(3-1)对筛选出的原始嘀嗒声信号进行波形编码，将2N1个嘀嗒声与2N1组N1位二进制比特信息构成编码映射关系，即每个原始嘀嗒声信号携带N1位比特信息，完成N1位比特信息的携带；
(3-2)对筛选出的原始嘀嗒声信号进行脉冲位置调制，将调制信息调制在相邻原始嘀嗒声信号之间的脉冲时间差，将调制信息q与2N2组N2位二进制比特信息构成编码映射关系，完成N2位比特信息的携带；其中，第i个和第i+1个原始嘀嗒声之间的时间差调制信息q＝1,2,…,2N2，Tc为每个嘀嗒声组合中相邻嘀嗒声之间的最大时间间隔，即调制时长，Tc占据2N2个最小时间差，即调制间隔/>共有2N2种可能的脉冲时间差。
(3-3)以筛选出的原始哨声信号为原型，构建仿生哨声信号，对仿生哨声信号进行波形参数调制：利用筛选出的原始哨声信号的时域包络，对单频信号和线性调频信号的包络进行调制，再进行波形参数调制，生成仿生哨声已调信号，表示为Al(t)表示原始哨声信号中t时刻第l个谐波的时域包络，l＝1,2,…,L，L为谐波个数，fc为原始哨声信号的中心频率，调频率/>δB为带宽的调制间隔，/>为持续时间的调制间隔，/>k＝0表示对单频信号进行调制，k≠0表示对线性调频信号进行调制，b和d为仿生哨声已调信号携带的调制信息，b和d的组合方案与2N3组N3位二进制比特信息构成编码映射关系，完成N3位比特信息的携带。
(3-4)对仿生哨声已调信号的时频形状进行编码，时频形状与2N4组N4位二进制比特信息构成编码映射关系，完成N4位比特信息的携带。
步骤(4)生成仿生通信探测一体化信号，由原始嘀嗒声的同步信号、嘀嗒声探测信号、嘀嗒声已调信号、仿生哨声探测信号、仿生哨声已调信号组成；其中，同步信号用于信号的帧同步，嘀嗒声探测信号和仿生哨声探测信号用于目标探测，嘀嗒声已调信号和仿生哨声已调信号用于进行通信。
步骤(5)主动目标探测：目标探测的接收端接收到通信探测一体化信号，采用时延测距法处理嘀嗒声探测信号，得到目标距离的估计；采用多普勒测速法和阵列测向法处理仿生哨声探测信号，分别得到目标速度和方位角信息的估计。具体如下：
(5-1)利用模糊度函数，估计发送信号与接收回波信号间的传播时间差，进行目标距离的估计：模糊度函数χ(τ,κ)是发送仿生通信探测一体化信号中的嘀嗒声探测信号u(t)及其回波信号u(κ(t-τ))的相关性函数，其中，κ为多普勒因子，τ为嘀嗒声探测信号模糊度函数峰值对应的时延，上标*表示共轭；被探测目标与回波信号接收端之间的距离估计值/> c为水下等效声速。
(5-2)采用分数阶傅里叶变换，估计多普勒因子，进行目标运动速度的估计。
对仿生哨声探测信号进行分数阶傅里叶变换，在最优阶傅里叶域，仿生哨声探测信号的能量得到聚集，产生峰值。此时，最优阶傅里叶域与时间轴的夹角为最佳旋转角
利用分数阶傅里叶变换最佳旋转角与仿生哨声已调信号调频率k之间的关系得到回波信号的调频率的估计值/> 其中，时间分辨率频率分辨率/>N为采样点数，fs为采样频率。
利用多普勒因子与调频率之间的关系，得到多普勒因子估计值 进而得到目标移动速度的估计值/>
(5-3)采用M元均匀分布线列阵，进行目标方位角的估计，具体是：
阵列输出信号其中，vp(t)是t时刻第p个仿生通信探测一体化信号中的哨声探测回波信号，p＝1,2,…,P，P是哨声探测回波信号的个数；D是阵元间距，θp是第p个仿生通信探测一体化信号中的哨声探测回波信号相对于阵列的入射方位角，nm(t)表示第m个阵元上的噪声，m＝1,2,…,M，M为阵列中阵元的个数，阵列输出协方差矩阵B＝E[x(t)xH(t)]，E[·]表示求均值运算，上标H表示转置共轭；
对每个阵元采用非相干信号子空间法进行目标方位角的估计：将宽带信号的频带划分成J个子带，对中心频率为fj的子带协方差矩阵B(fj)进行特征分解，B(fj)＝U(fj)Λ(fj)U(fj)H，j＝1,2,…,J；其中，特征向量U(fj)＝[US(fj)UN(fj)]，特征值US(fj)为信号特征向量，UN(fj)为噪声特征向量，ΛS(fj)为信号特征值，ΛN(fj)为噪声特征值；
当方向向量a(fj,θ)与信号子空间的欧氏距离最小时，在方位角θ处得到一个高的谱峰；利用每个子带的空间谱得到宽带的空间谱/>对宽带空间谱进行角度扫描，找到空间谱的峰值位置，峰值对应的信号入射角度为目标方位角的估计值/>
步骤(6)通信的接收端接收到通信探测一体化信号，进行仿生已调信号的组合解调和译码：
发射机和接收机共享叫声信号数据库，第i′个接收仿生已调信号ri′(t)与叫声信号数据库中的第j′个本地信号cj′(t)的互相关函数表示为τ′为相关峰值对应的时延；利用相关峰值所对应的本地信号的编号j′，与二进制序列进行逆映射，得到译码序列。具体如下：
(6-1)利用嘀嗒声已调信号与嘀嗒声本地信号计算的2N1个互相关函数值，得到相关峰值所在的位置，进行嘀嗒声本地信号编号的估计。利用嘀嗒声本地信号编号与二进制序列的逆映射，进行N1位比特信息的译码。
(6-2)利用相邻的第x个和第x′＝x+1个嘀嗒声已调信号相关峰值所在的位置tx和tx′，以及对应的已调信号的脉宽Tx和Tx′，得到两个嘀嗒声已调信号之间的时延差利用脉冲位置表达式/>得到q值，利用q和二进制序列的逆映射，进行N2位比特信息的解调。
(6-3)采用短时傅里叶变换，进行时频形状分析，确定时频形状对应的哨声本地信号；利用时频形状与二进制序列的逆映射，进行仿生哨声已调信号N4位比特信息的译码。
(6-4)利用(6-3)确定的哨声本地信号和仿生哨声已调信号计算的2N3个互相关函数值，进行哨声本地信号编号的估计；利用哨声本地信号编号与二进制序列的逆映射，进行N3位调制信息的解调。
本发明具有以下有益效果：
(1)本发明考虑了整体的原始哨声信号结构特点，而不是对哨声信号进行时域上的简单分割，一方面能让调制后的哨声信号携带更大的能量，降低通信的误码率，另一方面也是为了提高仿生性，不对哨声信号做分段处理，可以使信号的时频轮廓更平滑和连续。
(2)本发明考虑当前仿生隐蔽通信方法存在通信速率较低的问题，充分利用原始叫声信号的信号结构，即嘀嗒声和哨声的前后组合模式。分别对原始嘀嗒声信号和仿生哨声信号进行组合调制。提高了通信速率，并且使用原始嘀嗒声信号能进一步提高通信信号的隐蔽性。
(3)本发明基于叫声信号的良好探测性能，考虑主动声呐的隐蔽性问题，在仿生隐蔽水声通信方法设计的信号基础上，提出了一种仿生隐蔽水声通信探测一体化方法。采用原始嘀嗒声信号为距离探测信号，仿生哨声已调信号为速度和方位角探测信号。实现对目标距离、速度和方位角的准确估计，且能够保持仿生通信方法的通信性能。
附图说明
图1为本发明方法的流程图；
图2为调制信息q＝1,2,3,4时的脉冲位置调制示意图；
图3为仿生哨声已调信号示意图；
图4为仿生通信探测一体化信号帧结构示意图。
具体实施方式
以下结合附图并举实施例对本发明进一步详细说明。
一种仿生隐蔽水声通信探测一体化方法，具体流程如图1所示。
步骤(1)提取原始叫声信号：采用频带方差端点检测方法，提取原始嘀嗒声信号；采用双门限端点检测方法，提取原始哨声信号。
步骤(2)根据功率与相关性，筛选原始嘀嗒声信号和哨声信号：保留大于设定功率阈值，且大于自相关阈值、小于互相关阈值的原始嘀嗒声信号和原始哨声信号。
本实施例中原始嘀嗒声信号的功率阈值是0.1，原始哨声信号的功率阈值是0.6。原始嘀嗒声信号的自相关阈值是0.7，互相关阈值是0.3；原始哨声信号的自相关阈值是0.7，互相关阈值是0.1。
步骤(3)对筛选出的原始嘀嗒声信号和哨声信号进行组合信息编码和调制，具体如下：
(3-1)对筛选出的原始嘀嗒声信号进行波形编码，将2N1个嘀嗒声与2N1组N1位二进制比特信息构成编码映射关系，即每个原始嘀嗒声信号携带N1位比特信息，完成N1位比特信息的携带；
本实施例中嘀嗒声数据库中有128个高质量嘀嗒声，N1＝7。
(3-2)对筛选出的原始嘀嗒声信号进行脉冲位置调制，将调制信息调制在相邻原始嘀嗒声信号之间的脉冲时间差，将调制信息q与2N2组N2位二进制比特信息构成编码映射关系，完成N2位比特信息的携带；其中，第i个和第i+1个原始嘀嗒声之间的时间差调制信息q＝1,2,…,2N2，Tc为每个嘀嗒声组合中相邻嘀嗒声之间的最大时间间隔，即调制时长，Tc占据2N2个最小时间差，即调制间隔/>共有2N2种可能的脉冲时间差。
具体的调制方式原理如图2所示。本实施例中Tc＝200ms，向下取整等于153，N2＝7。
(3-3)以筛选出的原始哨声信号为原型，构建仿生哨声信号，对仿生哨声信号进行波形参数调制：利用筛选出的原始哨声信号的时域包络，对单频信号和线性调频信号的包络进行调制，再进行波形参数调制，生成仿生哨声已调信号，表示为Al(t)表示原始哨声信号中t时刻第l个谐波的时域包络，l＝1,2,…,L，L为谐波个数，fc为原始哨声信号的中心频率，调频率/>δB为带宽的调制间隔，/>为持续时间的调制间隔，/>k＝0表示对单频信号进行调制，k≠0表示对线性调频信号进行调制，b和d为仿生哨声已调信号携带的调制信息，b和d的组合方案与2N3组N3位二进制比特信息构成编码映射关系，完成N3位比特信息的携带。
如图3所示，本实施例中哨声信号的持续时间最大为500ms，取50ms时，有10(500ms/50ms)种信号持续时间方案(即b的取值个数为10)。哨声信号的带宽最大为4kHz，δB取200Hz时，有20(4000Hz/200Hz)种信号带宽方案(即d的取值个数为20)。一共有b×d＝200种组合，N3＝7，L＝3，fc＝1.5kHz。
(3-4)对仿生哨声已调信号的时频形状进行编码，时频形状与2N4组N4位二进制比特信息构成编码映射关系，完成N4位比特信息的携带。
本实施例中生成的仿生哨声的时频形状有三种：固定频率型、上扫频和下扫频，因此，N4＝2。
步骤(4)生成仿生通信探测一体化信号，由原始嘀嗒声的同步信号、嘀嗒声探测信号、嘀嗒声已调信号、仿生哨声探测信号、仿生哨声已调信号组成；其中，同步信号用于信号的帧同步，嘀嗒声探测信号和仿生哨声探测信号用于目标探测，嘀嗒声已调信号和仿生哨声已调信号用于进行通信。
仿生通信探测一体化信号帧结构示意图，如图4所示。
步骤(5)主动目标探测：目标探测的接收端接收到通信探测一体化信号，采用时延测距法处理嘀嗒声探测信号，得到目标距离的估计；采用多普勒测速法和阵列测向法处理仿生哨声探测信号，分别得到目标速度和方位角信息的估计。具体如下：
(5-1)利用模糊度函数，估计发送信号与接收回波信号间的传播时间差，进行目标距离的估计：模糊度函数χ(τ,κ)是发送仿生通信探测一体化信号中的嘀嗒声探测信号u(t)及其回波信号u(κ(t-τ))的相关性函数，其中，κ为多普勒因子，τ为嘀嗒声探测信号模糊度函数峰值对应的时延，上标*表示共轭；被探测目标与回波信号接收端之间的距离估计值/> c为水下等效声速。
本实施例中嘀嗒声探测信号u(t)的持续时间是0.02s，带宽为5kHz，c取1500m/s，κ为0.99，τ为2.7s，
(5-2)采用分数阶傅里叶变换，估计多普勒因子，进行目标运动速度的估计。
对仿生哨声探测信号进行分数阶傅里叶变换，在最优阶傅里叶域，仿生哨声探测信号的能量得到聚集，产生峰值。此时，最优阶傅里叶域与时间轴的夹角为最佳旋转角
利用分数阶傅里叶变换最佳旋转角与仿生哨声已调信号调频率k之间的关系得到回波信号的调频率的估计值/> 其中，时间分辨率频率分辨率/>N为采样点数，fs为采样频率。
利用多普勒因子与调频率之间的关系，得到多普勒因子估计值 进而得到目标移动速度的估计值/>
本实施例中仿生哨声探测信号的持续时间是0.1s，带宽为1kHz，中心频率为1.5kHz，fs＝44100Hz，N＝fs/Δf＝4410，Δf为10Hz，k＝10000，
(5-3)采用M元均匀分布线列阵，进行目标方位角的估计，具体是：
阵列输出信号其中，vp(t)是t时刻第p个仿生通信探测一体化信号中的哨声探测回波信号，p＝1,2,…,P，P是哨声探测回波信号的个数；D是阵元间距，θp是第p个仿生通信探测一体化信号中的哨声探测回波信号相对于阵列的入射方位角，nm(t)表示第m个阵元上的噪声，m＝1,2,…,M，M为阵列中阵元的个数，阵列输出协方差矩阵B＝E[x(t)xH(t)]，E[·]表示求均值运算，上标H表示转置共轭；
对每个阵元采用非相干信号子空间法进行目标方位角的估计：将宽带信号的频带划分成J个子带，对中心频率为fj的子带协方差矩阵B(fj)进行特征分解，B(fj)＝U(fj)Λ(fj)U(fj)H，j＝1,2,…,J；其中，特征向量U(fj)＝[US(fj)UN(fj)]，特征值US(fj)为信号特征向量，UN(fj)为噪声特征向量，ΛS(fj)为信号特征值，ΛN(fj)为噪声特征值；
当方向向量a(fj,θ)与信号子空间的欧氏距离最小时，在方位角θ处得到一个高的谱峰；利用每个子带的空间谱得到宽带的空间谱/>对宽带空间谱进行角度扫描，找到空间谱的峰值位置，峰值对应的信号入射角度为目标方位角的估计值/>
本实施例中仿生哨声探测信号的持续时间是0.1s，带宽为1kHz，中心频率为1.5kHz，有一个哨声探测回波信号，即P＝1。D＝0.5。M＝8，J＝5，P(θ)最大值是3.4452，对应的角度估计值
步骤(6)通信的接收端接收到通信探测一体化信号，进行仿生已调信号的组合解调和译码：
发射机和接收机共享叫声信号数据库，第i′个接收仿生已调信号ri′(t)与叫声信号数据库中的第j′个本地信号cj′(t)的互相关函数表示为τ′为相关峰值对应的时延；利用相关峰值所对应的本地信号的编号j′，与二进制序列进行逆映射，得到译码序列。具体如下：
(6-1)利用嘀嗒声已调信号与嘀嗒声本地信号计算的2N1个互相关函数值，得到相关峰值所在的位置，进行嘀嗒声本地信号编号的估计。利用嘀嗒声本地信号编号与二进制序列的逆映射，进行N1位比特信息的译码。
本实施例中编号为0的嘀嗒声本地信号，N1对应的比特信息是0000000，以此类推。
(6-2)利用相邻的第x个和第x′＝x+1个嘀嗒声已调信号相关峰值所在的位置tx和tx′，以及对应的已调信号的脉宽Tx和Tx′，得到两个嘀嗒声已调信号之间的时延差利用脉冲位置表达式/>得到q值，利用q和二进制序列的逆映射，进行N2位比特信息的解调。
本实施例中ΔTxx′＝10ms，q＝8，N2对应的比特信息是0001000。
(6-3)采用短时傅里叶变换，进行时频形状分析，确定时频形状对应的哨声本地信号；利用时频形状与二进制序列的逆映射，进行仿生哨声已调信号N4位比特信息的译码。
本实施例中固定频率型哨声信号，N4对应的比特信息是00；上扫频哨声信号，N4对应的比特信息是01；下扫频哨声信号，N4对应的比特信息是11。
(6-4)利用(6-3)确定的哨声本地信号和仿生哨声已调信号计算的2N3个互相关函数值，进行哨声本地信号编号的估计；利用哨声本地信号编号与二进制序列的逆映射，进行N3位调制信息的解调。
本实施例中编号为0的哨声本地信号，N3对应的比特信息是0000000，以此类推。
上述实施例所述的内容仅仅是对本发明实现形式的列举，本发明的保护范围不应限于实施例所陈述的具体形式，本发明的保护范围也应包含在本发明基础上所构思的同类发明方法。
Claims (4)
1.一种仿生隐蔽水声通信探测一体化方法，其特征在于，该方法具体是：
步骤(1)提取原始叫声信号：采用频带方差端点检测方法，提取原始嘀嗒声信号；采用双门限端点检测方法，提取原始哨声信号；
步骤(2)根据功率与相关性，筛选原始嘀嗒声信号和哨声信号：保留大于设定功率阈值，且大于自相关阈值、小于互相关阈值的原始嘀嗒声信号和原始哨声信号；
步骤(3)对筛选出的原始嘀嗒声信号和哨声信号进行组合信息编码和调制，编码和调制后，每个原始嘀嗒声信号分别携带N1和N2位比特信息，每个仿生哨声已调信号分别携带N3和N4位比特信息；
步骤(4)生成仿生通信探测一体化信号，由原始嘀嗒声的同步信号、嘀嗒声探测信号、嘀嗒声已调信号、仿生哨声探测信号、仿生哨声已调信号组成；其中，同步信号用于信号的帧同步，嘀嗒声探测信号和仿生哨声探测信号用于目标探测，嘀嗒声已调信号和仿生哨声已调信号用于进行通信；
步骤(5)主动目标探测：目标探测的接收端接收到通信探测一体化信号，采用时延测距法处理嘀嗒声探测信号，得到目标距离的估计；采用多普勒测速法和阵列测向法处理仿生哨声探测信号，分别得到目标速度和方位角信息的估计；
步骤(6)通信的接收端接收到通信探测一体化信号，进行仿生已调信号的组合解调和译码：
发射机和接收机共享叫声信号数据库，第i′个接收仿生已调信号ri′(t)与叫声信号数据库中的第j′个本地信号cj′(t)的互相关函数表示为τ′为相关峰值对应的时延；利用相关峰值所对应的本地信号的编号j′，与二进制序列进行逆映射，得到译码序列。
2.如权利要求1所述的一种仿生隐蔽水声通信探测一体化方法，其特征在于，步骤(3)具体过程如下：
(3-1)对筛选出的原始嘀嗒声信号进行波形编码，将2N1个嘀嗒声与2N1组N1位二进制比特信息构成编码映射关系，即每个原始嘀嗒声信号携带N1位比特信息，完成N1位比特信息的携带；
(3-2)对筛选出的原始嘀嗒声信号进行脉冲位置调制，将调制信息调制在相邻原始嘀嗒声信号之间的脉冲时间差，将调制信息q与2N2组N2位二进制比特信息构成编码映射关系，完成N2位比特信息的携带；其中，第i个和第i+1个原始嘀嗒声之间的时间差调制信息q＝1,2,…,2N2，Tc为每个嘀嗒声组合中相邻嘀嗒声之间的最大时间间隔，即调制时长，Tc占据2N2个最小时间差，即调制间隔/>共有2N2种可能的脉冲时间差；
(3-3)以筛选出的原始哨声信号为原型，构建仿生哨声信号，对仿生哨声信号进行波形参数调制：利用筛选出的原始哨声信号的时域包络，对单频信号和线性调频信号的包络进行调制，再进行波形参数调制，生成仿生哨声已调信号，表示为Al(t)表示原始哨声信号中t时刻第l个谐波的时域包络，l＝1,2,…,L，L为谐波个数，fc为原始哨声信号的中心频率，调频率/>δB为带宽的调制间隔，/>为持续时间的调制间隔，/>k＝0表示对单频信号进行调制，k≠0表示对线性调频信号进行调制，b和d为仿生哨声已调信号携带的调制信息，b和d的组合方案与2N3组N3位二进制比特信息构成编码映射关系，完成N3位比特信息的携带；
(3-4)对仿生哨声已调信号的时频形状进行编码，时频形状与2N4组N4位二进制比特信息构成编码映射关系，完成N4位比特信息的携带。
3.如权利要求2所述的一种仿生隐蔽水声通信探测一体化方法，其特征在于，步骤(5)具体过程如下：
(5-1)利用模糊度函数，估计发送信号与接收回波信号间的传播时间差，进行目标距离的估计：模糊度函数χ(τ,κ)是发送仿生通信探测一体化信号中的嘀嗒声探测信号u(t)及其回波信号u(κ(t-τ))的相关性函数，其中，κ为多普勒因子，τ为嘀嗒声探测信号模糊度函数峰值对应的时延，上标*表示共轭；被探测目标与回波信号接收端之间的距离估计值/>c为水下等效声速；
(5-2)采用分数阶傅里叶变换，估计多普勒因子，进行目标运动速度的估计；
对仿生哨声探测信号进行分数阶傅里叶变换，在最优阶傅里叶域，仿生哨声探测信号的能量得到聚集，产生峰值；此时，最优阶傅里叶域与时间轴的夹角为最佳旋转角
利用分数阶傅里叶变换最佳旋转角与仿生哨声已调信号调频率k之间的关系得到回波信号的调频率的估计值/>其中，时间分辨率频率分辨率/>N为采样点数，fs为采样频率；
利用多普勒因子与调频率之间的关系，得到多普勒因子估计值进而得到目标移动速度的估计值/>
(5-3)采用M元均匀分布线列阵，进行目标方位角的估计，具体是：
阵列输出信号其中，vp(t)是t时刻第p个仿生通信探测一体化信号中的哨声探测回波信号，p＝1,2,…,P，P是哨声探测回波信号的个数；D是阵元间距，θp是第p个仿生通信探测一体化信号中的哨声探测回波信号相对于阵列的入射方位角，nm(t)表示第m个阵元上的噪声，m＝1,2,…,M，M为阵列中阵元的个数，阵列输出协方差矩阵B＝E[x(t)xH(t)]，E[·]表示求均值运算，上标H表示转置共轭；
对每个阵元采用非相干信号子空间法进行目标方位角的估计：将宽带信号的频带划分成J个子带，对中心频率为fj的子带协方差矩阵B(fj)进行特征分解，B(fj)＝U(fj)Λ(fj)U(fj)H，j＝1,2,…,J；其中，特征向量U(fj)＝[US(fj)UN(fj)]，特征值US(fj)为信号特征向量，UN(fj)为噪声特征向量，ΛS(fj)为信号特征值，ΛN(fj)为噪声特征值；
当方向向量a(fj,θ)与信号子空间的欧氏距离最小时，在方位角θ处得到一个高的谱峰；利用每个子带的空间谱得到宽带的空间谱对宽带空间谱进行角度扫描，找到空间谱的峰值位置，峰值对应的信号入射角度为目标方位角的估计值/>
4.如权利要求3所述的一种仿生隐蔽水声通信探测一体化方法，其特征在于，步骤(6)具体过程如下：
(6-1)利用嘀嗒声已调信号与嘀嗒声本地信号计算的2N1个互相关函数值，得到相关峰值所在的位置，进行嘀嗒声本地信号编号的估计；利用嘀嗒声本地信号编号与二进制序列的逆映射，进行N1位比特信息的译码；
(6-2)利用相邻的第x个和第x′＝x+1个嘀嗒声已调信号相关峰值所在的位置tx和tx′，以及对应的已调信号的脉宽Tx和Tx′，得到两个嘀嗒声已调信号之间的时延差利用脉冲位置表达式/>得到q值，利用q和二进制序列的逆映射，进行N2位比特信息的解调；
(6-3)采用短时傅里叶变换，进行时频形状分析，确定时频形状对应的哨声本地信号；利用时频形状与二进制序列的逆映射，进行仿生哨声已调信号N4位比特信息的译码；
(6-4)利用(6-3)确定的哨声本地信号和仿生哨声已调信号计算的2N3个互相关函数值，进行哨声本地信号编号的估计；利用哨声本地信号编号与二进制序列的逆映射，进行N3位调制信息的解调。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310490069.6A CN116707661A (zh) | 2023-05-04 | 2023-05-04 | 一种仿生隐蔽水声通信探测一体化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310490069.6A CN116707661A (zh) | 2023-05-04 | 2023-05-04 | 一种仿生隐蔽水声通信探测一体化方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116707661A true CN116707661A (zh) | 2023-09-05 |
Family
ID=87832981
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310490069.6A Pending CN116707661A (zh) | 2023-05-04 | 2023-05-04 | 一种仿生隐蔽水声通信探测一体化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116707661A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117176536A (zh) * | 2023-09-15 | 2023-12-05 | 南京芯传汇电子科技有限公司 | 一种雷达探测与通信一体化波形设计与信号处理方法 |
-
2023
- 2023-05-04 CN CN202310490069.6A patent/CN116707661A/zh active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117176536A (zh) * | 2023-09-15 | 2023-12-05 | 南京芯传汇电子科技有限公司 | 一种雷达探测与通信一体化波形设计与信号处理方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | Biologically inspired covert underwater acoustic communication by mimicking dolphin whistles | |
| CN105227246A (zh) | 一种利用分段lfm信号仿海豚哨声信号的水声通信方法 | |
| CN109814094B (zh) | 一种多目标水声定位时延估计算法 | |
| Han et al. | Experimental demonstration of underwater acoustic communication using bionic signals | |
| CN116707661A (zh) | 一种仿生隐蔽水声通信探测一体化方法 | |
| CN114152916A (zh) | 一种基于先导信息的mimo雷达波形正交方法 | |
| Yu et al. | Estimating the delay-Doppler of target echo in a high clutter underwater environment using wideband linear chirp signals: Evaluation of performance with experimental data | |
| CN110535537A (zh) | 一种水下通信探测一体化方法 | |
| CN109765545A (zh) | 一种基于相关峰包络检测的水声定位时延估计算法 | |
| CN111736141A (zh) | 基于循环编码阵列的雷达通信一体化设计方法 | |
| CN111580110B (zh) | 一种基于浅海多径时延的复合码水声测距方法 | |
| Wei et al. | Doppler estimation based on dual-HFM signal and speed spectrum scanning | |
| CN109975842A (zh) | 一种基于小波变换的北斗卫星信号高精度盲捕获方法 | |
| CN110429954B (zh) | 一种基于混沌扩频的变换域水声通信方法 | |
| Jiang et al. | A sonar-embedded disguised communication strategy by combining sonar waveforms and whale call pulses for underwater sensor platforms | |
| Nusenu et al. | Dual-function MIMO radar-communications employing frequency-hopping chirp waveforms | |
| Han et al. | Signal design for underwater acoustic positioning systems based on orthogonal waveforms | |
| CN110927691A (zh) | 一种基于多时编码相位调制的低截获雷达信号设计方法 | |
| Peng et al. | JLHS: A joint linear frequency modulation and hyperbolic frequency modulation approach for speed measurement | |
| CN110109077A (zh) | 基于时间反演的mimo雷达相干角度估计方法 | |
| CN106910508B (zh) | 一种仿海洋打桩声源的隐蔽水声通信方法 | |
| CN110346802A (zh) | 基于计算水声信道参数的水下目标探测方法 | |
| Spindel | An underwater acoustic pulse compression system | |
| Ling et al. | Adaptive range-Doppler imaging and target parameter estimation in multistatic active sonar systems | |
| Chen et al. | On probing waveforms and adaptive receivers for active sonar |
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 |