CN107607934B - 一种时差、频差、频差变化率联合估计方法 - Google Patents
一种时差、频差、频差变化率联合估计方法 Download PDFInfo
- Publication number
- CN107607934B CN107607934B CN201710773846.2A CN201710773846A CN107607934B CN 107607934 B CN107607934 B CN 107607934B CN 201710773846 A CN201710773846 A CN 201710773846A CN 107607934 B CN107607934 B CN 107607934B
- Authority
- CN
- China
- Prior art keywords
- frequency difference
- estimation
- change rate
- difference
- time difference
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种时差、频差、频差变化率联合估计方法,其特征在于包括以下内容:1)两路接收机对目标辐射源信号进行时间同步采集,采集的数字信号分别为s1(n)和s2(n);2)对目标辐射源信号在时差范围和频差变化率范围内进行网格划分,对每个网格点计算二阶模糊函数H(τ,f),通过搜索每个网格点的二阶模糊函数H(τ,f)模值,得到时差和频差变化率的初始估计3)利用频差变化率初始估计通过互模糊函数CAF(τ,f)对频差扩展问题进行补偿,得到新的时差估计和频差估计;4)利用新的频差估计,通过频差变化率互相关函数,得到新的频差变化率估计。本发明适用于高动态条件下的时频差无源定位。
Description
技术领域
本发明是关于一种时差、频差、频差变化率联合估计方法,涉及时频差无源定位技术领域。
背景技术
无源定位是一种通过接收辐射源自身的辐射信号从而确定其位置的定位技术。由于该方法具有较强的抗干扰能力,较好的隐蔽性,所以近些年受到了各国的重视和广泛的应用,无源定位常用的方法是利用到达观测点的时间差(时差)和多普勒频率差(频差)等信息来确定目标辐射源的空间几何位置。
现有的时差和频差估计方法仅仅考虑时差和频差的观测量。在此条件下,为了提升参数的估计精度,需要增加观测时间。当观测时间增加时,频差将发生扩展。当频差发生扩展时,会降低频差的估计精度,从而导致后续无源定位的性能下降。为了解决频差的扩展,需要用频差变化率对频差的扩展进行补偿。由于时差、频差、频差变化率在实际的观测过程中都是未知的,因此需要完成对时差、频差、频差变化率的联合估计,解决联合估计过程中三维搜索带来的计算复杂度问题。
发明内容
针对上述问题,本发明的目的是提供一种高效及高精度的时差、频差、频差变化率联合估计方法。
为实现上述目的,本发明采取以下技术方案:一种时差、频差、频差变化率联合估计方法,其特征在于包括以下内容:1)两路接收机对目标辐射源信号进行时间同步采集,采集的数字信号分别为s1(n)和s2(n);2)对目标辐射源信号通过二阶模糊函数变换,得到时差和频差变化率的初始估计3)利用频差变化率初始估计通过互模糊函数对频差扩展问题进行补偿,得到新的时差估计和频差估计;4)利用新的频差估计通过频差变化率互相关函数,得到新的频差变化率估计。
进一步地,所述步骤1)数字信号s1(n)和s2(n)为:
s1(n)=s(n)+n1(n)-N/2≤n≤N/2
式中,s(n)表示目标辐射源信号,n1(n)和n2(n)分别表示接收机噪声信号,τ表示两路接收机采集到信号的到达时间差、fd表示两路接收机采集到信号的到达多普勒频率差、n表示采集数字信号的序号、N表示总的采集点数。
进一步地,所述步骤2)对目标辐射源信号通过二阶模糊函数变换,得到时差和频差变化率的初始估计的具体过程为:
2.1)对目标辐射源信号在时差范围和频差变化率范围内进行网格划分,假设时差的取值范围为τmin~τmax,步长为Δτ,频差变化率的取值范围为搜索步长为
2.2)对于网格内的任意网格点计算目标辐射源时差信号s1(n)和频差信号s2(n)的共轭积信号zτ(n):
式中,表示s1(n)的共轭信号,s2(n+τ)表示对s2(n)延迟τ之后的信号;
2.3)计算zτ(n)的二阶混合瞬时矩xτ(n):
式中,zτ(n+N/4)表示对zτ(n)延迟N/4的信号,表示对zτ(n)延迟-N/4并取共轭之后的信号;
2.4)计算二阶混合瞬时矩xτ(n)的二阶模糊函数
式中,j表示虚数单位,表示网格内的时差和频差变化率;
2.5)对取模;
2.6)对网格内的所有网格点均重复步骤2.2)~2.5),通过搜索所有网格点模值的极大值,得到时差和频差变化率的初始估计
进一步地,所述步骤3)利用频差变化率初始估计通过互模糊函数CAF(τ,f)对频差扩展问题进行补偿,得到新的时差估计和频差估计的具体过程为:
3.1)对目标辐射源信号在时差范围和频差范围内进行网格划分,假设时差的取值范围为τmin~τmax,步长为Δτ,频差的取值范围为fmin~fmax,搜索步长为Δf;
3.2)将频差变化率初始估计代入互模糊函数:
3.3)计算网格内的任意网格点(τ,f)的CAF(τ,f),并对CAF(τ,f)取模,通过搜索CAF(τ,f)模值最大值点的位置,得到新的时差估计以及频差估计
进一步地,所述步骤4)利用新的频差估计通过频差变化率互相关函数,得到新的频差变化率估计的具体过程为:
4.1)确定目标辐射源信号的频差变化率的取值范围为搜索步长为
4.2)对于频差变化率的取值范围的任意点分别计算其互相关函数:
其中,分别表示所述步骤3)得到的新的时差估计和频差估计;
4.3)搜索最大值点,最大值点对应的为频差变化率的估计值,从而得到新的频差变化率估计。
本发明由于采取以上技术方案,其具有以下优点:本发明可以在频差发生扩展条件下,扩展时频差的相干时间,提升时频差参数估计的精度,尤其是频差的精度;同时本发明还可以得到频差变化率的估计,为后续的目标定位提供更精确、更丰富的定位参数,显著提升对目标辐射源的定位精度。本发明适用于高动态条件下的时频差无源定位。
附图说明
图1是本发明联合估计方法的流程示意图;
图2是本发明实施例中时差、频差变化率网格化示意图;
图3是本发明实施例中时差、频差网格化示意图。
具体实施方式
以下结合附图来对本发明进行详细的描绘。然而应当理解,附图的提供仅为了更好地理解本发明,它们不应该理解成对本发明的限制。
如图1所示,本发明提供的时差、频差、频差变化率联合估计方法,包括以下内容:
1、两路接收机对目标辐射源信号进行时间同步采集,采集到的数字信号分别为s1(n)和s2(n):
s1(n)=s(n)+n1(n)-N/2≤n≤N/2
式中,s(n)表示目标辐射源信号,n1(n)和n2(n)分别表示接收机噪声信号,τ表示两路接收机采集到信号的到达时间差、fd表示两路接收机采集到信号的到达多普勒频率差、n表示采集数字信号的序号、N表示总的采集点数。
2、对目标辐射源信号通过二阶模糊函数变换,得到时差和频差变化率的初始估计,即:对目标辐射源信号在时差范围和频差变化率范围内进行网格划分,对每个网格点计算目标辐射源时差信号s1(n)和频差信号s2(n)的二阶模糊函数H(τ,f),通过搜索每个网格点的二阶模糊函数H(τ,f)模值,得到时差和频差变化率的初始估计,具体内容为:
1)如图2所示,对目标辐射源信号在时差范围和频差变化率范围内进行网格划分,假设时差的取值范围为τmin~τmax,步长为Δτ,频差变化率的取值范围为搜索步长为
2)对于网格内的任意网格点计算目标辐射源时差信号s1(n)和频差信号s2(n)的共轭积信号zτ(n):
式中,表示s1(n)的共轭信号,s2(n+τ)表示对s2(n)延迟τ之后的信号。
3)计算zτ(n)的二阶混合瞬时矩xτ(n):
式中,zτ(n+N/4)表示对zτ(n)延迟N/4的信号,表示对zτ(n)延迟-N/4并取共轭之后的信号。
4)计算二阶混合瞬时矩xτ(n)的二阶模糊函数
式中,j表示虚数单位,表示网格内的时差、频差变化率。
5)对取模。
6)对网格内的所有网格点均重复步骤2)~5),通过搜索所有网格点模值的极大值,得到时差和频差变化率的初始估计
3、利用得到频差变化率初始估计通过互模糊函数CAF(τ,f)对频差扩展问题进行补偿,得到新的时差估计和频差估计,具体过程为:
1)如图3所示,对目标辐射源信号在时差范围和频差范围内进行网格划分,假设时差的取值范围为τmin~τmax,步长为Δτ,频差的取值范围为fmin~fmax,搜索步长为Δf。
2)将频差变化率初始估计代入互模糊函数CAF:
式中,表示得到的频差变化率的初始估计。
3)计算网格内的任意网格点(τ,f)的CAF(τ,f),并取模,通过搜索其最大值点的位置,得到新的时差估计以及频差估计
4、利用新的频差估计,通过频差变化率的互相关函数,得到新的频差变化率估计,具体过程为:
1)确定目标辐射源信号的频差变化率的取值范围为搜索步长为
2)则对于频差变化率的取值范围的任意点分别计算其互相关函数:
其中,分别表示由步骤3得到的新的时差估计和频差估计。
3)搜索最大值点,最大值点对应的为频差变化率的估计值,从而得到新的频差变化率估计。
上述各实施例仅用于说明本发明,其中方法的各实施步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (3)
1.一种时差、频差、频差变化率联合估计方法,其特征在于包括以下内容:
1)两路接收机对目标辐射源信号进行时间同步采集,采集的数字信号分别为s1(n)和s2(n),其中数字信号s1(n)和s2(n)为:
s1(n)=s(n)+n1(n) -N/2≤n≤N/2
式中,s(n)表示目标辐射源信号,n1(n)和n2(n)分别表示接收机噪声信号,τ表示两路接收机采集到信号的到达时间差、fd表示两路接收机采集到信号的到达多普勒频率差、n表示采集数字信号的序号、N表示总的采集点数;
2)对目标辐射源信号通过二阶模糊函数变换,得到时差和频差变化率的初始估计具体过程为:
2.1)对目标辐射源信号在时差范围和频差变化率范围内进行网格划分,假设时差的取值范围为τmin~τmax,步长为Δτ,频差变化率的取值范围为搜索步长为
2.2)对于网格内的任意网格点计算目标辐射源时差信号s1(n)和频差信号s2(n)的共轭积信号zτ(n):
式中,表示s1(n)的共轭信号,s2(n+τ)表示对s2(n)延迟τ之后的信号;
2.3)计算zτ(n)的二阶混合瞬时矩xτ(n):
式中,zτ(n+N/4)表示对zτ(n)延迟N/4的信号,表示对zτ(n)延迟-N/4并取共轭之后的信号;
2.4)计算二阶混合瞬时矩xτ(n)的二阶模糊函数
式中,j表示虚数单位,表示网格内的时差和频差变化率;
2.5)对取模;
2.6)对网格内的所有网格点均重复步骤2.2)~2.5),通过搜索所有网格点模值的极大值,得到时差和频差变化率的初始估计
3)利用频差变化率初始估计通过互模糊函数对频差扩展问题进行补偿,得到新的时差估计和频差估计;
4)利用新的频差估计通过频差变化率互相关函数,得到新的频差变化率估计。
2.如权利要求1所述的一种时差、频差、频差变化率联合估计方法,其特征在于,所述步骤3)利用频差变化率初始估计通过互模糊函数CAF(τ,f)对频差扩展问题进行补偿,得到新的时差估计和频差估计的具体过程为:
3.1)对目标辐射源信号在时差范围和频差范围内进行网格划分,假设时差的取值范围为τmin~τmax,步长为Δτ,频差的取值范围为fmin~fmax,搜索步长为Δf;
3.2)将频差变化率初始估计代入互模糊函数:
3.3)计算网格内的任意网格点(τ,f)的CAF(τ,f),并对CAF(τ,f)取模,通过搜索CAF(τ,f)模值最大值点的位置,得到新的时差估计以及频差估计
3.如权利要求2所述的一种时差、频差、频差变化率联合估计方法,其特征在于,所述步骤4)利用新的频差估计通过频差变化率互相关函数,得到新的频差变化率估计的具体过程为:
4.1)确定目标辐射源信号的频差变化率的取值范围为搜索步长为
4.2)对于频差变化率的取值范围的任意点分别计算其互相关函数:
其中,分别表示所述步骤3)得到的新的时差估计和频差估计;
4.3)搜索最大值点,最大值点对应的为频差变化率的估计值,从而得到新的频差变化率估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710773846.2A CN107607934B (zh) | 2017-08-31 | 2017-08-31 | 一种时差、频差、频差变化率联合估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710773846.2A CN107607934B (zh) | 2017-08-31 | 2017-08-31 | 一种时差、频差、频差变化率联合估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107607934A CN107607934A (zh) | 2018-01-19 |
CN107607934B true CN107607934B (zh) | 2019-12-10 |
Family
ID=61056444
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710773846.2A Active CN107607934B (zh) | 2017-08-31 | 2017-08-31 | 一种时差、频差、频差变化率联合估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107607934B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108445445B (zh) * | 2018-02-06 | 2021-09-17 | 中国人民解放军战略支援部队信息工程大学 | 一种分布式无源定位方法及装置 |
CN108763158B (zh) * | 2018-06-01 | 2021-11-09 | 中国人民解放军战略支援部队信息工程大学 | 一种时频差联合计算方法和系统 |
CN109164438B (zh) * | 2018-09-30 | 2020-11-17 | 厦门大学 | 一种基于lfm相干脉冲串的到达时间和到达频率的联合估计方法 |
CN109633586B (zh) * | 2018-12-21 | 2022-07-12 | 电子科技大学 | 一种消除相位模糊度的时延估计方法 |
CN110109091B (zh) * | 2019-05-23 | 2021-11-09 | 中国人民解放军战略支援部队信息工程大学 | 一种针对高速目标的无源雷达参数估计方法及装置 |
JP7178960B2 (ja) * | 2019-06-05 | 2022-11-28 | 三菱電機株式会社 | 位置推定装置、位置推定システム、位置推定方法、およびプログラム |
CN110275134B (zh) * | 2019-06-27 | 2021-03-05 | 清华大学 | 一种基于虚拟到达频差的非共视连续信号无源定位方法 |
CN110703207B (zh) * | 2019-09-27 | 2023-02-24 | 天津大学 | 无源定位的低频多普勒频率差测量方法及装置 |
CN110632556B (zh) * | 2019-09-29 | 2022-08-05 | 电子科技大学 | 一种对静辐射源目标的微弱信号检测并定位的方法 |
CN111175734B (zh) * | 2019-11-28 | 2023-02-10 | 西安电子科技大学 | 一种基于虚拟等效合成孔径的辐射源定位测速方法及系统 |
CN111090109B (zh) * | 2019-12-27 | 2023-08-18 | 中国航天科工集团八五一一研究所 | 星载频差提取对频差快速变化的补偿方法 |
CN112307633B (zh) * | 2020-11-03 | 2022-10-04 | 清源智翔(重庆)科技有限公司 | 多站时频差变参量迭代定位方法及系统 |
CN116203501B (zh) * | 2023-04-27 | 2024-01-05 | 中国人民解放军战略支援部队航天工程大学 | 基于频域互模糊函数插值映射辐射源无源定位方法和设备 |
CN116593959B (zh) * | 2023-05-17 | 2024-03-15 | 中国人民解放军战略支援部队航天工程大学 | 基于载频搜索的互模糊函数映射辐射源定位方法和系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10268037A (ja) * | 1997-03-26 | 1998-10-09 | Nec Corp | 目標信号距離測定装置 |
CN101692163A (zh) * | 2009-09-24 | 2010-04-07 | 中国计量科学研究院 | 一种频率标准远程校准方法及其系统 |
CN103760522A (zh) * | 2014-01-09 | 2014-04-30 | 电子科技大学 | 用于时差估计与多站时钟误差校准的方法及系统 |
CN105607096A (zh) * | 2015-08-31 | 2016-05-25 | 中国电子科技集团公司第三十六研究所 | 一种双星时差频差定位方法和定位装置 |
CN105759241A (zh) * | 2015-12-18 | 2016-07-13 | 中国航天科工集团八五研究所 | 基于时差频差的直接定位方法 |
-
2017
- 2017-08-31 CN CN201710773846.2A patent/CN107607934B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10268037A (ja) * | 1997-03-26 | 1998-10-09 | Nec Corp | 目標信号距離測定装置 |
CN101692163A (zh) * | 2009-09-24 | 2010-04-07 | 中国计量科学研究院 | 一种频率标准远程校准方法及其系统 |
CN103760522A (zh) * | 2014-01-09 | 2014-04-30 | 电子科技大学 | 用于时差估计与多站时钟误差校准的方法及系统 |
CN105607096A (zh) * | 2015-08-31 | 2016-05-25 | 中国电子科技集团公司第三十六研究所 | 一种双星时差频差定位方法和定位装置 |
CN105759241A (zh) * | 2015-12-18 | 2016-07-13 | 中国航天科工集团八五研究所 | 基于时差频差的直接定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107607934A (zh) | 2018-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107607934B (zh) | 一种时差、频差、频差变化率联合估计方法 | |
CN103916201B (zh) | 一种天线信号初始相位差、时延和频率差估计装置与方法 | |
CN104049528B (zh) | 一种北斗授时方法及卫星导航接收机 | |
CN108845973B (zh) | 一种基于改进Quinn算法的多普勒频率估计方法 | |
CN102857961A (zh) | 一种存在频率偏移的通信信号时差测量方法 | |
CN103323822A (zh) | 一种估计通道误差的方法及装置 | |
CN109507704A (zh) | 一种基于互模糊函数的双星定位频差估计方法 | |
CN108196267B (zh) | 一种基于gnss cp技术的不间断时间传递方法 | |
CN103760522B (zh) | 用于时差估计与多站时钟误差校准的方法及系统 | |
CN109188362B (zh) | 一种麦克风阵列声源定位信号处理方法 | |
CN103856431B (zh) | 基于快速傅里叶变换的mpsk调制下开环载波同步算法 | |
CN111505584B (zh) | 基于相对速度和多普勒效应的音频信号到达时间修正方法 | |
CN110049440B (zh) | 一种水下无线传感网同步与定位方法及系统 | |
CN102183770A (zh) | 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法 | |
CN103580718A (zh) | 一种低信噪比下的快速时频同步方法 | |
CN109031261B (zh) | 一种时差估计方法及装置 | |
CN110161501A (zh) | 一种多时序sar图像的目标区域地表起伏信息提取方法 | |
CN104735713A (zh) | 一种适应于宽带无线定位系统中的精确信号传播时延估计方法及实现装置 | |
CN103823210B (zh) | 一种非合作式星地双基地sar时频同步方法 | |
CN106772510B (zh) | 一种基于载波相位测量的跳频测距方法 | |
CN103338024A (zh) | 天线组阵中时延的互补卡尔曼滤波装置与方法 | |
CN102624657B (zh) | 频差估计方法 | |
CN104007425B (zh) | 一种星间时差测量方法及系统 | |
CN110445740A (zh) | 基于重复序列的频偏估计方法及系统 | |
CN103259638B (zh) | 一种异地本振一致性误差下的基带时差估计方法 |
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 |