CN102608583A - 时频域混叠信号的时差和频差估计方法 - Google Patents
时频域混叠信号的时差和频差估计方法 Download PDFInfo
- Publication number
- CN102608583A CN102608583A CN2012100488923A CN201210048892A CN102608583A CN 102608583 A CN102608583 A CN 102608583A CN 2012100488923 A CN2012100488923 A CN 2012100488923A CN 201210048892 A CN201210048892 A CN 201210048892A CN 102608583 A CN102608583 A CN 102608583A
- Authority
- CN
- China
- Prior art keywords
- signal
- time
- frequency
- difference
- value
- 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
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
本发明公开了一种时频域混叠信号的时差和频差估计方法,信号检测与参数估计,包括:S1:计算主观测点接收的信号xi(t)和辅观测点接收的信号y(t)的互模糊函数值,并取模,所述xi(t)和y(t)包括若干时频域混叠的信号;S2:搜索互模糊函数最大模值,提取最大模值对应信号的时差和频差估计值S3:根据所述估计值抵消掉主观测点中所述最大模值对应的信号,得到抵消后的信号xi+1(t);S4:计算xi+1(t)与y(t)的互模糊函数值,并取模,转到步骤S2,直到混叠信号中最后一个信号的时差和频差估计值被估计出来。本发明减小了强信号的互模糊函数模值,从而解决了复杂电磁环境下,时频域混叠信号中弱信号时差和频差无法检测或估计误差大的问题,提高了弱信号的检测及时差和频差的估计性能。
Description
技术领域
本发明涉及信号检测与参数估计技术领域,特别涉及一种时频域混叠信号的时差和频差估计方法。
背景技术
通过接收辐射源目标发出的无线电信号测量目标的位置在航空、声纳、电子侦察及测控等诸多领域具有广泛的应用。其中通过估计信号到达两个观测点的时间差(时差)和频率差(频差)进行定位的方法,具有观测点资源利用率高和定位性能优越的特点,已成为目前研究的热点。
对时差和频差的参数估计,通常采用经典的互模糊函数算法,其原理如下:
定义互模糊函数:
表达式中,τ和f是时差和频差变量,T为积分时间,“*”表示共轭运算,u1(t)和u2(t)是两个具有相同参量的复包络信号,二者存在如下关系:
u2(t)=u1(t-D)e-j2πΔf(t-D) (2)
其中D,Δf分别为信号u1(t)和u2(t)之间的时差和频差。将(2)式代入(1)式并取模,可得:
注意到|ej2πΔf·(τ-D)|=1,并根据柯西-施瓦茨不等式可得:
|A(τ,f)|在τ=D,f=Δf处取得最大值,即在相应互模糊函数图中表现为一个显著的峰。
在实际情况中,随着无线电技术的快速发展,空间电磁环境日益复杂。因此,同时进入观测点接收机的通常是多个辐射源信号,甚至这些辐射源信号的频点相同或相近,导致面临时频域混叠信号处理及参数估计上的难题。采用经典的互模糊函数算法估计时差和频差,对于空间位置上分离的辐射源,可以充分利用辐射源的空间分离特性,形成多个相关峰,从而检测每个相关峰对应的时差和频差,这是利用互模糊函数进行时频域混叠信号的时差和频差联合估计的优势。但是时频域混叠信号中强信号会对弱信号形成严重干扰,即弱信号互模糊函数相关峰会被强信号相关峰旁瓣淹没,这将导致无法检测弱信号,或造成弱信号参数估计的误差很大。这一问题严重制约了互模糊函数方法在面临复杂电磁环境下的实际应用性能。
现有的采用互模糊函数算法估计时频域混叠信号的时差和频差参数的方法,只在多个信号的功率相当,或者对于弱信号的主瓣大于强信号最大副瓣时,可以有效对多个信号的时差和频差进行估计。而当强信号的副瓣强于弱信号的主瓣,即弱信号互模糊函数峰主瓣被淹没的情况,现有方法无法完成对弱信号时差和频差的估计,或估计结果误差大而无法使用。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是:在复杂电磁环境下,解决弱信号无法检测与参数估计误差大的问题。
(二)技术方案
为解决上述技术问题,本发明提供了一种时频域混叠信号的时差和频差估计方法,包括以下步骤:
S1:计算主观测点接收的信号xi(t)和辅观测点接收的信号y(t)的互模糊函数值,并取模,所述xi(t)和y(t)包括若干时频域混叠的信号;
S2:搜索互模糊函数最大模值,提取最大模值对应信号的时差和频差估计值
S4:计算xi+1(t)与y(t)的互模糊函数值,并取模,转到步骤S2,直到混叠信号中最后一个信号的时差和频差估计值被估计出来,其中i表示循环的次数,i=1……N,N≥1。
其中,步骤S3具体包括:
xi(t)减去y″(t),得到抵消后的信号xi+1(t)。
其中,所述步骤S2中还包括步骤,按以下公式计算全局门限值Th,
Th=λ·m+σ (5)
其中,m为当前所有模值的均值,σ为当前所有模值的方差,λ为计算全局门限的系数,若所述最大模值大于所述全局门限值则继续执行步骤S3和S4,否则结束。
其中,所述λ取值为3~5。
(三)有益效果
本发明利用辐射源信号的空间可分离特性,计算时频域混叠信号的互模糊函数并取模,通过最大模值搜索,获得强信号的时差和频差估计值,并将估计值通过参数反馈方法,用于强信号反馈抵消,以减小强信号的互模糊函数模值,从而解决了复杂电磁环境下,时频域混叠信号中弱信号时差和频差无法检测或估计误差大的问题,提高了弱信号的检测及时差和频差的估计性能。
附图说明
图1是本发明实施例的一种时频域混叠信号的时差和频差估计方法流程图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
本发明首先通过互模糊函数算法估计强信号的时差和频差,通过反馈强信号的时差和频差,抵消强信号对弱信号的干扰,完成对弱信号时差和频差估计。本实施例中,考虑存在两个时频混叠信号(S1(t)和S2(t))的典型情况,主观测点接收的信号x1(t)和辅观测点接收的信号y(t)的信号形式为:
x1(t)=S1(t)+S2(t)+w1(t) (6)
(6)式和(7)式中,τ1和τ2为两个信号分别到达两观测点的真实时差,f1和f2为两个信号分别到达两观测点的真实频差,本实施例中,设τ1=123us,τ2=-53us,f1=2104Hz,f2=5123Hz。S1(t)和S2(t)为主观测点接收到的两个辐射源的信号,设S1(t)为强信号,其功率值为0dbm,S2(t)为弱信号,功率值为-10dbm,S1(t)的互模糊函数相关峰的最大副瓣淹没S2(t)信号互模糊函数相关峰的主瓣。w1(t),w2(t)分别为独立于信号的高斯噪声。
具体流程如图1所示,包括:
步骤S101,计算主观测点接收的信号x1(t)和辅观测点接收的信号y(t)的互模糊函数值,并取模。
本步骤中,根据(1)式定义得到互模糊函数为:
根据(6)式和(7)式的信号定义形式,将x1(t)和y(t)代入(8)式可得:
对(9)式取模,得到互模糊函数模值。
步骤S102,搜索互模糊函数最大模值,提取其对应的时差和频差,即强信号S1(t)的时差和频差估计值本步骤中,搜索到互模糊函数最大模值为5×105,提取其对应的时差和频差为( ),可见时差估计值误差频差估计值误差
(11)
将x1(t)减去y″(t),得到x2(t)为:
(12)
通过反馈后将强信号抵消,因为δτ,δf很小,可以忽略不计,即强信号得到很好的抑制。强信号干扰抵消还可以通过多天线阵列波束合成方法来实现,但需要多副天线,而且系统复杂,计算量大,此处不再赘述。
步骤S104,计算x2(t)与y(t)的互模糊函数值,并取模。本步骤中,计算互模糊函数值并取模的方法同步骤S101。
转向步骤S102,搜索互模糊函数最大模值,提取其对应的时差和频差,即弱信号S2(t)的时差和频差估计值本步骤中,存在两个互模糊函数相关峰,其中之一为步骤S102中已得到的估计值搜索另一个相关峰的最大模值为9×104,提取最大模值对应的时差和频差为( ),即得到了弱信号S2(t)的时差和频差估计值。由于本实施是两个信号混叠的情况,在执行步骤S102后便结束。
上述步骤S102中,可根据(5)式计算门限值Th=1.16×105,最大模值5×105超过门限值,说明有信号存在,需要执行步骤S103抵消掉强信号。抵消后根据(5)式重新计算全局门限值为Th=9.8×103,由于本实施是两个信号混叠的情况,可以直接提取该弱信号的时差和频差估计值。
以上是两个信号叠加的情况,对于多个信号的叠加的情况,即:
x1(t)=S1(t)+…+Sk(t)+…Sn(t)+w1(t) (13)
(14)
循环地执行步骤S102~S104,即每次抵消掉最大模值对应的信号得到抵消后的信号xi(t)之后。继续计算xi(t)与y(t)的互模糊函数值,并取模,提取最大模值对应信号的时差和频差估计值按下面过程抵消掉最大模值对应的信号得到xi+1(t):
xi(t)减去y″(t),得到抵消后的信号xi+1(t)。
直到最后一个信号(最弱的信号)的时差和频差估计值被估计出来。其中i表示循环的次数,xi+1(t)表示第i次抵消后的信号,每一次执行步骤S102搜索互模糊函数最大模值,提取最大模值对应的信号也可以记为是第i个信号,提取的最大模值对应的信号的时差和频差估计值记为
优选地,按(5)式设定一个全局门限值Th=λ·m+σ,其中,m为进行互模糊函数计算后所有模值的均值,σ为进行互模糊函数计算后所有模值的方差(每一次抵消掉最大模值的信号后,总的模值数减少),λ为计算全局门限的系数,通常取3~5,如取5。
在一定的采样时间内,主观测点、辅观测点采样数据的互模糊函数模值是一个二维数组,假设该二维数组大小为M×N,则对这M×N个数据计算均值和方差,在步骤S102中利用(5)式Th=λ·m+σ计算得到全局门限值。如果检测有大于全局门限值的峰值,则说明还有信号存在,继续执行步骤S102~S104,否则结束。本发明通过检测大于全局门限的最大模值,得到强信号时频差,并对强信号进行反馈抵消,可以提高弱信号的检测能力,同时当检测到的模值不大于全局门限的最大模值时,不再执行,提高了方法执行效率。
以上实施方式仅用于说明本发明,而并非对本发明的限制,有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变化和变型,因此所有等同的技术方案也属于本发明的范畴,本发明的专利保护范围应由权利要求限定。
Claims (4)
3.如权利要求1或2所述的时频域混叠信号的时差和频差估计方法,其特征在于,所述步骤S2中还包括步骤,按以下公式计算全局门限值Th:
Th=λ·m+σ
其中,m为当前所有模值的均值,σ为当前所有模值的方差,λ为计算全局门限的系数,若所述最大模值大于所述全局门限值则继续执行步骤S3和S4,否则结束。
4.如权利要求3所述的时频域混叠信号的时差和频差估计方法,其特征在于,所述λ取值为3~5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210048892.3A CN102608583B (zh) | 2012-02-28 | 2012-02-28 | 时频域混叠信号的时差和频差估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210048892.3A CN102608583B (zh) | 2012-02-28 | 2012-02-28 | 时频域混叠信号的时差和频差估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102608583A true CN102608583A (zh) | 2012-07-25 |
CN102608583B CN102608583B (zh) | 2014-03-12 |
Family
ID=46526093
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210048892.3A Active CN102608583B (zh) | 2012-02-28 | 2012-02-28 | 时频域混叠信号的时差和频差估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102608583B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763158A (zh) * | 2018-06-01 | 2018-11-06 | 中国人民解放军战略支援部队信息工程大学 | 一种时频差联合计算方法和系统 |
CN112688716A (zh) * | 2020-11-06 | 2021-04-20 | 西安电子科技大学 | 一种时频混叠信号时频差估计方法 |
CN113567919A (zh) * | 2021-05-26 | 2021-10-29 | 北京科电航宇空间技术有限公司 | 一种通信信号的时频差参数联合估计gpu实现方法 |
CN116203501A (zh) * | 2023-04-27 | 2023-06-02 | 中国人民解放军战略支援部队航天工程大学 | 基于频域互模糊函数插值映射辐射源无源定位方法和设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1238086A (zh) * | 1996-11-27 | 1999-12-08 | 诺基亚电信公司 | 测量频差的方法和接收机 |
US20050032513A1 (en) * | 1999-12-14 | 2005-02-10 | Norman Charles P. | Strong signal cancellation to enhance processing of weak spread spectrum signal |
CN102238125A (zh) * | 2011-06-21 | 2011-11-09 | 西安电子科技大学 | 存在剩余时偏的ofdm系统整数倍频偏估计方法 |
-
2012
- 2012-02-28 CN CN201210048892.3A patent/CN102608583B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1238086A (zh) * | 1996-11-27 | 1999-12-08 | 诺基亚电信公司 | 测量频差的方法和接收机 |
US20050032513A1 (en) * | 1999-12-14 | 2005-02-10 | Norman Charles P. | Strong signal cancellation to enhance processing of weak spread spectrum signal |
CN102238125A (zh) * | 2011-06-21 | 2011-11-09 | 西安电子科技大学 | 存在剩余时偏的ofdm系统整数倍频偏估计方法 |
Non-Patent Citations (8)
Title |
---|
《信号处理》 20100930 魏石峰等 一种基于CAF的频谱混叠多信号检测算法 第1435页左栏第1行到第1436页右栏第4行 1-4 第26卷, 第9期 * |
《信息工程大学学报》 20070930 侯成郭 互模糊函数(CAF)在分布式卫星信号接收中的应用 全文 1-4 第8卷, 第03期 * |
《南京理工大学学报(自然科学版)》 20090831 朱伟强等 基于互模糊函数的时差频差估计算法及实现途径 全文 1-4 第33卷, 第04期 * |
《电波科学学报》 20061031 孙正波等 一种时差/频率差快速联合估计方法 全文 1-4 第21卷, 第05期 * |
侯成郭: "互模糊函数(CAF)在分布式卫星信号接收中的应用", 《信息工程大学学报》, vol. 8, no. 03, 30 September 2007 (2007-09-30) * |
孙正波等: "一种时差/频率差快速联合估计方法", 《电波科学学报》, vol. 21, no. 05, 31 October 2006 (2006-10-31) * |
朱伟强等: "基于互模糊函数的时差频差估计算法及实现途径", 《南京理工大学学报(自然科学版)》, vol. 33, no. 04, 31 August 2009 (2009-08-31) * |
魏石峰等: "一种基于CAF的频谱混叠多信号检测算法", 《信号处理》, vol. 26, no. 9, 30 September 2010 (2010-09-30) * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763158A (zh) * | 2018-06-01 | 2018-11-06 | 中国人民解放军战略支援部队信息工程大学 | 一种时频差联合计算方法和系统 |
CN108763158B (zh) * | 2018-06-01 | 2021-11-09 | 中国人民解放军战略支援部队信息工程大学 | 一种时频差联合计算方法和系统 |
CN112688716A (zh) * | 2020-11-06 | 2021-04-20 | 西安电子科技大学 | 一种时频混叠信号时频差估计方法 |
CN112688716B (zh) * | 2020-11-06 | 2021-12-24 | 西安电子科技大学 | 一种时频混叠信号时频差估计方法 |
CN113567919A (zh) * | 2021-05-26 | 2021-10-29 | 北京科电航宇空间技术有限公司 | 一种通信信号的时频差参数联合估计gpu实现方法 |
CN116203501A (zh) * | 2023-04-27 | 2023-06-02 | 中国人民解放军战略支援部队航天工程大学 | 基于频域互模糊函数插值映射辐射源无源定位方法和设备 |
CN116203501B (zh) * | 2023-04-27 | 2024-01-05 | 中国人民解放军战略支援部队航天工程大学 | 基于频域互模糊函数插值映射辐射源无源定位方法和设备 |
Also Published As
Publication number | Publication date |
---|---|
CN102608583B (zh) | 2014-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7180046B2 (ja) | 速度の曖昧さの解消を含むmimoレーダにおける速度検出のための方法及び装置 | |
CN101561490B (zh) | 基于码字结构的盲自适应gps干扰抑制方法 | |
CN102411140B (zh) | 一种基于多电台信号的无源雷达处理方法 | |
CN101291165B (zh) | 一种多天线系统的序列检测方法和装置 | |
He et al. | Cramer–Rao bound for MIMO radar target localization with phase errors | |
WO2010122370A1 (en) | Orientation and localization system | |
CN102135617A (zh) | 双基地多输入多输出雷达多目标定位方法 | |
US20120286994A1 (en) | Method and system for locating interferences affecting a satellite-based radionavigation signal | |
WO2006122316A2 (en) | Method for detecting navigation beacon signals using two antennas or equivalent thereof | |
CN104502906A (zh) | 基于rmdcft的空间超高速机动目标检测方法 | |
CN102608583B (zh) | 时频域混叠信号的时差和频差估计方法 | |
CN105137423A (zh) | 一种穿墙雷达对多运动目标实时探测、分离的方法 | |
CN115804170A (zh) | 移动设备的感测辅助定位 | |
US8879679B2 (en) | Frequency offset compensation improvement for a long term evolution (LTE) searcher | |
CN108337197B (zh) | 一种直射径提取方法及装置 | |
JP4183134B2 (ja) | スマートアンテナ及びそのビーム形成方法と装置 | |
CN104502924A (zh) | 一种gps信号处理方法和装置 | |
EP4014061A1 (en) | Radio station for client localization in multipath indoor environment | |
Xiao et al. | Efficient Target Search and Detection in RIS-aided Integrated Sensing and Communications System | |
CN113358931B (zh) | 一种基于互功率谱的时差计算方法 | |
WO2015162669A1 (ja) | 測位装置 | |
Ge et al. | V2X sidelink positioning in FR1: Scenarios, algorithms, and performance evaluation | |
CN101552630B (zh) | 一种基于航空通信信道下的波束成形方法 | |
Buvarp et al. | Probability-reduction of geolocation using reconfigurable intelligent surface reflections | |
WO2016202746A1 (en) | A method of processing offset carrier modulated ranging signals |
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 |