CN106301464B - 一种针对chirp信号的信噪比估计方法 - Google Patents
一种针对chirp信号的信噪比估计方法 Download PDFInfo
- Publication number
- CN106301464B CN106301464B CN201610711860.5A CN201610711860A CN106301464B CN 106301464 B CN106301464 B CN 106301464B CN 201610711860 A CN201610711860 A CN 201610711860A CN 106301464 B CN106301464 B CN 106301464B
- Authority
- CN
- China
- Prior art keywords
- signal
- infe
- sequence
- noise
- energy
- 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
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/69—Spread spectrum techniques
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/69—Spread spectrum techniques
- H04B2001/6912—Spread spectrum techniques using chirp
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种针对chirp信号的信噪比估计方法,包括:获取信噪比估计需要的同步参数估计的中间结果;根据同步参数估计的中间结果计算干扰信号能量;计算信号能量;计算噪声能量;计算chirp信号的信号功率;计算噪声功率;计算信噪比的线性估计值。本发明利用chirp信号通用的同步参数估计过程的中间处理结果,即经过同步参数估计算法处理过的频域中间结果,进一步计算得到信号功率和噪声功率,然后计算信噪比估计值;信噪比估计所需的输入数据和同步参数估计的部分处理过程的输出完全一致,信噪比估计所需的计算及实现复杂度大大降低。
Description
技术领域
本发明涉及信噪比计算技术领域,特别是涉及一种针对chirp信号的信噪比估计方法。
背景技术
在移动通信领域,信噪比因与误码率有直接对应关系,多倍通信测量技术采用作为衡量通信质量的重要指标。在通信过程中,一方面通过信噪比估计可以自适应的选择更有效的解调算法来提高解调性能;另一方面信噪比估计提供了功率控制、调制编码方案调整算法和新到分配算法所需的信道质量信息。目前,信噪比估计是移动通信领域测量通信信道的主要方式,信噪比估计的正确性直接影响着通信系统的性能。Chirp信号(线性调频信号)作为大时宽带宽积信号具有较强的抗干扰性,能同时完成时延和频偏的估计,因而被广泛应用于卫星移动通信系统中。
现有的信噪比估计方案多集中于对常规信号(例如PSK调制信号等)估计方案的描述(例如:徐卓异等在《电子设计工程》上发表的“信噪比估计方法研究”、李辉等在《无线电工程》上发表的“一种QPSK突发信号的信噪比估计方法”、专利“一种信噪比估计方法”(CN103916342A)、专利“一种信噪比忙估计的方法”(CN 103607363A))。而现有关于chirp信号的文献资料多是讨论怎么使用chirp信号进行同步参数估计或者是进行chirp信号本身参数估计的(例如:吴春航等在《电视技术》上发表的“卫星移动通信中利用chirp信号的时频同步算法”、王力男在《无线电通信技术》上发表的“卫星通信系统中chirp信号设计与捕获”、专利“基于chirp信号进行同步的方法、接收端装置及通信系统”(CN 105162571A)、专利“一种基于功率谱的线性调频信号参数估计方法”(CN 103063909B)、专利“Synchronization In Mobile Satellite Systems Using Dual-Chirp Waveform”(US6418158B1))。这两个领域的资料中都没有发现针对chirp信号的信噪比估计方法的相关描述。有些已有方案(例如“李辉等在《无线电工程》上发表的“一种QPSK突发信号的信噪比估计方法””)中的最小均方误差估计法经过修改可以用于chirp信号的信噪比估计。但是和上述文献中所描述的chirp信号同步参数估计方案对比来看,两个方法之间基本没有交集。也就是说在使用chirp信号进行同步参数估计的同时,还要有另一套独立的过程对chirp信号进行处理并进行信噪比估计。由此可见,这种方案实现复杂度较高,特别是在功率和尺寸都受限的手持移动终端中。
从chirp信号同步相关的文献资料中可以看到:由于chirp信号的特点,其同步参数估计方案大多都有如图1所示的处理结构。接收chirp信号与本地存储的上下扫频信号分别相乘,再把相乘的结果做FFT运算(快速傅里叶计算);然后再两个FFT序列结果中各自寻找能量峰值得到能量峰值所在的位置,再根据峰值位置与接收信号采样频率的关系求得其所对应的频率值。最后根据一定的策略判断是否满足同步标准,然后根据两个频率值计算当前的频偏和时偏值。也可以不进行同步状态判断直接进行频偏和时偏的计算,视具体的工作场景而定。
从上述过程可以知道经过FFT之后的频域信号中同时含有信号功率和噪声功率,理论上如果能有合适的方式进行处理应该可以得到信噪比的值。这种普遍的chirp信号同步参数估计处理器结构以及其中间结果特点都为融合信噪比估计方案提供了良好的基础。在此基础上如果能把chirp信号的同步参数估计和信噪比估计两个过程进行适当融合,这样的方案便可以降低计算、实现复杂度以及功耗但是目前还没有看到相关技术方案的提出。
发明内容
本发明的目的在于克服现有技术的不足,提供一种针对chirp信号的信噪比估计方法,信噪比估计所需的输入数据和同步参数估计的部分处理过程的输出完全一致,信噪比估计所需的计算及实现复杂度大大降低。
本发明的目的是通过以下技术方案来实现的:一种针对chirp信号的信噪比估计方法,包括:
S1.获取chirp信号同步过程中的接收信号、本地上扫频信号和本地下扫频信号;将所述接收信号与所述本地上扫频信号相乘,并对相乘结果进行快速傅里叶变换得到序列为mu;将所述接收信号与所述本地下扫频信号相乘,并对相乘结果进行快速傅里叶变换得到序列为md;
chirp信号同步过程中,计算序列mu的能量峰值Eu,max和能量峰值Eu,max在序列mu中的位置iu,max,其中1≤iu,max≤N,以及计算序列md的能量峰值Ed,max和能量峰值id,max在序列md中的位置iu,max,其中1≤id,max≤N;
S2.计算序列mu的干扰信号能量Eu,infe和序列md的干扰信号能量Ed,infe;
S3.计算序列mu的信号能量Eu,sig和序列md的信号能量Ed,sig;
S4.计算序列mu的噪声能量Eu,noise和序列md的噪声能量Ed,noise;
S5.计算chirp信号的信号功率Ps;
S6.计算噪声功率Pn;
S7.根据计算信号功率Ps和噪声功率Pn得到信噪比的线性估计值。
计算信号功率Ps除以噪声功率Pn得到信噪比的线性估计值。
所述步骤S1中,mu=[mu,1,...,mux,,...m,Nu,],其中N为快速傅里叶变换的长度;md=[md,1,...,md,x,...,md,N],其中N为快速傅里叶变换的长度。
所述步骤S2包括:
S21.定义chirp信号的采样频率为fsmp,chirp信号的扫频范围为[-fscan,fscan],计算得到其中表示向下取整,则序列md中干扰信号的频谱范围为[1,Ninfe]和[N-Ninfe+1,N];序列mu中干扰信号的频谱范围为[1,Ninfe]和[N-Ninfe+1,N];
S22.定义iu是区间[iu,max-16,iu,max+15]内的值,且1≤iu≤N,将序列mu对应位置的都变为0;
定义id是区间[id,max-16,id,max+15]内的值,且1≤id≤N,将序列md对应位置的都变为0;
S23.将序列mu中[1,Ninfe]和[N-Ninfe+1,N]两段中的所有值求模平方,再求所有模平方的和,最后将所有模平方的和除以2×Ninfe-32,得到干扰信号能量Eu,infe;
将序列md中[1,Ninfe]和[N-Ninfe+1,N]两段中的所有值求模平方,再求所有模平方的和,最后将所有模平方的和除以2×Ninfe-32,得到干扰信号能量Ed,infe。
所述信号能量Eu,sig的计算公式为:Eu,sig=Eu,max-Eu,infe。
所述信号能量Ed,sig的计算公式为:Ed,sig=Ed,max-Ed,infe。
所述噪声能量Eu,noise的计算方法为:将序列mu中[Ninfe+2,N-Ninfe-1]中的所有值求模平方,再求所有模平方的和,最后将所有模平方的和除以N-2×Ninfe-2,得到噪声能量Eu,noise。
所述噪声能量Ed,noise的计算方法为:将序列md中[Ninfe+2,N-Ninfe-1]中的所有值求模平方,再求所有模平方的和,最后将所有模平方的和除以N-2×Ninfe-2,得到噪声能量Ed,noise。
所述信号功率Ps的计算公式为:
式中,M为chirp信号的长度。
所述噪声功率Pn的计算公式为:
式中,M为chirp信号的长度。
所述信噪比的线性估计值为
本发明的有益效果是:本发明利用chirp信号通用的同步参数估计过程的中间处理结果,即经过同步参数估计算法处理过的频域中间结果,进一步计算得到信号功率和噪声功率,然后计算信噪比估计值;信噪比估计所需的输入数据和同步参数估计的部分处理过程的输出完全一致,信噪比估计所需的计算及实现复杂度大大降低。
附图说明
图1为现有chirp信号的同步参数估计方案的处理结构示意图;
图2为本发明的一个实施例的流程示意图;
图3为加入信噪比估计后chirp信号的同步参数估计方案的处理结构示意图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
如图2所示,一种针对chirp信号的信噪比估计方法,包括:
S1.获取信噪比估计需要的同步参数估计的中间结果:获取chirp信号同步过程中的接收信号、本地上扫频信号和本地下扫频信号;将所述接收信号与所述本地上扫频信号相乘,并对相乘结果进行快速傅里叶变换(FFT)得到序列为mu;将所述接收信号与所述本地下扫频信号相乘,并对相乘结果进行快速傅里叶变换得到序列为md。
mu=[mu,1,...,mu,x,...,mu,N],其中N为快速傅里叶变换的长度;md=[md,1,...,md,x,...,md,N],其中N为快速傅里叶变换的长度。
chirp信号同步过程中,计算序列mu的能量峰值Eu,max和能量峰值Eu,max在序列mu中的位置iu,max,其中1≤iu,max≤N,以及计算序列md的能量峰值Ed,max和能量峰值id,max在序列md中的位置iu,max,其中1≤id,max≤N。
S2.根据同步参数估计的中间结果计算干扰信号能量:计算序列mu的干扰信号能量Eu,infe和序列md的干扰信号能量Ed,infe。
所述步骤S2包括:
S21.定义chirp信号的采样频率为fsmp,chirp信号的扫频范围为[-fscan,fscan],计算得到其中表示向下取整,则序列md中干扰信号的频谱范围为[1,Ninfe]和[N-Ninfe+1,N];序列mu中干扰信号的频谱范围为[1,Ninfe]和[N-Ninfe+1,N]。
S22.定义iu是区间[iu,max-16,iu,max+15]内的值,且1≤iu≤N,将序列mu对应位置的都变为0;
定义id是区间[id,max-16,id,max+15]内的值,且1≤id≤N,将序列md对应位置的都变为0。
S23.将序列mu中[1,Ninfe]和[N-Ninfe+1,N]两段中的所有值求模平方,再求所有模平方的和,最后将所有模平方的和除以2×Ninfe-32,得到干扰信号能量Eu,infe;
将序列md中[1,Ninfe]和[N-Ninfe+1,N]两段中的所有值求模平方,再求所有模平方的和,最后将所有模平方的和除以2×Ninfe-32,得到干扰信号能量Ed,infe。
S3.计算信号能量:计算序列mu的信号能量Eu,sig和序列md的信号能量Ed,sig。
所述信号能量Eu,sig的计算公式为:Eu,sig=Eu,max-Eu,infe,其中,Eu,max为序列mu的能量峰值,Eu,infe为干扰信号能量。
所述信号能量Ed,sig的计算公式为:Ed,sig=Ed,max-Ed,infe,其中,Ed,max为序列md的能量峰值,Ed,infe为干扰信号能量。
S4.计算噪声能量:计算序列mu的噪声能量Eu,noise和序列md的噪声能量Ed,noise。
所述噪声能量Eu,noise的计算方法为:将序列mu中[Ninfe+2,N-Ninfe-1]中的所有值求模平方,再求所有模平方的和,最后将所有模平方的和除以N-2×Ninfe-2,得到噪声能量Eu,noise。
所述噪声能量Ed,noise的计算方法为:将序列md中[Ninfe+2,N-Ninfe-1]中的所有值求模平方,再求所有模平方的和,最后将所有模平方的和除以N-2×Ninfe-2,得到噪声能量Ed,noise。
S5.计算chirp信号的信号功率Ps。
所述信号功率Ps的计算公式为:
式中,M为chirp信号的长度S7.根据计算信号功率Ps和噪声功率Pn得到信噪比的线性估计值。
S6.计算噪声功率Pn。
所述噪声功率Pn的计算公式为:
式中,M为chirp信号的长度。
S7.计算信号功率Ps除以噪声功率Pn得到信噪比的线性估计值。
所述信噪比的线性估计值为
在一个实施例中:M=480,N=512,fsmp=64kHz,fscan=7.68kHz,则Ninfe=122,从而计算得到序列mu和序列md中干扰信号的频谱范围相同,均为[1,122]和[391,512],2×Ninfe-32=212,[Ninfe+2,N-Ninfe-1]为[124,389],N-2×Ninfe-2=266。
从图3可以看到信噪比估计过程与同步参数估计过程高度融合,没有任何冲突。只是需要添加若干针对信噪比计算得模块,但这些模块也基本只包含简单的加乘运算。
以上所述仅是本发明的优选实施方式,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
Claims (1)
1.一种针对线性调频chirp信号的信噪比估计方法,其特征在于,包括:
S1.获取chirp信号同步过程中的接收信号、本地上扫频信号和本地下扫频信号;将所述接收信号与所述本地上扫频信号相乘,并对相乘结果进行快速傅里叶变换得到序列为mu;
将所述接收信号与所述本地下扫频信号相乘,并对相乘结果进行快速傅里叶变换得到序列为md;
chirp信号同步过程中,计算序列mu的能量峰值Eu,max和能量峰值Eu,max在序列mu中的位置iu,max,其中1≤iu,max≤N,以及计算序列md的能量峰值Ed,max和能量峰值Ed,max在序列md中的位置id,max,其中1≤id,max≤N;
S2.计算序列mu的干扰信号能量Eu,infe和序列md的干扰信号能量Ed,infe;
S3.计算序列mu的信号能量Eu,sig和序列md的信号能量Ed,sig;
S4.计算序列mu的噪声能量Eu,noise和序列md的噪声能量Ed,noise;
S5.计算chirp信号的信号功率Ps;
S6.计算噪声功率Pn;
S7.根据计算信号功率Ps和噪声功率Pn得到信噪比的线性估计值,即
计算信号功率Ps除以噪声功率Pn得到信噪比的线性估计值;
其中,mu=[mu,1,...,mu,x,...,mu,N],md=[md,1,...,md,x,...,md,N],其中N为快速傅里叶变换的长度;
所述步骤S2还包括:
S21.定义chirp信号的采样频率为fsmp,chirp信号的扫频范围为[-fscan,fscan],计算得到其中表示向下取整,则序列md中干扰信号的频谱范围为[1,Ninfe]和[N-Ninfe+1,N];序列mu中干扰信号的频谱范围为[1,Ninfe]和[N-Ninfe+1,N];
S22.定义iu是区间[iu,max-16,iu,max+15]内的值,且1≤iu≤N,将序列mu对应位置的都变为0;
定义id是区间[id,max-16,id,max+15]内的值,且1≤id≤N,将序列md对应位置的都变为0;
S23.将序列mu中[1,Ninfe]和[N-Ninfe+1,N]两段中的所有值求模平方,再求所有模平方的和,最后将所有模平方的和除以(2×Ninfe-32),得到干扰信号能量Eu,infe;
将序列md中[1,Ninfe]和[N-Ninfe+1,N]两段中的所有值求模平方,再求所有模平方的和,最后将所有模平方的和除以(2×Ninfe-32),得到干扰信号能量Ed,infe;
所述信号能量Eu,sig的计算公式为:Eu,sig=Eu,max-Eu,infe;
所述信号能量Ed,sig的计算公式为:Ed,sig=Ed,max-Ed,infe;
所述噪声能量Eu,noise的计算方法为:将序列mu中[Ninfe+2,N-Ninfe-1]中的所有值求模平方,再求所有模平方的和,最后将所有模平方的和除以(N-2×Ninfe-2),得到噪声能量Eu,noise;
所述噪声能量Ed,noise的计算方法为:将序列md中[Ninfe+2,N-Ninfe-1]中的所有值求模平方,再求所有模平方的和,最后将所有模平方的和除以(N-2×Ninfe-2),得到噪声能量Ed,noise;
所述信号功率Ps的计算公式为:
所述噪声功率Pn的计算公式为:
式中,M为chirp信号的长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610711860.5A CN106301464B (zh) | 2016-08-23 | 2016-08-23 | 一种针对chirp信号的信噪比估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610711860.5A CN106301464B (zh) | 2016-08-23 | 2016-08-23 | 一种针对chirp信号的信噪比估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106301464A CN106301464A (zh) | 2017-01-04 |
CN106301464B true CN106301464B (zh) | 2018-10-12 |
Family
ID=57615020
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610711860.5A Active CN106301464B (zh) | 2016-08-23 | 2016-08-23 | 一种针对chirp信号的信噪比估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106301464B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111030762B (zh) * | 2019-12-23 | 2021-10-29 | 北京华力创通科技股份有限公司 | chirp信号功率测量方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6570394B1 (en) * | 1999-01-22 | 2003-05-27 | Thomas H. Williams | Tests for non-linear distortion using digital signal processing |
CN103116155A (zh) * | 2012-09-28 | 2013-05-22 | 北京理工大学 | 一种用于舰艇编队情况下的同型雷达同频干扰抑制方法 |
CN103888405A (zh) * | 2014-04-10 | 2014-06-25 | 厦门大学 | 一种全频谱载波调制方法 |
-
2016
- 2016-08-23 CN CN201610711860.5A patent/CN106301464B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106301464A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109510787B (zh) | 水声信道下线性调频信号参数估计方法及装置 | |
CN104601259B (zh) | 具有同相正交信号失衡估测及校正技术的无线通信接收器 | |
CN103837740A (zh) | 一种高精度数字瞬时测频方法及装置 | |
CN103428153B (zh) | 卫星移动通信中gmsk信号接收方法 | |
CN107404450A (zh) | 解调信号的方法及装置 | |
Harris et al. | The DST group ionospheric sounder replacement for JORN | |
CN108123723B (zh) | 一种基于无线电技术的定位系统及方法 | |
CN102645589A (zh) | 脉冲检测方法及系统 | |
CN104253659B (zh) | 一种频谱检测方法及其装置 | |
CN1256811A (zh) | 补偿调频接收机中喀音噪声的方法和设备 | |
CN109067684A (zh) | 一种低频2fsk通信解调方法、装置及计算机设备 | |
CN106301464B (zh) | 一种针对chirp信号的信噪比估计方法 | |
CN103916199B (zh) | 一种天线信号的时延和相位调整装置与方法 | |
CN106330342A (zh) | 一种低计算复杂度的水声通信多普勒因子估计方法 | |
CN108650203B (zh) | 一种基于侦察接收机的调制模式识别方法 | |
CN102540153A (zh) | 一种基于站间直达波干扰的阵列幅相误差校准方法 | |
CN108900445B (zh) | 一种信号符号率估计的方法及装置 | |
CN103529429A (zh) | 用于线性调频信号的频域中的脉冲压缩方法 | |
CN105933265B (zh) | 一种对qam信号进行非数据辅助的相位噪声盲估计方法 | |
US20190028322A1 (en) | Joint noncoherent demodulation and carrier frequency offset correction based on non-linear filtering | |
CN110113273B (zh) | 一种基于频谱对称的时频重叠信号盲分离方法 | |
CN110417696A (zh) | 一种基于随机共振的2ask信号解调方法 | |
CN102508060B (zh) | 一种抗外界脉冲干扰的阵列校准预处理方法 | |
CN102307164B (zh) | 数字频率估计方法及系统 | |
CN104270328A (zh) | 一种信噪比实时估计方法 |
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 | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: 610000 Building 1, No. 11, Gaopeng Avenue, Chengdu hi tech Zone, Chengdu, Sichuan Patentee after: Chengdu Guoheng Space Technology Engineering Co.,Ltd. Address before: No. 1 high tech Zone Gaopeng road in Chengdu city of Sichuan Province in 610041 Patentee before: CHENGDU GUOHENG SPACE TECHNOLOGY ENGINEERING Co.,Ltd. |