发明名称 一种基于温度补偿的系统时延误差校正方法
摘要 一种基于温度补偿的系统时延误差校正方法,在精密测距产品内部设置测温点并将温度遥测采集,通过对产品进行温度-伪距的测量得到先验值,建立n阶多项式回归模型进行拟合残差分析,当随着阶数的增加残差精度压缩或小于预设阀值,即可选取压缩点处/小于预设阀值作为该多项式回归模型的阶数。根据最小均方根准则估计n阶多项式回归模型时延校正系数,根据时延校正系数和实时采集的温度遥测进行精密测距产品温变下的时延校正,本发明通过在产品硬件上增加温度遥测采集作为校正输入量有效的解决了精密测距产品由于温度变化引起的测距系统误差,此外,本发明中的校正方法简单可靠,整体资源占用量较少。
申请公布号 CN105093239A 申请公布日期 2015.11.25
申请号 CN201510519287.3 申请日期 2015.08.21
申请人 西安空间无线电技术研究所 发明人 周晓平;蒙艳松;徐连军;郑先安;陈旭阳;朱向鹏;张中英;
分类号 G01S19/23(2010.01)I 主分类号 G01S19/23(2010.01)I
代理机构 中国航天科技专利中心 11009 代理人 臧春喜
主权项 一种基于温度补偿的系统时延误差校正方法,其特征在于步骤如下:(1)对待校正产品进行单向伪距值测量R<sub>i</sub>,i=1,2,...,M,并选取待校正产品的一处温度测量点采集该产品的遥测温度t<sub>i</sub>,i=1,2,...,M;(2)建立遥测温度t<sub>i</sub>与拟合得到的伪距值的n阶多项式回归模型,具体由公式:F<sub>i</sub>=a<sub>n</sub>(t<sub>i</sub>‑t<sub>0</sub>)<sup>n</sup>+a<sub>n‑1</sub>(t<sub>i</sub>‑t<sub>0</sub>)<sup>n‑1</sup>...+a<sub>2</sub>(t<sub>i</sub>‑t<sub>0</sub>)<sup>2</sup>+a<sub>1</sub>(t<sub>i</sub>‑t<sub>0</sub>)<sup>1</sup>+a<sub>0</sub>给出,式中,a<sub>n</sub>,a<sub>n‑1</sub>,…,a<sub>0</sub>为n阶多项式回归模型的拟合系数,t<sub>i</sub>为温度遥测量,t<sub>0</sub>为预先给定的基准温度,F<sub>i</sub>为拟合得到的伪距值;(3)从预先给定的n的取值范围中选取n的最小值,计算n阶多项式回归模型的拟合系数a<sub>n</sub>,a<sub>n‑1</sub>,…,a<sub>0</sub>,具体由公式:a=A/b给出,其中,A为(n+1)×(n+1)的矩阵,由公式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mi>n</mi><mo>+</mo><mi>n</mi></mrow></msup></mrow></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mi>n</mi><mo>+</mo><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msup></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msup></mrow></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mi>n</mi><mo>+</mo><mn>0</mn></mrow></msup></mrow></mtd></mtr><mtr><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn><mo>+</mo><mi>n</mi></mrow></msup></mrow></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn><mo>+</mo><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msup></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn><mo>+</mo><mn>1</mn></mrow></msup></mrow></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn><mo>+</mo><mn>0</mn></mrow></msup></mrow></mtd></mtr><mtr><mtd><mrow><mo>...</mo><mo>...</mo></mrow></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mn>1</mn><mo>+</mo><mi>n</mi></mrow></msup></mrow></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mn>1</mn><mo>+</mo><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msup></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mn>1</mn><mo>+</mo><mn>1</mn></mrow></msup></mrow></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mn>1</mn><mo>+</mo><mn>0</mn></mrow></msup></mrow></mtd></mtr><mtr><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mn>0</mn><mo>+</mo><mi>n</mi></mrow></msup></mrow></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mn>0</mn><mo>+</mo><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msup></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mn>0</mn><mo>+</mo><mn>1</mn></mrow></msup></mrow></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mn>0</mn><mo>+</mo><mn>0</mn></mrow></msup></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000786603830000011.GIF" wi="1482" he="666" /></maths>给出;a为(n+1)×1的矩阵,由公式:a=[a<sub>n</sub>,a<sub>n‑1</sub>,...,a<sub>1</sub>,a<sub>0</sub>]<sup>T</sup>给出;b为(n+1)×1的矩阵,由公式:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>b</mi><mo>=</mo><mo>&lsqb;</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mrow><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>n</mi></msup><msub><mi>R</mi><mi>i</mi></msub></mrow><mo>,</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mrow><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>R</mi><mi>i</mi></msub></mrow><mo>,</mo><mo>...</mo><mo>,</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mrow><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>1</mn></msup><msub><mi>R</mi><mi>i</mi></msub></mrow><mo>,</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mrow><msup><mrow><mo>(</mo><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>0</mn></msup><msub><mi>R</mi><mi>i</mi></msub></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000786603830000012.GIF" wi="1308" he="150" /></maths>给出,其中R<sub>i</sub>为温度遥测量t<sub>i</sub>时对应的原始测距值;(4)利用步骤(3)中确定的n阶多项式回归模型的拟合系数a<sub>n</sub>,a<sub>n‑1</sub>,…,a<sub>0</sub>,计算n阶多项式回归模型的拟合残差,n阶多项式回归模型的拟合残差平方和Sr;(5)利用步骤(4)中计算得到的n阶多项式回归模型的拟合残差,计算步骤(3)中确定的n阶多项式回归模型的拟合系数a<sub>n</sub>,a<sub>n‑1</sub>,…,a<sub>0</sub>的精度δ;(6)若步骤(5)中计算的精度δ小于等于预设的精度要求,则进入步骤(7),否则,将n的取值加1,返回步骤(3),直到选取满足预设精度要求的拟合阶数n;(7)计算遥测温度t<sub>i</sub>对应的原始测距值补偿量R<sub>bi</sub>;(8)利用步骤(7)中的原始测距值补偿量R<sub>bi</sub>对原始测距值进行补偿,具体由公式:R<sub>Fi</sub>=R<sub>bi</sub>+R<sub>i</sub>,i=1,2,...,M给出。
地址 710100 陕西省西安市长安区西街150号