一种主动声呐识别方法
【专利摘要】本发明公开了一种主动声呐识别方法，该方法是发射方在声呐发出的探测信号中嵌入水印，接收方通过对水印的探测来识别发送方信号，包括四个部分，数字水印预处理，水印生成，水印嵌入，水印探测。具体方法是在一个声呐波形的时频系数中嵌入水印，水印是由两个安全密钥产生。第一密钥是扩展码。第二密钥是用于选择和修改声呐的所选择的时间?频率单元的嵌入掩模，根据嵌入规则水印分布在声呐探测信号的分块DCT系数中，匹配滤波接收器产生经过信道模型的扩频水印的复制相关水印，探测器模型成功的检测需要访问这两个扩展码和嵌入掩模。本发明可以产生精度更高的检测概率，增加了可控性和识别性能。
【专利说明】
一种主动声呐识别方法
技术领域
[0001] 本发明涉及水声探测领域，尤其涉及一种主动声呐识别方法。
【背景技术】
[0002] 水声信道中充满了天然存在的或人为的来源不明的声发射，其中一个主要的发射 来源就是主动声呐。在水声对抗领域，利用虚假水声信号迷惑、诱骗对方已成为各国海军常 用的技术手段。随着水声对抗技术的发展，敌对方常常通过发射高度相似的信号实施干扰 或者诱骗。因此，迫切需要开展主动声呐信号的身份识别新技术的研究。因此出现通过对信 号嵌入水印来认证发送声呐波形。这种方法对实现有效的识别来自于敌方或是我方的声呐 探测信号至关重要，水下环境的复杂性，错综复杂的信息发送，有效的识别出正确信息对正 常探测发挥着重要作用。
[0003] 之前出现应用于探测的方法最早是正交频分复用(0FDM) "星座抖动"，其中的抖动 就是水印，根据抖动的特性来识别探测信号，这种方法容易受环境影响，错误识别率较高。 为解决这一问题数字水印技术应运而生，数字水印技术的关键就是如何提高信号识别的精 度。海洋水声探测中存在相似信号以及伪装信号或是复制信号，该技术可以有效的识别探 测信号的真实身份并且有效识别出来伪装或是复制信号，另一方面的应用就是在网络节点 为基础进行的通信导航，网络节点具有移动特性，他们发射的信号中嵌入发射信号的位置 信息对无人航行器的导航有重要的价值，基于现有通过目标特征进行鉴别分类这种方法是 不足的，有学者提出在主动声呐发射的信号中嵌入数字水印，通过检测接收信号中是否含 有水印鉴别其身份。针对主动声呐信号的鉴别提出分别用短时傅里叶变换(STFT)和离散余 弦变换(DCT)进行水印嵌入。但是算法比较复杂而且并没有充分考虑到海洋信道多途、衰减 以及多普勒频移等因素。又提出信道相关特性对水印声呐的影响以及利用嵌入水印提高探 测性能等方法，在已经提出的方法中根据水印嵌入并没有考虑水印嵌入的规则的改变，以 及嵌入带宽的改变，受水下环境影响较大。本文是基于在已经提出的理论基础上结合对水 下声信道的考虑，提出利用二维变换基本矩阵在时频域的分块DCT变换，通过改变系数嵌入 扩频水印，产生更高的检测性能，以及最佳的嵌入方式。本发明提出的方法全面考虑到水印 探测环境对探测性能的提高，该方法对这些信号的身份识别能力进一步加强，所以如何有 效的分辨出真伪信号以及同类别的声呐识别信号身份，成为现在亟需解决的问题。
【发明内容】
[0004] 发明目的：本发明针对声呐探测信号以及相似度极高的伪装信号的身份识别，能 有效的分辨这些信号，
[0005] 技术方案:一种主动声呐识别方法，包括如下步骤：
[0006] S1:将水印进行预处理
[0007] 步骤S1-1:将声呐探测信号生成x(t)，并且对其取样点为N产生离散信号X;
[0008] 步骤S1-2:对步骤S1-1所得到的离散信号进行分块，将T秒的时间间隔进行取样分 成Μ块，每一块包含着m = N/M;每一块长度是Tb = T/M，X被排列为m X Μ的矩阵；
[0009]步骤S1-3:将步骤S1-2分块后的信号利用二维DCT变换的基本酉矩阵Α进行DCT变 换产生DCT系数矩阵Θ ;其中A是mXM矩阵；
[0010] S2:生成以及嵌入水印序列
[0011] 步骤S2-1:生成水印序列，水印序列产生的基础是扩频矩阵C=[ci，c2,. . .，CM]e RmXM其中ce{±l}和二进制掩码Me{〇，l}两个秘钥的产生，扩频序列是结合信道模型产生， 掩码则是由嵌入位置决定；
[0012] 步骤S2-2:在步骤S2-1产生的扩频矩阵和二进制掩码两者的基础上进行Hadamard 乘积形成扩频水印W;
[0013] 步骤S2-3:扩频矩阵掩蔽水印矩阵，并且包含当在对应的位置的掩模具有一个逻 辑1，只对DCT系数对应于由水印所修改的1的位置的系数；
[0014] 步骤S2-4:生成的水印W如所述步骤S2-3嵌入DCT系数中，嵌入水印后的系数矩阵 为水印强度K为可控因子，根据SWR(信号和水印功率之比)来改变；
[0015] 步骤S2-5:根据产生的被水印的系数矩阵通过矩阵变换产生被水印的声呐信号矩 阵XW(3XW被重新排列为1 XN取样矢量；
[0016] 步骤S2_6:Xw(t)表示加入水印的声呐波形，由此可以得出水印w(t)是被水印波形 和没有水印波形之差；
[0017] S2:对水印进行检测
[0018] 步骤S3-1:接收含有水印的信号yw(t)，由步骤s2_6得到的xw(t)经过信道脉冲响h (t，T)所得到的接收响应；
[0019] 步骤S3-2:将所得到的yw(t)进行重新取样和重新排列为mXM矩阵Y w，并且对￥|进 行DCT系数转换，得到转换的DCT系数；
[0020] 步骤S3-3:水印是根据在步骤S2-6的基础上得到，水印通过被滤波器模型化的信 道产生补偿水印?(〇·根据接收到的水印的估计值，计算在接收器通过由相应的FIR滤波器 的滤波每个段，并在检测器用于副本相关，估计水印值表示为 h
[0021] 步骤S3-4:对#(〇进行取样并且排成mXM矩阵U，对矩阵U进行块DCT变换得到系数 矩阵# ;
[0022] 步骤S3-5:在上述得出和A，.检测标准定为|):和#_的;1^(^611；[118内积，这个相当于 整个对hadamard乘积之和，根据信号检测算法得到检验统计量；
[0023]步骤S3-6:在检验统计量的计算中，水印检测被规定为根据假设检验奈曼皮尔逊 引理对接收水印进行检测假设，根据公式产生检测统计量作为门限;利用检测概率，漏检概 率，误检率来评估性能。
[0024]有益效果:本发明采用一种基于数字水印技术的基础上进行信号鉴定的方法，主 要是包括三个主要部分水印预处理，水印序列的生成以及嵌入，水印探测三个主要的部分。 相对于现有技术，本发明直接在声呐的发射信号嵌入水印，探测方可根据嵌入水印的特征 进行相应的身份识别，并且嵌入的水印以及嵌入强度都是根据情况可控的，由于该方法嵌 入的方式的可变可选的，结合数字处理算法提高分辨精度。本发明在水声领域可以得到较 灵活的应用。更确切的说是声呐信号的鉴定基于复杂的水下环境，水声信道是充满了来源 不明的声发射，在水下自主航行器(AUV)、水下潜艇之间通过网络节点进行通信具有实时移 动特性，本发明可以有效的提高导航的精度或是身份识别的准确性。
【附图说明】
[0025] 图1是本发明水印生成及嵌入模块图；
[0026] 图2是本发明水印探测模块图；
[0027]图3是本发明水印嵌入框图图；
[0028]图4是本发明水印探测框图图。
【具体实施方式】
[0029]下面将结合附图，对本发明的实施案例进行详细的描述；
[0030]如图1-4所示，本发明的一种主动声呐识别方法，是发射方在声呐发出的探测信号 中嵌入水印，接收方通过对水印的探测来识别发送方信号，包括四个部分，数字水印预处 理，水印生成，水印嵌入，水印探测。具体方法是在一个声呐波形的时频系数中嵌入水印，水 印是由两个安全密钥产生。第一密钥是扩频码。第二密钥是用于选择和修改声呐的所选择 的时间-频率单元的嵌入掩模，根据嵌入规则水印分布在声呐探测信号的分块DCT系数中， 匹配滤波接收器产生经过信道模型的扩频水印的复制相关水印，探测器模型成功的检测需 要访问扩频码和嵌入掩模。这种方法可以产生精度更高的检测概率，增加了可控性和识别 性能。
[0031 ]本发明的一种主动声呐识别方法，包括如下步骤：
[0032] 水印进行预处理步骤：
[0033] 步骤S1-1:声呐探测信号的表示x = exp( j*t .*2*pi*f0+j*pi*k*t. ~2)其中f0表示 初始频率，k表示调频斜率，对信号进行取样得到
[0034] X= {x(n)，n = 1，2· · ·，N} (1)
[0035] N为采样点的个数，采样时间长度是T。
[0036] 步骤S1-2:对所得到的离散信号进行分块，分成Μ块，每一块长度为Tb = T/M，每一 块的采样点数为m，m = N/M，并且用矩阵X表示。
[0037]步骤S1-3:分块之后的信号矩阵X进行DCT变换，本发明采用的是二维DCT变换的基 本矩A，A表示成如下式
[0039] A是酉矩阵，变换对可以通过一下表示
[0040] 0 =AXeRmXM
[0041 ] χ=ΑτΘ eRmXM (3)
[0042] Θ是每一块的DCT系数
[0043] Θ =[θ!,θ2......ΘΜ] (4)
[0044] 水印产生及嵌入步骤：
[0045] 步骤S2-1:水印序列是由扩频矩ce {±1}，
[0046] C=[ci,C2,...,cM]eRmXM (5)
[0047] 二进制掩码Me {〇，1}两个秘钥结合产生的二进制掩码是在扩频的基础上产生，通 过对扩频序列的重塑，本发明采用的扩频序列是PN序列。分散到选择的时间-频率单元中的 水印位是由一个扩频率=矩阵的乘积c e {± 1}和一个二进制掩码Me {〇，1}控制。
[0048] 步骤S2-2:在已经产生的扩频矩阵和二进制掩码两者的基础上进行Hadamard乘积 形成扩频水印w。W是由C和Μ的hadamard相乘建模，其定义为两个矩阵的逐元素乘积，水印负 载是单一的比特we {〇，i}，w表示成如下：
[0049] ff=wC〇MeRmXM
[0050] =[ffi,ff2, , , ,Wm]
[0051] ffj e RmX1 (6)
[0052] 步骤S2-3:，当在对应的位置的掩模具有一个逻辑1，只对DCT对应于由水印所修改 的1的位置的系数，乘法由掩模部分地破坏了扩展序列的正交性但是这并没有实质性的影 响。
[0053]步骤S2-4:水印如上述S2-3所述嵌入DCT系数中，嵌入系数后的被水印的系数矩阵 可以表示如下
[0054] 0W= 0+kffeRmXM (7)
[0055] 其中K表示水印强度根据SWR(信号和水印功率之比）。
[0056] 步骤S2-5:根据产生的被水印的系数矩阵通过矩阵变换产生被水印的声呐信号
[0057] 矩阵 Xw，Xw=AT0w (8)
[0058] Xw被重新排列为1XN取样矢量。
[0059] 步骤S2_6:Xw(t)表示加入水印的声呐波形，由此可以得出水印w(t)是被水印波形 和没有水印波形之差，表示如下
[0060] w(t) =xw(t)-x(t) (9)
[0061 ] 水印检测步骤：
[0062]步骤S3-1:接收含有水印的信号yw(t)，它是由步骤s2_6得到的xw(t)经过信道脉冲 响h(t，W所得到的接受相应，其中信道脉冲响应表达式为
[0064] 其中P为路径指数，AP(t)表示时变路径相位，Tp(t) = Tp_aPt表示时变延时，aP表示 多普勒因子。通过脉冲响应的式子可以表示为：
[0066] 步骤S3-2:将所得到的yw(t)进行重新取样和重新排列为mXM矩阵Yw，并且对￥|进 行DCT系数转换，转换的DCT系数表示为
[0068] 其中$表示第i块接收水印信号。
[0069] 步骤S3-3:水印是根据在步骤s2_6的基础上得到，水印通过被滤波器模型化的信 道产生补偿水印何0。
[0071 ]表示为每一块的水印。FIR滤波器的系数
[0072] b = [bi，b2......bP] (14)
[0073] 接近时不变信道传递函数和由此延伸的脉冲响应H(z)
[0074] H(z) = l+biz_1+......+bPz-p (15)
[0075] 信号被分为M块每一个长度是Tb，信道被M个滤波器模型化，对于每一块，滤波器抽 头的滤波器组可以被布置作为Β矩阵的列
[0076] Β - [hv b2...,. ,bM), bi - [b^, bL . . . . .b^}, h(t) 0 hi (:16)
[0077] bi是表示对为时变信道脉冲响应/?,_(Γ)以上第i块的持续时间的估计值的滤波器抽 头
[0079]接收到的水印的估计值，计算在接收器通过由相应的FIR滤波器的滤波每个段，并 在检测器用于副本相关。估计水印值可以表示为
[0081 ] 步骤S3-4 :对#(〇进行取样并且排成m X Μ矩阵U，对矩阵U进行块DCT变换得到 二 ζ??-/ (.19)步骤 S3-5 :在上述得出 _ 和#， 检测标准定为(|和_#的frobenius内积，这个相当于整个对应的hadamard乘积；
[0085] 其中rj表示内积，其表达式如：
[0086] rj: ^wkOj^Wj > (22)
[0087] 如果水印有扩频码那么rj可以表示如下：
[0089]第一个和最后一个都是不相关的，因为扩频码和隐藏的噪声和DCT系数是没有关 系的，对于第二部分在达到m时有一个峰值，是每一块最大的扩频处理增益，如果声呐没有 被水印，w = 〇,和^最终T将会类似于高斯噪声分布，如果水印探测尝试了错误的密匙V要 么不属于所述码本或属于不同的或错误的来源，采用以下表示式子
[0091] η又会再一次接收到类噪声的信号。
[0092] 步骤S3-6:以下一个检验统计量的计算中，水印检测被规定为根据假设检验奈
[0093] 曼皮尔逊引理，特别是每一个假设都被定义为如下
[0095]如果假设Ho被拒绝有利于Hi，如果
[0097] L(.)表示似然比A被选择以实现一个指定的显着性水平α。
[0098] 上述Γ⑴决定这检测率PD，Pm，Pfp，水印检测率Pd表示检测出嵌入水印的信号含有 水印的概率，水印漏检率P M表示嵌入的水印信号中不含有水印的概率，水印误检率PFP表示 嵌入的水印信号没有被检测出来水印的概率。
[0099] 本发明的有益效果是直接在声呐的发射信号嵌入水印，探测方可根据嵌入水印的 特征进行相应的身份识别，并且嵌入的水印以及嵌入强度都是根据情况可控的，由于本发 明嵌入的方式的可变可选的，结合数字处理算法提高分辨精度。在水声领域可以得到较灵 活的应用。
【主权项】
1. 一种主动声响识别方法，其特征在于，包括如下步骤： S1:将水印进行预处理 步骤S1-1:将声响探测信号生成x(t)，并且对其取样点为N产生离散信号X; 步骤S1-2:对步骤S1-1所得到的离散信号进行分块，将T秒的时间间隔进行取样分成Μ 块，每一块包含着m = Ν/Μ;每一块长度是Tb = Τ/Μ，X被排列为m X Μ的矩阵； 步骤S1-3:将步骤S1-2分块后的信号利用二维DCT变换的基本酉矩阵A进行DCT变换产 生DCT系数矩阵Θ ;其中A是mXM矩阵； S2:生成W及嵌入水印序列 步骤S2-1 :生成水印序列，水印序列产生的基础是扩频矩阵C= ，C2，. . .，CM] e RmXM其 中ce{±l}和二进制掩码Me {0,1}两个秘钥的产生，扩频序列是结合信道模型产生，掩码 则是由嵌入位置决定； 步骤S2-2:在步骤S2-1产生的扩频矩阵和二进制掩码两者的基础上进行化damard乘积 形成扩频水印W; 步骤S2-3:扩频矩阵掩蔽水印矩阵，并且包含当在对应的位置的掩模具有一个逻辑1， 只对DCT系数对应于由水印所修改的1的位置的系数； 步骤S2-4:生成的水印W如所述步骤S2-3嵌入DCT系数中，嵌入水印后的系数矩阵为0W; 水印强度K为可控因子，根据SWR(信号和水印功率之比）来改变； 步骤S2-5:根据产生的被水印的系数矩阵通过矩阵变换产生被水印的声响信号矩阵Xw。 Xw被重新排列为1 X N取样矢量； 步骤S2-6: xw(t)表示加入水印的声响波形，由此可W得出水印w(t)是被水印波形和没 有水印波形之差； S2:对水印进行检测 步骤S3-1:接收含有水印的信号yw(t)，由步骤S2-6得到的xw(t)经过信道脉冲响h(t，T) 所得到的接收响应； 步骤S3-2:将所得到的yw(t)进行重新取样和重新排列为mXM矩阵Yw，并且对Yw进行DCT 系数转换，得到转换的DCT系数； 步骤S3-3:水印是根据在步骤S2-6的基础上得到，水印通过被滤波器模型化的信道产 生补偿水印成的，根据接收到的水印的估计值，计算在接收器通过由相应的FIR滤波器的滤 波每个段，并在检测器用于副本相关，估计水印值表示为別/); 步骤S3-4:对进行取样并且排成mXM矩阵U，对矩阵U进行块DCT变换得到系数矩阵 乐； 步骤S3-5:在上述得出.食和W，检测标准定为0和亦的frobenius内积，运个相当于整个 对hadamard乘积之和，根据信号检测算法得到检验统计量； 步骤S3-6:在检验统计量的计算中，水印检测被规定为根据假设检验奈曼皮尔逊引理 对接收水印进行检测假设，根据公式产生检测统计量作为口限；利用检测概率，漏检概率， 误检率来评估性能。
【文档编号】G01S7/537GK106093919SQ201610397288
【公开日】2016年11月9日
【申请日】2016年6月7日 公开号201610397288.X, CN 106093919 A, CN 106093919A, CN 201610397288, CN-A-106093919, CN106093919 A, CN106093919A, CN201610397288, CN201610397288.X
【发明人】王彪, 丁智慧, 杨奕飞, 戴跃伟
【申请人】江苏科技大学