发明名称 CPIII控制网平差方法
摘要 本发明涉及CPIII控制网平差方法。约束平差法得到的CPIII控制网不利于对轨道内部几何尺寸的评价与调整,是一种仅侧重于满足轨道外部几何尺寸要求的方法。本发明利用方向观测值和水平距离观测值的误差方程得到二者的改正数阵,通过误差方程的改正数阵、附加条件方程和平差准则获得附加条件阵,即建立平差模型,解算该模型完成平差计算。本发明既能够优先满足轨道内部几何尺寸要求,且能够满足轨道外部几何尺寸的要求,不受原始数据误差的影响,CPIII控制网不会受到扭曲,通过位置基准和方位基准将CPIII控制网纳入线路设计坐标系,使CPIII控制网能够满足高速铁路“三网合一”的要求。
申请公布号 CN106354008A 申请公布日期 2017.01.25
申请号 CN201610653683.X 申请日期 2016.08.10
申请人 中铁第一勘察设计院集团有限公司 发明人 任晓春;马洪磊;罗文彬
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 西安新思维专利商标事务所有限公司 61114 代理人 李罡
主权项 CPIII控制网平差方法,其特征在于:包括以下步骤:第一步:开列误差方程:参与平差的CPIII控制网观测数据分为方向观测值和水平距离观测值两种;方向观测值的误差方程为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>v</mi><mrow><mi>k</mi><mi>i</mi></mrow></msub><mo>=</mo><mo>-</mo><msub><mi>dz</mi><mi>k</mi></msub><mo>+</mo><msub><mi>a</mi><mrow><mi>k</mi><mi>i</mi></mrow></msub><msub><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi></msub><mo>+</mo><msub><mi>b</mi><mrow><mi>k</mi><mi>i</mi></mrow></msub><msub><mover><mi>y</mi><mo>^</mo></mover><mi>k</mi></msub><mo>+</mo><msub><mi>a</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mo>+</mo><msub><mi>b</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><msub><mover><mi>y</mi><mo>^</mo></mover><mi>i</mi></msub><mo>-</mo><msub><mi>l</mi><mrow><mi>k</mi><mi>i</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001074903420000011.GIF" wi="1139" he="63" /></maths>式中,<img file="FDA0001074903420000012.GIF" wi="1334" he="135" /><img file="DA00010749034235741.GIF" wi="401" he="142" /><img file="FDA0001074903420000013.GIF" wi="439" he="68" />式中,v<sub>ki</sub>是方向观测值的改正数;<img file="FDA0001074903420000014.GIF" wi="382" he="63" />分别为k、i点的近似坐标值;<img file="FDA0001074903420000015.GIF" wi="374" he="63" />分别为k、i点的坐标改正数;s<sub>ki</sub>是由点k和i的近似坐标计算的水平距离近似值;L<sub>ki</sub>为点k至点i的方向观测值;<img file="FDA0001074903420000016.GIF" wi="54" he="63" />为点k和i的近似坐标计算的近似方位角;<img file="FDA0001074903420000017.GIF" wi="46" he="63" />为定向角未知数的近似值;dz<sub>k</sub>为定向角近似值的改正数;ρ″≈206264.8;水平距离观测值的误差方程为:<maths num="0002"><math><![CDATA[<mrow><msub><mi>v</mi><mrow><mi>s</mi><mi>k</mi><mi>i</mi></mrow></msub><mo>=</mo><msub><mi>c</mi><mrow><mi>k</mi><mi>i</mi></mrow></msub><msub><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi></msub><mo>+</mo><msub><mi>d</mi><mrow><mi>k</mi><mi>i</mi></mrow></msub><msub><mover><mi>y</mi><mo>^</mo></mover><mi>k</mi></msub><mo>+</mo><msub><mi>c</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mo>+</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><msub><mover><mi>y</mi><mo>^</mo></mover><mi>i</mi></msub><mo>-</mo><msub><mi>l</mi><mrow><mi>s</mi><mi>k</mi><mi>i</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001074903420000018.GIF" wi="1027" he="63" /></maths>式中,<img file="FDA0001074903420000019.GIF" wi="1374" he="142" /><img file="FDA00010749034200000110.GIF" wi="310" he="63" />式中,v<sub>ski</sub>是水平距离观测值的改正数;<img file="FDA00010749034200000111.GIF" wi="373" he="71" />分别为k、i点的近似坐标值;<img file="FDA0001074903420000021.GIF" wi="358" he="63" />分别为k、i点的坐标改正数;s<sub>ki</sub>为点k和i的水平距离观测值;<img file="FDA0001074903420000026.GIF" wi="54" he="63" />为由两点近似坐标计算的水平距离近似值;方向观测值和水平距离观测值的误差方程由以下矩阵形式表示:<maths num="0003"><math><![CDATA[<mrow><mi>v</mi><mo>=</mo><mi>A</mi><mover><mi>x</mi><mo>^</mo></mover><mo>-</mo><mi>l</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001074903420000022.GIF" wi="510" he="61" /></maths>式中,v为方向观测值和水平距离观测值的改正数阵;A为误差方程系数阵;l为误差方程的常数项阵;第二步:建立平差模型:自由尺度拟稳平差法的平差模型由误差方程、附加条件方程和平差准则构成,如下:<maths num="0004"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>v</mi><mo>=</mo><mi>A</mi><mover><mi>x</mi><mo>^</mo></mover><mo>-</mo><mi>l</mi></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>G</mi><mi>T</mi></msup><mover><mi>x</mi><mo>^</mo></mover><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>v</mi><mi>T</mi></msup><mi>P</mi><mi>v</mi><mo>=</mo><mi>m</mi><mi>i</mi><mi>n</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001074903420000023.GIF" wi="422" he="220" /></maths>式中,v为观测值残差阵;A为误差方程系数阵;<img file="FDA0001074903420000025.GIF" wi="35" he="47" />为待估参数阵;l为误差方程常数项阵;G<sup>T</sup>为附加条件阵;P为观测值权阵;其中:<img file="FDA0001074903420000024.GIF" wi="1678" he="359" />式中,δx<sub>i</sub>、δy<sub>i</sub>为网点重心化后的坐标值;m<sub>1</sub>为网中联测的已知点的个数,包括CPI和CPII点;m<sub>2</sub>为网中待定点的个数,包括CPIII点和全站仪自由测站点;s为全区定向角未知数的个数;附加条件阵(5)中,第一行与已知点x坐标相对应位置的值为1,其余均为0;第二行与已知点y坐标相对应位置的值为1,其余均为0;第三行与待定点相对应位置的值均为0;第三步:模型解算:根据第二步的模型,可得出参数估值为: dX=(A<sup>T</sup>PA+GG<sup>T</sup>)<sup>‑1</sup>A<sup>T</sup>PL  (6)参数估值的协因数阵为: Q<sub>X</sub>=(N+GG<sup>T</sup>)<sup>‑1</sup>N(N+GG<sup>T</sup>)<sup>‑1</sup>  (7)式中,N=A<sup>T</sup>PA;验后单位权中误差为:<maths num="0005"><math><![CDATA[<mrow><msub><mover><mi>&delta;</mi><mo>^</mo></mover><mn>0</mn></msub><mo>=</mo><msqrt><mfrac><mrow><msup><mi>v</mi><mi>T</mi></msup><mi>P</mi><mi>v</mi></mrow><mrow><mi>n</mi><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>-</mo><mi>d</mi><mo>)</mo></mrow><mo>-</mo><mi>s</mi></mrow></mfrac></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001074903420000031.GIF" wi="598" he="158" /></maths>式中,n为总观测值数;t=m<sub>1</sub>+m<sub>2</sub>为全网总点数;d=3为必要起算基准数,包括2个位置基准和1个方位基准。
地址 710043 陕西省西安市西影路二号