CN102868653A - 一种基于双谱和稀疏矩阵的数字调制信号分类方法 - Google Patents
一种基于双谱和稀疏矩阵的数字调制信号分类方法 Download PDFInfo
- Publication number
- CN102868653A CN102868653A CN2012103308547A CN201210330854A CN102868653A CN 102868653 A CN102868653 A CN 102868653A CN 2012103308547 A CN2012103308547 A CN 2012103308547A CN 201210330854 A CN201210330854 A CN 201210330854A CN 102868653 A CN102868653 A CN 102868653A
- Authority
- CN
- China
- Prior art keywords
- digital modulation
- centerdot
- sparse matrix
- modulation signals
- lambda
- 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.)
- Granted
Links
Images
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)的双谱,得到其双谱表示,对每个序列进行以下处理:
ψ=[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.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.2:进行稀疏变换,计算其双谱的稀疏矩阵,得到稀疏矩阵yd,并构成稀疏矩阵字典,其中,d=1,2,3,4;d=1的时候表示2FSK,d=2的时候表示4FSK,d=3的时候表示QPSK,d=4的时候表示16QAM。
(2)、数字调制信号调制方式的分类
步骤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)的双谱,得到其双谱表示,对每个序列进行以下处理:
ψ=[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)同样的变换,得到变换结果
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 true CN102868653A (zh) | 2013-01-09 |
CN102868653B 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) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107292067A (zh) * | 2017-08-17 | 2017-10-24 | 湖南纬拓信息科技有限公司 | 一种基于压缩感知与双谱分析的齿轮故障诊断方法 |
CN109255313A (zh) * | 2018-08-30 | 2019-01-22 | 中国科学院国家授时中心 | 一种提升信号识别正确率的方法 |
CN109991578A (zh) * | 2019-03-25 | 2019-07-09 | 哈尔滨工程大学 | 基于盲压缩核字典学习的多分量雷达信号调制识别方法 |
Citations (3)
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 |
US20040264600A1 (en) * | 2003-06-30 | 2004-12-30 | Ping-Chieh Kao | Radio Receiver Supporting Multiple Modulation Formats with a Single Pair of ADCs |
CN101136893A (zh) * | 2007-10-10 | 2008-03-05 | 天津大学 | 基于全相位fft的通用解调方法 |
-
2012
- 2012-09-10 CN CN201210330854.7A patent/CN102868653B/zh not_active Expired - Fee Related
Patent Citations (3)
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 |
US20040264600A1 (en) * | 2003-06-30 | 2004-12-30 | Ping-Chieh Kao | Radio Receiver Supporting Multiple Modulation Formats with a Single Pair of ADCs |
CN101136893A (zh) * | 2007-10-10 | 2008-03-05 | 天津大学 | 基于全相位fft的通用解调方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107292067A (zh) * | 2017-08-17 | 2017-10-24 | 湖南纬拓信息科技有限公司 | 一种基于压缩感知与双谱分析的齿轮故障诊断方法 |
CN109255313A (zh) * | 2018-08-30 | 2019-01-22 | 中国科学院国家授时中心 | 一种提升信号识别正确率的方法 |
CN109991578A (zh) * | 2019-03-25 | 2019-07-09 | 哈尔滨工程大学 | 基于盲压缩核字典学习的多分量雷达信号调制识别方法 |
CN109991578B (zh) * | 2019-03-25 | 2022-05-20 | 哈尔滨工程大学 | 基于盲压缩核字典学习的多分量雷达信号调制识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102868653B (zh) | 2015-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107124381B (zh) | 一种数字通信信号调制方式自动识别方法 | |
CN101984612B (zh) | 基于压缩感知的非连续正交频分复用信道估计方法 | |
CN102984711A (zh) | 基于单比特压缩感知技术的多用户协作频谱感知方法 | |
CN104052702A (zh) | 一种复杂噪声下数字调制信号的识别方法 | |
CN106330805B (zh) | 一种信号调制方式自动识别方法及系统 | |
CN104333406B (zh) | 用于在多小区多用户多天线系统中的解码方法 | |
CN101764682A (zh) | 一种调制编码配置的选择方法和装置 | |
CN106357574A (zh) | 基于顺序统计量的bpsk/qpsk信号调制盲识别方法 | |
CN102868653A (zh) | 一种基于双谱和稀疏矩阵的数字调制信号分类方法 | |
CN106027184A (zh) | 一种基于lte-a的链路级到系统级仿真的接口方法 | |
CN104868962B (zh) | 基于压缩感知的频谱检测方法及装置 | |
CN102291188B (zh) | 一种认知无线电系统中的合作频谱感知方法 | |
CN104408027A (zh) | 一种基于广义协方差和张量分解的欠定盲辨识方法 | |
CN102868654B (zh) | 一种认知网络中数字调制信号的分类方法 | |
CN109309538A (zh) | 一种频谱感知方法、装置、设备、系统及存储介质 | |
CN101308651B (zh) | 音频暂态信号的检测方法 | |
CN104507106A (zh) | 8psk信号和pi/4-dqpsk信号的识别方法 | |
CN104022839B (zh) | 一种应用于合作频谱感知的信息融合判决方法 | |
CN101754343A (zh) | 一种信道传输质量的反馈方法、系统和设备 | |
CN116506273A (zh) | 新型mpsk调制信号识别分类方法 | |
CN104753842B (zh) | 基于峰值位置判别的信号调制方式识别方法 | |
CN105262506B (zh) | 一种时频重叠高斯调幅通信信号分离方法 | |
CN102412911A (zh) | 两极频谱检测方法 | |
CN106850096A (zh) | 空间相关性估算的方法及装置 | |
CN105245300B (zh) | 一种underlay频谱共享方式下用户个数估计方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150826 Termination date: 20180910 |