CN105897636A - 一种基于李雅普诺夫指数的cpfsk信号调制指数估计方法 - Google Patents

一种基于李雅普诺夫指数的cpfsk信号调制指数估计方法 Download PDF

Info

Publication number
CN105897636A
CN105897636A CN201610340841.6A CN201610340841A CN105897636A CN 105897636 A CN105897636 A CN 105897636A CN 201610340841 A CN201610340841 A CN 201610340841A CN 105897636 A CN105897636 A CN 105897636A
Authority
CN
China
Prior art keywords
signal
cpfsk
modulation index
exponent
estimation method
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
Application number
CN201610340841.6A
Other languages
English (en)
Other versions
CN105897636B (zh
Inventor
朱胜利
杜继萍
甘露
廖红舒
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201610340841.6A priority Critical patent/CN105897636B/zh
Publication of CN105897636A publication Critical patent/CN105897636A/zh
Application granted granted Critical
Publication of CN105897636B publication Critical patent/CN105897636B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0012—Modulated-carrier systems arrangements for identifying the type of modulation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/10—Frequency-modulated carrier systems, i.e. using frequency-shift keying
    • H04L27/12—Modulator circuits; Transmitter circuits
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/18—Phase-modulated carrier systems, i.e. using phase-shift keying
    • H04L27/183—Multiresolution systems
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/18—Phase-modulated carrier systems, i.e. using phase-shift keying
    • H04L27/20—Modulator circuits; Transmitter circuits
    • H04L27/2003—Modulator circuits; Transmitter circuits for continuous phase modulation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/32—Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明属于认知无线电、智能通信中信号参数估计技术领域，尤其涉及一种基于李雅普诺夫指数的连续相位移频键控(Continuous‑Phase Frequency Shift Keying,CPFSK)信号调制指数估计方法。本发明充分利用CPFSK信号的李雅普诺夫指数与其调制指数的成正比的数学关系，通过对CPFSK信号李雅普诺夫指数的计算估计其调制指数，适用于CPFSK信号的调制指数h的估计。通过仿真实验可以看出该算法结果直观，性能优异，同时拥有较低的运算复杂度。

Description

一种基于李雅普诺夫指数的CPFSK信号调制指数估计方法
技术领域
本发明属于认知无线电、智能通信中信号参数估计技术领域，尤其涉及一种基于李雅普诺夫指数的连续相位移频键控(Continuous-Phase Frequency Shift Keying,CPFSK)信号调制指数估计方法。
背景技术
CPFSK信号由于具有平滑的相位转移，因此具有比常规移频键控(FrequencyShift Keying,FSK)信号更好的频谱特性，使得CPFSK信号的带外功率小于常规FSK信号。
CPFSK调制的发送信号表示为其中，为振幅，fc为载频，φ0是载波的初始相位，表示为载波的时变相位，{In}表示幅度序列取值范围是±1,±3,...,±(M-1)，M为CPFSK信号调制进制，g(t)是一个幅度为1/2T且持续时间为T的矩形脉冲。因此，其中，h＝2fdT，h为调制指数，可以看出，CPFSK的瞬时相位是连续的。当表示为复数形式时，CPFSK信号可由式确定，信号李雅普诺夫指数表示为其中，r(t)为接收的时域信号，r′(ti)则为接收信号的微分。
认知无线电及智能通信中为了适应多种调制信号通信，需要对未知信号进行感知，获取其调制类型及其相应参数。CPFSK信号由于其较高的频谱利用效率，在现代通信中得到了广泛应用，而其中的调制指数h的估计一直是一个难点。
发明内容
针对现有技术的缺点，本发明提出一种基于李雅普诺夫指数的CPFSK信号调制指数估计方法。李雅普诺夫指数大量用于混沌信号处理研究中，可作为表征信号本质特征—相空间的一个重要指标。由于调制指数h不相等的CPFSK信号，其李雅普诺夫指数不一样。理论推导与仿真验证表明，CPFSK信号的李雅普诺夫指数与调制指数h成正比。本发明充分利用CPFSK信号的李雅普诺夫指数与其调制指数的成正比的数学关系，通过对CPFSK信号李雅普诺夫指数的计算估计其调制指数，适用于CPFSK信号的调制指数h的估计。
本发明使用的调制类型为2CPFSK、4CPFSK或8CPFSK中的任意一种。
一种基于李雅普诺夫指数的CPFSK信号调制指数估计方法，包括如下步骤：
S1、对接收信号进行采样，得到r(n)＝s(n)+n(n)，其中，s(n)为发送信号，调制指数h∈[0.2,1.5]，n(n)为加性高斯白噪声，n＝1,2,...,N，N是采样点总个数；
S2、对S1所述r(n)进行数字下变频处理，将接收信号r(n)从fc处搬移到零中频，得到使用Kaiser窗设计低通滤波器h(n)滤除所述r'(n)带外噪声得到新的信号r1(n)＝h(n)*r'(n)，fc为载波频率；
S3、对滤波后的信号r1(n)进行功率归一化得到
S4、求CPFSK信号的李雅普诺夫指数其中，r2′(n)＝fs(r2(n+1)-r2(n))，fs为采样率；
S5、更新李雅普诺夫指数，得到Le′＝Le*logRB，其中，RB为码元速率；
S6、求取CPFSK信号调制指数估计其中，thr为经验门限值，对于不同调制方式的信号thr取值不一样。
进一步地，S2所述滤波器参数设置如下：滤波器阶数为Nf，2CPFSK信号归一化通带截止频率ωp＝3RB/fs，4CPFSK信号归一化通带截止频率ωp＝7RB/fs，8CPFSK信号归一化通带截止频率ωp＝9RB/fs，归一化阻带截止频率ωs＝1.3ωp，其中，fs为采样率，RB为码元速率。
本发明的有益效果是：
本发明利用CPFSK信号的李雅普诺夫指数与其调制指数成正比的数学关系，通过求取CPFSK信号的李雅普诺夫指数，再根据先验已知的信号码速率，对李雅普诺夫指数进行修正得到新的参量，再根据不同调制阶数的CPFSK信号(本发明针对2CPFSK、4CPFSK、8CPFSK)设置不同的经验门限，有效的估计信号的调制指数。相对于传统的CPFSK信号调制指数估计方法，本发明具有过程简单、计算量小且复杂度低等优点。通过仿真实验可以看出该算法结果直观，性能优异，同时拥有较低的运算复杂度。
附图说明
图1是本发明流程图。
图2是h＝0.5时调制指数估计均方误差随着Es/N0变化的曲线图。
图3是h＝1时调制指数估计均方误差随着Es/N0变化的曲线图。
具体实施方式
下面结合附图，对本发明作进一步地详细描述。
图1是本发明基于李雅普诺夫指数的CPFSK信号调制指数估计方法的具体实施方式流程图。
图2是本发明基于李雅普诺夫指数的CPFSK信号调制指数估计方法对于以下参数设置的调制指数估计相对误差随着带内信噪比变化的曲线图。
为了描述清楚，首先对本发明中使用的参数进行定义：Es/N0表示信噪比，单位为dB；Es表示信号每个符号的能量；N0表示功率谱密度。
假设接收信号调制样式为2CPFSK、4CPFSK与8CPFSK，参数设置如下：取信噪比Es/N0从5dB到35dB，步长为2dB，采样率fs＝6MHz，符号速率为RB＝200ksps，调制指数为h＝0.5，载波频率fc＝0Hz，符号个数为6000个符号，在每个信噪比下进行100次蒙特卡洛实验。
图3是本发明基于李雅普诺夫指数的CPFSK信号调制指数估计方法对于以下参数设置的调制指数估计相对误差随着带内信噪比变化的曲线图。
假设接收信号调制样式为2CPFSK、4CPFSK与8CPFSK，参数设置如下：取信噪比Es/N0从5dB到35dB，步长为2dB，采样率fs＝4.5MHz，符号速率为RB＝150ksps，调制指数为h＝1，载波频率fc＝0Hz，符号个数为6000个符号，在每个信噪比下进行100次蒙特卡洛实验。
步骤1、对接受信号进行采样，得到r(n),n＝1,2,...,180000。
步骤2、使用Kaiser窗设计低通滤波器h(n)滤除接收信号带外噪声得到新的信号：r1(n)＝h(n)*r'(n),n＝1...180000，
对于图2与图3滤波器参数设置如下：2CPFSK信号滤波器阶数为Nf＝256，归一化通带截止频率ωp＝3RB/fs＝0.1，归一化阻带截至频率ωs＝1.3ωp＝0.13；4CPFSK信号滤波器阶数为Nf＝128，归一化通带截止频率ωp＝7RB/fs＝0.23，归一化阻带截至频率ωs＝1.3ωp＝0.303；8CPFSK信号滤波器阶数为Nf＝128，归一化通带截止频率ωp＝9RB/fs＝0.3，归一化阻带截至频率ωs＝1.3ωp＝0.39。
步骤3、r1(n)功率归一化得到r2(n),n＝1...180000，
步骤4、计算信号的李雅普诺夫指数其中，r2′(n)＝fs(r2(n)-r2(n-1)),n＝1...179999。图2中fs＝6Msps,图3中fs＝4.5Msps。
步骤5、得到新的参量Le′＝Le*logRB,图2中RB＝200ksps,图3中RB＝150ksps。
步骤6、计算CPFSK信号调制指数估计值为2CPFSK信号，thr＝0.972；4CPFSK信号，thr＝1.625；8CPFSK信号，thr＝2.270。

Claims (2)

1.一种基于李雅普诺夫指数的CPFSK信号调制指数估计方法，其特征在于，包括如下步骤：
S1、对接收信号进行采样，得到r(n)＝s(n)+n(n)，其中，s(n)为发送信号，调制指数h∈[0.2,1.5]，n(n)为加性高斯白噪声，n＝1,2,...,N，N是采样点总个数；
S2、对S1所述r(n)进行数字下变频处理，将接收信号r(n)从fc处搬移到零中频，得到使用Kaiser窗设计低通滤波器h(n)滤除所述r'(n)带外噪声得到新的信号r1(n)＝h(n)*r'(n)，fc为载波频率；
S3、对滤波后的信号r1(n)进行功率归一化得到
S4、求CPFSK信号的李雅普诺夫指数其中，r2′(n)＝fs(r2(n+1)-r2(n))，fs为采样率；
S5、更新李雅普诺夫指数，得到Le′＝Le*logRB，其中，RB为码元速率；
S6、求取CPFSK信号调制指数估计其中，thr为经验门限值，对于不同调制方式的信号thr取值不一样。
2.根据权利要求1所述的一种基于李雅普诺夫指数的CPFSK信号调制指数估计方法，其特征在于：S2所述滤波器参数设置如下：滤波器阶数为Nf，2CPFSK信号归一化通带截止频率ωp＝3RB/fs，4CPFSK信号归一化通带截止频率ωp＝7RB/fs，8CPFSK信号归一化通带截止频率ωp＝9RB/fs，归一化阻带截止频率ωs＝1.3ωp，其中，fs为采样率，RB为码元速率。
CN201610340841.6A 2016-05-19 2016-05-19 一种基于李雅普诺夫指数的cpfsk信号调制指数估计方法 Active CN105897636B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610340841.6A CN105897636B (zh) 2016-05-19 2016-05-19 一种基于李雅普诺夫指数的cpfsk信号调制指数估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610340841.6A CN105897636B (zh) 2016-05-19 2016-05-19 一种基于李雅普诺夫指数的cpfsk信号调制指数估计方法

Publications (2)

Publication Number Publication Date
CN105897636A true CN105897636A (zh) 2016-08-24
CN105897636B CN105897636B (zh) 2018-10-09

Family

ID=56716973

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610340841.6A Active CN105897636B (zh) 2016-05-19 2016-05-19 一种基于李雅普诺夫指数的cpfsk信号调制指数估计方法

Country Status (1)

Country Link
CN (1) CN105897636B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108960051A (zh) * 2018-05-28 2018-12-07 合肥工业大学 一种基于频率分析的自适应csi信号辅助滤波方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140064402A1 (en) * 2012-09-06 2014-03-06 Electronics And Telecommunications Research Institute Apparatus and method for modulation classification in wireless communication system
CN104022837A (zh) * 2014-05-28 2014-09-03 广州海格通信集团股份有限公司 一种基于特征提取的通信信号调制方式进行识别的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140064402A1 (en) * 2012-09-06 2014-03-06 Electronics And Telecommunications Research Institute Apparatus and method for modulation classification in wireless communication system
CN104022837A (zh) * 2014-05-28 2014-09-03 广州海格通信集团股份有限公司 一种基于特征提取的通信信号调制方式进行识别的方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
J.R. FONOLLOSA等: ""Estimation of the modulation index of CPM signals using higher-order statistics"", 《1993 IEEE INTERNATIONAL CONFERENCE ON ACOUSTICS, SPEECH, AND SIGNAL PROCESSING》 *
SAURABH U. PAWAR等: ""Modulation Recognition in Continuous Phase Modulation Using Approximate Entropy"", 《IEEE TRANSACTIONS ON INFORMATION FORENSICS AND SECURITY》 *
周晖: ""低信噪比下数字调制信号盲识别研究"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
彭耿等: ""角度调制信号调制指数盲估计研究"", 《电子学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108960051A (zh) * 2018-05-28 2018-12-07 合肥工业大学 一种基于频率分析的自适应csi信号辅助滤波方法

Also Published As

Publication number Publication date
CN105897636B (zh) 2018-10-09

Similar Documents

Publication Publication Date Title
CN106027432B (zh) 一种基于信号瞬时频率部分相关函数的cpfsk码速率估计方法
CN103199945B (zh) 一种低信噪比条件下认知无线电信号调制方式的识别方法
CN101478525B (zh) 一种多载波分离的方法及多载波分离装置
CN108365945A (zh) 基于两路索引调制的差分混沌移位键控调制解调器及方法
CN101860497B (zh) 一种利用改进球译码算法实现fbmc系统的均衡的方法
CN101977169A (zh) Ofdm信号的时域参数盲估计方法
CN106506428B (zh) 降低ufmc系统papr的低计算复杂度的pts方法
CN105516036B (zh) 一种cpfsk调制识别方法
CN103701492A (zh) 线性调频信号调制解调的水声跳频通信方法
CN105785324A (zh) 基于mgcstft的线性调频信号参数估计方法
CN113014523B (zh) 一种基于融合特征的复合调制信号盲识别方法
CN104980187A (zh) 一种信号去噪处理的方法及装置
CN105763499B (zh) 一种基于信号瞬时频率滑动平均的cpfsk调制识别方法
Wang et al. Frequency domain multi-carrier modulation based on prolate spheroidal wave functions
Jayan et al. Performance analysis of filtered OFDM for 5G
CN101447967A (zh) 一种用于ofdm/oqam系统中的导频序列结构及其构造方法
CN105897636A (zh) 一种基于李雅普诺夫指数的cpfsk信号调制指数估计方法
Hu et al. Preamble design with ICI cancellation for channel estimation in OFDM/OQAM system
CN102629859B (zh) 一种用于窄带信号的数字带通滤波器及滤波方法
CN107707499B (zh) 一种Alpha稳定分布噪声下OFDM信号调制参数估计方法
CN109039379B (zh) 电力线通信信号滤波方法
CN107566311B (zh) 基于资源块滤波rb f-ofdm系统的传输方法
CN107204952A (zh) 一种滤波ofdm系统的子带滤波处理方法
EP3033864A1 (en) Fine timing
Kosmowski et al. A novel OFDM sensing method based on CAF-max for hybrid detectors architecture

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant