本发明属于通信技术领域，特别是一种利用长肢领航鲸叫声谐波数量进行编码的隐蔽通信方法及装置。
背景技术：
水下隐蔽声波通信在军事应用中具有非常重要的作用。水下声波通信设备在进行通信时，必须向水中发射出声波信号，以传递通信信号给通信接收方。由于水下声波通信设备必须向水中发射出声波信号，而这些声波信号很容易被敌人的侦查系统探测和识别，所以水下隐蔽声波通信设备的隐蔽性与水下隐蔽声波通信设备所发出的声波信号有直接的关系。为了提高水声通信的隐蔽性，传统方法一般采用频率在跳变的线性调频信号作为载频波形来增加敌人的侦查系统解码通信信息的难度，但是，由于线性调频信号具有非常明显的特征，如它的频率线性变化、它的功率谱为矩形等，所以很容易引起敌人侦查系统的注意，一旦敌人的侦查系统注意到了这些信息，就能利用这些线性调频信号对水下声波通信设备进行定位或是摧毁。
技术实现要素：
为克服现有技术的不足，本发明旨在提高水声通信的隐蔽性。为此，本发明采用的技术方案是，利用长肢领航鲸叫声谐波数量进行编码的通信方法，从长肢领航鲸的叫声中，找出m个叫声，且m＝2n，m、n为正整数，利用m个不同的叫声编码n个比特的二进制数据，实现通信。
具体地：从m个不同的叫声中任意取出2n-n个叫声分为一组，n为一个小于n且大于零的整数，并将该组叫声命名为gn-n；
根据n-n的大小，选用对应的叫声分组gn-n，并让每一个叫声携带不同的n-n个比特通信信息，将通信过程中，每两个相邻叫声之间的时间间隔设置为t，并按照进行发送。
在接收端，对接收到的信号进行归一化的短时能量谱计算，通过设定一个合适的阈值tv1，去检测和定位每个叫声在所接收信号时间序列中的位置和持续时间长度；其中阈值tv1的取值范围为0.01<tv1<0.8；接下来，根据每个叫声的位置、持续时间长度和时频图谱去搜索该叫声中时频图谱所包含主波和谐波的数量，计算出该叫声所对应的主波与谐波数量之和的值，根据主波与谐波数量之和的值解码出该叫声所携带的通信信息。
更具体地，搜索过程按如下方法进行：第一步：在叫声的持续时间里，且在叫声的时频图中垂直于时间轴的方向等间隔地画q条直线，q为正整数，且0<q<10000000000，第二步：设定一个阈值tv2，阈值tv2的取值范围为0.01<tv2<0.8，找出位于每条所画直线中，超过阈值tv2的时频图峰值点，然后计算每条所画直线中对应时频图峰值点的数量，假设第l条直线上所对应的时频图峰值点的数量为pl，则其它直线上所对应的时频图峰值点的数量为p1,p2,…,pq，然后对所有峰值点的数量进行求和并取平均如下
利用长肢领航鲸叫声谐波数量进行编码的通信装置，包括鲸叫声波形存储器、鲸叫声分组器，通信参数设置器、通信编码器、换能器、水听器、通信解码器和输出器构成，在发送端，鲸叫声波形存储器、鲸叫声分组器，通信编码器、换能器依次相连，通信编码器还与通信参数设置器相连；在接收端，水听器、通信解码器和输出器依次相连。
具体地，鲸叫声波形存储器采用一个同步动态存储器sdram来实现，用于存储m个不同的叫声；
鲸叫声分组器用于：从鲸叫声波形存储器1出鲸叫声，并对这些鲸叫声进行分组，并根据所携带通信信息量的多少，分成多个不同的叫声组；更具体地，可以从m个不同的叫声中任意取出2n-n个叫声分为一组，n为一个小于n且大于零的整数，并将该组叫声命名为gn-n；例如，可以从m个不同的叫声中取出2n-1个叫声，分为一组，并将该组叫声命名为gn-1，那么这组叫声就能携带n-1bit的通信信息；鲸叫声分组器通过可编程逻辑器件或微处理器芯片来实现；
通信参数设置器用于设置n-n值的大小，用于设置通信过程中叫声之间的时间间隔t，且用于设置通信过程中叫声的声源级sl这三个参数，通信参数设置器可以通过可编程逻辑器件或微处理器芯片来实现；
通信编码器对通信参数设置器传过来的三个参数进行识别；然后，一方面，根据n-n的大小，从鲸叫声分组器2中取出对应的叫声分组gn-n，并让每一个叫声携带不同的n-nbit通信信息；另一方面，通信编码器根据叫声之间的时间间隔t，将通信过程中，每两个相邻叫声之间的时间间隔设置为t，其中，t的取值范围为0s<t<50s；再一方面，通信编码器根据叫声的声源级sl的值，来控制输出给换能器的叫声信号发射功率，其中sl的大小范围为0db<sl<300db；进一步地，通信编码器接收来待发送的通信信息，然后利用其所选择的叫声分组gn-n、时间间隔t、声源级sl来将待传输通信信息通过换能器转换成声波信号发送出去；进一步地，通信编码器可以通过可编程逻辑器件或微处理器芯片来实现；
水听器将声波转换为电信号后，输送给通信解码器；
通信解码器不停的对接收到的信号进行归一化的短时傅立叶变换，然后得到整个接收信号的时频图谱，当叫声被接收到时，通信解码器同时也对接收到的信号进行归一化的短时能量谱计算，通过设定一个合适的阈值tv1，去检测和定位每个叫声在所接收信号时间序列中的位置和持续时间长度；其中阈值tv1的取值范围为0.01<tv1<0.8；接下来，通信解码器根据每个叫声的位置、持续时间长度和时频图谱去搜索该叫声中时频图谱所包含主波和谐波的数量；更具体地，搜索过程可以按如下方法进行：第一步：在叫声的持续时间里，且在叫声的时频图中垂直于时间轴的方向等间隔地画q条直线，q为正整数，且0<q<10000000000；第二步：设定一个阈值tv2，其中阈值tv2的取值范围为0.01<tv2<0.8，找出位于每条所画直线中，超过阈值tv2的时频图峰值点，然后计算每条所画直线中对应时频图峰值点的数量，假设第l条直线上所对应的时频图峰值点的数量为pl，则其它直线上所对应的时频图峰值点的数量为p1,p2,…,pq，然后对所有峰值点的数量进行求和并取平均如下
进一步地，通信解码器就连续不断解码出的通信信息传输给通信信息输出器，以完成整个通信过程。
本发明的特点及有益效果是：
本发明利用原始的长肢领航鲸所发出的不同的叫声去携带不同的通信信息，此种通信方式能够提高水下隐蔽声波通信设备的隐蔽性。
附图说明：
图1示出本发明中所述的四种叫声的时频图。
图2示出本发明中所述的通信发送平台10。
图3示出本发明中所述的通信接收平台11。
图2中：1为鲸叫声波形存储器；2为鲸叫声分组器；3为通信参数设置器；4为通信编码器；5为待传输信息；6为换能器。
图3中：7为水听器；8为通信解码器；9为通信信息输出器。
具体实施方式
与传统方法不同，本发明提出了一种新的伪装隐蔽通信方法和装置，它利用原始的长肢领航鲸所发出的不同的叫声去携带不同的通信信息。由于海洋中鲸所发出的叫声一般被各国的侦查系统识别为噪声信号，所以，本发明利用原始的长肢领航鲸所发出的不同的叫声去携带不同的通信信息，从面实现伪装隐蔽通信。
本发明提供一种利用长肢领航鲸叫声谐波数量进行编码的隐蔽通信方法及装置。
根据现有文献记载，可以知道，长肢领航鲸的叫声特性非常丰富，其中一个很重要的特性是它们所发出的每个叫声一般由一个主波和多个附波(这些“附波”一般被称为“谐波”，在下面的所有描述中，本发明统一使用“谐波”这个名词)组成，且不同的叫声一般会具有不同数量的谐波。例如，在图1中，本发明给出了四种叫声的时频图，在图1(a)中包括1个主波和4个谐波，在图1(b)中包括1个主波和1个谐波，在图1(c)中包括1个主波和2个谐波，在图1(d)中包括1个主波和15个谐波。在不同种类的叫声中，主波与谐波的数量之和与其它种类的叫声不同。假设能从长肢领航鲸的叫声中，找出m(m为一个正整数)个叫声(并假设m＝2n(n为一个正整数))，在这m个叫声里，每个叫声的主波与谐波数量之和与其它叫声彼此不同；就像图1中的四个叫声一样，这四个叫声的主波与谐波数量之和与其它叫声彼此不同。由于m＝2n，所以，m个不同的叫声就可能编码nbit的二进制数据。例如，在图1中，有4个不同的叫声，所以图1中的四个叫声可以编码2bit的二进制数据。
进一步地，本发明中所述鲸均指长肢领航鲸；
更进一步地，被找出的m个不同的叫声被存储在鲸叫声波形存储器1里，该鲸叫声波形存储器1可以用一个同步动态存储器(sdram)来实现；
进一步地，鲸叫声分组器2从鲸叫声波形存储器1取出鲸叫声，并对这些鲸叫声进行分组，并根据所携带通信信息量的多少，分成多个不同的叫声组；更具体地，可以从m个不同的叫声中任意取出2n-n(n为一个小于n且大于零的整数)个叫声，分为一组，并将该组叫声命名为gn-n；例如，可以从m个不同的叫声中取出2n-1个叫声，分为一组，并将该组叫声命名为gn-1，那么这组叫声就能携带n-1bit的通信信息；再例如，可以从m个不同的叫声中取出2n-(n-1)个叫声，分为一组，并将该组叫声命名为g1，那么这组叫声就能携带1bit的通信信息；进一步地，鲸叫声分组器2可以通过可编程逻辑器件或微处理器芯片来实现；
进一步进，通信参数设置器3用于设置n-n值的大小，用于设置通信过程中叫声之间的时间间隔t，且用于设置通信过程中叫声的声源级sl这三个参数；进一步地，通信参数设置器3可以通过可编程逻辑器件或微处理器芯片来实现；
进一步地，通信编码器4对通信参数设置器3传过来的三个参数进行识别；然后，一方面，根据n-n的大小，从鲸叫声分组器2中取出对应的叫声分组gn-n，并让每一个叫声携带不同的n-nbit通信信息；例如，如果n-n＝2，那么，通信编码器4将从鲸叫声分组器2中取出对应的叫声分组g2，该叫声分组里就包含4个不同的叫声脉冲，然后，通信编码器4让4个不同的叫声脉冲中的第一个叫声携带“00”通信信息，让4个不同的叫声脉冲中的第二个叫声携带“01”通信信息，让4个不同的叫声脉冲中的第三个叫声携带“10”通信信息，让4个不同的叫声脉冲中的第四个叫声携带“11”通信信息；另一方面，通信编码器4根据叫声之间的时间间隔t，将通信过程中，每两个相邻叫声之间的时间间隔设置为t，其中，t的取值范围为0s<t<50s；再一方面，通信编码器4根据叫声的声源级sl的值，来控制输出给换能器6的叫声信号发射功率，其中sl的大小范围为0db<sl<300db；进一步地，通信编码器4接收来自待传输信息5传来的待发送的通信信息，然后利用其所选择的叫声分组gn-n、时间间隔t、声源级sl来将待传输信息5存放的通信信息通过换能器6转换成声波信号发送出去；进一步地，通信编码器4可以通过可编程逻辑器件或微处理器芯片来实现；
进一步地，换能器7发出的声波信号经过水下传输后，到达通信接收平台上的水听器7，水听器7将声波转换为电信号后，输送给通信解码器8，通信解码器8不停的对接收到的信号进行归一化的短时傅立叶变换，然后得到整个接收信号的时频图谱，当叫声被接收到时，就会出现类似于图1中叫声的时频图谱；进一步地，通信解码器8同时也对接收到的信号进行归一化的短时能量谱计算，通过设定一个合适的阈值tv1，去检测和定位每个叫声在所接收信号时间序列中的位置和持续时间长度；其中阈值tv1的取值范围为0.01<tv1<0.8；接下来，通信解码器8根据每个叫声的位置、持续时间长度和时频图谱去搜索该叫声中时频图谱所包含主波和谐波的数量；更具体地，搜索过程可以按如下方法进行：第一步：在叫声的持续时间里，且在叫声的时频图中垂直于时间轴的方向等间隔地画q条(q为正整数，且0<q<10000000000)直线，第二步：设定一个阈值tv2(其中阈值tv2的取值范围为0.01<tv2<0.8)，找出位于每条所画直线中，超过阈值tv2的时频图峰值点，然后计算每条所画直线中对应时频图峰值点的数量，假设第l条直线上所对应的时频图峰值点的数量为pl，则其它直线上所对应的时频图峰值点的数量为p1,p2,…,pq，然后对所有峰值点的数量进行求和并取平均如下
进一步地，通信解码器8就连续不断解码出的通信信息传输给通信信息输出器9，以完成整个通信过程。