CN109450829B - 数字调制信号码速率估计方法和装置 - Google Patents
数字调制信号码速率估计方法和装置 Download PDFInfo
- Publication number
- CN109450829B CN109450829B CN201811350243.2A CN201811350243A CN109450829B CN 109450829 B CN109450829 B CN 109450829B CN 201811350243 A CN201811350243 A CN 201811350243A CN 109450829 B CN109450829 B CN 109450829B
- Authority
- CN
- China
- Prior art keywords
- frequency
- sequence
- instantaneous
- code rate
- signal
- 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.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0262—Arrangements for detecting the data rate of an incoming signal
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明公开了数字调制信号码速率估计方法和装置,方法包括:对以预设采样率fs进行采样后的信号提取瞬时幅度A(l)和瞬时频率f(l);进行Haar小波变换后再分别进行离散傅里叶变换;提取脉冲尖峰频率位置并对所得位置进行排序;对序列做差分,分别对序列中出现次数最多的数值及其对应的票数分别记录;按照预设条件改变小波变换尺度因子,并重复步骤直到小波变换的变换尺度a大于数据长度;分别确定瞬时幅度A(l)和瞬时频率f(l)中所有数值中得票最多的数值及其对应的票数;比较两个数值的得票数,最终确定码速率。本方法不依赖于先验知识,在较低信噪比的情况下识别正确率高;提出的“票数统计”策略能够有效解决小波变换尺度选择的问题,具有自适应性。
Description
技术领域
本发明涉及一种数字调制信号码速率估计方法和装置,属于通信技术领域。
背景技术
通信信号调制识别目的就是对截获的一段通信信号,依据较少的先验信息,判断出信号的调制方式。随着通信技术的飞速发展,通信信号的体制和调制样式变得更加复杂多样,信号环境日趋密集,这使得通信信号的识别变得越加困难。通信信号中包含两大类调制信号:数字调制信号和模拟调制信号,利用判决理论时,模拟调制信号和数字调制信号的区分门限往往难以设置。由于模拟调制信号没***速率之说,其估计结果为任意无规律值,因此利用码速率估计的方法可以有效区分模拟和数字调制,码速率是区分两种调制方式的重要特征。
码速率对于数字调制信号特征提取和分析以及识别和解调有着重要的意义,在电子侦察中,数字调制信号通过正交变频和低通滤波后得到系带信号,为了分析与识别信号,需要知道数字调制信号的码速率。因此正确估计数字调制信号的码速率是十分必要的。
估计数字调制信号码速率的方法有多种方法,其核心思想是根据调制信号码元之间存在突变,根据突变之间的间隔来估计码速率。对于数字调制基带信号,可以运用延迟相乘的方法来获得信号的码速率。经过延迟相乘之后,输出信号的双边功率谱密度主要包含三项,第一项是直流分量,第二项是码速率及高次谐波,第三项为连续谱。可以根据第二项对调制信号的码速率进行估计。但是时延长度必须大于码元宽度的一半,这在实际的应用中很难保证。
论文《无线电信号的调制识别研究》(作者王生兵)公开了利用承载数字信息的幅度调制波形p(n)进行码速率估计。对于ASK信号,p(n)就是信号的瞬时幅度;对于PSK信号,p(n)就是信号的瞬时相位;对于FSK信号,p(n)就是信号的瞬时频率。该方法的流程可概括如下:
1.对p(n)求L阶差分序列,并求取绝对值,得到p′(n)序列,L选为15;
2.计算p′(n)序列的均值mp,以3mp为门限对p′(n)序列进行判决;
3.当p′(n)>3mp,判决器输出1;否则,判决器输出0;得到阶跃位置序列l(n)
4.以阶跃带的中心作为符号转换点,得到符号转换位置序列lc(n),在对lc(n)取一阶差分,得到差分序列lcd(n),其长度为Nd;
该方法的不足之处在于码速率估计方法中的门限通过经验取值,缺乏理论依据,门限的取值对不同参数的调制信号不具有适应性;
除此之外,因小波变换对瞬变信号有着很好的检测性能,基于小波变换的码速率估计是另外一种思路。文献【基于VSSR法的数字与模拟调制信号的区分-戚玉鹏】对该方法进行了详细的描述,其码速率估计方法可概括为如下:
1.对采集的信号s(t)进行小波变换,得到A(t);
2.对A(t)小波变换后取包络值,得到B(t);
3.对B(t)求其傅里叶变换,得到C(f);
4.求得|C(f)|中的相邻两尖峰之间的距离d,该距离d与|C(f)|的点个数的比值,即为符号率与信号的采样率的比值,从而可求得信号的符号率。
基于小波变换的方法,对于小波变换的尺度选择没有给出说明,如果对采集信号的码速率完全未知,很难选择合适的小波变换的尺度,因而导致码速率估计错误;
此外现***速率估计方法对各类数字调制信号分析的不够全面,一般是对某种特定调制方式的信号进行码速率估计,缺乏一种统一的码速率估计方法且可能依赖于先验知识,如需要预先知道调制方式(幅度、相位和频率)、信号的载频和载波相位等,实际情况很难准确获知这些参数;
发明内容
本发明的目的是为了克服目前码速率估计技术的不足,提出一种无需先验知识的基于小波变换的数字调制信号码速率估计方法,其自适应程度高。
为解决上述技术问题,本发明采用了以下的技术方案:
数字调制信号码速率估计方法,包括:
步骤1:对以预设采样率fs及总采样点数L进行采样后的信号提取瞬时幅度A(l)和瞬时频率f(l);
步骤2:按照设定的小波变换尺度分别对瞬时幅度A(l)和瞬时频率f(l)进行小波变换,小波变换母函数选用Haar小波,得到CA(l)和Cf(l);
步骤3:对CA(l)进行离散傅里叶变换,得到YA(fk);对频率f(l)进行离散傅里叶变换,得到Yf(fk);
步骤4:提取YA(fk)和Yf(fk)的脉冲尖峰频率位置并对所得位置进行排序,分别得到序列xA(N1)和序列xf(N2),N1和N2为正整数;
步骤5:对序列xA(N1)做差分,得到DxA(N1)序列,对DxA(N1)中的数值进行统计:选出DxA(N1)中重复次数最多的数值,记为XAi,相应的重复次数为NAi;对序列xf(N2)做差分,得到Dxf(N2)序列,对Dxf(N2)中的数值进行统计:选出Dxf(N2)中重复次数最多的数值,记为Xfj,相应的重复次数为Nfj;其中下角标i和j都是整数且初始值为零;
步骤6:按照预设条件改变小波变换尺度因子m,并令i等于i加1,j等于j加1,并重复步骤2~5,直到小波变换的变换尺度a大于数据长度,并记录此时的i和j均等于M;
步骤7:挑选出NA0,NA1,NA2,…NAM的最大值,并记录其下标为ma,ma为0到M中的某一整数,把XAma记为RSA,NAma记为票数PA;挑选出Nf0,Nf1,Nf2,…NfM的最大值,并记录其下标为mf,mf为0到M中的某一整数,把Xfmf记为RSf,Nfmf记为票数Pf;
步骤8:若PA大于Pf,则码速率RS为RSA与采样率fs的乘积与总采样点数L的比值;否则,则码速率RS为RSf与采样率fs的乘积与总采样点数L的比值。
进一步地,提取瞬时幅度A(l)和瞬时频率f(l)包括:
对于连续波调制,已调实信号的表达式为:
对信号进行希尔伯特变换,得到信号的正交分量,表达式如下:
得到信号s(t)的解析表达式z(t)=s(t)+jv(t)。
再进一步地,提取提取瞬时幅度A(l),方法如下:
对信号进行采样后,得到离散形式的瞬时幅度A(l),l=1,2,…L。
进一步地,提取瞬时频率f(l)的方法:瞬时相位提取,瞬时相位表达式如下:
去相位折叠,首先计算修正相位序列C(l),l=1,2,…L:
则无折叠的相位为:
瞬时频率提取,表达式为
其中fs为采样率。
进一步地,小波变换尺度因子为m,m取不小于零的整数,则小波变换的变换尺度a的表达式为:a=2m。
再进一步地,设定初始小波变换尺度因子m为零。
进一步地,改变小波变换尺度因子时令小波变换尺度因子m等于m加1。
进一步地,脉冲尖峰频率位置的提取方法如下:
|Y(fk)|在码速率对应的频率及码速率倍频位置会出现脉冲尖峰,如果满足以下表达式:
|Y(fk)|>|Y(fk-1)|&|Y(fk)|>|Y(fk+1)| (9),
则表明在频率fk处存在脉冲尖峰。
在另一方面本发明提供了数字调制信号码速率估计装置,包括:
采样单元:用于根据预设采样率fs及总采样点数L对信号进行采样;
瞬时幅度和瞬时频率提取单元:用于对采用后的信号进行瞬时幅度A(l)和瞬时频率f(l)的提取;
小波变换单元:用于按照设定的小波变换尺度分别对瞬时幅度A(l)和瞬时频率f(l)进行小波变换,小波变换母函数选用Haar小波,得到CA(l)和Cf(l);
离散傅里叶变换单元:用于对CA(l)进行离散傅里叶变换,得到YA(fk);对频率f(l)进行离散傅里叶变换,得到Yf(fk);脉冲尖峰频率位置提取单元:用于提取YA(fk)和Yf(fk)的脉冲尖峰频率位置并对所得位置进行排序,分别得到序列xA(N1)和序列xf(N1),N1和N2为正整数;
差分单元:用于对序列xA(N1)做差分,得到DxA(N1)序列以及对序列xf(N2)做差分,得到Dxf(N2)序列;
票数统计单元:用于对DxA(N1)中的数值进行统计:选出DxA(N1)中重复次数最多的数值,记为XAm,相应的重复次数为NAm;并且对Dxf(N2)中的数值进行统计:选出Dxf(N2)中重复次数最多的数值,记为Xfm,相应的重复次数为Nfm;
其中下角标m是每次循环都不同的变量;
小波变换尺度因子设定单元:用于设定的小波变换尺度并且按照预设条件改变小波变换尺度因子直到小波变换的变换尺度a大于数据长度,并记录M等于此时的小波变换尺度因子m,m取不小于零的整数;
码速率计算单元:用于挑选出NA0,NA1,NA2,…NAM的最大值,并记录其下标为ma,ma为0到M中的某一整数,把XAma记为RSA,NAma记为票数PA;
挑选出Nf0,Nf1,Nf2,…NfM的最大值,并记录其下标为mf,mf为0到M中的某一整数,把Xfmf记为RSf,Nfmf记为票数Pf;
若PA大于Pf,则码速率RS为RSA与采样率fs的乘积与总采样点数L的比值;否则,则码速率RS为RSf与采样率fs的乘积与总采样点数L的比值。
本发明所达到的有益效果:
1.数字调制信号的种类多,对常见信号都进行了全面分析和码速率估计;
2.不依赖于先验知识,包括信号调制方式、信噪比、载频、载波相位和码元同步等,在对信号完全未知的情况下进行码速率估计,并且在较低信噪比的情况下识别正确率高;
3.本发明提出的“票数统计”策略能够有效解决小波变换尺度选择的问题;
4.具有自适应性,对于数字调制信号的不同码速率都能够自动选择合适的尺度;算法规则简单且收敛,码速率的估计可以在很短时间内完成,适合在线分析。
附图说明
图1是本发明方法流程示意图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
图1为本发明方法流程示意图,图1示出了本发明数字调制信号码速率估计方法,
步骤1:对以预设采样率fs及总采样点数L进行采样后的信号提取瞬时幅度A(l)和瞬时频率f(l),l为采样后的离散自变量;因为信号未知的情况下,无法确定码元调制在幅度上还是频率上,因此本发明方法分幅度和频率处理:
具体方法包括:
对于连续波调制,已调实信号的一般表达式为:
(1).对信号进行希尔伯特变换,得到信号的正交分量。
希尔伯特变换实际是对原信号进行了π/2相移,使它成为自身的正交对。可得信号s(t)的解析表达式z(t)=s(t)+jv(t)。
(2).瞬时幅度提取
信号包络瞬时幅度估计:
(3).瞬时相位提取
瞬时相位估计:
(5).去相位折叠
由于实际相位计算时是以模2π来计算的,因此在求得瞬时相位之后需要对相位去折叠,恢复出没有折叠的相位,为此首先计算修正相位序列C(l),l=1,2,…L:
则无折叠的相位为:
(6).去载频线性相位
(7).瞬时频率提取
瞬时频率可以由瞬时相位的导数计算求得,对数字信号可直接求差分,即:
步骤2:按照设定的小波变换尺度分别对瞬时幅度A(l)和瞬时频率f(l)进行小波变换,由于Haar小波简单、计算快和最常用,因此小波变换母函数选用Haar小波,得到CA(l)和Cf(l);
小波变换尺度因子为m,m取不小于零的整数,则小波变换的变换尺度a的表达式为:a=2m。
从图1可以看出设定初始小波变换尺度时设小波变换尺度因子为m为零。
步骤3:对CA(l)进行离散傅里叶变换,得到YA(fk);对频率f(l)进行离散傅里叶变换,得到Yf(fk);
步骤4:提取YA(fk)和Yf(fk)的脉冲尖峰频率位置并对所得位置进行排序,分别得到序列xA(N1)和序列xf(N2),N1和N2为正整数;
YA(fk)和Yf(fk)在码速率对应的频率及码速率倍频位置会出现脉冲尖峰,记录尖峰对应频率位置;在对所得位置进行排序时升序和降序均可。
步骤5:对序列xA(N1)做差分,得到DxA(N1)序列,对DxA(N1)中的数值进行统计:选出DxA(N1)中重复次数最多的数值,记为XAi,相应的重复次数为NAi;对序列xf(N2)做差分,得到Dxf(N2)序列,对Dxf(N2)中的数值进行统计:选出Dxf(N2)中重复次数最多的数值,记为Xfj,相应的重复次数为Nfj;其中下角标i和j都是整数且初始值为零;
步骤6:按照预设条件改变小波变换尺度因子m,并令i等于i加1,j等于j加1,并重复步骤2~5,直到小波变换的变换尺度a大于数据长度,并记录此时的i和j均等于M;
优选地,改变小波变换尺度因子时令小波变换尺度因子m等于m加1。
步骤7:挑选出NA0,NA1,NA2,…NAM的最大值,并记录其下标为ma,ma为0到M中的某一整数,把XAma记为RSA,NAma记为票数PA;挑选出Nf0,Nf1,Nf2,…NfM的最大值,并记录其下标为mf,mf为0到M中的某一整数,把Xfmf记为RSf,Nfmf记为票数Pf;
步骤8:若PA大于Pf,则码速率RS为RSA与采样率fs的乘积与总采样点数L的比值;否则,则码速率RS为RSf与采样率fs的乘积与总采样点数L的比值。
小波变换原理
称ψ(t)为一个基本小波或者母小波。将母函数经过平移和伸缩后得:
称其为一个小波序列。其中a为伸缩因子,b为平移因子。对于任意的函数x(t)∈L2(R)的连续小波变换为:
对应采样后信号的小波变换为:
本发明中,x(l)对应于瞬时幅度A(l)和瞬时频率f(l)。
由于基小波在小波变换中对被分析的信号起着观测窗的作用,所以ψ(t)还应该满足一般函数的约束条件:
(2).Haar母小波
Haar母小波的数学表达式为
脉冲尖峰频率位置的提取方法
|Y(fk)|在码速率对应的频率及码速率倍频位置会出现脉冲尖峰,如果满足以下规则:
|Y(fk)|>|Y(fk-1)|&|Y(fk)|>|Y(fk+1)| (17)
则表明在频率fk处存在脉冲尖峰,‘&’表示“与门”。为了消除噪声起伏带来的影响,在算法实现中,仅考虑3dB带宽对应频率内的脉冲尖峰。
在另一方面本发明提供数字调制信号码速率估计装置,包括:
采样单元:用于根据预设采样率fs及总采样点数L对信号进行采样;
瞬时幅度和瞬时频率提取单元:用于对采用后的信号进行瞬时幅度A(l)和瞬时频率f(l)的提取;
小波变换单元:用于按照设定的小波变换尺度分别对瞬时幅度A(l)和瞬时频率f(l)进行小波变换,小波变换母函数选用Haar小波,得到CA(l)和Cf(l);
离散傅里叶变换单元:用于对CA(l)进行离散傅里叶变换,得到YA(fk);对频率f(l)进行离散傅里叶变换,得到Yf(fk);脉冲尖峰频率位置提取单元:用于提取YA(fk)和Yf(fk)的脉冲尖峰频率位置并对所得位置进行排序,分别得到序列xA(N1)和序列xf(N1),N1和N2为正整数;
差分单元:用于对序列xA(N1)做差分,得到DxA(N1)序列以及对序列xf(N2)做差分,得到Dxf(N2)序列;
票数统计单元:用于对DxA(N1)中的数值进行统计:选出DxA(N1)中重复次数最多的数值,记为XAm,相应的重复次数为NAm;并且对Dxf(N2)中的数值进行统计:选出Dxf(N2)中重复次数最多的数值,记为Xfm,相应的重复次数为Nfm;
其中下角标m是每次循环都不同的变量;
小波变换尺度因子设定单元:用于设定的小波变换尺度并且按照预设条件改变小波变换尺度因子直到小波变换的变换尺度a大于数据长度,并记录M等于此时的小波变换尺度因子m,m取不小于零的整数;
码速率计算单元:用于挑选出NA0,NA1,NA2,…NAM的最大值,并记录其下标为ma,ma为0到M中的某一整数,把XAma记为RSA,NAma记为票数PA;
挑选出Nf0,Nf1,Nf2,…NfM的最大值,并记录其下标为mf,mf为0到M中的某一整数,把Xfmf记为RSf,Nfmf记为票数Pf;
若PA大于Pf,则码速率RS为RSA与采样率fs的乘积与总采样点数L的比值;否则,则码速率RS为RSf与采样率fs的乘积与总采样点数L的比值。
进一步地,所述小波变换单元设定初始小波变换尺度因子m为零,小波变换的变换尺度a的表达式为:a=2m。
在具体实施例中实际通信***中常见的数字调制信号2ASK,4ASK,BPSK,QPSK,8PSK,OQPSK,PI/4-DQPSK,2FSK,4FSK,16QAM,64QAM为例,进行码速率的自动估计计算机仿真,使用随机序列作为数字调制信号的信源。用均值为0的正态分布的随机数发生器产生加性高斯白噪声AWGN序列。信噪比SNR变化范围从5dB到15dB,步进为1dB。每个识别样本用20000个点,每个识别样本仿真400次。仿真平台操作***是WIN7,采用软件MATLAB进行编程。仿真表明,当信噪比达到6dB时,码速率的估计准确率高达99%。
现代通信***中最常用的是数字调制技术,因此,数字调制方式的识别以及关键参数码速率的估计对通信信号的分析以及特征提取有着至关重要的作用。本发明提出的基于小波变换的数字调制信号码速率估计方法无需任何先验知识,并且根据“投票”策略能够自适应确定小波变换的尺度,算法规则简单且收敛。对常见数字调制信号的仿真结果表明,在信噪比较低的情况下,码速率的估计正确率仍能达到99%以上
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (10)
1.数字调制信号码速率估计方法,其特征是,包括:
步骤1:对以预设采样率fs及总采样点数L进行采样后的信号提取瞬时幅度A(l)和瞬时频率f(l);
步骤2:按照设定的小波变换尺度分别对瞬时幅度A(l)和瞬时频率f(l)进行小波变换,小波变换母函数选用Haar小波,得到CA(l)和Cf(l);
步骤3:对CA(l)进行离散傅里叶变换,得到YA(fk);对频率f(l)进行离散傅里叶变换,得到Yf(fk);
步骤4:提取YA(fk)和Yf(fk)的脉冲尖峰频率位置并对所得位置进行排序,分别得到序列xA(N1)和序列xf(N2),N1和N2为正整数;
步骤5:对序列xA(N1)做差分,得到DxA(N1)序列,对DxA(N1)中的数值进行统计:选出DxA(N1)中重复次数最多的数值,记为XAi,相应的重复次数为NAi;对序列xf(N2)做差分,得到Dxf(N2)序列,对Dxf(N2)中的数值进行统计:选出Dxf(N2)中重复次数最多的数值,记为Xfj,相应的重复次数为Nfj;其中下角标i和j都是整数且初始值为零;
步骤6:按照预设条件改变小波变换尺度因子m,并令i等于i加1,j等于j加1,并重复步骤2~5,直到小波变换的变换尺度a大于数据长度,并记录此时的i和j均等于M;
步骤7:挑选出NA0,NA1,NA2,…NAM的最大值,并记录其下标为ma,ma为0到M中的某一整数,把XAma记为RSA,NAma记为票数PA;挑选出Nf0,Nf1,Nf2,…NfM的最大值,并记录其下标为mf,mf为0到M中的某一整数,把Xfmf记为RSf,Nfmf记为票数Pf;
步骤8:若PA大于Pf,则码速率RS为RSA与采样率fs的乘积与总采样点数L的比值;否则,则码速率RS为RSf与采样率fs的乘积与总采样点数L的比值。
5.根据权利要求1所述的数字调制信号码速率估计方法,其特征是,小波变换尺度因子为m,m取不小于零的整数,则小波变换的变换尺度a的表达式为:a=2m。
6.根据权利要求5所述的数字调制信号码速率估计方法,其特征是,设定初始小波变换尺度因子m为零。
7.根据权利要求5所述的数字调制信号码速率估计方法,其特征是,改变小波变换尺度因子时令小波变换尺度因子m等于m加1。
8.根据权利要求1所述的数字调制信号码速率估计方法,其特征是,脉冲尖峰频率位置的提取方法如下:
|Y(fk)|在码速率对应的频率及码速率倍频位置会出现脉冲尖峰,如果满足以下表达式:
|Y(fk)|>|Y(fk-1)|&|Y(fk)|>|Y(fk+1)| (9),
则表明在频率fk处存在脉冲尖峰。
9.数字调制信号码速率估计装置,其特征在于,包括:
采样单元:用于根据预设采样率fs及总采样点数L对信号进行采样;
瞬时幅度和瞬时频率提取单元:用于对采用后的信号进行瞬时幅度A(l)和瞬时频率f(l)的提取;
小波变换单元:用于按照设定的小波变换尺度分别对瞬时幅度A(l)和瞬时频率f(l)进行小波变换,小波变换母函数选用Haar小波,得到CA(l)和Cf(l);
离散傅里叶变换单元:用于对CA(l)进行离散傅里叶变换,得到YA(fk);对频率f(l)进行离散傅里叶变换,得到Yf(fk);脉冲尖峰频率位置提取单元:用于提取YA(fk)和Yf(fk)的脉冲尖峰频率位置并对所得位置进行排序,分别得到序列xA(N1)和序列xf(N1),N1和N2为正整数;
差分单元:用于对序列xA(N1)做差分,得到DxA(N1)序列以及对序列xf(N2)做差分,得到Dxf(N2)序列;
票数统计单元:用于对DxA(N1)中的数值进行统计:选出DxA(N1)中重复次数最多的数值,记为XAm,相应的重复次数为NAm;并且对Dxf(N2)中的数值进行统计:选出Dxf(N2)中重复次数最多的数值,记为Xfm,相应的重复次数为Nfm;
其中下角标m是每次循环都不同的变量;
小波变换尺度因子设定单元:用于设定的小波变换尺度并且按照预设条件改变小波变换尺度因子直到小波变换的变换尺度a大于数据长度,并记录M等于此时的小波变换尺度因子m,m取不小于零的整数;
码速率计算单元:用于挑选出NA0,NA1,NA2,…NAM的最大值,并记录其下标为ma,ma为0到M中的某一整数,把XAma记为RSA,NAma记为票数PA;
挑选出Nf0,Nf1,Nf2,…NfM的最大值,并记录其下标为mf,mf为0到M中的某一整数,把Xfmf记为RSf,Nfmf记为票数Pf;
若PA大于Pf,则码速率RS为RSA与采样率fs的乘积与总采样点数L的比值;否则,则码速率RS为RSf与采样率fs的乘积与总采样点数L的比值。
10.根据权利要求9所述的数字调制信号码速率估计装置,其特征在于,包括:所述小波变换单元设定初始小波变换尺度因子m为零,小波变换的变换尺度a的表达式为:a=2m。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811350243.2A CN109450829B (zh) | 2018-11-14 | 2018-11-14 | 数字调制信号码速率估计方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811350243.2A CN109450829B (zh) | 2018-11-14 | 2018-11-14 | 数字调制信号码速率估计方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109450829A CN109450829A (zh) | 2019-03-08 |
CN109450829B true CN109450829B (zh) | 2020-12-25 |
Family
ID=65551756
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811350243.2A Active CN109450829B (zh) | 2018-11-14 | 2018-11-14 | 数字调制信号码速率估计方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109450829B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110690931B (zh) * | 2019-10-11 | 2020-12-22 | 中国科学院软件研究所 | 一种基于多小波基联合的数字信号自适应码率估计方法和装置 |
CN113259027B (zh) * | 2021-05-12 | 2022-09-06 | 深圳华创电科技术有限公司 | 一种基于Haar变换的敌我识别信号码速率计算方法 |
CN116032709B (zh) * | 2022-12-06 | 2024-04-12 | 中国电子科技集团公司第三十研究所 | 无先验知识fsk信号盲解调和调制特征解析方法及装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103441964B (zh) * | 2013-08-30 | 2016-05-18 | 清华大学 | 一种基于随机共振的通信信号参数估计方法与装置 |
CN103490824A (zh) * | 2013-09-04 | 2014-01-01 | 中国电子科技集团公司第四十一研究所 | 一种evm分析的参考信号恢复方法 |
CN103841600B (zh) * | 2014-03-13 | 2017-02-22 | 中国人民解放军信息工程大学 | 信号码元速率估计方法及传感器网络电磁频谱监测*** |
KR101611534B1 (ko) * | 2015-08-28 | 2016-04-11 | 한화탈레스 주식회사 | 심볼율 추정 방법 |
CN106209701B (zh) * | 2016-07-01 | 2019-08-02 | 中国人民解放军信息工程大学 | Alpha稳定分布噪声环境下MFSK信号符号速率估计方法及装置 |
CN108270700A (zh) * | 2016-12-30 | 2018-07-10 | 中国航天科工集团八五研究所 | 一种改进的数字信号码元速率特征提取算法 |
CN108282428A (zh) * | 2018-01-23 | 2018-07-13 | 大连理工大学 | 一种无需先验知识的常用通信信号调制方式的自动识别方法 |
-
2018
- 2018-11-14 CN CN201811350243.2A patent/CN109450829B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109450829A (zh) | 2019-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wong et al. | Specific emitter identification using convolutional neural network-based IQ imbalance estimators | |
CN109450829B (zh) | 数字调制信号码速率估计方法和装置 | |
CN107948107B (zh) | 一种基于联合特征的数字调制信号分类方法 | |
CN105429719B (zh) | 基于功率谱和多尺度小波变换分析强干扰信号检测方法 | |
CN110113278B (zh) | 一种基于全数字接收机的调制方式识别方法 | |
CN110690931B (zh) | 一种基于多小波基联合的数字信号自适应码率估计方法和装置 | |
CN112511477A (zh) | 一种基于星座图和深度学习的混合卫星通信调制识别方法及*** | |
CN111935046B (zh) | 一种低复杂度的频移键控信号符号率估计方法 | |
CN113037663B (zh) | 一种改进的适用于非恒包络信号的码元速率估计方法 | |
CN104363194A (zh) | 基于波形变换的psk调制识别方法 | |
CN105785324A (zh) | 基于mgcstft的线性调频信号参数估计方法 | |
CN114422311B (zh) | 联合深度神经网络和专家先验特征的信号调制识别方法及*** | |
CN116257752A (zh) | 一种信号调制样式识别方法 | |
Majhi et al. | Novel blind modulation classification of circular and linearly modulated signals using cyclic cumulants | |
Wong et al. | Emitter identification using CNN IQ imbalance estimators | |
Haq et al. | Recognition of digital modulated signals based on statistical parameters | |
CN108737302A (zh) | 低信噪比条件下随机共振联合小波变换的符号速率估计方法及其装置 | |
CN109274626B (zh) | 一种基于星座图正交扫描特征的调制识别方法 | |
Hassanpour et al. | A robust algorithm based on wavelet transform for recognition of binary digital modulations | |
CN115460048B (zh) | 一种基于码元速率的msk调制识别方法、介质及装置 | |
CN116032709B (zh) | 无先验知识fsk信号盲解调和调制特征解析方法及装置 | |
CN109660475B (zh) | 一种非合作相位编码水声通信信号自主辨识方法 | |
CN116506273A (zh) | 新型mpsk调制信号识别分类方法 | |
CN104270328B (zh) | 一种信噪比实时估计方法 | |
CN115664905A (zh) | 基于多域物理层指纹特征的Wi-Fi设备识别***及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |