5800导线平差程序
5800通用线路测量程序及编程原理
5800通用线路测量程序及编程原理关于坐标计算内核1、程序采用数据保存在扩充存储的方式，更换或补充数据不用修改程序。
2、线路采用线路号来区分，程序内部和使用时都用自然里程参加计算。
3、平曲线不用区分匝道和主线，所有平曲线都转化成交点资料，此交点可以缓和曲线不对称、不完整，主线通常由交点起算、匝道均由曲线起点起算。
每条平曲线交点格式为9个数据。
4、直线、圆曲线、缓和曲线分别用数学模型不同的计算子程序，保证了计算速度。
直线、圆曲线子程序为精确公式，缓和曲线计算子程序为原创的复化段数可调的复化高斯2节点公式，计算精度可自由控制。
（现暂改为固定4次复化）************************************************************************************** *************数据库编制规则基本规则：圆曲线半径右转为正，坡度上坡为正，角度单位为度，边桩转角右转为正,线路号为非负整数,数据库中的直圆里程应该在起点和缓圆各输入一个互差0.0001的里程。
圆直点同样处理。
一、目录部分：格式为：起点里程、终点里程、平曲线资料起始位置、竖曲线资料起始位置、坡度资料起始位置起点里程、终点里程、平曲线资料起始位置、竖曲线资料起始位置、坡度资料起始位置……各条线路依次连续排列。
二、平曲线部分：格式为：圆半径起点里程、缓圆里程、圆缓里程、终点里程、切线长、交点X、交点Y、切线方位角、圆半径起点里程、缓圆里程、圆缓里程、终点里程、切线长、交点X、交点Y、切线方位角、圆半径……终点里程、圆半径各交点依次连续排列。
主线按图纸设计数据整理，匝道的线元数据组合成交点组合时，每条曲线由直线、缓和曲线、圆曲线、缓和曲线、直线组成，对不完整缓曲要对数据进行简单处理，不完整缓曲可以划归任一侧，不完整缓曲所在的起（终）点里程改为负，交点坐标输入为起点坐标，切线长输入0即可。
fx5800平曲线竖曲线程序
PQX(5800平曲线程序) 算左、中、右任意角度及结构物坐标1、Lb1 0：Deg：Prog“A”2、E－D→J3、Arg (cos（J)+ i sin(J))→J4、Abs (J)÷J→H5、S 2÷(24R)-S4÷(2268R3)程序中有两个括弧+ i(S÷2-S3÷(240R2 ))→P6、JH÷2→M7、90÷(Rπ) →N8、（R+ReP (P)）tan(M)+ImP (P)→T9、M÷N+S →L10、C-T →U11、A-T(cos(D)+ i sin(D))→B12、If K<U:Then C-K→I:0→F:A-I(cos(D)+i sin(D))→X:Goto1:Else IfK<U+S:Then K-U→W:Goto2:Else If K<U+L-S:Then GotoA:Else If K<U+L:Then U+L-K →W:Goto2:Else GotoB :IfEnd: IfEnd: IfEnd: IfEnd13、Lbl 2:NW 2÷S→F:W-W5÷(40R2S2)+W9÷(3456R6S6)程序中有两个括弧+ i(W3÷(6RS)-W 7÷(336R4S4)程序中有两个括弧→I: If K<U+S:Then Goto3:Else GOto4: IfEnd14、Lbl 3:B+Abs(I)(cos(D+H Arg(I))+ i sin(D+H Arg(I))→X:Goto115、Lbl A:K-U-S→W:2NW+NS→F:R sin(F)+Imp (P)+i(R-R cos(F)+Rep (P)) →I:Goto316、Lbl B:K-U-L→W: 0→F:A+(W+T)(cos(E)+ i sin(E))→X:E→D:Goto117、Lbl4:E-HArg(I)→V:E→D:-H→H:A+T(cos(E)+i sin(E))-Abs(I)(cos(V)+i sin(V))→X:Goto118、Lbl 1:D+HF+Z→Q:“X+Y i=”:X+G(cos(Q)+i sin(Q))+Ｏ（cos(Q-Y)+i sin(Q-Y)）→X ▲19、Goto 0A(数据库)PQX子程序1、“DK”？K: “L-R”?G: “XJ”?Z:? Ｏ:?Y2、If K<141216.781(第二曲线起点)：Then 1300(半径) →R: 165 (缓和曲线长) →S:73688.363+61082.739i(交点坐标X+Y i) →A: 140083.413(交点桩号)→C: 297 51 48.67（与前一交点方位角）→D: 236 24 44.84(与后一交点方位角) →E: Return:If End例题：1、”DK”?K:”L-R”?G:”XJ”?Z:O:?Y2、If K<424.235:Then 1000→R:120→S:43812.504+48756.851i→A:-33.942→C：103°3°30°→D:91°8°38°→E:Return:IfEnd3、If K<1663.769:Then 1000→R:120→S:43797.960+49485.331i→A:693.821→C：91°8°38°→D:114°48°17°→E:Return:IfEnd4、If K<3972.554:Then 2000→R: 0→S:43277.777+50610.861i→A:1927.535→C：114°48°17°→D:99°46°44°→E:Return:IfEnd5、If K<6129.653:Then 3000→R:0→S:42846.791+53111.511i→A:4462.027→C：99°46°44°→D:118°18°43°→E:Return:IfEnd6、If K<7223.231:Then 1000→R:120→S:41916.431+54838.510i→A:6415.133→C：118°18°43°→D:80°12°52°→E:Return:IfEndSQX(5800竖曲线程序)1、Lbl 0:Prog“B”2、(B-C)R÷200→T3、If K<J-T:Then H-0.01B(J-K) →X:Else If K<J+T:Then K-J+T→D:H-0.01B(T-D)-D2÷(2R) →X:ElseH+0.01C(K-J) →X:IfEnd:IfEnd：4、X-0→X:Prog“Z”:“Z+Xi=”：X+Z+Xi→Z▲“Y=”：X+Y→Y▲4、Goto 0B（数据库）1、“DK”？K：If K≦140815（下一个竖曲线起点）：Then -1000(半径凸 +，凹 -)→R:-2.2(前一纵坡坡度) →B:2.9(后一纵坡坡度)→C: 140200（边坡点桩号）→J:219.755（边坡点高程）→H: Return:IfEnd例子：“DK”？K：If K≦2241.25：Then 8000→R: -0.4 →B: -2→C: 1920→J:31→H: Return:IfEnd If K≦2857：Then-5000→R: -2→B:0.35→C: 2300→J:23.4→H: Return:IfEndIf K≦2939.44：Then18000→R: 0.35→B:-0.35→C: 2650→J:24.625→H: Return:IfEnd If K≦3296.165：Then-10000→R: -0.35→B:0.861→C: 3000→J:23.4→H: Return:IfEnd If K≦3560：Then5000→R:0.861→B:-1.692→C: 3360→J:26.5→H: Return:IfEndIf K≦4117.9：Then-5000→R: -1.692→B:0.714→C: 3620→J:22.1→H: Return:IfEndIf K≦4955.317：Then12000→R:0.714→B:-0.321→C:4180→J:26.1→H: Return:IfEndIf K≦5255.341：Then-7000→R: -0.321→B:1.527→C: 5020→J:23.4→H: Return:IfEnd If K≦5594.478：Then4500→R: 1.527→B:-1.347→C:5320→J:27.98→H: Return:IfEnd If K≦5876.658：Then-7000→R: -1.347→B:0.525→C: 5660→J:23.4→H: Return:IfEnd If K≦6000：Then15000→R:0.525→B:-0.32→C: 5940→J:24.87→H: Return:IfEndZ (5800超高程序)1、Prog“C”:“B”?V2、0.04S÷(0.02+C) →D:A-T→E:E+D→F:E+S→G:G+L→I:I+S→M:M-D→J3、If K<E:Then -0.02V→Z:Z→Y::Goto1:Else If K<F:Then Goto2:Else If K<G:Then H((0.02+C)V(K-E)÷S-0.02V→Z:-Z→Y:Goto1:Else If K<I:Then HVC→Z:-Z→Y:Goto1：Else If K<J:Then H((0.02+C)V(M-K)÷S-0.02V→Z:-Z→Y:Goto1:Else If K<M:Then Goto3:Else -0.02V→Z:Z→Y:Goto1:IfEnd: IfEnd: IfEnd: IfEnd: IfEnd: IfEnd4、Lbl 2:IfH<0（零）:Then0.02VH→Z:(0.02+C)V(K-E)/S-0.02V→Y:Else(0.02+C)V(K-E)/S-0.02V→Z:-0.02VH→Y: IfEnd: Goto15、Lbl 3:IfH<0:Then0.02VH→Z:(0.02+C)V(M-K)/S-0.02V→Y:Else(0.02+C)V(M-K)/S-0.02V→Z:-0.02VH→Y: IfEnd: Goto1: Lbl 1C程序If K≦141216.781：Then 140083.413→A:0.03→C:-1→H:130→S:1229.28→L:855.675→T: Return:IfEndIf K≦142472.319：Then 141700.107→A:0.03→C:1→H:130→S:608.144→L:483.326→T: Return:IfEndIf K≦143839.524：Then 142976.487→A:0.02→C:-1→H:105→S:634.484→L:504.168→T: Return:IfEndIf K≦145570.280：Then 144509.787→A:0.02→C: 1→H:105→S:928.816→L:670.263→T: Return:IfEndIf K≦148062.26：Then 146191.812→A:0.02→C:-1→H:105→S:822.148→L:621.531→T: Return:IfEndPQXZFS平曲线正反算( 主程序)Lbl 4：“1.SZ=>XY”: “2.SZ=>XY”:?N:?S:Prog “SUB0”1/P→C:(P-R)/(2HPR) →D:180/兀→E:N=1=>Goto1:Goto2:Lbl 1：？Z:Abs(S-O) →W: Prog “SUB1”: “XS=”:X▲“YS=”:Y▲F-90→F: “FS=”:F▲DMS▲Goto4Lbl 2:?X:?Y:X→I:Y→J: Prog “SUB2”:O+W→S: “S=”:S▲“Z=”:Z▲Goto4正算子程序SUB10.1739274226→A:0.3260725774→B:0.0694318442→K:0.3300094782→L:1-L→F:1-K→M:U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C +MWD))) →X:V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C +MWD))) →Y:G+QEW(C+WD)+90→F:X+Zcos(F) →X:Y+Zsin(F) →Y反算子程序SUB2G-90→T:Abs((Y-V)cosT-(X-U)sinT) →W:0 （零）→Z: Lbl 0: Prog “SUB1”:T+QEW(C+WD) →L:(J-Y)cosL-(I-X)sinL →Z:IfAbs(Z)<1**10-6：Then0 →Z: Prog “SUB1”:(J-Y)/sinF →Z:ElseW+Z→W:Goto0:IfEnd数据库子程序SUB0Goto1Lbl 1：If S<***(线元终点里程)：Then***→G(线元起点方位角)：***→O（线元起点里程）：***→U（线元起点X）：***→V（线元起点Y）：***→P（线元起点曲率半径）：***→R（线元终点曲率半径）：***→H（线元起点至终点长度）：0或1、-1→Q：Return：IfEndLbl 1：If S<***(线元终点里程)：Then***→G(线元起点方位角)：***→O（线元起点里程）：***→U（线元起点X）：***→V（线元起点Y）：***→P（线元起点曲率半径）：***→R（线元终点曲率半径）：***→H（线元起点至终点长度）：0或1、-1→Q：Return：IfEnd ………………使用说明：1、一道路中线的前进方向（即里程增大的方向）区分左右，当线元往左偏时，Q=-1；当线元往右偏时Q=1；当线元为直线时Q=02、当所求点位于中线时，Z=0；当所求点位于中线左侧时Z取负值；位于右侧时Z取正直3、当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替4、当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧半径5、当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；当与圆曲线相接时，曲率半径等于圆曲线的半径。
5800线元法修正版程序编写
5800线元法修正版程序编写1.2版本，添加了遗漏的标点符号，修改了竖曲线与隧道，让隧道测量步骤简化一些。
此程序根据4850线元法程序改编，本人第一次编写程序内容可能有些不尽人意敬请谅解.程序名称不可更改，包括符号，也不要擅自删减子程序，如需更改必须与各子程序相呼应，否则将无法调动程序进行计算。
1~16页为程序说明与应用步骤具体解释，17~26页为程序编写，所有S开头后面有横杠的为各项数据库，如换线路更改数据库即可完成各项计算。
注意：程序编写时有下划线字符为一次性输入字符，不可拆分。
置仪点X坐标输入时必须输入小数，如果无小数将默认调动导线点数据库（例如3843.000那么请输入3843.0001并不影响计算）程序名（0-0）主程序目录（U=1）:数据库输入原始数据计算中边桩*（U=2）: 数据库输入原始数据反算线外点垂桩距状态（线路任一点坐标反算桩号及偏距）*（U=3）:计算平面位置检查记录表（U=4）:人工输入原始数据计算中边桩（数据库计算线元方位角）* （U=5）:人工输入原始数据推算开挖、填筑桩状态（不常用）（U=6）:数据库输入原始数据推算开挖、填筑桩状态*（U=7/7.1）:计算斜交斜做涵洞单跨桥、斜交正做涵洞七点放样状态并能避免曲线内单跨桥涵进出水口长度不一致，同时适用于曲线桥扇形布桩的盖梁放样，曲线桥扇形布桩的支座放样。
*（U=8）: 数据库输入原始数据计算锥坡放样（U=9）: 后方交会（U=10）: 数据库输入原始数据求线外两点交路线的桩号与夹角（U=11）: 圆心辐射法隧道断面超欠挖（适合各种隧道线型）*（U=12）: 手动输入原始数据求多边型面积、亩数*（U=13）: 数据库输入原始数据求设计高程（竖曲线）*（U=14）: 坐标正反算小程序（U=1）数据库输入原始数据计算中边桩（Y）置仪点X坐标（也可输入导线点编号，输入导线点编号即可使用导线点数据库）（O）置仪点Y坐标（当X坐标输入导线点编号后则此行不显示，因数据库已装入X和Y坐标）（B）步长（步长就是间隔长度，当要计算大量有相同距离间隔且相同中边桩，步长就是好东西了，输入数量每次计算自动加减桩号距离）（Q）所求桩号（计算第二个桩开始Q已经自动加入步长不用手动输入，但仍然会显示给你看，并可修改成别的桩号）（V）中桩到边桩的方位角，即与中桩前进方向方位角的夹角（注意：左转为负数，右转为正数）（D）中桩到边桩的距离（中桩输入0）计算完成会显示四个成果：DMS放样方位角（置仪点与所求桩号的放样方位角）I放样距离（置仪点与所求桩号的放样距离）所求桩号的X坐标所求桩号的Y坐标按—号继续下一步放样计算（U=3）数据库输入原始数据计算中边桩（偏差检测、资料编写）（B）步长（步长就是间隔长度，当要计算大量有相同距离间隔且相同中边桩，步长就是好东西了，输入数量每次计算自动加减桩号距离）(MM)所允许的最大误差（毫米为单位，比如高速的中线最大误差是20MM，就输入20）（Q）所求桩号（计算第二个桩开始Q已经自动加入步长不用手动输入，但仍然会显示给你看，并可修改成别的桩号）（V）中桩到边桩的方位角，即与中桩前进方向方位角的夹角（注意：左转为负数，右转为正数）（D）中桩到边桩的距离（中桩输入0）计算完成会显示四行成果：DMS放样方位角（置仪点与所求桩号的放样方位角）I放样距离（置仪点与所求桩号的放样距离）所求桩号的X坐标所求桩号的Y坐标按（—）号显示偏差检测成果：所求桩号的X坐标所求桩号的Y坐标所求桩号的双轴偏差X轴偏差Y轴偏差按-号下一步计算（U=4）人工输入原始数据计算中边桩（线元角计算）（Y）置仪点X坐标（O）置仪点Y坐标（B）步长（H+）起算点桩号(E)起算点X坐标(F) 起算点Y坐标(A)起算点方位角(+-C)起算点曲率（即为半径分之一，也就是1/半径，直线为0，左转为负，右转为正）(+-Z)终算点曲率（即为半径分之一，也就是1/半径，直线为0，左转为负，右转为正）(G)终算点桩号(Q)所求桩号(既是起算点桩号，因程序内字符转换冲突只好再输入) (V)左右夹角（左转为负，右转为正，中桩输入0）(D)中边桩距离（中桩输入0）计算完成会显示四个成果：（此处可忽略，参照U=1计算成果）计算完成按（AC/on）一次，然后按（EXIT）将内容（Prog 0-0）清除输入字母(R)，再按（EXE）显示为线元方位角数字（注：必须用度，分，秒抄记，因为计算结果为百进位，度分秒为60进位）注：此计算成果谨适用于数据库线元角输入（即：数据库曲线段线元角ZH~HY、HY~YH、YH~HZ段）（U=5）:人工输入原始数据推算开挖、填筑桩状态（不常用）（Y）置仪点X坐标（O）置仪点Y坐标（H+）起算点桩号(E)起算点X坐标(F) 起算点Y坐标(A)起算点方位角(+-C)起算点曲率（即为半径分之一，也就是1/半径，直线为0，左转为负，右转为正）(+-Z)终算点曲率（同上）（G）终算点桩号（W）实测大地高程（S）变坡点设计高程(I) =变坡点设计高程与实测大地高程之高差（正为填方、负为挖方）(P) 填方或挖方的坡比（若1:0.75则输入0.75填挖坡比全为正）（B）中桩到变坡点的宽度（D）=推算的中桩到开挖、填筑桩距离（D?）采用的中桩到开挖、填筑桩距离（Q）所求开挖、填筑桩号(V)左右夹角（注：只能输入线路右边90和左边-90因为开挖、填筑桩总是垂直的）计算完成显示四个结果：（同U=6显示结果相同，按—号下一步计算）（U=6）:数据库输入原始数据推算开挖、填筑桩状态（Y）置仪点X坐标（O）置仪点Y坐标（W）实测大地高程（S）变坡点设计高程（距离地面最近的那个平台或路面高程）（I）=变坡点设计高程与实测大地高程之高差（自动计算所得结果，此处不需输入。
卡西欧FX5800P编程计算器应用及程序
程序示例
U“Z=0 Y≠0”=0
K0+=0 H0=3014.6 ≥=1/74.09 D=2.9 E=2.46 R=1.671
【1】 直线段
实测坐标（m）
圆拱段超欠(m) 直墙超欠(m) 底板超欠(m) 至拱肩高差(m)
111.139 93.763 64.030 94.543 116.748
18°52ˊ49.81" 91°24ˊ13.56" 88°26ˊ53.49" 49°40ˊ03.20" 60°41ˊ09.73"
移 动 值（m）
DV=、0.254 D∧=、0.700 DV=、0.049 DV=、0.140 DV=、0.168
据目录（个数）是利用“ ”，按 FUNCTION->7->2->1->EXE 就可以看到个数了。
在程序中调用统计储存器中的数据时，是利用 List 调用的。具体格式如下： List X［A］+Y［A］i->C
式中 A 是点号，C 是任意储存单元。具体示例详见“程序部分”->“极程序放样” （2）、利用统计计算（SD/REG）
2、以 形式计算。以坐标换算为例：
坐标重合点 A 旧坐标 98042.831->SHIFT->RCL->A 11430.122->SHIFT->RCL->B A+Bi->E 新坐标值均为零，就不必另行输值，新旧坐标方位角差值 69°56"23"->W。 要求计算施工坐标 X=27，Y=29.1 的 B 点的测量坐标。 计算赋值：27->U，29.1->V。U+Vi->F
<=、0.663 =>、0.023 <=、0.218 <=、0.561 =>、0.298
5800平曲线和竖曲线计算程序。以验证过
5800平曲线和竖曲线计算程序。
以验证过5800平曲线程序1.主程序Lbl 4："1.SZ => XY"："2.XY => SZ"：？N：？S：Prog“SUB0”↙1÷P→C: (P-R)÷(2HPR) →D:180÷∏→E:N=1 => Goto1：Goto2：↙Lbl 1:? Z：Abs(S-O) →W：Prog "SUB1"："XS="：X◢"YS="：Y◢F-90→F：“FS=”：F▲DMS◢D:X+DCos(F)→X: "X=":X◢Y+DSin(F)→Y: "Y="：Y◢Goto4↙Lbl 2：？X：？Y：X→I：Y→J：Prog“SUB2”：O+W→S：“S=”：S◢“Z=”：Z◢Goto4↙2. 正算子程序(SUB1)0.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L：1-L→F：1-K→M：U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+B cos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD))) →X：V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QE MW(C+MW D))) →Y：G+QEW(C+WD)+90→F：X+Zcos（F）→X：Y+Zsin（F）→Y2. 反算子程序(SUB2)G-90→T：Abs((Y-V)cosT-(X-U)sinT) →W：0→Z：Lbl 0：Prog "SUB1"：T+QEW(C+WD) →L：(J-Y)cos（L）-(I-X)sin（L）→Z：IF Abs （Z）<10-6：T hen Goto1：Else W+Z→W：Goto 0：IfEndLb1 1：0→Z：Prog "SUB1" ：(J-Y)÷sin（F）→ZSUB0 数据库子程序Goto1↙同时保存多个曲线时的指针Lbl 1：IF S<***（线元终点里程）:Then***→G（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→H（线元起点至终点长度）：0或1、-1→Q：Return:IfEnd↙Lbl 1：IF S<***（线元终点里程）:Then***→G（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→H（线元起点至终点长度）：0或1、-1→Q：Return:IfEnd ……………..为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个Lbl，每增加个工程增加一个文件。
fx-5800P_单一闭附合图根导线近似平差程序PM3-1N
基于统计串列输入观测数据的单一闭附合图根导线近似平差计算程序(PM3-1N)(2) 程序PM3-1N，占用内存1848字节。
"CLOSED Or CONNECTING MAPPING TRAVERSE PM3-1N"\显示程序标题Norm 1n¼N:N-1¼DimZ定义额外变量"SIDE NUM=":N-1显示测边数"UNKNOWN POINT NUM=":N-2\显示未知点数"CLOSED(0) Or CONNECTING(Else)="?Z输0为闭合导线，其余数为附合导线Deg:Fix 3:FreqOn基本设置"XA(m),0 To BEAR A¼B(Deg)="?A输入已知点A的x坐标或输入≤0的数值If A>0:Then "YA(m)="?B:Else "BEAR A¼B(Deg)="?R:IfEndA>0时为输入A点y坐标，否则为输入A→B的方位角"XB(m)="?C:"YB(m)="?D输入已知点B的坐标C+Di¼U:U¼V B坐标的复数形式If A>0:Then Pol(C-A,D-B):Cls计算A→B的方位角J<0J+360¼J:J¼R"DIST A¼B(m)=":I显示A→B的水平距离"BEAR A¼B(DMS)=":R DMS显示A→B的方位角IfEndIf Z0:Then "XC(m)="?E:"YC(m)="?F附合导线，继续输入已知点C的坐标E+Fi¼V C坐标的复数形式"XD(m),0 To BEAR C¼D(Deg)="?G输入已知点D的x坐标或输入≤0的数值If G>0:Then "YD(m)="?H:Else "BEAR C¼D(Deg)="?S:IfEndG>0时为输入D点y坐标，否则为输入C→D的方位角If G>0:Then Pol(G-E,H-F):Cls计算C→D的方位角J<0J+360¼J:J¼S"DIST C¼D(m)=":I\显示C→D的水平距离"BEAR C¼D(DMS)=":S`DMS\显示C→D的方位角IfEndElse If R>180:Then R-180¼S:Else R+180¼S:IfEnd闭合导线计算A→B的反方位角"BEAR B¼A(DMS)=":S`DMS\显示B→A方向的方位角IfEndΣx¼M导线边长之和For 1¼I To NIf I=1:Then R+List Y[I]¼L:Else L+List Y[I]¼L:IfEnd导线边方位角If L>180:Then L-180¼L:Else L+180¼L:IfEndL>360L-360¼L:L<0L+360¼L判断方位角是否位于0°~360°之间Next3600(L-S)¼T以秒为单位的方位角闭合差"ANGLE CLOSE ERROR(S)=":T\显示方位角闭合差60(N)¼W以秒为单位的方位角闭合差限差Abs (T)>W"ANGLE CLOSE ERROR OVER!"显示角度闭合差超限-T÷N÷3600¼P:0¼O角度改正数，存储坐标增量累加和复数清零For 1¼I To N分配角度改正数并计算导线边的坐标增量If I=1:Then R+List Y[I]+P¼L改正后的导线边方位角Else L+List Y[I]+P¼L:IfEndIf L>180:Then L-180¼L:Else L+180¼L:IfEndL>360L-360¼L:L<0L+360¼L判断方位角是否位于0°~360°之间If I N:Then List X[I]∠L¼Z[I]:Z[I]+O¼O:IfEnd累加复数形式坐标增量Next3600(L-S)¼T以秒为单位的方位角闭合差检核计算"CHECK ANGLE CLOSE ERROR(S)=":T\显示方位角闭合差检核结果O-(V-U)¼O计算附合导线的坐标增量闭合差复数形式Int(M÷Abs(O))¼K计算导线全长相对闭合差"DELTA X(m)=":ReP(O)\显示导线X坐标增量闭合差"DELTA Y(m)=":ImP(O)\显示导线Y坐标增量闭合差"RELAT CLOSE ERROR=1÷":K\显示导线全长相对闭合差K<2000"RELAT CLOSE ERROR OVER!"显示全长相对闭合差超限For 1¼I To N-1分配坐标闭合差并计算未知点的坐标Z[I]-OList X[I]÷M¼X改正后的导线边坐标增量复数If I=1:Then U+X¼Y:Else Y+X¼Y:IfEnd未知点坐标复数形式Norm 1:"POINT n=":I\显示未知点号Fix 3:"Xp(m)=":Rep(Y)\显示未知点的X坐标"YP(m)=":ImP(Y)\显示未知点的Y坐标NextY-V¼Q计算导线的坐标检核结果"CHECK X(m)=":ReP(Q)\显示X坐标检核计算结果"CHECK Y(m)=":ImP(Q)\显示Y坐标检核计算结果"PM3-1N END"(3) 程序说明运行程序前，先按键进入REG模式，按键多次删除统计串列的数据，按坐标推算路线、从起点向终点方向将水平距离观测值顺序输入统计串列List X，将水平夹角观测值顺序输入统计串列List Y，当观测的水平夹角位于坐标推算路线左边时，输入正角度；当观测的水平夹角位于坐标推算路线右边时，输入负角度。
fx-5800p坐标正、反算;面积;土方量;单一水准路线平差;复数形式单一闭、附合与无定向导线近似平差程序。
fx-5800p坐标正、反算；面积计算；土方量计算；单一水准路线近似平差计算；复数形式单一闭、附合与无定向导线近似平差计算程序。
一.ZBZS(坐标正算程序)LbI 0(零):〝X1(零)=〞?S:〝Y1(零)=〞?G:〝FWJ=〞?F:〝HD=〞?D↙〝X2=〞:S+Dcos(F)→X◢〝Y2=〞:G+Dsin(F)→X◢Goto 0(零)↙本程序中输入及显示说明输入部分：X1=? 输入第1点X坐标；Y1=? 输入第1点Y坐标；FWJ=? 输入第1点至第2点方位角（即水平角）；HD=? 输入第1点至第2点水平距离。
显示部分：X2= 显示第2点X坐标；Y2= 显示第2点Y坐标。
二.ZBFS(坐标反算程序)LbI 0(零):〝X1=〞?Z:〝Y1=〞?W:〝X2=〞?X:〝Y2=〞?Y↙X-Z→A:Y-W→B↙√(A2+B2)→D:cos-1(A÷D)→J:If B≥0(零):Then J→J:EIse 360-J →J:IfEnd↙〝FWJ=〞:J◣DMS◢〝HD=〞:D◢Goto 0(零)↙本程序中输入及显示说明输入部分：X1=? 输入第1点的X坐标；Y1=? 输入第1点的Y坐标；X2=? 输入第2点的X坐标；Y2=? 输入第2点的Y坐标；显示部分：FWJ= 显示第1点至第2点的方位角（即水平角）； HD= 显示第1点至第2点的水平距离。
三.MIAN JI(面积计算程序)LbI 0(零):〝DIAN SHU=〞?N:〝X=〞?A:〝Y=〞?B↙0(零)→K:N-1→N:A→C:B→D↙LbI 1:〝X=〞?X:〝Y=〞?Y↙(X-C)(Y+D)÷2→Q:K+Q→K:X→C:Y→D↙Dsz N: Goto 1↙〝MIAN JI=〞:Abs(K+(A-C)(B+D)÷2→M◢Goto 0(零) ↙本程序中输入及显示说明输入部分：DIAN SHU=? 输入总点数；X=? 输入第1点的X坐标；Y=? 输入第1点的Y坐标；连续输入第2、3、4、………、各点的X、Y坐标。
卡西欧fx-5800p计算器测量程序
卡西欧fx-5800p计算器测量程序（作者：水江华）1.坐标反算程序"ZBFS QH1-4"◢Deg：Fix 3←┘"X+YiS(m）="?A←┘Lb1 0:"X+YiE(m)，0=>END="?B←┘Abs（B）=0=>Goto E←┘Arg（B-A)→JJ<0=>J+360→J"DIST(m）="：Abs（B-A)◢"α（DMS)=":J▶DMS◢Goto 0←┘Lb1 E:"ZBFS QH1-4=>END"2.坐标正算程序"ZBZS QH1-3"◢Deg：Fix 3←┘"XS(m)= "？A:"YS(m)= "？B←┘Lb1 1:"DIST(m)= "？L←┘If L≤0:Then Goto 2:IfEnd←┘"α（Deg)="？R←┘A+Lcos（R）→C←┘B+Lsin（R）→D←┘Fix 3："XE(m)= ":C◢"YE(m)= ":D◢Goto 1←┘Lb1 2:"ZBZS QH1-3=>END"3.单一闭、附合导线近似平差程序"BFHDXPC QH1-8"◢Deg:Norm 1:FreqOn←┘n→N←┘If List Y[N]=0:Then N-1→DimZ:Goto 0←┘Else N→DimZ:Goto 1: IfEnd←┘Lb1 0: "CLOSE Or CONECT TRA VE"◢"ANGLE NUM=":N◢"SIDE NUM=":N-1◢"UNKNOWN P NUM=":N-2◢"CLOSE(0)，CONNECT(≠0)="?→Z←┘Fix 3："XA(m)，≤0=>αA→B(Deg)="？A←┘If A>0：Then "YA(m)= "？B：Else "αA→B(Deg)="？R：IfEnd←┘"XB(m)= "？C:"YB(m)= "？D←┘C+Di→U：U→V←┘If A>0：Then Pol（C-A，D-B）：Cls←┘J<0=>J+360→J：J→R←┘"DIST A→B(m)="：I◢"αA→B(DMS)="：R▶DMS◢IfEnd←┘If Z≠0：Then "XC(m)= "？E:"YC(m)= "？F←┘E+Fi→V←┘"XD(m)，≤0=>αC→D(Deg)="？G←┘If G>0：Then "YD(m)= "？H←┘Pol（G-E，H-F）：Cls←┘J<0=>J+360→J：J→S←┘"DIST C→D(m)="：I◢"αC→D(DMS)="：S▶DMS◢Else "αC→D(Deg)="？S：IfEnd←┘Else If R>180：Then R-180→S：Else R+180→S：IfEnd←┘IfEnd←┘∑y→M←┘For 1→I To N←┘If I=1：Then R+List X[I]→L：Else L+List X[I]→L：IfEnd←┘If L>180：Then L-180→L：Else L+180→L：IfEnd←┘Nent←┘3600(L-S)→T←┘"α CLOSE ERROR(S)="：T◢40→W←┘Abs(T)>W=>"α CLOSE ERROR OVRE!"←┘-T÷N÷3600→P：0→O←┘For 1→I To N←┘If I=1：Then R+List X[I]+P→L←┘Else L+List X[I] +P→L：IfEnd←┘If L>180：Then L-180→L：Else L+180→L：IfEnd←┘If I<N：Then List Y[I]<L→Z[I] ：Z[I]+O→O：IfEnd←┘Nent←┘3600(L-S)→T←┘"CHECK α CLOSE ERROR(S)="：T◢U+O-V→O←┘Int(M÷Abs(O))→K←┘"DELTA X(m)= "：ReP(O)◢"DELTA Y(m)= "：ImP(O)◢"RELAT CLOSE ERROR=1÷"：K◢K<4000=>"RELAT CLOSE ERROR OVER! "←┘For 1→I To N-1←┘Z[I]-OList Y[I]÷M→X←┘If I=1：Then U+X→Y：Else Y+X→Y：IfEnd←┘Norm 1："POINT n="：I◢Fix 3："X(m)= "：ReP(Y)◢"Y(m)= "：ImP(Y)◢Nent←┘Y-V→Q←┘"CHECK X (m)= "：ReP(Q)◢"CHECK Y (m)= "：ImP(Q)◢Goto E←┘Lb1 1："NO DIRECT TRA VE"◢"ANGLE NUM="：N-1◢"SIDE NUM="：N◢"UNKNOWN P NUM="：N-1◢"XA(m)= "？A:"YA(m)= "？B←┘"XB(m)= "？C:"YB(m)= "？D←┘A+Bi→U:C+Di→V←┘0→H:U+List Y[1]<H→Z[1]←┘For 1→I To N-1←┘H+List X[1]→H←┘If H>180：Then H-180→H：Else H+180→L：IfEnd←┘以下未完程序由读者完成。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
5800计算器导线平差程序
一、程序用途及使用范围
本程序适用于一般导线复测平差计算，利用左角复测复合导线、闭合导线的平差时可直接使用。
复合导线平差时输入起始和终止边两个方位角，如果将终止边方位角输入为0，程序会自动转入闭合导线平差界面。
如果想用观测右角平差时，只需将程序中带下划线的地方按使用说明稍加改动即可。
二、源程序清单：1、JDPCA（文件名称）
0→M:0→I:0→K:“AFWJ”?A:“BFWJ”?B:“CZS”?→N:“JDRXBHCA”:24√¯(N)→P◣“∑(ZJ)”?→C:If B≠0 Then A-B-180N+C→D: EIseC-180(N-2)→D:IfEnd: “JDBHCA.f=”:D ▶DMS◣P<Abs(D)＝>Goto 4:“JDGZ=”:-D÷N→F:F▶DMS◣LbI 1:“CJ”?→P:B=0 And P=0＝>Goto3:“JGH=”:P+F→J:J▶DMS ◣A+J-180→E:If E>0 And E<360:ThenE:IfEnd:If E>360:ThenE-360→E:IfEnd: :If E<360:ThenE+360→E:IfEnd: “FWJ=”:E▶DMS◣“L”?→L:L=0 And B≠0＝>Goto3:“XO=”:Lcos(E)→X◣“YO=”:Lsin(E)→Y◣M+L→M:I+X→I:K+Y→K:E→A:Goto1:LbI3:“∑(L)=”:M◣“∑(XO)=”:I◣“∑(YO)=”:K◣Prog“DXPCA”：LbI4
DXPCA（文件名称，可单独运行）
“AX”?→A:“A Y”?→Z:If B≠0:Then “BX”?→C:“BY”?→D:“FX=”:I-C+A→F◣“FY=”:K-D+Z→W◣EIse “FX=”:I
→F◣“FY=”:K→W◣IfEnd:“DXBHCA.F=”:√¯((F)²+（W）²)→P◣“RXDXBHCAFM”?→R:“DXBHCAFM=”:M÷P→E ◣E<R＝>Goto1:F÷M→N:W÷M→O:LbI0:“L”?→L:“XG=”:-LN→Q◣“YG=”:-LO→V◣“XO”?→R:“YO”?→S:“XGH=”:Q+R→T◣“YGH=”:S+V→U◣“X=”:A+T→X ◣“Y=”:Z+U→Y◣X→A:Y→Z:Goto0:LbI1
三、程序符号含义
AFWJ?；BFWJ-?-----分别为起始边和终止边方位角，（闭合导线
平差时BFWJ输入0）。
JDRXBHCA=-----（24√¯(N)，三级导线角度容许闭合差（秒），根据导线级别改常数，如二级导线将24改为16，一级改10. CZS---------测站个数。
∑(ZJ)？-------各站测量角度总和。
JDBHCA.f=----------角度闭合差。
CJ?-----------各站角度实测值。
JDGZ=-------毎站角度该正值。
JGH=----------平差修正后角值。
FWJ=----------平差后两站间方位角。
L?--------两站间实测距离。
XO=,YO=---------坐标增量计算值。
RXDXBHCAFM?-------容许导线闭合差分母。
一级为15000，
二级为10000，三级5000.
DXBHCA.F----------计算的导线闭合差
∑(L)=∑(XO)=∑(YO)=------分别为边长X、Y坐标增量值总和。
XF=,YF=--------分别为X,Y的闭合差
AX,A Y,BX,BY-------分别为起点和终点的设计坐标值。
XG=，YG=-------分别为X、Y坐标增量改正值。
XGH=，YGH=-------分别为改正后X、Y坐标增量值。
X=,Y=----------平差后各站坐标值。
JD.ERR;DX.ERR-----分别为角度闭合差和导线闭合差不符合求。
四、程序使用方法说明
1、首先确定侧角方向，如果是观测的左角，可以直接启动程序进行计算，如果观测的是右角，将程序中带下划线的地方按下式修改：-D/N改成D/N; A+J-180改成A-J+180。
复合导线平差将A-B-180N+C改成A-B+180N-C;闭合导线平差将C-180(N-2)改成180（N-2）-C。
2、确定导线级别，按程序符号含义的规定修改程序参数。
3、按程序提示输入参数进行平差计算，本程序分两部分进行：
①、输入起始和终止方位角、测站个数，实测角度总和进行角度平差，如果不符合要求，给出错误提示，程序终止运行，符合要求时，输入实测角值和边长，程序给出平差改正后的角值和坐标方位角并且计算出坐标增量值，当复合导线计算边长输入0，闭合导线测角输入0时，本部分计算结束，程序会自动转入下一部分计算。
②、程序会自动给出边长、X增量、Y增量总和，接着提示输入起始点和终止点的坐标，程序自动进行导线平差，精度达不到规范要求时，程序给出错误提示，终止运行；符合要求时，提示输入边长，给出改正值，提示输入X、Y坐标增量值，程序给出改正后坐标增量和平差后的测站坐标值。
继续按提示输入直至平差结束。
五、计算实例。