CN102868653B - 一种基于双谱和稀疏矩阵的数字调制信号分类方法 - Google Patents
一种基于双谱和稀疏矩阵的数字调制信号分类方法 Download PDFInfo
- Publication number
- CN102868653B CN102868653B CN201210330854.7A CN201210330854A CN102868653B CN 102868653 B CN102868653 B CN 102868653B CN 201210330854 A CN201210330854 A CN 201210330854A CN 102868653 B CN102868653 B CN 102868653B
- Authority
- CN
- China
- Prior art keywords
- digital modulation
- modulation signals
- lambda
- sparse matrix
- 4fsk
- 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.)
- Expired - Fee Related
Links
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
本发明公开了一种基于双谱和稀疏矩阵的数字调制信号分类方法,提出了两步法来对数字调制信号进行分类。第一步已知调制类型的数字调制信号进行双谱估计,利用压缩感知的方法对数字调制信号的双谱进行压缩,得到其稀疏表示,根据这个稀疏表示制作字典;第二步对输入待分类的数字调制信号进行第一步相同的变换,得到的变换结果与字典进行对比,计算残差均值,将残差均值最小的那个类型确定为信号的调制方式。通过与以前的基于经典特征的支持向量机(SVM)分类算法相比较,本发明提出的分类方法在分类正确率上有了明显的提高。
Description
技术领域
本发明属于认知无线电技术领域,更为具体地讲,涉及一种基于双谱和稀疏矩阵的数字调制信号分类方法。
背景技术
无线电频谱是一种宝贵的自然资源,它的分配通常是由无线电法规部门确定的。目前世界各国采用的是基于固定频谱分配的原则,随着无线通信需求的不断增长,人们对通信速率的需求也越来越高。根据香农理论,通信速率越高,通信系统所需占用的频谱带宽也越来越宽,从而导致适用于无线通信的频谱资源变得日益紧张。另一方面调查显示,无线电通信的频谱资源的利用率很低。
在解决日益稀缺的频谱问题上,认知无线电(Cognitive Radio,简称CR)是一种很有前景的技术。认知无线电被视为以适应在真实环境中无线电频谱的使用不断变化的智能无线通信系统。这种频谱随机接入策略是依赖于频谱感知,频谱感知是认知无线电技术的突出特征。为提高频谱感知正确率,已进行了很多关于频谱感知技术的研究,然而,实现这样的认知无线电感知技术,需要知道输入数字调制信号的调制类型。
自动调制分类(AMC)技术在非合作和动态通信环境,防止有害干扰授权用户,提高频谱的利用率中起着重要的作用。一般情况下,自动调制分类分为两个类型:决策理论和调制分类。决策理论通常涉及数字调制信号的统计特征和似然比;调制分类则主要依赖于特征的选取。数字调制信号分类的特征选取没有一个特定的方式,当已经确定特征的时候,只有通过仿真计算才能确定是否是最优的,数字调制信号的分类准确率有待提高。
发明内容
本发明的目的在于克服现有技术的不足,提供一种分类准确率更高的基于双谱和稀疏矩阵的数字调制信号分类方法。
为实现上述目的,本发明基于双谱和稀疏矩阵的数字调制信号分类方法,其特征在于,包括以下步骤:
(1)、稀疏矩阵字典的建立
1.1)、分别计算确知的四种数字调制信号2FSK、4FSK、QPSK、16QAM的训练序列x2FSK(t),x4FSK(t),xQPSK(t),x16QAM(t)的双谱,得到其双谱表示,对每个序列进行以下处理:
将长度为N的训练序列分割成M段,计算每个分段的双谱
其中L1是预设整数,Δ0是双谱域的频率采样间隔, K=N/M,xi(·)表示数字调制信号的第i个部分;
训练序列的双谱采用均值定义,表示为
这样四种数字调制信号2FSK、4FSK、QPSK、16QAM的训练序列x2FSK(t),x4FSK(t),xQPSK(t),x16QAM(t)的双谱为:以及
1.2)、对四种数字调制信号的双谱分别进行稀疏表示,得到各自的稀疏矩阵ψ:
ψ=[v1,v2,...vn,]∈Rn
Rn是n维实数域,vj是第j个数据,j=1,2…,K;
由于维度为n相对较高,计算复杂度高,需要进行降维处理:
用矩阵Ф同时对两边进行相乘,得到:
Rf表示f(f<n)维的实数域;
这样得到四种数字调制信号2FSK、4FSK、QPSK、16QAM的稀疏矩阵为:
y2FSK,y4FSK,yQPSK,y16QAM,并构成稀疏矩阵字典;
(2)、数字调制信号调制方式的分类
2.1)对输入待分类的数字调制信号进行步骤(1)同样的变换,得到变换结果
2.2)、计算输入待分类数字调制信号变换后结果与字典中每个稀疏矩阵的残差:
d=1的时候yd=y2FSK,d=2的时候yd=y4FSK,
d=3的时候yd=yQPSK,d=4的时候yd=y16QAM,
为了提高分类的稳定性,采用残差均值,得到其残差均值为:
其中,l表示对同一待分类数字调制信号的测试次数;
2.3)、鉴别待分类数字调试信号的调制方式:对哪种调制方式的稀疏矩阵的残差均值E[rd]最小,就确定为该调制类别。
本发明的发明目的是这样实现的:
本发明基于双谱和稀疏矩阵的数字调制信号分类方法,提出了两步法来对数字调制信号进行分类。第一步已知调制类型的数字调制信号进行双谱估计,利用压缩感知的方法对数字调制信号的双谱进行压缩,得到其稀疏表示,根据这个稀疏表示制作字典;第二步对输入待分类的数字调制信号进行第一步相同的变换,得到的变换结果与字典进行对比,计算残差均值,将残差均值最小的那个类型确定为信号的调制方式。通过与以前的基于经典特征的支持向量机(SVM)分类算法相比较,本发明提出的分类方法在分类正确率上有了明显的提高。
附图说明
图1是本发明基于双谱和稀疏矩阵的数字调制信号分类方法一种具体实施方式的原理框图;
图2是四种调制方式数字信号的双谱图;
图3是四种调制方式数字信号的不同信噪比下的分类正确率。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
实施例
图1是本发明基于双谱和稀疏矩阵的数字调制信号分类方法一种具体实施方式的原理框图;
如图1所示,本发明基于双谱和稀疏矩阵的数字调制信号分类方法首先是建立稀疏矩阵字典,用这个稀疏矩阵字典去判断是哪种调制方式。
(1)、稀疏矩阵字典的建立
步骤ST1.1:稀疏矩阵字典的建立采用确知调制类型的数字调制信号2FSK、4FSK、QPSK和16QAM分别计算它们的双谱,得到的双谱表示为: 以及
步骤ST1.2:进行稀疏变换,计算其双谱的稀疏矩阵,得到稀疏矩阵yd,并构成稀疏矩阵字典,其中,d=1,2,3,4;d=1的时候表示2FSK,d=2的时候表示4FSK,d=3的时候表示QPSK,d=4的时候表示16QAM。
(2)、数字调制信号调制方式的分类
步骤ST2.1:对于输入加噪的待分类数字调制信号,同样通过双谱变换和稀疏变换,得到其变换结果
步骤ST2.2:计算输入待分类数字调制信号变换与字典每个稀疏矩阵的残差
为了提高稳定性,采用残差均值,得到其残差均值为
与那种调制方式的字典残差最小,就确定待分类数字调制信号的调制方式为此种调制方式
图2是四种调制方式(2FSK,4FSK,QPSK,16QAM)的双谱图。
如图2所示,双谱是制作稀疏矩阵字典的一个关键中间环节,同时,可以看出,四种数字调制信号的双谱有较大差异,计算其双谱作为得到稀疏矩阵字典的前提。
图3是四种调制方式数字信号的不同信噪比下的分类正确率。
如图3所示,采用本发明方法,在不同信噪比条件下,分类正确率在信噪比高于0时都高于95%,具有较高的分类正确率。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (1)
1.一种基于双谱和稀疏矩阵的数字调制信号分类方法,其特征在于,包括以下步骤:
(1)、稀疏矩阵字典的建立
1.1)、分别计算确知的四种数字调制信号2FSK、4FSK、QPSK、16QAM的训练序列x2FSK(t),x4FSK(t),xQPSK(t),x16QAM(t)的双谱,得到其双谱表示,对每个序列进行以下处理:
将长度为N的训练序列分割成M段,计算每个分段的双谱
其中L1是预设整数,Δ0是双谱域的频率采样间隔,K=N/M,xi(·)表示数字调制信号的第i个部分;
训练序列的双谱采用均值定义,表示为
这样四种数字调制信号2FSK、4FSK、QPSK、16QAM的训练序列x2FSK(t),x4FSK(t),xQPSK(t),x16QAM(t)的双谱为:以及
1.2)、对四种数字调制信号的双谱分别进行稀疏表示,得到各自的稀疏矩阵ψ:
ψ=[v1,v2,...vn,]∈Rn
Rn是n维实数域,vj是第j个数据,j=1,2…,K;
由于维度为n相对较高,计算复杂度高,需要进行降维处理:
用矩阵Φ同时对两边进行相乘,得到:
Rf表示f(f<n)维的实数域;
这样得到四种数字调制信号2FSK、4FSK、QPSK、16QAM降维后的稀疏矩阵为:
y2FSK,y4FSK,yQPSK,y16QAM,并构成稀疏矩阵字典;
(2)、数字调制信号调制方式的分类
2.1)对输入待分类的数字调制信号进行步骤(1)同样的变换,得到变换结果
2.2)、计算输入待分类数字调制信号变换后结果与字典中每个稀疏矩阵的残差:
d=1的时候yd=y2FSK,d=2的时候yd=y4FSK,
d=3的时候yd=yQPSK,d=4的时候yd=y16QAM,
为了提高分类的稳定性,采用残差均值,得到其残差均值为:
其中,l表示对同一待分类数字调制信号的测试次数;
2.3)、鉴别待分类数字调试信号的调制方式:对哪种调制方式的稀疏矩阵的残差均值E[rd]最小,就确定为该调制类别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210330854.7A CN102868653B (zh) | 2012-09-10 | 2012-09-10 | 一种基于双谱和稀疏矩阵的数字调制信号分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210330854.7A CN102868653B (zh) | 2012-09-10 | 2012-09-10 | 一种基于双谱和稀疏矩阵的数字调制信号分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102868653A CN102868653A (zh) | 2013-01-09 |
CN102868653B true CN102868653B (zh) | 2015-08-26 |
Family
ID=47447253
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210330854.7A Expired - Fee Related CN102868653B (zh) | 2012-09-10 | 2012-09-10 | 一种基于双谱和稀疏矩阵的数字调制信号分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102868653B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107292067B (zh) * | 2017-08-17 | 2021-02-02 | 湖南纬拓信息科技有限公司 | 一种基于压缩感知与双谱分析的齿轮故障诊断方法 |
CN109255313A (zh) * | 2018-08-30 | 2019-01-22 | 中国科学院国家授时中心 | 一种提升信号识别正确率的方法 |
CN109991578B (zh) * | 2019-03-25 | 2022-05-20 | 哈尔滨工程大学 | 基于盲压缩核字典学习的多分量雷达信号调制识别方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1178601A1 (en) * | 2000-08-04 | 2002-02-06 | Motorola, Inc. | Frequency modulation using a digital filter for baseband waveshaping |
CN101136893A (zh) * | 2007-10-10 | 2008-03-05 | 天津大学 | 基于全相位fft的通用解调方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7403556B2 (en) * | 2003-06-30 | 2008-07-22 | Via Technologies Inc. | Radio receiver supporting multiple modulation formats with a single pair of ADCs |
-
2012
- 2012-09-10 CN CN201210330854.7A patent/CN102868653B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1178601A1 (en) * | 2000-08-04 | 2002-02-06 | Motorola, Inc. | Frequency modulation using a digital filter for baseband waveshaping |
CN101136893A (zh) * | 2007-10-10 | 2008-03-05 | 天津大学 | 基于全相位fft的通用解调方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102868653A (zh) | 2013-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101827060B (zh) | 基于分数阶傅里叶变换的自适应调制解调方法 | |
CN101984612B (zh) | 基于压缩感知的非连续正交频分复用信道估计方法 | |
CN103685096B (zh) | 一种基于最优导频的mimo-ofdm系统信道估计方法 | |
CN102984711A (zh) | 基于单比特压缩感知技术的多用户协作频谱感知方法 | |
CN103259760B (zh) | 基于多维星座图的脉冲波形调制方法 | |
CN102868653B (zh) | 一种基于双谱和稀疏矩阵的数字调制信号分类方法 | |
CN102724147A (zh) | 一种水声正交频分复用的信道估计方法 | |
CN102035609A (zh) | 基于复数连续全反馈神经网络的信号盲检测方法 | |
CN106209716B (zh) | 一种降低大规模mu-mimo-ofdm系统峰均功率比的方法 | |
CN106506415A (zh) | 一种多用户mimo‑ofdm系统信道估计的方法 | |
CN107528649A (zh) | 一种基于Bagging的认知无线电网络频谱的预测方法 | |
CN107018106A (zh) | 一种用于蓝牙接收机的码元波形频率偏移补偿的方法 | |
CN103595679A (zh) | 降低lte上行单载波频分多址信号峰均比的方法 | |
CN106027184A (zh) | 一种基于lte-a的链路级到系统级仿真的接口方法 | |
CN102901855B (zh) | 一种特高压直流电晕电流信号去噪方法 | |
CN103152308A (zh) | 正交频分复用系统中的频偏、直流和失衡的联合估计方法 | |
CN107094043B (zh) | 基于块迭代法的改进后的mmse低复杂度信号检测方法 | |
CN103888146B (zh) | 一种数据压缩的方法、装置和通信设备 | |
CN102291188B (zh) | 一种认知无线电系统中的合作频谱感知方法 | |
CN106209717A (zh) | 一种提升功放能效的自适应极化‑qam联合调制方法 | |
CN103220016A (zh) | 一种脉冲超宽带系统正交稀疏字典的产生系统及方法 | |
CN102111362B (zh) | 降低连续相位调制接收机复杂度的方法 | |
CN102035769B (zh) | 基于复数离散全反馈神经网络的相移键控信号盲检测方法 | |
CN106656879A (zh) | 一种高速高阶变步长自适应均衡方法 | |
CN105262506B (zh) | 一种时频重叠高斯调幅通信信号分离方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150826 Termination date: 20180910 |
|
CF01 | Termination of patent right due to non-payment of annual fee |