[发明专利]一种亚音识别压制系统及其识别方法

权利要求书

1.一种亚音识别方法，其特征在于，包括如下步骤：

S1：对收到的中频信号采样数据进行解调，得到基带信号，解调采用快速计算公式：

其中，Freq为频率；i为采样数据的序号；I为同相位分量，代表向量在横轴上的投影；Q为90度相移分量，代表向量在纵轴上的投影；

S2：对基带信号进行带通滤波，带通范围为30Hz～300Hz；

S3：根据数字亚音码的码速率134.4bit/s对基带信号解码；

S4：在亚音识别模块读取数字亚音存储表，根据标准MPT1381，存储的数字亚音码后添加固定码与校验码作为标准码，数字亚音码传输时每帧为23位，帧结构为9位信息码+3位固定码(001)+12位校验码，校验码采用(23，12)Golay编码；

S5：将S3解码后的二进制数据依次与S4中得到的二进制码进行对比，若匹配，则输出匹配的数字亚音码，识别结束；若不匹配，则重复步骤S4与步骤S5，若与亚音识别模块存储表中的数字亚音码均无法匹配，说明数字亚音码不存在，转到S6进行模拟亚音频的识别；

S6：对带通滤波后的信号进行抽取，采样率≥800Hz；

S7：对抽样后的信号，采用Goertzel算法，计算存储的每个模拟亚音频频点的频谱幅值、最大值与均值，若最大值与均值的比值大于给定的门限，则最大值所对应的亚音频频点即为信号的亚音频；输出该模拟亚音品，识别结束；否则，信号不存在亚音频，识别结束。

2.根据权利要求1所述的亚音识别方法，其特征在于，S7中用Goertzel算法计算指定点的频谱值，步骤如下：

S71：对每个采样点进行递归运算(n＝0,1,…,N)

其初始条件为：v(-1)＝0，v(-2)＝0；

S72：当N个采样点计算完成后，计算频谱值

其中，x(n)为抽样后的信号，fs为采样率，f为给定的频率点。