FX-4850计算器曲线计算程序2011.721
ZHY程序FX5800
辛普森计算器公路测量常用程序修改版一、程序功能
LbI 0:＂X=＂?A: ＂Y=＂?B：＂Z=＂?C：＂1/（R1）＂?D：＂1/（R2）＂?E：＂QT=＂?F: ＂QT=＂?G: ＂V＂?V：＂W＂?W LbI 1：＂L＂?L：＂O＂?O：＂H＂?H：IF L>G：Then Goto 1：IFEnd
（E-D）/（G-F）→P：L-F→Q：P*Q→I
C+（I+2*D）*Q*（90/∏）→J
C+（I/4 +2*D）*Q*（45/（2*∏））→M
C+（3/4*I+2*D）*Q*（135/（2*∏））→N
C+(I÷2+2*D)*Q*(45÷∏) →K
＂X=＂:A+Q÷12*(COS（C）+4*（COS（M）+C0S（N））+2*COS （K）+COS（J）)+O*COS（J+H）→X◢
＂Y=＂:B+Q÷12*(SIN（C）+4*（SIN（M）+SIS（N））+2*SIN （K）+SIN（J）)+O*SIN（J+H）→Y◢
X-V→I：Y-W→J
PoL（I，J）
”J=”:J▲DMS▲
If J <0:Then “J=”:J+360°→J:J▲DMS ▲ if End Goto1
F5800计算器曲线计算JD程序（ZH~HZ）坐标计算放样公式
主程序：
30→Dim Z: "XJD="?B: "YJD="?D: "KJD="?M: "A="?A: "ZJ="?
Q: "R="?R: "LS1="?C: "LS2="?T: "F="?F:prog "A":Lb1 0: "XO="?G: "YO="?H: "LC="?L: "JJ="?Z: "PJ="?S:If L≤M-Z[11]+C:Then M- Z[11]COS(F) →N:D- Z[11]Sin(F) →E:Prog "ZHHY":Prog "FWJ": Goto 0:Else If L≤M-Z[11]+ Z[13]-T:
Then M- Z[11] →K:B- Z[11] COS(F) →N:D- Z[11]Sin(F) →E: Prog "HYYH":Prog "FWJ": Goto 0: Else If L≤M-Z[11]+ Z[13]:Then M- Z[11]+ Z[13] →K:F+AQ+180→V:B-Z[12] COS(V) →N:D- Z[12] Sin(V) →E: Prog "YHHZ":Prog "FWJ": Goto 0: Else "End" ▲
If End: If End: If End
说明：以上主程序需要四个副程序，缺少任何一个副程序不能运行！
"XJD="?: "YJD="?: "KJD="?为曲线交点坐标及桩号
"A="?为曲线左转取-1；右转取1
"ZJ="?为曲线转角值；"R="?为曲线半径
"LS1="?： "LS2="?为前后段缓和曲线长度
"F="?曲线起始方位角； " X0="?、" Y0="?:全站仪测站点坐标
" LC="?为曲线ZH-HZ段内所求任意点桩号；
"JJ="?为切线顺时针夹角；"PJ="?为左右侧所求点距设计线距离；。
Casio 4850计算器数据库型全线高程计算程序
CASIO4850计算器公路测量计算程序程序设计:魏加训2009.3.10Casio 4850计算器数据库型全线任意位置高程计算程序一、主程序：HLbl 0: X “BM-H”：U “HS”：“SXG=”: U+X◢输入水准点高程及后视，显示视线高程Lbl 1：{ P K }: “QS”K “K××+×××”: Prog“SJK”↙（换路修改此处为对应数据库名即可）C= C÷100: D= D÷100:↙D>C=> W=1:≠>W = -1⊿K<A=> I= -C:≠> I=D⊿↙Abs (A-K):L> R Abs(D-C)÷2 =>H=Z+IL: "H(Z-SJ)=": H◢Goto2↙≠>H= Z+IL+W(R Abs(D-C)÷2- Abs (A-K))2÷2÷R: "H(Z-SJ)=": H◢Lbl 2: {ONJP}:O “JGC(h)” : N“Z”:J “HP”:N=Abs N : B=H-O+N J÷100: "H(C) =" : U+X-P :Pause 0 : “H(S)=”: B◢“⊿h=”: B-U-X+P◢输入前视后显示实测高、设计高及高差（+填，-挖）Goto 1二、竖曲线数据库格式：K<下一竖曲线起点里程=>C=本竖曲线变坡点之前坡度(%前分子有正负): D=本竖曲线变坡点之后坡度(%前分子有正负): R=本竖曲线半径: A=本竖曲线变坡点桩号: Z=本竖曲线变坡点高程: Goto 1⊿↙三、数据库样例: SJKK<下一竖曲线起点里程=> C=本竖曲线变坡点之前坡度(%前分子有正负):D=本竖曲线变坡点之后坡度(%前分子有正负):R=本竖曲线半径: A=本竖曲线变坡点桩号: Z=本竖曲线变坡点高程: Goto 1⊿↙K<下一竖曲线起点里程=> C=本竖曲线变坡点之前坡度(%前分子有正负):D=本竖曲线变坡点之后坡度(%前分子有正负):R=本竖曲线半径: A=本竖曲线变坡点桩号: Z=本竖曲线变坡点高程: Goto 1⊿↙…………………………………………………………………………………………………………………..(注：如一条路有多个竖曲线要素则按顺序排列依次添加入数据库S.Z中)K≤下一竖曲线起点里程=> C=本竖曲线变坡点之前坡度(%前分子有正负): D=本竖曲线变坡点之后坡度(%前分子有正负): R=本竖曲线半径: A=本竖曲线变坡点桩号: Z=本竖曲线变坡点高程: Goto 1⊿↙Lbl 1↙四、输入与显示说明BM-H ？输入水准点的高程HS？输入水准点上的后视读数SXG= ×××显示视线高程“QS”？输入所求点前视读数K××+×××? 输入所求的桩号H(Z-SJ)= ×××显示所求桩号的路面中桩设计高程JGC(h) ? 输入结构层厚度或从设计路面高程点处垂直向下的厚度Z ？输入所求点距中线的边距(注:不分正负号均可，有中央分隔带的输入时需自行扣除中央分隔带半宽)HP？输入所求桩号处的横坡(%前分子有正负).H(C) = 显示计算得出的所求点的实测高程H(S) = 显示计算得出的所求点的设计高程⊿h= 显示设所求点的计与实测的高差（+填，-挖）注：（①当Z和JGC(h)为0时，则H=×××为路面中桩设计高程等于前面显示的H(S) ；②当Z为0，JGC(h)不为0时，则H=×××为路面垂直向下至某层的中桩设计高程；③当Z、JGC(h)、HP三者同时不为0时，则H=×××为路面垂直向下至某层的边桩对应Z位置的设计高程；④当Z不为0、JGC(h)为0、HP不为0时，则H=×××为路面边桩对应Z位置的设计高程；）附：数据库中用到的竖曲线基本计算公式：竖曲线起点桩号=本竖曲线变坡点桩号+本竖曲线之切线长T竖曲线切线长T=RAbs（i2-i1）÷2。
4850全线连线计算程序
CASIO4850全线坐标正反算程序主程序：ZHANG LIANGLbI Ø ：Deg:T"WJM":C=2Ø+4Ø(T-1):T=Ø:N"1.SR,2.JS":N=2=> Goto A⊿M"JZ-XX":O"QD-ZH":U"X-Ø ":V"Y-Ø ":G"FWJ ":Lbl B:{HRQP}: P"QD-R"：H"XX-CD":R"ZD-R":Q").-1,(.+1,--.Ø":W=H:Z =Ø:Prog"S1":Prog"S4":U=X:V=Y:O=O+H:G=F-9Ø:P=R:"…SR…NEXT…!”:Isz T: Dsz M:Goto B⊿Lbl A：Fixm:{XYSZN}:T=Ø:N"1.JS-XY,2.XY-FS": N=2=>Goto 2⊿S"JS-ZH=":Z"D=":Lbl 1:S＞Z［8T+1+C］=>Isz T:Goto 1⊿Prog "S3":W=S-O: Prog "S1": "JS-X=":X：Pause Ø:"JS-Y=":Y◢Goto AΔLbl 2: X"X…":Y"Y…":I=X:J=Y:Lbl 3: Prog "S3":W=Abs((Y-V)cos(G-9Ø)-(X-U)sin(G-9Ø:W＞H=> Isz T:Goto 3⊿Prog "S2":"FS-ZH…":S=O+W:Pause Ø:"D…":Z ◢Goto A⊿子程序S1A=.1184634425:B=.2393143352:N=.2844444444:K=. Ø4691ØØ77:L=.23Ø7653449:D=Q r W2(P-R)÷2HPR:E=Q r WP-1:Z[1]=G+KE+K2D:Z[2]=G+LE+L2D:Z[3]=G+.5E+.25D:Z[4]=G +(1-L)2D+E(1-L:Z[5]=G +(1-K)2D+E(1-K:F=G+E+D+9Ø：X=U+ZcosF+W(AcosZ[1]+BcosZ[2]+NcosZ[3]+BcosZ[4]+AcosZ[5:Y=V+ZsinF+W(AsinZ[1]+BsinZ[2]+NsinZ[3]+Bsin Z[4]+Asin Z[5：子程序S2W=Abs((Y-V)cos(G-9Ø)-(X-U)sin(G-9Ø:Z=Ø:Lbl 4:Prog "S1":L=G-9Ø+E+D:Z=(J-Y)cosL-(I-X)sinL:AbsZ<.ØØ1=>Goto5:≠>W=W+Z:Goto 4⊿Lbl 5:Z=Ø:Prog "S1":Z=(J-Y)÷sinF:子程序S3H=Z［8T+C+6:O=Z［8T+C+1］-H:U=Z［8T+C+2:V=Z［8T+C+3:G=Z［8T+C+4:P=Z［8T+C+5:R=Z［8T+C+7: Q=Z［8T+C+8:子程序S4Z［8T+C+1］=O+H: Z［8T+C+2］=U: Z［8T+C+3］=V:Z［8T+C+4］=G: Z［8T+C+5］=P:Z［8T+C+6］=H:Z［8T+C+7］=R: Z［8T+C+8］=Q: 结束程序编制说明：1.本程序是为CASIO-fx4850p编制的。
CASIO fx—4850P计算器高程控制网按间接平差计算程序
则 法 方 程 式为 ＣＸ ＋Ｌ ＝０ Ｌ
解 上 式 得 ｘ＝ 一 ｃ Ｌ Ｌ （） ３
Ｍｅ ： ｅ ： Ｎ ： ｉ ： ＝Ｔ ＋Ｎ： ＝０ ５ Ｎ ＋１） ｌ Ｄ ｇ Ｔ Ｍ Ｆｘ Ａ ｍ Ｗ ． Ｎ（ ：
式 中 ： ～ 观 测 值 改 正 数 向 量 ，Ａ 一 误 差 方 程
＝
单位 权 中误 差 ： ＝ ￣（ Ｐ Ｖ Ｕ ／ ［ Ｖ ］÷（ｎ—ｎ ）待 ｒ ）
Ｌ ｌ ： ［ Ｍ ＋ ］＝ ＝ｚ Ｉ ＝ｚ Ｐ 一 ［ ＋ Ｍ ＋ ｂ Ｆ ｚ Ａ＋ Ｊ Ｐ ［］ ［ ］ ｚ Ａ ２
Ｊｌ Ｌ ｌ Ｇ： ≠Ｍｊ Ｇ ｔ Ｊ＝０ Ｉ ｂ Ｊ ｏｏＤ ： ≠Ａｊ Ｇ ｔ Ｃ ｏｏ
将程 序 中划 线部 分 内容删 除 ， 节 省 程序 容 量 ２ ５字 可 ３
其 法方 程 式 为 ：
Ａ Ｐ Ｘ ＋Ａ Ｐ Ａ Ｌ＝０ （） ２
令
Ｃ＝ Ａ Ｐ Ａ， Ｌ Ｌ＝Ａ Ｐ Ｌ
节， 中问结果在计算结束后从 内存 中提出查看 。 ３ １ 主 程序 ．
ＧＣＷ ＪＰＣ Ｊ
Ｚ［ ＋ Ｍ ＋Ｉ Ｉ Ａ ６ ］ ≠Ｍ： ｏｏＱ Ｉ ： ｔ ＝０／ ￣Ｇ
Ｌ ｌ ＩｚＩ “ Ｌ Ｎ）＝” Ｚ『 ＋Ｎ＋８ ＋Ｉ Ｉ ｂ ｓ ： Ｌ （ Ｘ： Ａ Ｍ ］ ≠Ｎ
ｊ Ｇｏ Ｉ ｔＸ ｏ ＝０： ｂ ＩｚＩ “ Ｗ ）＝” ＺＩ ＋ Ｌ ｌＹ：ｓ ： Ｃ（ Ａ
求点高程中误差 Ｍｕｉ （）＝±Ｕ￣ｃ ／
２ 程序 功 能 ． 程序可进行高程控制 网按间接平 差计算 ， 计算数
据在扩大储存 变量 ｚＥ］中按 一维存 放 ， 了节省 内 ｊ 为 存单元 ， 解算 中采用 了内存单元叠加 的办法 ， 即由误
4850计算器计算程序
ZY
切线元素 各弧长所对 园心角(度） 1.137681652 4.957400286 8.77711892 12.59683755 16.41655619 20.23627482 24.05599346 25.8028 切线方位角 施工坐标X 施工坐标Y施工转测量坐标X 施工转测量坐标Y (度） 99.17336673 7635.349 3210.668 100.3110484 104.130767 107.9504856 111.7702043 115.5899229 119.4096415 123.2293602 124.9761667 5.957 25.925 45.777 65.427 84.786 103.768 122.289 130.583 0.059 1.122 3.513 7.221 12.230 18.518 26.056 29.911 7634.341 7630.108 7624.583 7617.790 7609.758 7600.525 7590.131 7585.003 3216.539 3236.082 3255.300 3274.107 3292.419 3310.156 3327.239 3334.812
5.9569 25.957 45.957 65.957 85.957 105.957 125.957 135.103
坐标计算模板 主编：黄绍敏 曲线测设元素与主点桩号计算（范本）
要素计算值(m)
切线方位角 99.17336673 JD点X坐标 7624.394 JD点Y坐标 3278.506 偏向(左-右+)) 1
曲线长Li
方位角(度） 9.1734 10.3110 14.1308 17.9505 21.7702 25.5899 29.4096 33.2294 34.9762
最新CASIOFX4850公路计算程序汇总
C A S I O F X4850公路计算程序0前言：传统公路测量中，使用的仪器设备和方法都很落后，需带着数学用表、曲线用表、计算盘、计算尺和算盘等一类的工具，完成外业测量工作。
计算器的出现，改变了这一局面。
高速公路建设中，长大曲线比比皆是，传统中对公路中线的测设方法，被极坐标法彻底的否定与取代，但大量的计算工作，只能带着提前计算好的线路逐桩坐标、高程资料，进行外业测量工作，机动性很差，现场查找也不方便。
这些问题都能在CASIO系列可编程计算器上得到很好的解决，对CASIO系列可编程计算器如何使用，直接影响到测量成果的质量和工作效率，本文将对CASIO系列可编程计算器快捷的计算方法进行分析与介绍。
1：以知线外任意点坐标，求对应线路里程在缓和曲线上，要计算任意里程的法线方向及任意宽度的边线坐标，非常简单。
但要计算任意一个已知坐标点，是对应哪一个里程法线方向上的点，就有一些困难。
很难推导一个这样的计算公式。
唯一的方法“渐进”，如果手工计算这可不是一个好方法。
但在有CASIO系列可编程计算器，如：FX-4500的情况下就变的非常简单了。
亦可用于直线和圆曲线的计算。
首先在缓和曲线上任选一点A为起始点，计算该点的坐标和切线方位角，通过坐标反算求起始点A与计算点B的方位角和距离，B点肯定对应A点切线方向上有一个垂足C点，把三点看成一个直角三角形，通过解直角三角形计算A C的距离，当该距离大于某一数值，如0。
001m，A点里程加AC的距离等于C点的里程，回到开始重新进入新一轮的计算，如果AC的距离小于某一规定值，则计算C点的里程与BC的距离即可。
求对应线路里程程序：主程序QLC (已知坐标求里程)Lb1 0：｛LＤE｝：Prog XH：Goto 0子程序：XH (循环)L1 Lb1 1L2 Norm：Prog LYYD：L3 PO1(D-X,E-Y)：W≤0=>W=W+360⊿L4 Z=W-I：A=V×cos Z：L=L+AL5 Abs A≥0.001=>Goto1：≠=>B=V×sinZ：Fix 3：“FXJL=”◢L6 L:Fix3:“DYLC=”◢程序中字母代表D 任意点X坐标，E 任意点Y坐标，DYLC 对应里程，FXJL 中线法线距离。
CASIO-fx4850多交点道路平曲线程序
№
PQX FANG YANG(平曲线放样主程序)（544步）
1
Deg:T“WENJIAN MING”: Q=1Ø5 (T-1) :LbI A:｛J｝: J“1SURU-6YIZB-7FY-8FS”
2
=7 =>Goto F⊿J=8=> Goto 7⊿J=6=> Goto 5⊿LbI 9:﹛ZL﹜: Z“JIGE JIAODIAN”:
8
“SQX QD=”:W =Z［7T+Q+6］-S◢“SQX ZD=”:D =Z［7T+Q+6］+S◢
9
Prog“A”:Isz T:Y=T:X=T:Dsz Z: Goto 8: O“LU KUANG”N“PO DU”: Fixm: Prog“AI”:
1Ø
LbI G:﹛VE﹜: V“YIQI H”E“MIAN XIAQU h”:
6
D= G:﹛YXUHGD﹜: Prog“AF”:“SHURU ZB”: Z［7T+Q+13］=Y“X3”:Z［7T+Q+14］
7
=X“Y3”:Prog“AL”:Prog“AG”:Prog“AH”:Prog“A”:Isz T:Y=T:X=T:Dsz Z: Goto 8:
8
LbI 5:fixm:﹛NOKL﹜: Prog“H”:Prog“AI”:
2
M=6 => Goto 5⊿M=8 => Goto 3⊿LbI 6:Fixm :｛ZL｝: Z“JIZU XIANXIN”: Z＞14 => Goto 6⊿
3
Z［12Ø+Q］=Z: T=3: L“1-JISUAN 2-SHURU”=2 => Goto 4⊿Fixm :｛WGHD｝: Prog“D”:
4850计算程序
ZDZBZBJS（匝道中桩、边桩坐标计算）LbI 0：{EG}：A″XA″：B″YA″：C″CA″：D″1/RA″：E″1/RB″：F″DKA″：G″DKB″LbI 1：{HLR}：H″DKI″：L″DL″：R″DR″：H>G＝＞Goto 3⊿P=（E-D）÷Abs（G-F）：Q= Abs（H-F）：I=P×Q：T=D+IJ=C+（I+2D）Q×90÷π◢M=C+（I÷4+2D）×Q×45÷（2π）：N=C+（3I÷4+2D）Q×135÷（2π）K=C+（I÷2+2D）Q×45÷πX=A+Q÷12×（Cos C+4（Cos M+ Cos N）+2 Cos K+ Cos J）◢Y=B+Q÷12×（Sin C+4（Sin M+ Sin N）+2 Sin K+ Sin J）◢U″XL″=X+L Cos（J-90）◢V″YL″=Y+L Sin（J-90）◢W″XR″=X+R Cos（J+90）◢Z″YR″=Y+R Sin（J+90）◢A=X：B=Y：C=J：D=T：F=H：Goto 1LbI 3：A=X：B=Y：D=E：F=G：C=J：Goto 0说明：XA？、YA？——ZH或HZ的X、Y坐标；CA？——切线方位角；1/RA？——起算点的曲率：直线段代0，曲线段代+1/R（右转）或-1/R（左转）；1/RB？——代终点曲率：即±1/R；DKA？——代起算点里程；DKB？——代待求段的终点里程；DKI？——代所求点里程；DL？——代左边桩距中桩的距离；DR？——代右边桩距中桩的距离；——是回车键符号。
（本程序适合于任何线形的计算）（2）ZXZBZJS（直线段中桩、边桩坐标计算）G″KA″A″XA″B″YA″D″XB″E″YB″：POl（D-A，E-B）：J＜0＝＞J=J+360：⊿F=J◢LbI 0：{KSP}：H=K-G：L=Abs H：Rec（L，F）：X=A+I◢Y=B+J◢Rec（S，F-90）：M″XL″=X+I◢N″YL″=Y+J◢Rec（P，F+90）：V″XR″=X+I◢W″YR″=Y+J◢Goto 0说明：KA——起算点里程 XA、YA——起算点坐标 XB、YB——后视点坐标K ——待求点里程 X、Y——待求点中桩坐标 S——左侧桩距中桩距离XL、YL——左侧桩坐标 P——右侧桩距中桩距离 XR、YR——右侧桩坐标（3）DBXMJJS（多边形面积计算）G=0：A=X：B=Y：C=A：D=B：LbI 0：{XY}：X＜0＝＞G=G+AD-BC：A= AbsG/2：A″A=″◢A=0.0015A：A″AA=″◢⊿G=G+XD-YC：C=X：D=Y：Goto 0说明：X、Y——待求点坐标 A——面积（m2） AA——面积（亩）测点坐标输入完毕后在X中输入-1。
Fx4850线元法测量程序
Fx4850最精简的万能坐标计算程序(线元法）1、JS-XY(计算主程序)Norm↙Z“KM”：D：N“PJ”:Prog“D”↙P=(S-R)÷T:Q= Z-C:I=PQ↙J=F+90Q(I+2R)÷π:J≥360=>J=J-360: ≠>J<0=>J=J+360⊿↙L=F+45Q(I÷4+2R)÷(2π):○=F+135Q(3I÷4+2R)÷(2π):K=F+45Q(I÷2+2R)÷π↙U=A+Q÷12×(cosF+4(cosL+cos○)+2 cos K+COS J):V =B+Q÷12×(sin F+4(sin L+ sin ○)+2 sin K+ sin J): W=J+N:X=U+D cosW:Y=V+D sin W:Fix 3: “XD=”:X▲“YD=”:Y↙本程序可以计算任何线型，包括直线、缓和曲线、圆曲线、不等半径曲线等，由其是匝道各曲线计算。
注：0为数字（零），○为字母（噢）KM-----计算点桩号（输入）D-------计算边桩距离左负，右正，中0（输入）PJ------路线右偏角，正交边桩为90（输入）XD------计算所得X坐标（计算得出）YD-----计算所得Y坐标（计算得出）FW---所求点桩号的切线方位角（计算得出）2、子程序：D（数据源）Z<803.961=>A=46670.1：B=13269.591:C=700：F=84°50′23″：R=0：S=0：T=103.961：Goto 0: ≠>Z<1085.902=> A=46679.45：B=13373.13:C=803.961：F=83°33′23″：R=0：S=0：T=281.941：Goto 0: ≠>Z<1363.366=> A=46711.091：B=13653.29:C=1085.902:F=83°33′23″：R=-1÷1000：S=-1÷1000：T=277.464：Goto 0:LbI 0:Fixm:A:B:C:F:R:S:T: ↙Z<本段线元的终点桩号=>A=起点X坐标：B=起点Y坐标:C=线元起点桩号：F=起点切线方位角：R=线元起点半径：S=线元终点半径：T=本段线元长度：Goto 0: ≠>Z<本段线元的终点桩号=>A=起点X坐标：B=起点Y坐标:C=线元起点桩号：F=起点切线方位角：R=线元起点半径：S=线元终点半径：T=本段线元长度：Goto 0: ≠>Z<？=>A=？：B=？:C=？：F=？°？′？″：R=？：S=？：T=？：Goto 0:Z<？=>A=？：B=？:C=？：F=？°？′？″：R=？：S=？：T=？：Goto 0:LbI 0:Fixm:A:B:C:F:R:S:T: ↙说明：字程序Z〈后面为每段线元的终点桩号，可以把你标段所有线元段都输入：其程序中的字母：A、B为起点坐标；C为起点桩号；F为起点方位角：R、S为起终点半径（以倒数输入，如R60，1÷60）分左右转如左转为-1÷60,右转为1÷60，直线为0；T为线元长度。
卡西欧4850计算器程序
卡西欧fx-4850P万能放线程序主程序：ZBJS↙◢“◢：H◢↙注：↙↙↙↙注：此程序可把全线路要素输入数据库，能执行直线、高程放线正反算。
卵形曲线。
卡西欧fx-4850P直线、各种曲线任意放中线(匝道)1：主程序：TYQXJS“1，KS=＞XY”：“2，XY=＞KS”：N：U“X0”：V“Y0”：o“K0”：G“FWJ”：H“LS”：P“R1”：R“R2”：Q“Y+1；0；Z-1”：C=P-1：D=（P-R）÷（2HPR）：E=180÷π：N=1=＞Goto 1:≠＞Goto 2↙◢“Goto 1↙2：正算子程序：S U B1A=0.1739274226:B=0.3260725774:Z=0.0694318442:L=0.33 00094782:F=1-L:M=1-Z:X=U+W(Acos(G+QEZW(C+ZWD))+Bcos (G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+M WD))):Y=V+W(Asin(G+QEZW(C+ZWD))+Bsin(G+QELW(C+LWD)) +Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))):F=G+QEW(C +WD)+90:X=X+Scos F :Y=Y+SsinF↙3坐标反算子程序 SUB2T=G-90：W=Abs((Y-V)cosT-(X-U)sinT) :S=0:Lbl 0: Prog “SUB1”：L=T+QEW（C+WD）:S=（J-Y）cos L-(I-X)sin L : Abs S＜10-6=＞Goto 1: ≠＞W=W+S:Goto 0↙Lbl 1:S=0:Prog“SUB1”:S=(J-Y)÷sin F↙注：o→代表：欧 0→代表零使用说明：1规定：⑴本程序可执行直线、圆曲线、缓和曲线（完整或非完整型）、复曲线、卵形曲线的任意里程中边桩坐标的正反算。
4850曲线计算程序
ZDXLbI 0Prog “FWJ”I”S”▲LbI 1｛ALP｝︰A”J”︰L︰E=X+Lcos (A+J) ︰F=Y+Lsin (A+J) ︰”X”︰E▲“Y”︰F▲P=1=＞Goto 0︰≠=＞Goto 1︰⊿FWJFix 4｛XYMN｝︰X”X0”︰Y”Y 0”︰M”X1”︰N”Y1”︰P ol(M-X,N-Y)FixmJ＜0=>J=J+360︰⊿”A”︰J→DMS▲“S”︰IHYQXLbI 3Fixm﹛K﹜︰Prog “DAT”Prog “QJS”“X”:X▲“Y”:Y▲Goto 3SQXFixmLbl 6Prog “DAT1”Prog “SJS”“H”:X▲Goto 6SDDMDefm 1LbI 0Prog “DAT1”Prog “SJS”Z[1]=XProg “DAT”Prog “QJS”（此处计算坐标必须是隧道中心点坐标）A=X:B=Y:C=D:F=Z[1]｛Z｝︰Z”DM”FixmLbI 1｛VGEH｝︰G”X”︰H”Y”︰E”H”︰V”JING”M=(G-A)cos C+(H-B)sin C︰L=(H-B)cos C-(G-A)sin C︰L=Abs L︰N=E-F︰”K-GAI”︰M▲Prog “DAT2”{W}:W=1=＞Goto 0︰≠=＞Goto 1︰⊿JSKDefm 3LbI 0﹛K﹜︰Prog “DAT”G=Z[3]Prog “QJS”A=X:B=Y:C=D:FixmZ[3]=G:G=Z[1]:H=Z[2]:｛GH｝︰G”X”︰H”Y”Z[1]=G:Z[2]=H:M=(G-A)cos C+(H-B)sin C:”K-GAI”︰M▲K=K+M:”LICH”:K▲L=(H-B)cos C-(G-A)sin C:”BIAN”︰L▲Goto 0QJSS=B²÷24R︰M=C²÷24R︰N=B÷2-B^3÷240R²︰O=C÷2-C^3÷240R²︰T=(R+M-(R+S)cos J)÷sin J+N︰U=(R+S-(R+M)cos J)÷sin J+O︰L=JRπ÷180-B÷2-C÷2+B+C︰H=K-I︰V=L-C︰H≤V=＞W=T︰Z=A+180︰≠=＞W=U︰Z=A+PJ︰⊿E=X+UcosZ︰F=Y+Usin ZH>0=＞H<B=＞Prog”H”︰⊿⊿H≥B=＞H≤V=＞Prog “y”︰⊿⊿H>V=＞H<L=＞H= L-H︰Prog “H”︰H=K-I︰Y=-Y︰D=2A-D+PJ︰⊿⊿H≤0=＞X=H︰Y=0︰D=A︰⊿H≥L＝＞X=L-H︰Y=0︰D=A+PJ︰⊿H≤V=＞O=A︰≠＝＞O=A+PJ +180︰⊿M=E+Xcos O-YsinO︰N=F+Xsin O+Ycos O﹛GQ﹜︰G”B”︰G=0=＞Q= 0︰≠=＞Q”J”︰⊿X=M+Gcos (D+Q) ︰Y=N+Gsin (D+Q) YX=N+Rsin (((2H-B)÷2R)×180÷π) ︰Y=P (S+R(1-cos (90(2H-B)÷πR))) ︰D=A+P (90(2H-B)÷Rπ)HX=H-H^5÷40R²B²+H^9÷3456R^4B^4︰Y=(H^3÷6RB-H^7÷336R^3B^3+H^11÷42240R^5B^5)P︰D=A+P(90H²÷RBπ)SJSW=C-B︰L=WR︰W＜0=＞G=-1︰L= -L︰≠＝＞G=1︰⊿T=L÷2︰L=WR︰E=T²÷2÷R︰H=K-I︰H≤0=＞F= BH︰⊿H>0=＞F=CH︰⊿Abs H<T=＞F=F+(T-A bsH)²÷2÷RG︰⊿X=F+AJDFixmO=0︰V≠1=＞Goto 2︰⊿Pol(N-P,L-Q) ︰J<S=>Goto 3︰⊿J>T=>Goto 3︰⊿o=1︰U=I-R︰U<0=>”R-WAI”︰≠=>”R- NEI”︰⊿U=Abs U︰U▲Goto 3Lbl 2U=R²-(N-P)²︰U<0=> “WUJIE”Goto 3︰⊿U=√U+Q ︰Pol(N-P,U-Q) ︰J<S=>Goto 3 ︰⊿J>T=>Goto 3 ︰⊿o=1 ︰U=L-U ︰U<0=>”H-WAI “︰≠=>”H-NEI”︰⊿U=Abs U ︰U▲Lbl 3DATK<135200＝＞X=1408︰Y=0︰C=280：B=280︰R=6000︰J=23º51º50º︰A=0︰I=131624.36︰P=-1︰Goto 0︰⊿K<137000＝＞X=4475.753︰Y= -1357.128︰A=336º08º10º︰J=20º29º59º︰B=310︰C=310︰R=5000︰I=13 5290.64︰P=1︰Goto0︰⊿K<158950＝＞X=158082.391︰Y= 10000︰A=11º16º41º︰J=11º16º41º︰B=260︰C=260︰R=7000︰I=15726 5.73︰P=-1︰Goto0︰⊿K<166710＝＞X=16096.745︰Y= 10000︰A=0︰J=13º34º14º︰B=280︰C=280︰R=6000︰I=164 591.784︰P=1︰Goto0︰⊿K≥166710＝＞X=14549.512︰Y= 7357.2034︰A=176º15º7.8º︰J=13º42º56.5º︰B=260︰C=260︰R=7000︰I=1671 80.249︰P=1︰Goto0︰⊿Lbl 0DAT1{K}:K<159649=>A=16.692 ︰B=0.003 ︰C=-0.003 ︰R=20000 ︰I=157950 ︰Goto 4 ︰⊿K<161960=>A=8.772︰B=-0.006 ︰C=0︰R=20000 ︰I=161900 ︰G oto 4 ︰⊿K<164558=>A=8.772︰B=0 ︰C=0︰R=20000 ︰I=164500 ︰Goto 4 ︰⊿K<165700=>A=15.152 ︰R=20000︰B=0.0058 ︰C=-0.002︰I=16 5600 ︰Goto 4 ︰⊿K<166750=>A=12.939︰B=-0.002 ︰C=0.003︰R=20000 ︰I=166700 ︰Goto 4 ︰⊿K<168287=>A=17.439︰B=0.003 ︰C=-0.0057︰R=20000 ︰I=16 8200 ︰Goto 4 ︰⊿Lbl 4DAT2Z=330=> Goto A︰⊿Z=320=> Goto B︰⊿Lbl AP=2.27︰Q=0︰R=7.01︰S=0︰T=60︰Prog “JD”:o=1=>Goto Z︰⊿P=2.27︰Q=0︰R=7.01︰S=60︰T=107º53º55º︰Prog “JD” :o=1=>Goto Z︰⊿P=1.217︰Q=3.273︰R=3.57︰S=107º53º55º︰T=157º35º24º︰Prog “JD”:o=1=>Goto Z︰⊿Lbl BP=2.27︰Q=0︰R=7.01︰S=0︰T=60︰Prog“JD”:o=1=>Goto Z︰⊿P=2.27︰Q=0︰R=7.01︰S=60︰T=107º53º55º︰Prog “JD” :o1=>Goto Z︰⊿P=1.217︰Q=3.273︰R=3.57︰S=107º53º55º︰T=157º35º24º︰Prog “JD”:o=1=>Goto Z︰⊿Labl Z程序使用说明此程序适应于CASIO 4850计算器ZDX程序（支导线计算）X0，Y0测站点坐标X1，Y1后视点坐标A 反算方位角S 反算距离J 观测左角L 观测距离X，Y 计算坐标FWJ程序（方位角计算）X0，Y0测站点坐标X1，Y1后视点坐标A 反算方位角S 反算距离HYQX程序（基本型曲线计算）K 计算里程B 计算边距J 右偏角（90为右侧，270为左侧，支持任意角）X，Y 计算坐标SQX程序（竖曲线计算）K 计算里程H 中桩高程SDDM程序（隧道断面计算）K 计算断面所在里程，半径越小，里程越需要准确，7000半径里程偏差十米，误差约1cm，直线无误差。