一、QXFY 辛甫森公式放样程序1. “X0”? U：“Y0”?V2. “XA”? A：“Y A”? B：“CA”? C：“1÷RA”?D：“1÷RB”?E：“KA”?F：“KB”? G3. Lb1 1：“KI”?H：“JJ”?L：“Y+Z－”?R4. If H>G Or H<F : Then Goto 1: IfEnd6. (E-D) ÷Abs(G-F) P: Abs(H –F)→Q：P×Q→Z：D+Z→T7. C+（Z+2D）Q×90÷π→W :“QXFWJ”：8. C+(Z÷4+2D)Q×22.5÷π→M：C+（3Z÷4+2D）Q×67.5÷π→N：C+(Z÷2+2D)Q×45÷π→K9. A+Q（cosC+4(cosM+cosN)+2cosK+cosW）÷12+Rcos(W+L) →X：“X＝”：10. B+Q(sin(C)+4(sin(M)+sin(N))+2sin(K)+sin(W))÷12+Rsin (W+L)→Y：“Y＝”：Y11. Pol((X-U)，（Y-V）)If J≤0 Then J+360→J: IfEnd : “FWJ”:J12. “JU”: I13. Goto 1程序显示说明：须输入参数：X0:置镜点X (对应变量U)Y0:置镜点Y (对应变量V)XA:曲线起算点X (对应变量A)Y A:曲线起算点Y (对应变量B)CA：曲线起算点切线方位角(对应变量C) 1÷RA：1÷半径，即起算点曲率，右偏为正左偏为负(对应变量D)1÷RB：1÷半径，即曲终点曲率，右偏为正左偏为负(对应变量E)KA:起算点里程(对应变量F)KB:曲终点里程(对应变量G)KI:待求点里程(对应变量H)JJ:夹角(与前进方向切线向右的夹角)(对应变量L)D“Y + Z－：偏距，右正左负(对应变量R)计算结果: W: 待求点切线方位角(对应变量W) X:计算点x (对应变量X)Y:计算点y (对应变量Y)FWJ:放样方位角(对应变量J)S:放样距离(对应变量I)置X0Y0KB本程序依据复化辛甫森公式计算原理改进而成，特点是把曲线按曲率变化点分成若干计算单元单独计算，不论直线（曲率为0）、圆曲线（曲率为1/R）、卵形曲线，只要按曲率变化点分解弄清变化点曲率半径，右偏左偏曲线由曲率的正负号分别，既可求得该计算单元内任意里程中桩也可求斜交或正交的边桩。

考虑到程序中调用变量较多与其它程序变量冲突，造成不方便来回交叉使用，便把反算程序嵌套其中，作为一个完整独立放样和坐标计算程序，不含子程序当输入待求里程超出该曲线单元范围时程序会拒绝运算，并提示重新输入里程，避免计算错误结果；如要求另一曲线单元内里程须重新进入程序输入相应曲线元参数计算由于应用中待求里程不一定按从小到大顺序便去掉了把计算里程替代为起算里程的语句，造成待求里程与起算里程距离较远时精度会略有影响，但对于一般工程中较大半径曲线，影响非常小仅在毫米上可以忽略不予考虑。

