本发明主要针对fm复合调制信号，主要涉及一种基于级联svm和全数字接收机的复合信号调制识别方法。
背景技术：
无线通信的快速发展衍生出越来越多复杂的信号类型，尤其随着物联网时代的来临，海量连接的逐步实现，各种广域低功耗的连接信号越来越多。传统的单一调制信号已经无法满足现代通信的要求，复合调制信号由于一般同时具有两种调制方式的优点并且两者的缺点也得到互补，已被广泛应用于新型通信信号中。因此针对复合调制信号和单一调制的混合信号开展调制识别方法和参数估计方法研究，对现在复杂通信环境下的调制识别技术发展具有重要意义。
目前对于复合调制信号的识别研究主要是基于局部信号特征和分类树的识别方法。该方法实现简单，实时性好，但其分类性能易受阈值影响且存在误差累积现象。支持向量机(svm)算法是深度学习算法的一种，该算法可以从标记数据中进行鲁棒且快速的特征学习，通过提取一组数据的支持向量来实现分类，很好的避免了决策树方法的缺点，因此被广泛应用于模式识别等领域。
由此可见，将深度学习技术与通信信号调制识别技术结合起来，来实现对复合调制信号的模式识别有着非常广泛的应用前景。
技术实现要素：
本发明的目的是解决测控系统中fm复合调制信号的识别问题，提出了一种基于级联svm和全数字接收机的复合信号调制识别方法。通过全数字接收机的鉴相误差提取fm复合信号的内调制信息，使用级联svm分类器对内调制信号的频谱和平方谱整体图形特征进行分类识别。理论分析和仿真实验表明，该算法与传统的基于局部信号特征和决策树分类的识别方法相比，有更好识别性能，特别是在低信噪比条件下。
本发明解决其技术问题所采用的技术方案包括以下步骤：
步骤1、将接收到的采样信号输入全数字接收机得到鉴相误差；
步骤2、获取鉴相误差的频谱图形用于区分mpsk-fm、mt-fm和2ask-fm信号；
步骤3、对被判别成mpsk信号的数据进行平方谱图形获取，用于区分bpsk-fm和qpsk-fm信号；
步骤4、通过级联svm分类器实现对上述信号的分类。
步骤1具体实现如下：
1-1.将接收到的采样信号与数控振荡器的两路输出相乘，通过fir低通滤波器滤除和频分量，得到正交的两路信号si(n)和sq(n)：
其中kf和ts均为常数，δf为nco初始频率与信号载频的差值，n为采样点数，θ(n)为nco初始相位。
1-2.将si(n)和sq(n)通过鉴相器得到鉴相误差
鉴相误差经过二阶环路滤波器滤波后反馈回数控振荡器直至环路达到稳定状态。此时有δf→0且
其中，kf和ts均为常数；因此鉴相后得到的
步骤2具体实现如下：
2-1.对经过接收机处理后得到的内调制信号进行fft变换，得到不同信号的频谱图形。
2-2.2ask信号的功率谱密度表达式为：
式中，ps(f)为单极性基带信号功率谱。因此2ask信号的功率谱和频谱在其载频f＝fc处存在单根谱线。
将多音调制信号看成多个不同载波的2ask信号的叠加，因此有:
由式(6)可知多音调制信号的频谱在其载频f＝fk处存在多根独立谱线。
mpsk信号的功率谱密度表达式与2ask信号类似，但是不同的是其基带码元为双极性码，其频带宽度较宽，通常表现为宽带信号；因此2ask信号的频谱图形中仅有一个峰值且其带宽极窄，多音信号的频谱峰值个数大于等于2，而mpsk信号的频谱图形具有单个峰值但带宽很宽的特点。
步骤3具体实现如下：
3-1.对mpsk信号，进行平方处理：
再对平方后的数据进行fft，获取平方谱图形特征。
3-2.对于bpsk信号，a(n)，b(n)为相同的双极性码，因此将式(7)改写为：
因此bpsk信号的平方谱在直流和2倍载频处存在峰值，而对于qpsk信号，a(n)，b(n)并不相同，乘积在正负1之间跳变，因此其功率谱不包含2倍载频分量，仅在直流处存在峰值。
步骤4具体实现如下：
4-1.在分类时首先将两类psk信号归为一类，并标记为mpsk。利用频谱图形数据训练svm分类器1，实现对2ask、mtone和mpsk信号的区分。之后对识别为mpsk信号的数据进行平方谱图形数据获取，然后输入svm分类器2，利用平方谱图形区分bpsk和qpsk信号。
4-2.训练时通过向信号文件中添加不同功率的噪声形成具有不同信噪比的混合信号，将此混合信号作为训练集，通过交叉验证方法得到svm分类器的最佳参数c和γ。
4-3.最后将未知信号的内调制频谱图形输入svm分类器，完成对未知信号的分类。
本发明有益效果如下：
1.本发明针对传统方法难以进行调制识别的复合信号，通过接收机和级联支持向量机可以实现了调制方式识别。
2.发明通过接收机输出的鉴相误差数据，提取内调制信号的频域图形特征，利用级联svm算法实现复合调制信号(2ask-fm、mtone_fm、bpsk-fm和qpsk-fm)的调制方式识别，识别率高于传统的基于局部信号特征和决策树的分类方法。
3.本发明提供了解决复合信号调制识别问题的新思路，避免了传统算法对阈值的依赖，具有更高的识别率和更强的鲁棒性。
附图说明
图1本发明中使用的全数字接收机原理图；
图2基于鉴相误差得到的内调制信号的频谱图；
图3mpsk信号的平方谱图形；
图4基于级联svm的调制识别流程图；
图5svm的关键参数c和γ对识别性能的影响；
图6不同算法在不同信噪比下的识别率曲线。
具体实施方式
下面结合附图对本发明的具体实施例做进一步的说明。
如图1-6所示，一种基于级联svm和全数字接收机的复合信号调制识别方法，具体实现步骤如下：
步骤1、将接收到的采样信号输入如图1所示的全数字接收机得到鉴相误差
步骤2、获取鉴相误差的频谱图形如图2所示，使用该频谱图形特征区分mpsk-fm、mt-fm和2ask-fm信号；
步骤3、对被判别成mpsk信号的数据进行平方谱图形获取，其平方谱特征如图3所示，使用该平方谱特征区分bpsk-fm和qpsk-fm信号；
步骤4、通过如图4所示的级联svm分类器实现对上述bpsk-fm和qpsk-fm信号的分类。
步骤1具体实现如下：
1-1.接收到的采样信号与数控振荡器的两路输出相乘，通过fir低通滤波器滤除和频分量，得到正交的两路信号si(n)和sq(n)如式(1)和式(2)所示：
1-2.将si(n)和sq(n)通过鉴相器得到鉴相误差如式(3)所示，鉴相误差经过二阶环路滤波器滤波后反馈回数控振荡器直至环路达到稳定状态。此时有δf→0且
其中kf和ts均为常数，因此鉴相后得到的
步骤2具体实现如下：
2-1.对经过接收机处理后得到的内调制信号进行fft变换，得到不同信号的频谱图形。
2-2.式(5)为2ask信号的功率谱密度表达式，式中ps(f)为单极性基带信号功率谱。如图2(a)所示，其频谱在f＝fc处存在单根谱线。
多音调制信号可以看成多个不同载波的2ask信号的叠加，其表达式为式(6)。因此其频谱如图2(b)所示，其频谱在f＝fk处存在多根独立谱线。
mpsk信号的功率谱密度表达式与2ask信号类似，但是不同的是其基带码元为双极性码，如图2(c)和2(d)所示，其频谱表现为宽带信号。
步骤3具体实现如下：
3-1.对mpsk信号，利用式(7)对其进行平方处理，再对平方后的数据进行fft，获取平方谱图形特征。
3-2.对于bpsk信号，a(n)，b(n)为相同的双极性码，因此式(7)改写为：
因此bpsk信号的平方谱在直流和2倍载频处存在峰值，去直流后其平方谱如图3中a所示。而对于qpsk信号，a(n)，b(n)并不相同，乘积在正负1之间跳变，因此其功率谱不包含2倍载频分量，仅在直流处存在峰值。去直流后其平方谱如图3中b所示。
步骤4具体实现如下：
4-1.分类的整体流程如图4所示：在分类时首先将两类psk信号归为一类，并标记为mpsk。利用频谱图形数据训练svm分类器1，实现对2ask、mtone和mpsk信号的区分。之后对识别为mpsk信号的数据进行平方谱图形数据获取，然后输入svm分类器2，利用平方谱图形区分bpsk和qpsk信号。
4-2.训练时通过向信号文件中和添加不同功率的噪声形成具有不同信噪比的混合信号，将此混合信号作为训练集，通过交叉验证方法得到svm分类器的最佳参数c和γ如图5所示。
4-3.最后将未知信号的内调制频谱图形输入svm分类器，完成对未知信号的分类。
实施例：
仿真时使用中国电子科技集团某研究所提供的复合调制信号文件(信噪比很大，仿真时视为不含噪声)。已知信号采样频率fs＝2mhz，载波频率fc＝500khz，码元速率fb＝40kbps，fm调频系数为α＝0.005。副载波频率为f1＝40khz，mt-fm调制信号的副载波为f1＝20khz，f2＝30khz，f3＝40khz，f4＝50khz，f5＝60khz。
仿真时首先向信号文件中分别添加不同功率的高斯白噪声，形成信噪比分别为10db和0db的带噪声信号，然后每4096点为一组，产生不同信号类型、不同信噪比的信号各100组作为训练集，使用交叉验证方法寻找参数c和γ的最优值。即将样本分为多份，对于不同的参数按比例将部分数据作为训练集，另一部分作为测试集，并对识别率进行统计，将识别率最高的参数c和γ作为最佳参数。
使用最佳参数，保持训练集不变，在相同条件下产生不同信噪比的四种信号各200组作为测试集进行识别率测试。测试结果如图6所示，对于2ask-fm信号和mt-fm信号，本专利方法的识别性能优于决策树分类算法，在0db信噪比时识别率达到95％。对于bpsk-fm和qpsk-fm信号，决策树方法的识别性能较差，在8db时识别率仅为92％，而本专利方法可以在3db时识别率即可达到95％以上。
传统的svm分类方案对四种信号使用一对一支持向量机时需要的分类器个数为m(m-1)/2＝6个，而本专利提出的级联svm方案需要的分类器个数仅为m1(m1-1)/2 m2(m2-1)/2＝4个，且本专利方案中每个向量机要训练的特征数仅为传统方案的一半，因此在训练和识别速度上更有优势。综上所述可以看出，本算法与基于局部信号特征和决策树的分类方法相比，避免了对阈值设置的依赖，提高了低信噪比下的识别性能，识别率有了很大提高。
最后需要注意的是，公布实施例的目的在于帮助进一步理解本发明，但是本领域的技术人员可以理解：在不脱离本发明及所附的权利要求的精神和范围内，各种替换和修改都是可能的。因此，本发明不应局限于实施例所公开的内容，本发明要求保护的范围以权利要求书界定的范围为准。
1.基于级联svm和全数字接收机的复合信号调制识别方法，通过全数字接收机的鉴相误差提取fm复合信号的内调制信息，使用级联svm分类器对内调制信号的频谱和平方谱整体图形特征进行分类识别，其特征在于包括以下步骤：
步骤1、将接收到的采样信号输入全数字接收机得到鉴相误差；
步骤2、获取鉴相误差的频谱图形用于区分mpsk-fm、mt-fm和2ask-fm信号；
步骤3、对被判别成mpsk信号的数据进行平方谱图形获取，用于区分bpsk-fm和qpsk-fm信号；
步骤4、通过级联svm分类器实现对bpsk-fm和qpsk-fm信号的分类。
2.根据权利要求1所述的基于级联svm和全数字接收机的复合信号调制识别方法，其特征在在于所述的步骤1具体实现如下：
1-1.将接收到的采样信号与数控振荡器的两路输出相乘，通过fir低通滤波器滤除和频分量，得到正交的两路信号si(n)和sq(n)：
其中kf和ts均为常数，δf为nco初始频率与信号载频的差值，n为采样点数，θ(n)为nco初始相位；
1-2.将si(n)和sq(n)通过鉴相器得到鉴相误差
鉴相误差经过二阶环路滤波器滤波后反馈回数控振荡器直至环路达到稳定状态；此时有δf→0且
其中，kf和ts均为常数；因此鉴相后得到的
3.根据权利要求2所述的基于级联svm和全数字接收机的复合信号调制识别方法，其特征在在于步骤2具体实现如下：
2-1.对经过接收机处理后得到的内调制信号进行fft变换，得到不同信号的频谱图形；
2-2.2ask信号的功率谱密度表达式为：
式中，ps(f)为单极性基带信号功率谱；因此2ask信号的功率谱和频谱在其载频f＝fc处存在单根谱线；
将多音调制信号看成多个不同载波的2ask信号的叠加，因此有:
由式(6)可知多音调制信号的频谱在其载频f＝fk处存在多根独立谱线；
根据频谱图形对信号进行判断：2ask信号的频谱图形中仅有一个峰值且其带宽极窄，多音信号的频谱峰值个数大于等于2，而mpsk信号的频谱图形具有单个峰值但带宽很宽的特点。
4.根据权利要求3所述的基于级联svm和全数字接收机的复合信号调制识别方法，其特征在在于步骤3具体实现如下：
3-1.对mpsk信号，进行平方处理：
再对平方后的数据进行fft，获取平方谱图形特征；
3-2.对于bpsk信号，a(n)，b(n)为相同的双极性码，因此将式(7)改写为：
因此bpsk信号的平方谱在直流和2倍载频处存在峰值，而对于qpsk信号，a(n)，b(n)并不相同，乘积在正负1之间跳变，因此其功率谱不包含2倍载频分量，仅在直流处存在峰值。
5.根据权利要求4所述的基于级联svm和全数字接收机的复合信号调制识别方法，其特征在在于步骤4具体实现如下：
4-1.在分类时首先将两类psk信号归为一类，并标记为mpsk；利用频谱图形数据训练svm分类器1，实现对2ask、mtone和mpsk信号的区分；之后对识别为mpsk信号的数据进行平方谱图形数据获取，然后输入svm分类器2，利用平方谱图形区分bpsk和qpsk信号；
4-2.训练时通过向信号文件中添加不同功率的噪声形成具有不同信噪比的混合信号，将此混合信号作为训练集，通过交叉验证方法得到svm分类器的最佳参数c和γ；
4-3.最后将未知信号的内调制频谱图形输入svm分类器，完成对未知信号的分类。
技术总结