【技术实现步骤摘要】
数字调制信号码速率估计方法和装置
本专利技术涉及一种数字调制信号码速率估计方法和装置,属于通信
技术介绍
通信信号调制识别目的就是对截获的一段通信信号,依据较少的先验信息,判断出信号的调制方式。随着通信技术的飞速发展,通信信号的体制和调制样式变得更加复杂多样,信号环境日趋密集,这使得通信信号的识别变得越加困难。通信信号中包含两大类调制信号:数字调制信号和模拟调制信号,利用判决理论时,模拟调制信号和数字调制信号的区分门限往往难以设置。由于模拟调制信号没有码速率之说,其估计结果为任意无规律值,因此利用码速率估计的方法可以有效区分模拟和数字调制,码速率是区分两种调制方式的重要特征。码速率对于数字调制信号特征提取和分析以及识别和解调有着重要的意义,在电子侦察中,数字调制信号通过正交变频和低通滤波后得到系带信号,为了分析与识别信号,需要知道数字调制信号的码速率。因此正确估计数字调制信号的码速率是十分必要的。估计数字调制信号码速率的方法有多种方法,其核心思想是根据调制信号码元之间存在突变,根据突变之间的间隔来估计码速率。对于数字调制基带信号,可以运用延迟相乘的方法来获得信号的码速率。经过延迟相乘之后,输出信号的双边功率谱密度主要包含三项,第一项是直流分量,第二项是码速率及高次谐波,第三项为连续谱。可以根据第二项对调制信号的码速率进行估计。但是时延长度必须大于码元宽度的一半,这在实际的应用中很难保证。论文《无线电信号的调制识别研究》(作者王生兵)公开了利用承载数字信息的幅度调制波形p(n)进行码速率估计。对于ASK信号,p(n)就是信号的瞬时幅度;对于PSK信 ...
【技术保护点】
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, ...
【技术特征摘要】
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
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。