线路测量坐标正、反算原理及卡西欧5800计算器程序说明_secret.doc

线路测量坐标正、反算计算原理 及卡西欧fx-5800P程序说明 一、计算原理 在各测量书中对于坐标正算的相关计算式均有说明，故在此不做详解，仅对正算过程中需要用到的原理及公式做一汇总。对于坐标反算，虽然都采用无限趋近原理进行计算，但计算方式各有差别，本文仅对其中一种自认为相对简单易懂并便于操作的原理进行详解。 1.1 坐标转换[1] 如图1，设XP、YP为P点在国家控制网坐标系中的坐标；xP、yP为P点在工程独立控制网坐标系中的坐标。XO、YO为工程独立坐标系原点o在国家坐标系中的坐标，Δα为两坐标系纵坐标轴的交角。如果一条边在国家坐标系中的坐标方位角为A，而在工程独立坐标系中的坐标方位角为α，则：  (1-1) 当由工程独立坐标系坐标换算至国家坐标系坐标时，换算公式为：  (1-2) 当由国家坐标系坐标换算至工程独立坐标系坐标时，也可使用式(1-2)，此时应将X、Y与x、y互换，且。 1.2 坐标方位角关系计算 1.2.1 正、反坐标方位角[2] 一条直线的坐标方位角与直线的前进方向有关，沿直线前进方向的坐标方位角称为正坐标方位角，与其相反方向的坐标方位角称为反坐标方位角。如图2，由于轴子午线之间是互相平行的，因此同一直线的正、反坐标方位角相差180°，即：  (2-1) 当时，取“+”号；当时，取“-”号。 1.2.2 坐标方位角的推算[3] 1.2.2.1 转折角为右角 如图3(a)，为已知边坐标方位角，为推算边的坐标方位角，为该两边所夹的右角，则：  (2-2) 1.2.2.2 转折角为左角 如图3(b)，为已知边坐标方位角，为推算边的坐标方位角，  (2-3) 无论用右角还是左角推算，如遇出现负数的情形，应加上。当推算出的坐标方位角超过，则应减去。 如果已知导线起始边的坐标方位角和导线各右角或左角时，可直接用以下公式计算第n边的坐标方位角： 按右角： (2-4) 按左角： (2-5) 1.3基本曲线要素及坐标计算 各类测量书均对曲线要素及坐标计算有详细推到过程，此处仅列出计算公式，便于程序编制说明。本编制程序适用于圆曲线两端缓和曲线不等长的基本曲线，即基本形式为：第一直线段→第一缓和曲线段→圆曲线段→第二缓和曲线段→第二直线段。曲线形式及各要素如图4所示，由于圆曲线两端的缓和曲线不等长，故曲线中点与圆曲线的中点不一致，一般取曲线交点与圆心的连线与圆曲线相交之点M作为曲线中点。 1.3.1 曲线要素计算[4] 切线角： (3-1) 内移值： (3-2) 切线增值： (3-3) 切线长： (3-4) 圆曲线长： (3-5) 曲线长： (3-6) 或者： (3-7)  (3-8)  (3-9) 1.3.2 切线支距坐标计算[5] 缓和曲线： (3-10) 式中：为计算点至直缓点ZH或缓直点HZ的曲线长，为缓和曲线长度。圆 曲 线： (3-11) 式中：为计算点至缓圆点HY或YH点的曲线长，为该点圆曲线部分的长度，即该点里程与HY或YH点里程之差；为弧长对应圆心角。 1.4 中桩坐标计算[6] 参见图5，分段对各曲线组成段进行坐标计算。 1.4.1 HZ点(包括路线起点)至ZH点之间的中桩坐标计算 此段为直线段，即第一直线段，计算点坐标按下式计算：  (4-1) 式中：为路线导线至的坐标方位角；为计算点至点的距离，及计算点里程与点里程之差。 实际计算过程中，通过换算的反方位角即至的坐标方位角来计算第一直线段中桩坐标。计算式如下：  (4-2) 式中：为计算点至的距离，及计算点里程与里程之差。 1.4.2 ZH点至YH点之间的中桩坐标计算 此段包括第一缓和曲线及圆曲线，可按式(3-10)和式(3-11)先分别计算出切线支距坐标x、y，然后通过坐标变换将其转换为测量坐标X、Y。坐标变换公式为：  (4-3) 当曲线为左转角时，应以代入；点的坐标采用(4-2)计算，在此不做计算说明。 1.4.3 YH点至HZ点之间的中桩坐标计算 此段为第二缓和曲线，仍按(3-10)计算切线支距坐标，再按下式转换为测量坐标：  (4-4) 当曲线为右转角时，应以代入。 实际编程计算过程中，应换算成与转角的关系式，即：  (4-5) 当为右偏角时取“+”号，当为左偏角时取“-”号。 1.5 中桩切线方位角计算[7] 如图5，第一直线段切线方位角即，第二直线段切线方位角即。第一缓和曲线段、圆曲线段及第二缓和曲线段切线方位角计算分别如下： 1.5.1 第一缓和曲线段切线方位角计算  (5-1) 式中：右转角时，左转角时；，其中为计算点至点的里程差。 1.5.2 圆曲线段切线方位角计算  (5-2) 式中：同前；按式(3-1)计算。 1.5.3 第二缓和曲线段切线方位