数字调制信号码速率估计方法和装置制造方法及图纸

技术编号:20550072 阅读:18 留言:0更新日期:2019-03-09 22:21
本发明专利技术公开了数字调制信号码速率估计方法和装置,方法包括:对以预设采样率fs进行采样后的信号提取瞬时幅度A(l)和瞬时频率f(l);进行Haar小波变换后再分别进行离散傅里叶变换;提取脉冲尖峰频率位置并对所得位置进行排序;对序列做差分,分别对序列中出现次数最多的数值及其对应的票数分别记录;按照预设条件改变小波变换尺度因子,并重复步骤直到小波变换的变换尺度a大于数据长度;分别确定瞬时幅度A(l)和瞬时频率f(l)中所有数值中得票最多的数值及其对应的票数;比较两个数值的得票数,最终确定码速率。本方法不依赖于先验知识,在较低信噪比的情况下识别正确率高;提出的“票数统计”策略能够有效解决小波变换尺度选择的问题,具有自适应性。

【技术实现步骤摘要】
数字调制信号码速率估计方法和装置
本专利技术涉及一种数字调制信号码速率估计方法和装置,属于通信

技术介绍
通信信号调制识别目的就是对截获的一段通信信号,依据较少的先验信息,判断出信号的调制方式。随着通信技术的飞速发展,通信信号的体制和调制样式变得更加复杂多样,信号环境日趋密集,这使得通信信号的识别变得越加困难。通信信号中包含两大类调制信号:数字调制信号和模拟调制信号,利用判决理论时,模拟调制信号和数字调制信号的区分门限往往难以设置。由于模拟调制信号没有码速率之说,其估计结果为任意无规律值,因此利用码速率估计的方法可以有效区分模拟和数字调制,码速率是区分两种调制方式的重要特征。码速率对于数字调制信号特征提取和分析以及识别和解调有着重要的意义,在电子侦察中,数字调制信号通过正交变频和低通滤波后得到系带信号,为了分析与识别信号,需要知道数字调制信号的码速率。因此正确估计数字调制信号的码速率是十分必要的。估计数字调制信号码速率的方法有多种方法,其核心思想是根据调制信号码元之间存在突变,根据突变之间的间隔来估计码速率。对于数字调制基带信号,可以运用延迟相乘的方法来获得信号的码速率。经过延迟相乘之后,输出信号的双边功率谱密度主要包含三项,第一项是直流分量,第二项是码速率及高次谐波,第三项为连续谱。可以根据第二项对调制信号的码速率进行估计。但是时延长度必须大于码元宽度的一半,这在实际的应用中很难保证。论文《无线电信号的调制识别研究》(作者王生兵)公开了利用承载数字信息的幅度调制波形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;5.令Nm=min[lcd(n)],计算序列ni(n)=round(lcd(n)/Nm),得到平均符号间距的估计6.采集信号的码速率估计为:其中fb为输入数据的速率。该方法的不足之处在于码速率估计方法中的门限通过经验取值,缺乏理论依据,门限的取值对不同参数的调制信号不具有适应性;除此之外,因小波变换对瞬变信号有着很好的检测性能,基于小波变换的码速率估计是另外一种思路。文献【基于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)包括:对于连续波调制,已调实信号的表达式为:式中,fc表示载波的频率,A(t)表示已调信号的瞬时幅度;f(t)表示频率;为相位;对信号进行希尔伯特变换,得到信号的正交分量,表达式如下:得到信号s(t)的解析表达式z(t)=s(t)+jv(t)。再进一步地,提取提取瞬时幅度A(l),方法如下:对信号进行采样后,得到离散形式的瞬时幅度A(l),l=1,2,…L。进一步地,提取瞬时频率f(l)的方法:瞬时相位提取,瞬时相位表达式如下:对信号采样后,得到离散形式的瞬时相位序列,l=1,2,…L;去相位折叠,首先计算修正相位序列C(l),l=1,2,…L:则无折叠的相位为:去载频线性相位:采用相性规划法估计出线性相位分量并使误差最小,求出C1和C2两个常数,表示为:其中为非线性相位分量,为无折叠的相位;瞬时频率提取,表达式为其中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本文档来自技高网
...

【技术保护点】
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的比值。...

【技术特征摘要】
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的比值。2.根据权利要求1所述的数字调制信号码速率估计方法,其特征是,提取瞬时幅度A(l)和瞬时频率f(l)包括:对于连续波调制,已调实信号的表达式为:式中,fc表示载波的频率,A(t)表示已调信号的瞬时幅度;f(t)表示频率;为相位;对信号进行希尔伯特变换,得到信号的正交分量,表达式如下:得到信号s(t)的解析表达式z(t)=s(t)+jv(t)。3.根据权利要求2所述的数字调制信号码速率估计方法,其特征是,提取瞬时幅度A(l),方法如下:对信号进行采样后,得到离散形式的瞬时幅度A(l),l=1,2,…L。4.根据权利要求2所述的数字调制信号码速率估计方法,其特征是,提取瞬时频率f(l)的方法:瞬时相位提取,瞬时相位表达式如下:对信号采样后,得到离散形式的瞬时相位序列,l=1,2,…L;去相位折叠,首先计算修正相位序列C(l),l=1,2,…L:则无折叠的相位为:去载频线性相位:采用相性规划法估计出线性相位分量并使误差最小,求出C1和C2两个常数,表示为:其中为非线性相位分量,为无折叠的相位;瞬时频率提取,表达式为其中fs为采样率。5.根据权利要求1...

【专利技术属性】
技术研发人员:王玉桃王胜张昊
申请(专利权)人:南京长峰航天电子科技有限公司
类型:发明
国别省市:江苏,32

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1