CN106354008A - Cpiii控制网平差方法 - Google Patents
Cpiii控制网平差方法 Download PDFInfo
- Publication number
- CN106354008A CN106354008A CN201610653683.XA CN201610653683A CN106354008A CN 106354008 A CN106354008 A CN 106354008A CN 201610653683 A CN201610653683 A CN 201610653683A CN 106354008 A CN106354008 A CN 106354008A
- Authority
- CN
- China
- Prior art keywords
- point
- observation
- adjustment
- cpiii
- coordinate
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及CPIII控制网平差方法。约束平差法得到的CPIII控制网不利于对轨道内部几何尺寸的评价与调整,是一种仅侧重于满足轨道外部几何尺寸要求的方法。本发明利用方向观测值和水平距离观测值的误差方程得到二者的改正数阵,通过误差方程的改正数阵、附加条件方程和平差准则获得附加条件阵,即建立平差模型,解算该模型完成平差计算。本发明既能够优先满足轨道内部几何尺寸要求,且能够满足轨道外部几何尺寸的要求,不受原始数据误差的影响,CPIII控制网不会受到扭曲,通过位置基准和方位基准将CPIII控制网纳入线路设计坐标系,使CPIII控制网能够满足高速铁路“三网合一”的要求。
Description
技术领域
本发明涉及高速铁路轨道控制网(平面网)数据处理技术领域,具体涉及一种CPIII控制网平差方法。
背景技术
为了保证列车在高速行驶条件下的安全性和舒适性,高速铁路轨道必须具有非常高的平顺性和精确的几何线形参数,包括轨道内部几何尺寸(轨道自身的几何尺寸)与外部几何尺寸(轨道与周围建筑物的相对尺寸),如轨距、轨向、高低、水平、扭曲、与设计高程及中线的偏差等。其中,内部几何尺寸用于描述轨道的几何形状,主要影响轨道的平顺性;外部几何尺寸用于描述轨道在空间中的位置和标高。如何确保高速铁路轨道内、外部结构的合理性和稳定性是高速铁路建设和维护面临的关键问题之一。
无砟轨道的内、外部几何尺寸的测量称为轨道的相对定位和绝对定位测量。目前,无砟轨道的相对定位通常是利用绝对定位方式来实现的,因此绝对定位精度必须满足相对定位的精度要求,即轨道平顺性的要求。由于CPIII控制网(本文特指CPIII平面网)是轨道绝对定位的控制基准,因此轨道的相对定位对绝对定位的要求就转嫁成了对CPIII控制网相对精度的要求。由此可见,CPIII控制网必须同时满足轨道绝对定位和相对定位的精度要求。由于高速铁路对轨道内部几何尺寸的要求严于对轨道外部几何尺寸的要求,因此作者认为,CPIII控制网应优先满足轨道内部几何尺寸的要求。德国采用“控制网置平”方法得到的PS4(相当于我国的CPIII)便优先满足了轨道内部几何尺寸的要求,同时兼顾了轨道外部几何尺寸的要求。而我国采用约束平差方法得到的CPIII控制网则侧重于满足轨道外部几何尺寸的精度要求,并未优先考虑轨道内部几何尺寸的精度要求。
然而,由于德国对“控制网置平”方法采取了保密措施,国内仅有极少数人声称掌握了其平差原理,导致该算法未能在我国普及。目前,我国仍普遍采用约束平差方法进行CPIII控制网处理。约束平差法的实质是将CPIII控制网强制拉入到线路设计坐标系中,这不仅会造成CPIII控制网的扭曲,还会造成控制网的尺度与现场实际尺度不符。由此可见,约束平差法得到的CPIII控制网不利于对轨道内部几何尺寸的评价与调整,是一种侧重于满足轨道外部几何尺寸要求的方法。因此,有必要研究能够优先满足轨道内部几何尺寸要求(且能够满足轨道外部几何尺寸的要求)的CPIII控制网平差新方法,形成具有我国自主知识产权的CPIII控制网平差新方法,这对高速铁路建设具有重要的应用价值,对进一步提高我国高铁精密工程测量技术水平具有重要意义。
发明内容
本发明的目的是提供一种CPIII控制网平差方法,能够优先满足轨道内部几何尺寸要求,同时又能够满足轨道外部几何尺寸要求,克服约束平差用于CPIII控制网数据处理的缺陷。
本发明所采用的技术方案为:
CPIII控制网平差方法,其特征在于:
包括以下步骤:
第一步:开列误差方程:
参与平差的CPIII控制网观测数据分为方向观测值和水平距离观测值两种;
方向观测值的误差方程为:
式中,
式中,vki是方向观测值的改正数;分别为k、i点的近似坐标值;分别为k、i点的坐标改正数;ski是由点k和i的近似坐标计算的水平距离近似值;Lki为点k至点i的方向观测值;为点k和i的近似坐标计算的近似方位角;为定向角未知数的近似值;dzk为定向角近似值的改正数;ρ″≈206264.8;
水平距离观测值的误差方程为:
式中,
式中,vski是水平距离观测值的改正数;分别为k、i点的近似坐标值;分别为k、i点的坐标改正数;ski为点k和i的水平距离观测值;为由两点近似坐标计算的水平距离近似值;
方向观测值和水平距离观测值的误差方程由以下矩阵形式表示:
式中,v为方向观测值和水平距离观测值的改正数阵;A为误差方程系数阵;l为误差方程的常数项阵;
第二步:建立平差模型:
自由尺度拟稳平差法的平差模型由误差方程、附加条件方程和平差准则构成,如下:
式中,v为观测值残差阵;A为误差方程系数阵;为待估参数阵;l为误差方程常数项阵;GT为附加条件阵;P为观测值权阵;
其中:
式中,δxi、δyi为网点重心化后的坐标值;m1为网中联测的已知点的个数,包括CPI和CPII点;m2为网中待定点的个数,包括CPIII点和全站仪自由测站点;s为全区定向角未知数的个数;
附加条件阵(5)中,第一行与已知点x坐标相对应位置的值为1,其余均为0;第二行与已知点y坐标相对应位置的值为1,其余均为0;第三行与待定点相对应位置的值均为0;
第三步:模型解算:
根据第二步的模型,可得出参数估值为:
dX=(ATPA+GGT)-1ATPL (6)
参数估值的协因数阵为:
QX=(N+GGT)-1N(N+GGT)-1 (7)
式中,N=ATPA;
验后单位权中误差为:
式中,n为总观测值数;t=m1+m2为全网总点数;d=3为必要起算基准数,包括2个位置基准和1个方位基准。
平差后得到已知点的新坐标,新坐标与原坐标之差控制在5mm的范围内,否则剔除该已知点或将该已知点视为待定点参与平差。
本发明具有以下优点:
1、该平差方法不受原始数据误差的影响,保留了CPIII控制网的原始测量精度和原始尺度;
2、该平差方法采用了拟合的思想,未将CPIII控制网强制拉入到高斯投影面中,从而使高精度的CPIII控制网不会受到扭曲,有利于轨道内部几何尺寸的评价与调整;
3、该平差方法提供的平差基准不含尺度基准,其实质是,平差时,以外业采集的水平距离作为尺度基准。而外业采集的水平距离所在高程面与该处的轨道高程面近似相同,故能够近似最佳的满足轨道内部几何尺寸的要求,即优先满了轨道内部几何尺寸的要求;
4、该平差方法通过位置基准和方位基准将CPIII控制网纳入到了线路设计坐标系中,使CPIII控制网能够满足高速铁路“三网合一”的要求,从而也能够满足轨道外部几何尺寸的要求。
具体实施方式
下面结合具体实施方式对本发明进行详细的说明。
本发明涉及的CPIII控制网平差方法,由以下步骤实现:
第一步:开列误差方程
参与平差的CPIII控制网(平面网)观测数据可分为方向观测值和水平距离观测值两种。方向观测值的误差方程为:
式中,
式中,vki是方向观测值的改正数;分别为k、i点的近似坐标值;分别为k、i点的坐标改正数;ski是由点k和i的近似坐标计算的水平距离近似值;Lki为点k至点i的方向观测值;为点k和i的近似坐标计算的近似方位角;为定向角未知数的近似值;dzk为定向角近似值的改正数;ρ″≈206264.8。
水平距离观测值的误差方程为:
式中,
式中,vski是水平距离观测值的改正数;分别为k、i点的近似坐标值;分别为k、i点的坐标改正数;ski为点k和i的水平距离观测值;为由两点近似坐标计算的水平距离近似值。
方向观测值和水平距离观测值的误差方程可由以下矩阵形式表示:
式中,v为方向观测值和水平距离观测值的改正数阵;A为误差方程系数阵;l为误差方程的常数项阵。
第二步:建立平差模型
自由尺度拟稳平差法的平差模型由误差方程、附加条件方程(也称基准方程)和平差准则构成,如下:
式中,v为观测值残差阵;A为误差方程系数阵;为待估参数阵;l为误差方程常数项阵;GT为附加条件阵;P为观测值权阵。
其中
式中,δxi、δyi为网点重心化后的坐标值;m1为网中联测的已知点(包括CPI和CPII点)的个数;m2为网中待定点(包括CPIII点和全站仪自由测站点)的个数;s为全区定向角未知数的个数。上述附加条件矩阵(5)中,第一行与已知点x坐标相对应位置的值为1,其余均为0;第二行与已知点y坐标相对应位置的值为1,其余均为0;第三行与待定点相对应位置的值均为0。
第三步:模型解算
根据上述模型,可得出参数(控制点近似坐标的改正数)估值为:
dX=(ATPA+GGT)-1ATPL (14)
参数估值的协因数阵为:
QX=(N+GGT)-1N(N+GGT)-1 (15)
式中,N=ATPA。
验后单位权中误差为:
式中,n为总观测值数;t=m1+m2为全网总点数;d=3为必要起算基准数(包括2个位置基准和1个方位基准)。
本发明涉及的CPIII控制网平差方法实施时需注意以下几点:
1、不需要对全站仪原始观测数据做两化改正处理;
2、平差后会得到已知点的新坐标,新坐标与原坐标之差建议控制在5mm的范围内,否则应剔除该已知点或将该已知点视为待定点参与平差;
3、当已知点的新坐标与原坐标之差普遍偏大时,应分析已知点的稳定性,不得采用变动的或不兼容的已知点坐标参与平差,否则将难以满足轨道外部几何尺寸的要求。建议将已知点新坐标与原坐标的较差作为控制指标来评价CPIII控制网的绝对精度。
4、由于相邻点的相对中误差与方位基准有关,当CPIII控制网采用该平差方法时,不应采用该指标来评价CPIII控制网的相对精度。建议采用距离中误差和角度中误差来评价CPIII控制网的相对精度。
实例验证及精度分析:
本发明收集国内某高铁一段长约5km的CPIII控制网的平面观测数据,验证该平差方法的可行性。利用该平差方法对进行处理,所得主要精度信息统计情况如下表:
表1该CPIII网主要精度信息统计表
由上表1可知,利用自由尺度拟稳平差方法得到的CPIII控制网总体精度良好,其中距离观测中误差和方向观测中误差两项精度指标均远小于《高速铁路工程测量规范》中的要求;相邻点的相对点位中误差的最大值为0.6mm,满足规范中限差为1mm的要求。
本发明的内容不限于实施例所列举,本领域普通技术人员通过阅读本发明说明书而对本发明技术方案采取的任何等效的变换,均为本发明的权利要求所涵盖。
Claims (2)
1.CPIII控制网平差方法,其特征在于:
包括以下步骤:
第一步:开列误差方程:
参与平差的CPIII控制网观测数据分为方向观测值和水平距离观测值两种;
方向观测值的误差方程为:
式中,
式中,vki是方向观测值的改正数;分别为k、i点的近似坐标值;分别为k、i点的坐标改正数;ski是由点k和i的近似坐标计算的水平距离近似值;Lki为点k至点i的方向观测值;为点k和i的近似坐标计算的近似方位角;为定向角未知数的近似值;dzk为定向角近似值的改正数;ρ″≈206264.8;
水平距离观测值的误差方程为:
式中,
式中,vski是水平距离观测值的改正数;分别为k、i点的近似坐标值;分别为k、i点的坐标改正数;ski为点k和i的水平距离观测值;为由两点近似坐标计算的水平距离近似值;
方向观测值和水平距离观测值的误差方程由以下矩阵形式表示:
式中,v为方向观测值和水平距离观测值的改正数阵;A为误差方程系数阵;l为误差方程的常数项阵;
第二步:建立平差模型:
自由尺度拟稳平差法的平差模型由误差方程、附加条件方程和平差准则构成,如下:
式中,v为观测值残差阵;A为误差方程系数阵;为待估参数阵;l为误差方程常数项阵;GT为附加条件阵;P为观测值权阵;
其中:
式中,δxi、δyi为网点重心化后的坐标值;m1为网中联测的已知点的个数,包括CPI和CPII点;m2为网中待定点的个数,包括CPIII点和全站仪自由测站点;s为全区定向角未知数的个数;
附加条件阵(5)中,第一行与已知点x坐标相对应位置的值为1,其余均为0;第二行与已知点y坐标相对应位置的值为1,其余均为0;第三行与待定点相对应位置的值均为0;
第三步:模型解算:
根据第二步的模型,可得出参数估值为:
dX=(ATPA+GGT)-1ATPL (6)
参数估值的协因数阵为:
QX=(N+GGT)-1N(N+GGT)-1 (7)
式中,N=ATPA;
验后单位权中误差为:
式中,n为总观测值数;t=m1+m2为全网总点数;d=3为必要起算基准数,包括2个位置基准和1个方位基准。
2.根据权利要求1所述的CPIII控制网平差方法,其特征在于:
平差后得到已知点的新坐标,新坐标与原坐标之差控制在5mm的范围内,否则剔除该已知点或将该已知点视为待定点参与平差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610653683.XA CN106354008B (zh) | 2016-08-10 | 2016-08-10 | Cpiii控制网平差方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610653683.XA CN106354008B (zh) | 2016-08-10 | 2016-08-10 | Cpiii控制网平差方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106354008A true CN106354008A (zh) | 2017-01-25 |
CN106354008B CN106354008B (zh) | 2019-05-03 |
Family
ID=57845026
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610653683.XA Active CN106354008B (zh) | 2016-08-10 | 2016-08-10 | Cpiii控制网平差方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106354008B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106980598A (zh) * | 2017-04-11 | 2017-07-25 | 中铁第勘察设计院集团有限公司 | 自由测站边角交会网的三维平差严密方法 |
CN108563129A (zh) * | 2018-05-30 | 2018-09-21 | 江苏师范大学 | 高速铁路cpiii高程控制网条件方程的快速建立方法 |
CN109916430A (zh) * | 2019-02-28 | 2019-06-21 | 中铁二院工程集团有限责任公司 | 一种自由测站边角交会网测站观测数据质量检查与评价方法 |
CN111473737A (zh) * | 2020-04-13 | 2020-07-31 | 中国三峡建设管理有限公司 | 一种水电站地下工程自由设站精密控制网的测设方法 |
CN111891176A (zh) * | 2020-06-23 | 2020-11-06 | 中铁第一勘察设计院集团有限公司 | 基于轨道几何状态测量仪的轨道测量数据整体平差方法 |
CN112906198A (zh) * | 2021-01-21 | 2021-06-04 | 西南交通大学 | 任意平面测量控制网非线性最小二乘平差数据处理方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1120690A (ja) * | 1997-07-04 | 1999-01-26 | Railway Technical Res Inst | 鉄道車両の車体傾斜制御方法 |
CN101914881A (zh) * | 2010-07-27 | 2010-12-15 | 唐粮 | 一种高速铁路基桩控制网cpiii的快速测量方法 |
CN102168396A (zh) * | 2011-03-18 | 2011-08-31 | 中铁第一勘察设计院集团有限公司 | 轨道基准网数据采集与数据处理外业实时一体化测量方法 |
CN102425085A (zh) * | 2011-10-20 | 2012-04-25 | 中铁第一勘察设计院集团有限公司 | 高速铁路运营维护阶段cpiii平面网复测方法 |
CN103759714A (zh) * | 2014-01-20 | 2014-04-30 | 国家测绘地理信息局卫星测绘应用中心 | 一种三线阵卫星影像区域网平差方法 |
-
2016
- 2016-08-10 CN CN201610653683.XA patent/CN106354008B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1120690A (ja) * | 1997-07-04 | 1999-01-26 | Railway Technical Res Inst | 鉄道車両の車体傾斜制御方法 |
CN101914881A (zh) * | 2010-07-27 | 2010-12-15 | 唐粮 | 一种高速铁路基桩控制网cpiii的快速测量方法 |
CN102168396A (zh) * | 2011-03-18 | 2011-08-31 | 中铁第一勘察设计院集团有限公司 | 轨道基准网数据采集与数据处理外业实时一体化测量方法 |
CN102425085A (zh) * | 2011-10-20 | 2012-04-25 | 中铁第一勘察设计院集团有限公司 | 高速铁路运营维护阶段cpiii平面网复测方法 |
CN103759714A (zh) * | 2014-01-20 | 2014-04-30 | 国家测绘地理信息局卫星测绘应用中心 | 一种三线阵卫星影像区域网平差方法 |
Non-Patent Citations (3)
Title |
---|
任晓春,周东卫: "高速铁路运营阶段CPIII复测方法研究与应用", 《铁道工程学报》 * |
张玉柱,陈龙: "CPIII轨道控制网平差", 《铁路工程》 * |
段泽洋: "顾及CPIII点精度的严密平差方法研究", 《测绘技术装备》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106980598A (zh) * | 2017-04-11 | 2017-07-25 | 中铁第勘察设计院集团有限公司 | 自由测站边角交会网的三维平差严密方法 |
CN106980598B (zh) * | 2017-04-11 | 2019-05-03 | 中铁第一勘察设计院集团有限公司 | 自由测站边角交会网的三维平差严密方法 |
CN108563129A (zh) * | 2018-05-30 | 2018-09-21 | 江苏师范大学 | 高速铁路cpiii高程控制网条件方程的快速建立方法 |
CN109916430A (zh) * | 2019-02-28 | 2019-06-21 | 中铁二院工程集团有限责任公司 | 一种自由测站边角交会网测站观测数据质量检查与评价方法 |
CN109916430B (zh) * | 2019-02-28 | 2023-01-13 | 中铁二院工程集团有限责任公司 | 一种自由测站边角交会网测站观测数据质量检查与评价方法 |
CN111473737A (zh) * | 2020-04-13 | 2020-07-31 | 中国三峡建设管理有限公司 | 一种水电站地下工程自由设站精密控制网的测设方法 |
CN111891176A (zh) * | 2020-06-23 | 2020-11-06 | 中铁第一勘察设计院集团有限公司 | 基于轨道几何状态测量仪的轨道测量数据整体平差方法 |
CN112906198A (zh) * | 2021-01-21 | 2021-06-04 | 西南交通大学 | 任意平面测量控制网非线性最小二乘平差数据处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106354008B (zh) | 2019-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106354008A (zh) | Cpiii控制网平差方法 | |
CN104236546B (zh) | 一种卫星星光折射导航误差确定与补偿方法 | |
CN104061932B (zh) | 一种利用引力矢量和梯度张量进行导航定位的方法 | |
CN108562288A (zh) | 一种激光捷联惯组系统级在线自标定系统及方法 | |
CN104680008B (zh) | 一种基于多参考站的网络rtk区域大气误差建模方法 | |
CN104359492B (zh) | 惯性导航和轮速计组成的航迹推算定位系统误差估算算法 | |
CN107192381A (zh) | 一种任意无序空间结构的测量放样方法 | |
CN103994765B (zh) | 一种惯性传感器的定位方法 | |
CN108375985A (zh) | 一种土地三维规划设计平台及其设计方法 | |
CN108307301A (zh) | 基于rssi测距和轨迹相似性的室内定位方法 | |
CN110793508A (zh) | 三维控制网中控制点的观测数据处理方法与装置 | |
CN104613966B (zh) | 一种地籍测量事后数据处理方法 | |
CN108490461A (zh) | 一种卫星导航信号采集回放仪定位精度一致性的测试方法 | |
CN107300377A (zh) | 一种绕飞轨迹下的旋翼无人机三维目标定位方法 | |
CN107066747A (zh) | 一种视觉测量网络组网规划方法 | |
CN107179533A (zh) | 一种多参数机载LiDAR系统误差自动检校方法 | |
CN104931071B (zh) | 一种基于选权迭代的星敏感器在轨几何定标方法及系统 | |
CN104121927A (zh) | 一种适用于低精度无方位基准单轴转位设备的惯性测量单元标定方法 | |
CN103646136A (zh) | 基于gis的飞行校验方法 | |
CN109855623A (zh) | 基于Legendre多项式和BP神经网络的地磁模型在线逼近方法 | |
CN108508463A (zh) | 基于Fourier-Hermite正交多项式扩展椭球集员滤波方法 | |
CN110006448A (zh) | 星载激光测高仪指向角系统误差在轨标定精度的评估方法 | |
CN110779503B (zh) | 一种三维精密控制网测量方法 | |
CN109856619B (zh) | 一种雷达测向相对系统误差修正方法 | |
CN109282802B (zh) | 一种基于单盘位测量的自由设站方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |