CN104678382B - 一种适用于扩频体制通信测控***下的远程高精度测距方法 - Google Patents
一种适用于扩频体制通信测控***下的远程高精度测距方法 Download PDFInfo
- Publication number
- CN104678382B CN104678382B CN201310636701.XA CN201310636701A CN104678382B CN 104678382 B CN104678382 B CN 104678382B CN 201310636701 A CN201310636701 A CN 201310636701A CN 104678382 B CN104678382 B CN 104678382B
- Authority
- CN
- China
- Prior art keywords
- distance
- centerdot
- precision
- ranging
- moment
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S11/00—Systems for determining distance or velocity not using reflection or reradiation
- G01S11/02—Systems for determining distance or velocity not using reflection or reradiation using radio waves
- G01S11/08—Systems for determining distance or velocity not using reflection or reradiation using radio waves using synchronised clocks
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开一种适用于扩频体制通信测控***下的远程高精度测距方法。采用的技术方案是:1)采集数据,2)分析数据,3)建立模型,4)卡尔曼滤波,5)输出数据;由于本发明采用了优化的卡尔曼滤波方法,创新性地应用于远程测距***,充分利用了卡尔曼滤波的预测估计特性,从而使测距***的精度提高了一个数量级,因而也解决了远程测距方法***复杂精度不高的问题,取得了良好的效果。本发明具有不增加***复杂度,易于***集成,易于软件实现的特点。主要应用于无线数据链领域。
Description
技术领域
本发明属于无线数据链领域,涉及一种相关数据链的高精度无线电导航定位方法,尤其是一种扩频体制通信测控***的远程高精度测距方法。
背景技术
无线通信中,由于信号传输环境复杂多变,需要多种导航定位方法保证***可靠运行。无线电导航定位作为GPS等卫星定位信号失灵时的一种备用导航模式,应用范围非常广泛,而高精度测距是保障无线电导航定位的必要手段。目前测控***中测距方法主要有侧音测距及伪码测距。
侧音测距是利用正弦信号经过往返距离所对应的相位变化进行距离测量。其中低侧音保证测距距离,高侧音保证测距精度。但在距离增大时,纯侧音测距***由于设备复杂、抗干扰能力差等特点,很难实现远距离测量。
如图1所示,伪码测距在扩频体制通信测控***中虽然具有实现方便、灵活、抗干扰能力强等特点,但由于普通伪码测距的测量距离取决于伪码长度,测距精度受限于时钟速率,因此在要求中远程测距就要求伪码尽可能长,这与通信测控***要求快速捕获与跟踪形成矛盾;同时测距精度也受限于***的时钟频率。
因此,如何保证远程通信测控***的快速捕获与跟踪性能,既不大幅度提高***复杂度(***时钟频率提高会提高***复杂度),又能保证远程距离测量的高精度,成为亟待解决的问题。
发明内容
本发明公开了一种适用于扩频体制通信测控***下的远程高精度测距方法,旨在提供一种适用于远程通信测控***下的高精度测距方法,通过对上一时刻测量值与当前时刻估算值进行滤波运算,最终达到高精度测距的目的。
本发明采用如下的技术方案:
本方法分下述几步进行:
1)采集数据:
接收端机通过AD采样器量化接收扩频信号,通过相关器进行解扩,提取信号作为测距与通信解算模块的输入;
如图3所示,接收端机通过AD采样器量化接收扩频信号,通过相关器进行解扩,把隐没于噪声下的有用信号提取出来,作为测距与通信解算模块的输入。相关器设计过程中,伪随机码越长,捕获时间越长,在码片速率受限的***中,普通伪码测距实现远程高精度测距所带来的时间开销和***复杂度是不可忍受的。因此,普通的伪码测距方法测量距离受限于伪码长度,测距精度受限于码片速率与***时钟。
2)分析数据:
对被测量对象的运动建立动态模型,并将最优估计方法应用于动态模型,对测距误差进行滤波估计处理;
分析上述运算过程,我们发现***的同步误差与测距误差是该方法的关键量,远程高精度测距的过程即是减少这两个量的过程;由此,我们可以创造性地应用扩频吗同步+卡尔曼滤波的方法,也即对被测量对象的运动建立动态模型,并将最优估计方法应用于动态模型,对测距误差进行滤波估计处理,就可以有效减小误差,提高测距精度;
3)建立模型:
以测量者为基准点,被测对象与测量者之间的距离作为被测对象在坐标系中的坐标值,建立***模型;
伪码测距方法的测量值是基于实时测量的测距值,即某一时刻被测对象与测量者之间的距离变化情况;故此可以等效建立一维坐标系,即以测量者为基准点,被测对象与测量者之间的距离为被测对象在该坐标系中的坐标值,状态方程为:
其中x(t)为距离值,为距离的一次变率,为距离的二次变率,为距离的三次变率,为当前加速度的均值,τ为时间常数,ω(t)为高斯白噪声;
4)卡尔曼滤波:
设定***参数和初始数值,进行卡尔曼滤波;卡尔曼滤波的更新方程:
其中,Xi为i时刻的被估计状态;Φi,i-1为i-1时刻至i时刻的一步转移矩阵;Xi-1为i-1时刻的被估计状态;Γi-1为***噪声驱动阵;Hi为量测阵;Zi为i时刻的被观测状态;Vi为量测噪声序列;Wi-1为***激励噪声序列,{ui-1}、{Bi-1}与{Yi}为已知的确定性输入序列;
5)输出数据:
通过上述步骤,对测距误差实时进行滤波估计与矫正,逐步减小测距误差,输出测距数据;达到远程高精度测距的目的。
本发明依托远程通信测控***而开发,该***采用直接扩频体制,要求捕获和同步时间不大于10ms,测距误差不大于0.1m。本发明解决了上述***要求。
本发明的积极效果:
由于本发明采用了优化的卡尔曼滤波方法,创新性地应用于远程测距***,充分利用了卡尔曼滤波的预测估计特性,从而使测距***的精度提高了一个数量级,因而也解决了远程测距方法***复杂精度不高的问题,取得了良好的效果。
本方法基于以下思想:普通的伪码测距方法具有在扩频体制下易于实现、资源利用率高的优点,但是测量距离受限于伪码长度,测距精度受限于码片速率与***时钟;而卡尔曼滤波理论成熟,可以通过对上一时刻测量值与当前时刻估算值进行滤波运算,减小***误差;结合上述两者优点,提出了基于扩频体制远程通信测控***的高精度测距方法,并成功地将该方法应用于相关的数据链***中。
本方法与目前其他方法相比,具有不增加***复杂度,易于***集成,易于软件实现的特点。
附图说明
图1是普通的伪码测距方法示意图。
图2是本方法加以改进后的伪码测距方法示意图。
图3是整个方法的***模型示意图。
图4针对本方法测量误差仿真的效果图。
具体实施方式
下面,结合附图和具体实施例,对发明作进一步的说明。
实施例1
如图3所示,接收端机通过AD采样器量化接收扩频信号,通过相关器进行解扩,把隐没于噪声下的有用信号提取出来,作为测距与通信解算模块的输入。相关器设计过程中,伪随机码越长,捕获时间越长,在码片速率受限的***中,普通伪码测距实现远程高精度测距所带来的时间开销和***复杂度是不可忍受的。因此,普通的伪码测距方法测量距离受限于伪码长度,测距精度受限于码片速率与***时钟。
普通的伪码测距方法测量距离受限于伪码长度,测距精度受限于码片速率与***时钟。为了能够保证在通信测控***下伪码测距的远程高精度,必须对普通的伪码测距方法进行改进,如图2所示。在此基础上,我们对被测量对象的运动建立动态模型,并将最优估计方法应用于动态模型,对测距误差进行滤波估计处理,就可以有效减小误差,提高测距精度。
因此,关键是动态模型的建立与最优估计方法的选取。由上可知,伪码测距方法的测量值是基于实时测量的测距值,即某一时刻被测量对象与测量者之间的距离变化情况。故此可以等效建立一维坐标系,即以测量者为基准点,则被测量对象与测量者之间的距离即为被测对象在该坐标系中的坐标值。在此基础上建立初步的动态模型。状态方程为:
其中x(t)为距离值,为距离的一次变率,为距离的二次变率,为距离的三次变率,为当前加速度的均值,τ为时间常数,ω(t)为高斯白噪声。整个***的框架结构如图3所示:
其次,考虑到计算的复杂度,收敛时间,以及收敛结果的无偏性,在本发明中选取卡尔曼滤波方法,该方法的更新方程
其中,Xi为i时刻的被估计状态;Φi,i-1为i-1时刻至i时刻的一步转移矩阵;Xi-1为i-1时刻的被估计状态;Γi-1为***噪声驱动阵;Hi为量测阵;Zi为i时刻的被观测状态;Vi为量测噪声序列;Wi-1为***激励噪声序列,{ui-1}、{Bi-1}与{Yi}为已知的确定性输入序列。
通常情况下,估计值与观测值是一一对准的。测距应用中,距离二次变率变化范围和强度均比较小,即状态突变很小。实验表明,此范围内的状态变化不会影响卡尔曼滤波效果,也不会引入时延。所以,滤波前后的距离值是一一对准的。
根据被测量对象与测量者之间的相对距离关系,可以模拟一段距离曲线,作为真实值。在此基础上添加高斯白噪声,模拟原始观测值,作为卡尔曼滤波的输入,比较滤波前后之间的测距误差。为便于仿真,参数选取距离二次变率方差为1,观测噪声方差R=1.22,加速度时间相关常数为1,时间相关常数为0.5,仿真结果如图4所示。
由此可见,通过仿真可以证明该普通发明的远程测距精度不大于0.1m,完全满足***要求。
以上所述仅为本发明的较佳可行实施例,并非因此局限本发明的专利范围,故凡是运用本发明说明书及附图内容所作的等效结构变化,均包含于本发明的保护范围。
Claims (3)
1.一种适用于扩频体制通信测控***下的远程高精度测距方法,包括:
1)采集数据:接收端机通过AD采样器量化接收扩频信号,通过相关器进行解扩,提取信号作为测距与通信解算模块的输入;
2)分析数据:对被测量对象的运动建立动态模型,并将最优估计方法应用于动态模型,对测距误差进行滤波估计处理;
3)建立模型:以测量者为基准点,被测对象与测量者之间的距离作为被测对象在坐标系中的坐标值,建立***模型;
4)卡尔曼滤波:设定***参数和初始数值,进行卡尔曼滤波;
5)输出数据:通过上述步骤,对测距误差实时进行滤波估计与矫正,逐步减小测距误差,输出测距数据。
2.根据权利要求1所述的适用于扩频体制通信测控***下的远程高精度测距方法,其特征在于,步骤3)建立模型中,所述***模型的状态方程为:
其中x(t)为距离值,为距离的一次变率,为距离的二次变率,为距离的三次变率,为当前加速度的均值,τ为时间常数,ω(t)为高斯白噪声。
3.根据权利要求1所述的适用于扩频体制通信测控***下的远程高精度测距方法,其特征在于,步骤4)卡尔曼滤波中的更新方程:
其中,Xi为i时刻的被估计状态;Φi,i-1为i-1时刻至i时刻的一步转移矩阵;Xi-1为i-1时刻的被估计状态;Γi-1为***噪声驱动阵;Hi为量测阵;Zi为i时刻的被观测状态;Vi为量测噪声序列;Wi-1为***激励噪声序列,{ui-1}、{Bi-1}与{Yi}为已知的确定性输入序列。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310636701.XA CN104678382B (zh) | 2013-11-29 | 2013-11-29 | 一种适用于扩频体制通信测控***下的远程高精度测距方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310636701.XA CN104678382B (zh) | 2013-11-29 | 2013-11-29 | 一种适用于扩频体制通信测控***下的远程高精度测距方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104678382A CN104678382A (zh) | 2015-06-03 |
CN104678382B true CN104678382B (zh) | 2017-07-14 |
Family
ID=53313709
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310636701.XA Active CN104678382B (zh) | 2013-11-29 | 2013-11-29 | 一种适用于扩频体制通信测控***下的远程高精度测距方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104678382B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102721951A (zh) * | 2012-05-04 | 2012-10-10 | 西安电子科技大学 | 一种高机动目标跟踪方法 |
CN102749633A (zh) * | 2012-06-29 | 2012-10-24 | 北京航空航天大学 | 一种卫星导航接收机的动态定位解算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2002315595B2 (en) * | 2001-06-04 | 2007-06-07 | Novatel Inc. | An inertial/GPS navigation system |
-
2013
- 2013-11-29 CN CN201310636701.XA patent/CN104678382B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102721951A (zh) * | 2012-05-04 | 2012-10-10 | 西安电子科技大学 | 一种高机动目标跟踪方法 |
CN102749633A (zh) * | 2012-06-29 | 2012-10-24 | 北京航空航天大学 | 一种卫星导航接收机的动态定位解算方法 |
Non-Patent Citations (5)
Title |
---|
On the Design of Direct Sequence Spread-Spectrum Signaling for Range Estimation;Brian Bingham et al.;《OCEANS 2007》;20071004;第1-7页 * |
Research on INS/Terrestrial Pseudo-range Integrated Navigation;Zhang Kang et al.;《2010 International Conference on Computer Application and System Modeling》;20101024;第9卷;第V9-262-V9-265页 * |
伪码扩频航天测控体制设计;任松等;《信息与电子工程》;20070630;第5卷(第3期);第165-169页 * |
航天扩频测控通信***中伪码测距方法及精度分析;唐军等;《电讯技术》;20061231(第4期);第91-95页 * |
高度机动目标的"当前"统计Jerk 模型;乔向东等;《***工程与电子技术》;20021231;第24卷(第6期);第53-56页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104678382A (zh) | 2015-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103033828B (zh) | 一种高灵敏度北斗辅助授时装置和授时接收机及授时方法 | |
US8170082B2 (en) | Cross-talk mitigation in global navigation satellite systems | |
CN106291619B (zh) | 一种导航卫星星间链路无线长周期扩频码信号高性能捕获方法 | |
CN100553168C (zh) | 一种长周期扩频码的频率折叠时频并行搜索方法 | |
CN104280746A (zh) | 一种惯性辅助gps的深组合半实物仿真*** | |
CN102183770A (zh) | 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法 | |
CN101902288B (zh) | 一种直序扩频二元相移键控调制器时延的测量方法 | |
CN108011653B (zh) | 基于自适应的快速捕获跟踪***及方法 | |
CN101839987A (zh) | 一种自适应gps软件接收机的实现方法 | |
CN105099498A (zh) | 扩频信号捕获***和方法 | |
CN104581926A (zh) | 一种低轨卫星通信的上行链路准同步时间精确测量方法 | |
CN102279402B (zh) | 一种gnss数据/导频混合信号的频域相干联合捕获方法 | |
CN106054214A (zh) | 一种可削弱多路径效应的gnss数据处理方法 | |
CN101567709A (zh) | 一种减少多径对接收机天线定位精度影响的方法与装置 | |
CN103869337A (zh) | 跟踪全球导航卫星***信号的方法和设备及其接收机 | |
CN103809192A (zh) | 一种gnss接收机的动态校正算法 | |
CN105242287B (zh) | 基于gpu和imu的卫星导航软件接收机及其导航方法 | |
CN106104301A (zh) | 用于无线电接收机的处理器 | |
CN103809193B (zh) | 提高微弱gnss信号处理增益的捕获*** | |
CN102243309A (zh) | Gnss互相关干扰抑制方法和装置 | |
CN104735713A (zh) | 一种适应于宽带无线定位***中的精确信号传播时延估计方法及实现装置 | |
CN102565825A (zh) | 接收信号可靠度判定装置、方法及码相位误差算出方法 | |
CN102565816B (zh) | 一种载波平滑伪距分组平滑方法和装置 | |
CN109541651A (zh) | 一种gpsl1l2双频信号联合捕获方法 | |
CN104678382B (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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20180404 Address after: 300308 Tianjin District of Dongli City Airport Economic Zone bonded Road No. 357 Patentee after: TIANJIN JINHANG INSTITUTE OF COMPUTING TECHNOLOGY Address before: 300308 Tianjin District of Dongli City Airport Economic Zone bonded Road No. 357 Patentee before: NO.8357 Research Institute of the Third Academy of China Aerospace Science & Industry Corp. |
|
TR01 | Transfer of patent right |