CN109450487A - 一种基于能量检测的中频信号幅度估计方法 - Google Patents
一种基于能量检测的中频信号幅度估计方法 Download PDFInfo
- Publication number
- CN109450487A CN109450487A CN201811644642.XA CN201811644642A CN109450487A CN 109450487 A CN109450487 A CN 109450487A CN 201811644642 A CN201811644642 A CN 201811644642A CN 109450487 A CN109450487 A CN 109450487A
- Authority
- CN
- China
- Prior art keywords
- signal
- energy
- value
- freuqncy signal
- range coefficient
- 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
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
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/7073—Synchronisation aspects
- H04B1/7085—Synchronisation aspects using a code tracking loop, e.g. a delay-locked loop
-
- 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
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/7073—Synchronisation aspects
- H04B1/7075—Synchronisation aspects with code phase acquisition
-
- 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
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/7097—Interference-related aspects
- H04B1/7103—Interference-related aspects the interference being multiple access interference
- H04B1/7107—Subtractive interference cancellation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radio Relay Systems (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
本发明公开了基于能量检测的中频信号幅度估计方法。该方法包括在卫星接收端将接收的混合信号与本地产生信号相减,得到差值信号,根据已经获取的载波信息和伪码信息,对差值信号进行解调和解扩运算,得到基带信号,计算基带信号的能量,以最强接收信号的幅度为基准,对幅度系数值进行归一化设置,再以归一化能量的1/Q为步进间隔进行累加变化,每改变一次值计算一次对应的能量值,统计所有幅度系数值对应的能量值,输出能量最小值对应的幅度系数,即为该路中频信号最接近的接收信号幅度的估计值。通过该方法可以减少和消除多址干扰的影响,会使得弱信号支路的码跟踪环路不会出现误差。
Description
技术领域
本发明涉及卫星测控技术领域,尤其涉及一种基于能量检测的中频信号幅度估计方法。
背景技术
为了对卫星实施定轨,测控系统需要采用多站测距体制,多个地面测控站发送的上行测距信号采用码分多址体制。由于多站之间扩频伪码的互相关特性,在卫星接收机收到来自多个地面测控站的多路上行信号的混合信号后,会使得弱信号支路的码跟踪环路受到其它站的多址干扰的影响,从而使码环跟踪出现误差,导致测距值出现偏差。当系统对测距精度要求较高时,多址干扰引起的测距误差不能忽视,需要采取相应的措施。
发明内容
本发明主要解决的技术问题是提供一种基于能量检测的中频信号幅度估计方法,解决卫星通信现有技术中的普勒频移补偿工作体制复杂、补偿精度低、实时性不强等问题。
为解决上述技术问题,本发明采用的一个技术方案是提供一种基于能量检测的中频信号幅度估计方法,包括以下步骤:第一步:在卫星接收端将接收的混合信号与本地产生信号G'M(t)相减,得到差值信号第二步:根据已经获取的载波信息cosM(ωt+θ)和伪码信息PNM(t),对差值信号进行解调和解扩运算,得到解调解扩后的基带信号SM(t)=(AM-αM)D(t),αM为幅度系数;第三步:计算解调后基带信号的能量:
EM=∫SM(t)2dt=∫(AM-αM)2dt;第四步:以最强接收信号的幅度为基准,对幅度系数值αM进行归一化设置,αM以归一化能量的1/Q为步进间隔,以该步进间隔进行累加变化,αM每改变一次值计算一次对应的能量值EM;第五步:重复上述第一步至第四步,直至幅度系数值αM=1为止;第六步:统计所有幅度系数值αM对应的能量值,输出能量最小值EMmin对应的幅度系数αMmin,即为第M路中频信号最接近的接收信号幅度的估计值。
在本发明基于能量检测的中频信号幅度估计方法另一实施例中,所述本地产生信号G'M(t)=αMD(t)PNM(t)cosM(ωt+θ),αM为幅度系数。
在本发明基于能量检测的中频信号幅度估计方法另一实施例中,多次重复所述第一步至第六步,对第M路中频信号进行多次幅度估计,然后再对多次幅度估计的结果进行平均处理。
在本发明基于能量检测的中频信号幅度估计方法另一实施例中,在卫星接收端,在干扰消除之前先利用所述中频信号幅度估计方法对各路中频信号Ri(t)的幅度进行估计,得到对应的幅度系数αimin,然后从混合信号中消除不需要的其他地址的中频信号,而保留所需要的单路中频信号。
本发明的有益效果是:本发明公开了基于能量检测的中频信号幅度估计方法。该方法包括在卫星接收端将接收的混合信号与本地产生信号相减,得到差值信号,根据已经获取的载波信息和伪码信息,对差值信号进行解调和解扩运算,得到基带信号,计算基带信号的能量,以最强接收信号的幅度为基准,对幅度系数值进行归一化设置,再以归一化能量的1/Q为步进间隔进行累加变化,每改变一次值计算一次对应的能量值,统计所有幅度系数值对应的能量值,输出能量最小值对应的幅度系数,即为该路中频信号最接近的接收信号幅度的估计值。通过该方法可以减少和消除多址干扰的影响,会使得弱信号支路的码跟踪环路不会出现误差。
附图说明
图1是根据本发明基于能量检测的中频信号幅度估计方法一实施例的原理图;
图2是根据本发明基于能量检测的中频信号幅度估计方法另一实施例的原理图;
图3是根据本发明基于能量检测的中频信号幅度估计方法另一实施例的流程图;
图4是根据本发明基于能量检测的中频信号幅度估计方法另一实施例的原理图。
具体实施方式
为了便于理解本发明,下面结合附图和具体实施例,对本发明进行更详细的说明。附图中给出了本发明的较佳的实施例。但是,本发明可以以许多不同的形式来实现,并不限于本说明书所描述的实施例。相反地,提供这些实施例的目的是使对本发明的公开内容的理解更加透彻全面。
需要说明的是,除非另有定义,本说明书所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的,不是用于限制本发明。本说明书所使用的术语“和/或”包括一个或多个相关的所列项目的任意的和所有的组合。
在卫星接收机端收到来自地面多个测控站的多路上行信号的混合信号后,对各路信号下变频至中频,然后进行捕获、跟踪、解调,通过载波环和码环的跟踪可以得载波相位和码相位信息,同时解调出符号信息。
设中频信号Ri(t)的幅度为Ai,gi(t)是Ri(t)的幅度归一化信号,则混合信号表示为:
通过各路信号解扩解调获取载波相位、伪码相位和符号信息后,只需要估计出每路信号的幅度,就可以在接收机内恢复出各路中频信号,在此基础上对多址干扰信号进行消除。例如,要想恢复出第i路中频信号Ri(t),在接收端用混合信号减去不含Ri(t)的其他多路信号即可,例如,要获得R1(t),则有:
与此类似,其他路的中频信号恢复如图1所示,其中包括:
优选的,结合图2,对于每一路接收的中频信号进行准确的幅度估计,例如对其中的第M路中频信号进行幅度估计。多路混合信号信号减G'M(t):
其中,RM(t)表示第M路中频信号,利用已经获取的载波信息cosM(ωt+θ)和伪码信息PNM(t)进行解调和解扩,特别是解扩中所采用的扩频码PNM(t)是针对于第M路中频信号RM(t)进行解扩,其他路信号因为解扩效应而被视为噪声滤除,因此解扩之后得到的基带信号只与其中的第M路中频信号RM(t)有关,即:
SM(t)=(AM-αM)D(t)
对基带信号SM(t)进行积分求能量,能量计算可控制在一个数据符号D(t)所持续的区间内,这样可以避开符号跳变的影响,所以数据D(t)不作考虑。根据E=∫SM(t)2dt,可得:
E=∫SM(t)2dt=∫(AM-αM)2dt
优选的,由此进一步根据计算出的能量E来控制幅度系数进行调整。由上式可知:幅度系数αM和接收信号RM(t)的幅度AM越相近,能量E越小。当最小能量Emin≈0时,αM≈AM,说明该路本地信号幅度已经接近该路接收信号的实际幅度值,输出当前的αM值作为该路信号的幅度估计值。
进一步的,基于上述说明,如图3所示,对每一路中频信号幅度估计方法如下:
第一步S101:在卫星接收端将接收的混合信号与本地产生信号G'M(t)相减,得到差值信号
第二步S102:根据已经获取的载波信息cosM(ωt+θ)和伪码信息PNM(t),对差值信号进行解调和解扩运算,得到解调解扩后的基带信号SM(t)=(AM-αM)D(t);
第三步S103:计算解调后基带信号的能量:
EM=∫SM(t)2dt=∫(AM-αM)2dt;
第四步S104:以最强接收信号的幅度为基准,对幅度系数值αM进行归一化设置,αM以归一化能量的1/Q(例如Q=1000)为步进间隔,以该步进间隔进行累加变化,αM每改变一次值计算一次对应的能量值EM;
第五步S105:重复上述第一步至第四步,直至幅度系数值αM=1为止;
第六步S106:统计所有幅度系数值αM对应的能量值,输出能量最小值EMmin对应的幅度系数αMmin,即为第M路中频信号最接近的接收信号幅度的估计值。
优选的,还可以多次重复上述第一步至第六步,对第M路中频信号进行多次幅度估计,然后再对多次幅度估计的结果进行平均处理,依此提高估计的准确度。
优选的,可以进一步根据图3所示实施例,对混合信号中的其他路中频信号的幅度进行估计,得到对应的幅度系数。
优选的,结合图4,在对每一路中频信号的幅度进行估计以后,可以应用到图1所示的卫星接收端,在干扰消除之前先进行中频信号的幅度估计,然后利用幅度估计结果直接从多址信号的混合信号中消除不需要的其他地址的中频信号,而保留所需要的单路中频信号。如图4所示,经过幅度估计后得到各路中频信号的幅度系数α1min、α2min、…α2min,然后应用到每一路中频信号进行多址消除,例如其中: 由此可以从混合信号中增强对每一路信号提取和接收。
由此可见,本发明公开了基于能量检测的中频信号幅度估计方法。该方法包括在卫星接收端将接收的混合信号与本地产生信号相减,得到差值信号,根据已经获取的载波信息和伪码信息,对差值信号进行解调和解扩运算,得到基带信号,计算基带信号的能量,以最强接收信号的幅度为基准,对幅度系数值进行归一化设置,再以归一化能量的1/Q为步进间隔进行累加变化,每改变一次值计算一次对应的能量值,统计所有幅度系数值对应的能量值,输出能量最小值对应的幅度系数,即为该路中频信号最接近的接收信号幅度的估计值。通过该方法可以减少和消除多址干扰的影响,会使得弱信号支路的码跟踪环路不会出现误差。本发明运算量小,降低了对硬件设备的要求。经性能对比测试,使用本方法后能明显改善多址干扰信号引起的码环鉴相偏差,具有重要的实用价值。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (4)
1.一种基于能量检测的中频信号幅度估计方法,其特征在于,包括以下步骤:
第一步:在卫星接收端将接收的混合信号与本地产生信号G'M(t)相减,得到差值信号
第二步:根据已经获取的载波信息cosM(ωt+θ)和伪码信息PNM(t),对差值信号进行解调和解扩运算,得到解调解扩后的基带信号SM(t)=(AM-αM)D(t),αM为幅度系数;
第三步:计算解调后基带信号的能量:
EM=∫SM(t)2dt=∫(AM-αM)2dt;
第四步:以最强接收信号的幅度为基准,对幅度系数值αM进行归一化设置,αM以归一化能量的1/Q为步进间隔,以该步进间隔进行累加变化,αM每改变一次值计算一次对应的能量值EM;
第五步:重复上述第一步至第四步,直至幅度系数值αM=1为止;
第六步:统计所有幅度系数值αM对应的能量值,输出能量最小值EMmin对应的幅度系数αMmin,即为第M路中频信号最接近的接收信号幅度的估计值。
2.根据权利要求1所述的基于能量检测的中频信号幅度估计方法,其特征在于,所述本地产生信号G'M(t)=αMD(t)PNM(t)cosM(ωt+θ),αM为幅度系数。
3.根据权利要求2所述的基于能量检测的中频信号幅度估计方法,其特征在于,多次重复所述第一步至第六步,对第M路中频信号进行多次幅度估计,然后再对多次幅度估计的结果进行平均处理。
4.根据权利要求3所述的基于能量检测的中频信号幅度估计方法,其特征在于,在卫星接收端,在干扰消除之前先利用所述中频信号幅度估计方法对各路中频信号Ri(t)的幅度进行估计,得到对应的幅度系数αimin,然后从混合信号中消除不需要的其他地址的中频信号,而保留所需要的单路中频信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811644642.XA CN109450487B (zh) | 2018-12-29 | 2018-12-29 | 一种基于能量检测的中频信号幅度估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811644642.XA CN109450487B (zh) | 2018-12-29 | 2018-12-29 | 一种基于能量检测的中频信号幅度估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109450487A true CN109450487A (zh) | 2019-03-08 |
CN109450487B CN109450487B (zh) | 2020-02-07 |
Family
ID=65539804
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811644642.XA Active CN109450487B (zh) | 2018-12-29 | 2018-12-29 | 一种基于能量检测的中频信号幅度估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109450487B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5003395A (en) * | 1988-10-26 | 1991-03-26 | U.S. Philips Corporation | Satellite receiver having an independent active energy dispersal signal generator |
US20020051433A1 (en) * | 1999-12-23 | 2002-05-02 | Institut National De La Recherche Scientifique | Interference suppression in CDMA systems |
US20030031234A1 (en) * | 2001-05-17 | 2003-02-13 | Smee John Edward | System and method for received signal prediction in wireless communications systems |
CN103139118A (zh) * | 2013-02-26 | 2013-06-05 | 航天恒星科技有限公司 | 幅度估计方法及针对数据链传输的干扰消除方法和装置 |
CN103916148A (zh) * | 2013-01-05 | 2014-07-09 | 华为技术有限公司 | 一种自适应射频干扰抵消装置、方法以及接收机 |
CN105917622A (zh) * | 2014-01-10 | 2016-08-31 | 迪尔公司 | 用于接收复合信号的方法和接收器 |
CN107682053A (zh) * | 2017-11-08 | 2018-02-09 | 南京天际易达通信技术有限公司 | 一种卫星通信多普勒频移补偿方法、装置及卫星通信系统 |
-
2018
- 2018-12-29 CN CN201811644642.XA patent/CN109450487B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5003395A (en) * | 1988-10-26 | 1991-03-26 | U.S. Philips Corporation | Satellite receiver having an independent active energy dispersal signal generator |
US20020051433A1 (en) * | 1999-12-23 | 2002-05-02 | Institut National De La Recherche Scientifique | Interference suppression in CDMA systems |
US20030031234A1 (en) * | 2001-05-17 | 2003-02-13 | Smee John Edward | System and method for received signal prediction in wireless communications systems |
CN103916148A (zh) * | 2013-01-05 | 2014-07-09 | 华为技术有限公司 | 一种自适应射频干扰抵消装置、方法以及接收机 |
CN103139118A (zh) * | 2013-02-26 | 2013-06-05 | 航天恒星科技有限公司 | 幅度估计方法及针对数据链传输的干扰消除方法和装置 |
CN105917622A (zh) * | 2014-01-10 | 2016-08-31 | 迪尔公司 | 用于接收复合信号的方法和接收器 |
CN107682053A (zh) * | 2017-11-08 | 2018-02-09 | 南京天际易达通信技术有限公司 | 一种卫星通信多普勒频移补偿方法、装置及卫星通信系统 |
Also Published As
Publication number | Publication date |
---|---|
CN109450487B (zh) | 2020-02-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1173778B1 (en) | Signal detector employing correlation analysis of non-uniform and disjoint sample segments | |
US6219391B1 (en) | Wireless communication apparatus and wireless communication method | |
CN101266292B (zh) | 一种gnss反射信号的频域处理装置和方法 | |
CN105607096B (zh) | 一种双星时差频差定位方法和定位装置 | |
US8804808B1 (en) | Dynamic equalization systems and methods for use with a receiver for a multipath channel | |
CN102338878B (zh) | 一种多径抑制boc码跟踪方法以及码跟踪环 | |
CN101132191A (zh) | 一种gnss接收机基带信号处理方法 | |
CN101247138B (zh) | 一种多相关值定宽拟合伪码跟踪鉴相方法 | |
CN102798870A (zh) | 一种卫星导航系统兼容性分析方法 | |
CN105445767B (zh) | 基于平均模糊函数的boc信号参数盲估计方法 | |
CN114578387A (zh) | 多径干扰消除方法及装置、估计方法及装置、接收机 | |
CN103698777A (zh) | 一种基于子空间投影的gnss接收机互相关干扰抑制方法 | |
CN107621643B (zh) | 一种适用于导航信号质量评估的相关域参数精确解算方法 | |
Gentner et al. | OFDM indoor positioning based on TDOAs: Performance analysis and experimental results | |
CN105334521A (zh) | 多用户条件下提高码分测距精度的电路 | |
Wang et al. | Bias and CRB analysis of LoS-based and RSS-based ranging methods | |
CN109450487A (zh) | 一种基于能量检测的中频信号幅度估计方法 | |
CN101154983B (zh) | 一种基于单倍采样的卫星通信抗干扰技术的实现方法 | |
Pan et al. | Analysis of an improved acquisition method for high-dynamic BOC signal | |
WO1999003220A2 (en) | Signal quality measurement | |
CN101945474A (zh) | 无线传感网中基于rssi的定位方法 | |
CN112213751B (zh) | 一种接收的随机跳时脉冲伪卫星信号中脉冲时延估计方法 | |
US11303017B2 (en) | Method and device for estimating the squint of an antenna and method and device for following the pointing of an antenna implementing such a method and device, based on a harmonic analysis | |
CN115079211B (zh) | 一种可用作频率协调依据的卫星导航信号性能评估方法 | |
WO2003036851A1 (en) | Amplitude-offset invariant template detection for pulse position estimation |
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 |