主权项 |
一种基于温度补偿的系统时延误差校正方法,其特征在于步骤如下:(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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>Σ</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>[</mo><munderover><mo>Σ</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>Σ</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>Σ</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>Σ</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>]</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给出。 |