CN105093239B - 一种基于温度补偿的系统时延误差校正方法 - Google Patents
一种基于温度补偿的系统时延误差校正方法 Download PDFInfo
- Publication number
- CN105093239B CN105093239B CN201510519287.3A CN201510519287A CN105093239B CN 105093239 B CN105093239 B CN 105093239B CN 201510519287 A CN201510519287 A CN 201510519287A CN 105093239 B CN105093239 B CN 105093239B
- Authority
- CN
- China
- Prior art keywords
- sigma
- temperature
- regression model
- formula
- distance measurement
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/23—Testing, monitoring, correcting or calibrating of receiver elements
- G01S19/235—Calibration of receiver components
Abstract
一种基于温度补偿的系统时延误差校正方法,在精密测距产品内部设置测温点并将温度遥测采集,通过对产品进行温度‑伪距的测量得到先验值,建立n阶多项式回归模型进行拟合残差分析,当随着阶数的增加残差精度压缩或小于预设阀值,即可选取压缩点处/小于预设阀值作为该多项式回归模型的阶数。根据最小均方根准则估计n阶多项式回归模型时延校正系数,根据时延校正系数和实时采集的温度遥测进行精密测距产品温变下的时延校正,本发明通过在产品硬件上增加温度遥测采集作为校正输入量有效的解决了精密测距产品由于温度变化引起的测距系统误差,此外,本发明中的校正方法简单可靠,整体资源占用量较少。
Description
技术领域
本发明涉及一种系统时延误差校正方法,特别是特别涉及一种星载产品在大范围工作温度变化条件下系统时延测量要求较高的校正方法,通过在产品硬件上增加温度遥测采集作为校正输入量校正精密测距产品由于温度变化引起的测距系统误差,达到星载导航接收机工程化的目的,属于精密测距领域。
背景技术
星间精密测距产品作为精密测距系统的核心单机,接收扩频测距信号进行伪距的精密测量,为了保证精密测距系统的精度,伪距的系统时延变化要求小于0.5ns。卫星在轨运行的整个生命周期内,星间精密测距产品的环境温度会不断发生变化。环境温度的变化会导致星间精密测距产品测距系统的零值时延发生改变,进而直接影响到测距的系统误差(经测试,-25℃~+60℃范围内变化大于1ns),该误差对于高精度测距系统是不可接受,因此为消除环境温度变化对测距误差的影响,需要采取一种时延误差校正方法解决该问题。
目前未在公开的文献上查询到关于温度下时延变化的处理方法。
国内外对于时延随环境变化的控制方法主要有以下几种:
《双向时间同步系统的设备时延校准技术研究》,李星、耿淑敏、李垣陵,中国空间科学,2011年8月中通过通过在测量系统中增加用于测量发射/接收通道时延变化的系统测量回路,将测量结果用于校正由于环境变化引起的系统测量误差,该方法的主要缺点是为了校正由于环境、老化等因素引起的时延变化增加了测量时延变化系统测量回路,整个系统的设计硬件、软件资源额外增加。
《影响卫星导航系统测距设备时延特性的因素浅析》,高帅、宫磊、董继峰,第五届中国卫星导航学术年会论文集,2014年5月中的测距系统设备中选取温度时延随温度变化较小惰性器件,同时对环境影响系统时延较大设备温控。
产品设计中对测距系统有源设备选取温度时延随温度变化较小惰性器件,此外在建立测量系统中时延随温度变化较大的设备进行温度控制,对裸露在外部无法进行温控设备(无源设备)选择时延随温度变化系数较小的设备,该方法的主要缺点对系统中时延随温度变化较大设备需要采取专用的温控措施,对裸露在外部无法进行温控设备在产品设计时需要对产品的器件温度特性进行筛选,约束了产品的器件选型范围,增加了产品的设计复杂难度。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供了一种基于温度补偿的系统时延误差校正方法,通过在产品硬件上增加温度遥测采集作为校正输入量有效的解决了精密测距产品由于温度变化引起的测距系统误差,此外,本发明中的校正方法简单可靠,整体资源占用量较少。
本发明的技术解决方案是:一种基于温度补偿的系统时延误差校正方法,步骤如下:
(1)对待校正产品进行单向伪距值测量Ri,i=1,2,...,M,并选取待校正产品的一处温度测量点采集该产品的遥测温度ti,i=1,2,...,M;
(2)建立遥测温度ti与拟合得到的伪距值的n阶多项式回归模型,具体由公式:
Fi=an(ti-t0)n+an-1(ti-t0)n-1...+a2(ti-t0)2+a1(ti-t0)1+a0
给出,式中,an,an-1,…,a0为n阶多项式回归模型的拟合系数,ti为温度遥测量,t0为预先给定的基准温度,Fi为拟合得到的伪距值;
(3)从预先给定的n的取值范围中选取n的最小值,计算n阶多项式回归模型的拟合系数an,an-1,…,a0,具体由公式:
a=A/b
给出,其中,A为(n+1)×(n+1)的矩阵,由公式:
给出;a为(n+1)×1的矩阵,由公式:
a=[an,an-1,...,a1,a0]T
给出;
b为(n+1)×1的矩阵,由公式:
给出,其中Ri为温度遥测量ti时对应的原始测距值;
(4)利用步骤(3)中确定的n阶多项式回归模型的拟合系数an,an-1,…,a0,计算n阶多项式回归模型的拟合残差,n阶多项式回归模型的拟合残差平方和Sr;
(5)利用步骤(4)中计算得到的n阶多项式回归模型的拟合残差,计算步骤(3)中确定的n阶多项式回归模型的拟合系数an,an-1,…,a0的精度δ;
(6)若步骤(5)中计算的精度δ小于等于预设的精度要求,则进入步骤(7),否则,将n的取值加1,返回步骤(3),直到选取满足预设精度要求的拟合阶数n;
(7)计算遥测温度ti对应的原始测距值补偿量Rbi;
(8)利用步骤(7)中的原始测距值补偿量Rbi对原始测距值进行补偿,具体由公式:
RFi=Rbi+Ri,i=1,2,...,M
给出。
所述步骤(4)中n阶多项式回归模型的拟合残差平方和Sr具体由公式:
给出,其中Ri为温度遥测量ti时对应的原始测距值。
所述步骤(5)中计算步骤(3)中确定的n阶多项式回归模型的拟合系数an,an-1,…,a0的精度δ,具体由公式:
δ=(Sr/(M-(n+1)))0.5
给出。
所述步骤(7)中计算遥测温度ti对应的原始测距值补偿量Rbi,具体由公式:
Rbi=-(an(ti-t0)n+an-1(ti-t0)n-1...+a2(ti-t0)2+a1(ti-t0)1),i=1,2,...,M
给出。
本发明与现有技术相比的有益效果是:
(1)本发明提出了一种星载产品在大范围工作温度变化条件下系统时延测量要求较高的校正方法,通过在产品硬件上增加温度遥测采集作为校正输入量,有效解决了精密测距产品由于温度变化引起的测距系统误差的问题;
(2)本发明通过采集温度-原始测距值作为初始校正模型的输入量,采用拟合残差作为校正模型阶数选取标准,可保证选取合理的多项式回归阶数n满足温度下时延测量的精度需求,该校正方法简单可靠、整体资源占用量少;
(3)本发明通过评估拟合回归后残差进行多项式阶数的选取,实现了整个校正模型精度可量化评估,阶数选取合理;
(4)本发明给出了精确的误差回归模型,通过采集的温度-原始测距值对应关系即可实现精确校正。
附图说明
图1为本发明的流程图。
具体实施方式
下面结合附图对本发明的具体实施方式进行进一步的详细描述。
一种基于温度补偿的系统时延误差校正方法,其特征在于步骤如下:
(1)对待校正产品进行单向伪距值测量Ri,i=1,2,...,M,并选取待校正产品的一处温度测量点采集该产品的遥测温度ti,i=1,2,...,M(实际实现过程中应该采集的为与产品环境温度相匹配的电压遥测,其遥测采集系统变化范围应适应产品实际温度范围,遥测采集系统分辨率小于2℃);
(2)建立遥测温度ti与拟合得到的伪距值的n阶多项式回归模型,具体由公式:
Fi=an(ti-t0)n+an-1(ti-t0)n-1...+a2(ti-t0)2+a1(ti-t0)1+a0
给出,式中,an,an-1,…,a0为n阶多项式回归模型的拟合系数,ti为温度遥测量,t0为基准温度(t0选取测距终端工作温度范围内一点即t0∈[TMIN,TMAX],对补偿后的测距均值需求为Rmean,则选取与Rmean误差最小的单向伪距测量值Ri对应的t0),Fi为拟合得到的伪距值;
(3)从预先给定的n的取值范围中选取n的最小值,计算n阶多项式回归模型的拟合系数an,an-1,…,a0,具体过程为:
根据最小均方差准则,残差平方和计算如下:
其中,Ri为温度遥测量ti时对应的原始测距值;
由上式可得
Sr=∑i=1:M(Ri-an(ti-t0)n+an-1(ti-t0)n-1...+a2(ti-t0)2+a1(ti-t0)1+a0)2
通过上式拟合系数an,an-1……a0,分别求偏导:
…..
令上式等于0,进一步整理可得:
(∑i=1:M(ti-t0)n+n)an+(∑i=1:M(ti-t0)n-1+n)an-1+...+(∑i=1:M(ti-t0)1+n)a1+(∑i=1:M(ti-t0)0+n)a0=∑i=1:M((ti-t0)nRi)
(∑i=1:M(ti-t0)n+n-1)an+(∑i=1:M(ti-t0)n-1+n-1)an-1+...+(∑i=1:M(ti-t0)1+n-1)a1+(∑i=1:M(ti-t0)0+n-1)a0=∑i=1:M((ti-t0)n-1Ri)
……
(∑i=1:M(ti-t0)n+1)an+(∑i=1:M(ti-t0)n-1+1)an-1+...+(∑i=1:M(ti-t0)1+1)a1+(∑i=1:M(ti-t0)0+1)a0=∑i=1:M((ti-t0)1Ri)
(∑i=1:M(ti-t0)n+0)an+(∑i=1:M(ti-t0)n-1+0)an-1+...+(∑i=1:M(ti-t0)1+0)a1+(∑i=1:M(ti-t0)0+0)a0=∑i=1:M((ti-t0)0Ri)
从上式可看出为an,an-1……a0的线性方程,整理可得:
a=A/b
给出,其中,A为(n+1)×(n+1)的矩阵,由公式:
给出;a为(n+1)×1的矩阵,由公式:
a=[an,an-1,...,a1,a0]T
给出;
b为(n+1)×1的矩阵,由公式:
给出,其中Ri为温度遥测量ti时对应的原始测距值;
(4)利用步骤(3)中确定的n阶多项式回归模型的拟合系数an,an-1,…,a0,计算n阶多项式回归模型的拟合残差,n阶多项式回归模型的拟合残差平方和Sr具体由公式:
给出,其中Ri为温度遥测量ti时对应的原始测距值;
(5)利用步骤(4)中计算得到的n阶多项式回归模型的拟合残差,计算步骤(3)中确定的n阶多项式回归模型的拟合系数an,an-1,…,a0的精度δ,具体由公式:
δ=(Sr/(M-(n+1)))0.5
给出;
(6)若步骤(5)中计算的精度δ小于等于预设的精度要求,则进入步骤(7),否则,将n的取值加1,返回步骤(3),直到选取满足预设精度要求的拟合阶数n;
(7)计算遥测温度ti对应的原始测距值补偿量Rbi,具体由公式:
Rbi=-(an(ti-t0)n+an-1(ti-t0)n-1...+a2(ti-t0)2+a1(ti-t0)1),i=1,2,...,M
给出;
(8)利用步骤(7)中的原始测距值补偿量Rbi对原始测距值进行补偿,具体由公式:
RFi=Rbi+Ri,i=1,2,...,M
给出。
具体实施例
利用两台精密测距产品进行双向测距,其中一台处于温箱外(标识为精密测距产品A)恒温条件下,另一台(标识为精密测距产品B)处于温箱内的变温条件下。精密测距产品A的接收信号测距时延变化反应了收发信机发射时延随温度的变化规律,精密测距产品B的接收信号测距时延变化反应了收发信机接收时延随温度的变化规律。环境温度分别在4个温度档位上循环切换,包括:30℃、20℃、10℃、0℃。
利用本发明中的方法,温度补偿前后发射时延和接收时延随温度变化数据具体表1和表2所示,其中,表1为发射时延0℃~30℃范围内温度校正前后偏差,表2为接收时延0℃~30℃范围内温度校正前后偏差。
表1
产品环境温度(℃) | 发射时延原始均值偏差(ns) | 发射时延温度校正后偏差(ns) |
0 | 0.43 | -0.03 |
10 | 0.15 | -0.04 |
20 | -0.18 | -0.02 |
30 | -0.42 | 0.01 |
表2
产品环境温度(℃) | 接收时延原始均值偏差(ns) | 接收时延温度校正后偏差(ns) |
0 | -0.28 | 0.03 |
10 | -0.17 | 0.04 |
20 | 0.155 | -0.035 |
30 | 0.27 | -0.07 |
从表1可看出,发射时延在30摄氏度温度范围内时延均值变化约0.85ns,经采用温度补偿后(该补偿模型选用1阶),时延均值在30摄氏度温度范围内变化为0.05ns;从表2可看出,接收时延在30摄氏度温度范围内时延均值变化约0.55ns,经采用温度补偿后(该补偿模型选用1阶),时延均值在30摄氏度温度范围内变化为0.11ns。从该数据可以看出,采用本专利所述方法可以使由于温度变化引起系统误差降低至系统要求范围内。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (3)
1.一种基于温度补偿的系统时延误差校正方法,其特征在于步骤如下:
(1)对待校正产品进行单向伪距值测量Ri,i=1,2,...,M,并选取待校正产品的一处温度测量点采集该产品的遥测温度ti,i=1,2,...,M;
(2)建立遥测温度ti与拟合得到的伪距值的n阶多项式回归模型,具体由公式:
Fi=an(ti-t0)n+an-1(ti-t0)n-1...+a2(ti-t0)2+a1(ti-t0)1+a0
给出,式中,an,an-1,…,a0为n阶多项式回归模型的拟合系数,ti为温度遥测量,t0为预先给定的基准温度,Fi为拟合得到的伪距值;
(3)从预先给定的n的取值范围中选取n的最小值,计算n阶多项式回归模型的拟合系数an,an-1,…,a0,具体由公式:
a=A/b
给出,其中,A为(n+1)×(n+1)的矩阵,由公式:
给出;a为(n+1)×1的矩阵,由公式:
a=[an,an-1,...,a1,a0]T
给出;
b为(n+1)×1的矩阵,由公式:
给出,其中Ri为温度遥测量ti时对应的原始测距值;
(4)利用步骤(3)中确定的n阶多项式回归模型的拟合系数an,an-1,…,a0,计算n阶多项式回归模型的拟合残差,n阶多项式回归模型的拟合残差平方和Sr;
(5)利用步骤(4)中计算得到的n阶多项式回归模型的拟合残差,计算步骤(3)中确定的n阶多项式回归模型的拟合系数an,an-1,…,a0的精度δ;
(6)若步骤(5)中计算的精度δ小于等于预设的精度要求,则进入步骤(7),否则,将n的取值加1,返回步骤(3),直到选取满足预设精度要求的拟合阶数n;
(7)计算遥测温度ti对应的原始测距值补偿量Rbi;
所述步骤(7)中计算遥测温度ti对应的原始测距值补偿量Rbi,具体由公式:
Rbi=-(an(ti-t0)n+an-1(ti-t0)n-1...+a2(ti-t0)2+a1(ti-t0)1),i=1,2,...,M
给出;
(8)利用步骤(7)中的原始测距值补偿量Rbi对原始测距值进行补偿,具体由公式:
RFi=Rbi+Ri,i=1,2,...,M
给出。
2.根据权利要求1所述的一种基于温度补偿的系统时延误差校正方法,其特征在于:所述步骤(4)中n阶多项式回归模型的拟合残差平方和Sr具体由公式:
给出,其中Ri为温度遥测量ti时对应的原始测距值。
3.根据权利要求1所述的一种基于温度补偿的系统时延误差校正方法,其特征在于:所述步骤(5)中计算步骤(3)中确定的n阶多项式回归模型的拟合系数an,an-1,…,a0的精度δ,具体由公式:
δ=(Sr/(M-(n+1)))0.5
给出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510519287.3A CN105093239B (zh) | 2015-08-21 | 2015-08-21 | 一种基于温度补偿的系统时延误差校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510519287.3A CN105093239B (zh) | 2015-08-21 | 2015-08-21 | 一种基于温度补偿的系统时延误差校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105093239A CN105093239A (zh) | 2015-11-25 |
CN105093239B true CN105093239B (zh) | 2017-07-28 |
Family
ID=54574095
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510519287.3A Active CN105093239B (zh) | 2015-08-21 | 2015-08-21 | 一种基于温度补偿的系统时延误差校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105093239B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109541523B (zh) * | 2018-11-30 | 2020-09-18 | 西安空间无线电技术研究所 | 一种干涉仪温变误差的热真空标校系统及标校方法 |
CN109669130A (zh) * | 2018-12-29 | 2019-04-23 | 深圳市新威尔电子有限公司 | 电池内阻精确计算方法、计算装置以及可读存储介质 |
CN109634093B (zh) * | 2019-01-16 | 2021-07-27 | 和芯星通科技(北京)有限公司 | 一种基于gnss接收机的授时方法及gnss接收机 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN86207227U (zh) * | 1986-10-07 | 1987-07-22 | 李振山 | 温度仪表冷端补偿非线性误差校正器 |
CN102650527A (zh) * | 2012-05-25 | 2012-08-29 | 北京航空航天大学 | 一种基于时间序列分析消噪的光纤陀螺温度补偿方法 |
CN102750443A (zh) * | 2012-06-04 | 2012-10-24 | 河海大学 | 一种碾压混凝土坝层面性态综合评价方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4981359A (en) * | 1989-06-19 | 1991-01-01 | Litton Systems, Inc. | Ring laser gyroscope dither drive system and method |
US6959583B2 (en) * | 2002-04-30 | 2005-11-01 | Honeywell International Inc. | Passive temperature compensation technique for MEMS devices |
CN1762091A (zh) * | 2003-03-17 | 2006-04-19 | 精工爱普生株式会社 | 振荡器的特性自动补偿装置、特性自动补偿方法、特性自动补偿程序以及定位信号接收机 |
KR100668910B1 (ko) * | 2006-02-15 | 2007-01-12 | 삼성전자주식회사 | 전세계 위치 확인 시스템을 이용한 이동 단말의 기준시계공유 장치 및 방법 |
US8976065B2 (en) * | 2008-10-15 | 2015-03-10 | Nokia Corporation | Temperature sensor for oscillator and for satellite positioning circuitry |
US8643444B2 (en) * | 2012-06-04 | 2014-02-04 | Broadcom Corporation | Common reference crystal systems |
-
2015
- 2015-08-21 CN CN201510519287.3A patent/CN105093239B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN86207227U (zh) * | 1986-10-07 | 1987-07-22 | 李振山 | 温度仪表冷端补偿非线性误差校正器 |
CN102650527A (zh) * | 2012-05-25 | 2012-08-29 | 北京航空航天大学 | 一种基于时间序列分析消噪的光纤陀螺温度补偿方法 |
CN102750443A (zh) * | 2012-06-04 | 2012-10-24 | 河海大学 | 一种碾压混凝土坝层面性态综合评价方法 |
Non-Patent Citations (2)
Title |
---|
《光纤陀螺温度补偿技术研究》;左瑞芹;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技II辑》;20061115;第25-37页 * |
《硅微机械振动陀螺零偏温度补偿研究》;程龙 等;《传感技术学报》;20080331;第21卷(第3期);第484-485页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105093239A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105093239B (zh) | 一种基于温度补偿的系统时延误差校正方法 | |
CN105911858B (zh) | 一种亚纳秒级远程时间比对/传递系统误差测试方法 | |
CN104215977B (zh) | 一种基于卫星导航系统的精度评估方法及系统 | |
CN108196267B (zh) | 一种基于gnss cp技术的不间断时间传递方法 | |
CN103176400A (zh) | 智能电表时钟校准方法 | |
CN104460304A (zh) | 一种具有自动校正功的高分辨率时间间隔测量装置 | |
CN108459331B (zh) | 多模卫星导航接收机的时延绝对校准方法 | |
CN103698602A (zh) | 一种大动态高精度同步连续频率测量方法 | |
CN105258779A (zh) | 电子秤自动校准方法 | |
CN103983957A (zh) | 一种多普勒偏移测量方法及其装置 | |
CN113258288A (zh) | 一种相控阵天线波束控制装置及控制方法 | |
Zhou et al. | A new method of D-TDOA time measurement based on RTT | |
WO2017070974A1 (zh) | 无线射频接收机灵敏度快速检测方法及装置 | |
CN114647178B (zh) | 基于北斗和地面基准传递的原子钟自动校准方法和系统 | |
CN104155664A (zh) | 一种星载接收机校时功能测试系统及方法 | |
CN114063125A (zh) | 一种基于分时双向测距的精密轨道确定和时间同步方法 | |
CN111157081B (zh) | 一种电子式燃气表的校准方法 | |
CN107817480B (zh) | 一种时间比对测量装置及方法 | |
US11525812B2 (en) | Gas concentration measuring device and method of calibrating same | |
CN112946693B (zh) | 一种卫星导航系统的系统时间偏差的确定方法及系统 | |
CN211180593U (zh) | 时频检测设备 | |
CN112198537B (zh) | 一种基于差分的罗兰高精度定位解算方法 | |
CN103792021A (zh) | Gps授时的石英晶体温度计 | |
CN109635421B (zh) | 一种基于威布尔模型的一般压力表检测周期动态优化方法 | |
CN103698044B (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 |