二、坐标反算ZBFS1. “X0”？U：V“Y0”？V2. Lbl 0：“X”？X：“Y”？Y：Pol((X－U),(Y－V))3. If J≤0：Then J+360→J：IfEnd ：“FWJ”：J4. “JU”：I5. Goto 0程序显示说明：须输入参数：X0？: 起算点XY0？: 起算点YX ？: 计算点XY？: 计算点Y计算结果: FWJ: 反算方位角JU: 反算距离X ？Y？: 重新输入计算点X Y三、坐标正算ZBZS1. “X0”?X：“Y0”? Y：“S”? S：“T”? T2. “X1”:X+Scos(T)→M“Y1”:Y+Ssin (T)→N程序显示说明：须输入参数：X0？: 起算点XY0？: 起算点YS？: 正算距离T？: 方位角计算结果: X1: 计算点XY1: 计算点Y四、竖曲线SQX1. “R＝”? A：“T＝”? B：“I1＝”? C：“I2＝”? D：“K0＝”?E：“H0＝”? F2. If C>D : Then －1→T : Else 1→T :IfEnd3. Lbl 0：“K＝”?K：Abs(E－K）→S：4. If S>B :Then Goto 0 : IfEnd5. If E≥K : Then Goto 1: Else Goto 2 :IfEnd6. Lbl 1：F－SC+T（B－S）2÷2÷A→G:“G”: G7. Goto 08. Lbl 2：F+SD+T（B－S）2÷2÷A→G:“G”: G9. Goto 0程序显示说明：须输入参数：R＝？: 竖曲线半径T＝？: 切线长I1＝？: 变坡点前坡度值，上坡为正，下坡为负（如-2.5%则输-0.025）I2＝？: 变坡点后坡度值，上坡为正，下坡为负K0＝？：变坡点里程H0＝？：变坡点坡线高程K＝？：待求点里程计算结果: G: 待求点曲线高程K＝？：重新输入待求点里程注：计算里程仅在竖曲线范围内凹型竖曲线凸型竖曲线五、直线段坐标换算ZXZH (X,Y)→(A,B) 1. “X0”?M ：“Y0”? N ：“K ”? Z ：“T ”? T2. Lbl 0 ：“X ”?X ：“Y ”?Y3． (X －M) cos(Ｔ)＋(Y －N)sin(Ｔ)＋Ｚ→Ａ: “LC ”: A4． (Y －N)cos(Ｔ)－(X －M)sin(Ｔ)→Ｂ: “KD ”: B 5． Goto 0程序显示说明：须输入参数：X0? Y0?－直线起算点坐标（Ｘ，Ｙ） (对应变量M N)K?－起算点桩号 (对应变量Z) Ｔ?－起算点至ＪＤ的方位角 (对应变量T) Ｘ? Ｙ?－待求点坐标 (对应变量X Y)计算结果: LC －待求点对应桩号 (对应变量A )KD －距离中线距离，左负，右正 (对应变量B )Ｘ? Ｙ? －重新输入待求坐标点α切Y X 圆心O X Y α切六、圆曲线坐标换算YHS (X,Y)→(A,B)1．“X0”?M ：“Y0”? N ：“KQZ ”? Z ：“T ”? T:“R ”? R : “K ”? K 2. Lbl 0 :“X ”?X ：“Y ”? Y3．Ｐol(Ｘ－Ｍ,Ｙ－Ｎ） 注：计算(Ｘ,Ｙ)点至圆心的距离 4. If J ≤0：Then J+360→J ：IfEnd 5. Ｚ＋ＫπＲ(Ｊ－Ｔ)÷180→Ａ:“LC ”: A 注：计算(Ｘ,Ｙ)点的桩号 6. Ｋ（Ｒ－I ）→Ｂ: “KD ”: B 注：计算(Ｘ,Ｙ)点的距中线距离 7．Goto 0程序显示说明：须输入参数：X0? Y0?－圆心坐标（Ｘ,Ｙ） (对应变量M N)KQZ?－曲线中点（QZ ）的桩号 (对应变量Z) Ｔ?－圆心至QZ 点的方位角 (对应变量T) Ｒ?－圆曲线半径 (对应变量R) Ｋ?－根据ＪＤ点处的偏角α确定，左偏取－１，右偏取＋１(对应变量K)Ｘ? Ｙ?－待求点坐标 (对应变量X Y)计算结果： LC －待求点对应桩号(对应变量A)KD －距离中线距离，左负，右正 (对应变量B ) Ｘ? Ｙ? －重新输入待求坐标点α切Y X 圆心O X Y α切七、缓和曲线段坐标换算HHS (X,Y)→(A,B)1．“XQ”?M：“YQ”? N：“Z”? Z：“S”? S:“R”? R: “L0”? L:“T”? T:“K”? K2. Lbl 0 ：“X”?X：“Y”? Y3．Z →A ：M→E ：N→F ：T→O4．Lbl 2：Ｐol（X－E，Y－F）5．If J≤0：Then J+360→J：IfEnd6．U=J－O7 . A+S I cosU→A8．If Abs(I cos (U)）≤0.001：Then A→A: “LC”: A9. S I sin (U)→B : “KD”: B10. Goto 0 : IfEnd11．RL→C：Abs(A－Z）→D12．If D＞L：Then“D＞L”13. Goto 0 : IfEnd14．D－D∧5÷（40C 2）+D∧9÷（3456C∧4）→V15．D∧3÷（6C）－D∧7÷（336C∧3）+D∧11÷（42240C∧5）→W 16．Ｐol（V，W）17．T+K J→P：T+90 K D2÷(Cπ)→O18. M+I cos( P)→E：N+I Sin(P)→F19. Goto 2程序显示说明须输入参数：X0? Y0?－为ZH（或HZ）点坐标（X，Y）(对应变量M N) Z?－为ZH（或HZ）点的桩号(对应变量Z)S?－由ZH点计算时, S取+1，否则-1 (对应变量S)R?－圆曲线半径(对应变量R)L0?－缓和曲线长度（L0）(对应变量L)T?－为ZH（或HZ）点至本曲线JD点的方位角(对应变量T)K?－根据JD点处的偏角α确定，左偏取-1，右偏取+1；当由HZ点计算时，左偏取+1，右偏取-1 (对应变量K)Ｘ? Ｙ?－待求点坐标(对应变量X Y)运算结果: LC－待求点对应桩号(对应变量A)KD－距离中线距离，左负，右正(对应变量B)Ｘ? Ｙ? －重新输入待求坐标点8-1边长交汇，在两点不通视或不方便架仪器时可任意架设仪器，依次照准两已知点并观测记录边长（夹角不用测）由本程序推算出测站坐标及后视方位角通过配盘完成设站，在外业观测中提高测量效率，特别对于没有程序功能的全站仪方便实用注意：运算结果为P在AB点右侧点，观测务必由A到B顺时针观测，不然结果为镜像点位置为AB连线左侧点，应用中必须判断明确！！！1. “XA”？U: “Y A” ?V: “XB”？X: “YB”？Y2. Pol((U－X),(V－Y))J：“JB-A” JIf J≤0：Then J+360→J：IfEnd3.4. Lbl 1 : “SA”？E: “SB”？F5. cos－1((I2+F2－E2) ÷(2IF)) →B：J－B→A: “XP”：7. X+FcosA→C8. Y+FsinA →D:“YP”：D9. G=A－180：“JPB”：If G≤0 ：Then G+360→G：IfEnd9. Goto 1程序显示依次说明须输入参数：XA？YA？－输入观测第一点坐标对应变量U V XB？YB？－输入观测第二点坐标对应变量X YJB-A 计算得出第二点指向第一点方位角SBA－计算得出第二点到第一点距离（以上两项可视需要可不显示）SA？－输入由未知点P点到第一点A观测距离对应变量ESB？－输入由未知点P点到第二点B观测距离对应变量F 运算结果: XP YP－解得未知P点坐标对应变量C D JPB－由P点后视第二已知点推算方位角，（配盘用）对应变量J九-1、ZBZH程序（坐标新旧系转换）（已验算）″X1″?A：″Y1″?B：″X0″?D：″Y0″?F：″JJ″? T Lb1 0 :″XP″?C : ″YP″?EA＋(C－D)cos（T）＋(E－F)sin（T）→G:″XP1″:GB－(C－D) sin（T）＋(E－F) cos（T）→H″YP1″：HGoto 0输入说明：X1？Y1？为（公共点）新坐标系中坐标X0？Y0?为（公共点）原坐标系中坐标T? 为新坐标系X轴在原坐标中方位角（即坐标轴旋转角）XP? YP?为待求点在原坐标系中坐标计算结果：XP1=、YP1=为待求点在新坐标系中坐标。