卡西欧4850计算器程序

合集下载

4850计算器程序

卡西欧4850计算器计算程序一、XLZBJS(线路坐标计算)主程序：P＝L2÷(24R)-L4÷(2688R3)M＝L÷2－L3÷(240R2)B＝90L÷(πR)T＝(R＋P)tan(AbsA÷2) ＋MC“ZH”＝D－TG“HY”＝C＋LH“YH”＝G＋(AbsA－2B) πR÷180I“HZ”＝H＋LLbl 0：{K}K≤G＝＞Goto 1：≠＝＞Goto 2△Lbl 1：（K－C）≤0＝＞Goto 5△X＝K－C－(K－C)5÷(40R2L2)Y＝(K－C)3÷(6LR) －(K－C)7÷(336R3L3)J＝√（X2＋Y2）O＝tan-1(Y÷X)“X0＝”：Q＝N＋T cos (F＋180) ＋J cos (AbsAO÷A ＋F)◢“Y0＝”：S＝E＋T sin (F＋180) ＋J sin (AbsAO÷A ＋F)◢W＝F＋3AbsAO÷AGoto 4△Lbl 2：K≥H＝＞Goto 3△X＝R sin ((K－G) ×180÷(πR) ＋B) ＋MY＝R(1－cos ((K－G) ×180÷(πR) ＋B)) ＋PJ＝√（X2＋Y2）O＝tan-1(Y÷X)“X0＝”：Q＝N＋T cos (F＋180) ＋J cos (AbsAO÷A ＋F)◢“Y0＝”：S＝E＋T sin (F＋180) ＋J sin (AbsAO÷A ＋F)◢W＝F＋AbsAB÷A＋AbsA(K－G) ×180÷(πR) ÷AGoto 4△Lbl 3：K≥I＝＞Goto 6△X＝K－I－(K－I)5÷(40R2L2)Y＝(I－K)3÷(6LR) －(I－K)7÷(336R3L3)J＝√（X2＋Y2）O＝tan-1(Y÷X)“X0＝”：Q＝N＋T cos (F＋A) ＋J cos (F＋A＋180－AbsAO÷A )◢“Y0＝”：S＝E＋T sin (F＋A) ＋J sin (F＋A＋180－AbsAO÷A )◢W＝F＋A－3AbsAO÷AGoto 4△Lbl 4：“FWJ”：W◢{Z}U＝Q＋Z cos (W－90)“XL＝”：U◢“YL＝”：V＝S＋Z sin (W－90)◢“XR＝”：U＝Q＋Z cos (W＋90)◢“YR＝”：V＝S＋Z sin (W＋90)◢Goto 0△Lbl 5：“X0＝”：Q＝N＋(D－K) cos (F＋180 )◢“Y0＝”：S＝E＋(D－K) sin (F＋180 )◢W＝FGoto 4△Lbl 6：“X0＝”：Q＝N＋(T＋K－I) cos (F＋A )◢“Y0＝”：S＝E＋(T＋K－I) sin (F＋A )◢W＝F＋AGoto 4△“END”子程序（运算程序）D＝交点里程N＝交点的X坐标E＝交点的Y坐标A＝曲线转向角R＝曲线半径L＝缓和曲线长F＝第一切线方位角（ZH—JD）Prog“XLZBJS”◢“END”说明：计算时只需进入子程序（运算程序），无需进入主程序。

4850计算器正反算程序

4850计算器正反算程序一、程序功能及原理1.功能说明：本程序由一个主程序(TYQXjs)和两个子程——正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。

另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。

本程序也可以在CASIO fx-4500P计算器及CASIO fx-4850P计算器上运行。

2．计算原理：利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据。

二、源程序1.主程序(TYQXjs)Defm 2:"1.SZ => XY"："2.XY => SZ"：N：Prog”DAT1”:U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"RN"：Q：C=1÷P：D=(P-R)÷(2HPR)：E=180÷π：N=1=>Goto 1：≠>Goto 2Δ←Lbl 1：{SZT}：SZ:T”ANG”：W=Abs(S-O)：Prog "SUB1"："XS=":X=X◢"YS=":Y=Y◢Goto 3←┘Lbl 2：{XY}：XY：I=X：J=Y：Prog "SUB2"："S=":S=O+W ◢"Z=":Z=Z◢Goto 2Lbl3:J=0:I=Pol(X-Z[1],Y-Z[2]):F=J:F<0=>F=F+360:≠>F>360=>F =F-36ΔΔ”DIST=”: I◢”FW=”:F→DMS◢Goto 12. 正算子程序(SUB1)A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469 100770:L=0.2307653449:M=0.5:F=1-L:M=1-K:X=U+W(Acos( G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Ncos(G+QE MW(C+MWD))+Bcos(G+QE(1-L)W(C+(1-L)WD))+Acos(G+ QE(1-K)W(C+(1-K)WD))):Y=V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD) )+Nsin(G+QEMW(C+MWD))+Bsin(G+QE(1-L)W(C+(1-L)WD))+Asin(G+QE(1-K)W(C+(1-K)WD)))：F=G+QEW(C+WD)+T：X=X+ZcosF：Y=Y+ZsinF3. 反算子程序(SUB2)T=G-90：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog "SUB1"：L=T+QEW(C+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ←┘Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF4.曲线元要素数据库：DAT1J=1=>U=xxxx.xxx:V=xxxx.xxx:O=xxxx.xxx:G=xxxx.xxx:H=x xxx.xxx:P=xxxx.xxx:R=xxxx.xxx:Q=xxΔJ=2=>U=xxxx.xxx:V =xxxx.xxx:O=xxxx.xxx:G=xxxx.xxx:H=xxxx.xxx:P=xxxx.xxx: R=xxxx.xxx:Q=xxΔ……….J=n=>U………. Δ(注：如有多个曲线元要素继续添加入数据库DAT1中)三、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

CASIO4850坐标正反算程序说明书

CASIO fx4850线路正反算测量程序（全站仪伴侣）CL-测量Norm :Lbl 1:Cls:{U}:U"1 ZS 2 FS"=1 =>Prog "ZS":≠>U=2 =>Prog "FS":≠>Goto 1△△D-读Lbl 1:K[B]=Z[G+B]:Dsz B:Goto 1 FS-反算Lbl 1:Cls:{KVW}:V"X"W"Y":Lbl 2:Prog "P":Pol(V-X,W-Y):J=J-A:I=Rec(I,J):K=K+I:AbsI≥E-3=>Goto 2:≠>Prog "WZ":Goto 1△HQ-缓曲N=I-I^5÷40S2+I^9÷3456S^4-I^13÷599040S^6:M=I^3÷6S-I^7÷336S^3+I^11÷42240S^5:Z=90IAbsI÷πS:J=90-(90-tan-1(M÷N))AbsI÷I:I=√(M2+N2)HY-缓圆I=O:Prog "HQ":Prog "XY":I=K-L-O HZ-缓直I=P: Prog "Y": Prog "XY":I=Q:S=RQ:Prog “HQ”:J=Z-J:Prog "XY":I=K-L-O-P-QML-目录Cls:Fix 3:C=Z[100]:D=Z[101] :Lbl 1:B=6:G=C-6D:Prog "D":K≥L =>K≤M =>C=N:D=O:E=P:F=Q: Goto 2△△Dsz D:Goto 1:Lbl 2P-平Prog"ML": Prog"PM"PM-平面Lbl 1:B=9:G=C-8D-1:Prog "D":K≤T =>Prog "PQ":Goto 2△Dsz D:Goto 1:Lbl 2 PQ-平曲Prog "ZH":K≤L =>Prog "ZX":≠>K≤L+O =>Prog "HQ":≠>Prog "HY":K≤L+O+P =>Prog "Y":≠> Prog "HZ":K≤L+O+P+Q => Prog "HQ" :≠>Prog "ZX"△△△△Prog "XY"PY-偏移Lbl 1:I=0:J=90:{IJ}:I"L"≠0 =>Z=J"<":Prog "XY":Prog "ZB":Goto 1△PZ-平纵Prog "ML":Prog "PM":Prog "ZM"SQ-竖曲Z=N+(K-M)L:K≤M-T =>Z:≠>Z=Z- (K-M+T)2÷2O△WZ-位置"KJ":K:Pause1:J▲XY-XYJ=A+J: A=A+Z:X=X+Rec(I,J):Y=Y+JY-圆J=90I÷πR: I=2RSinJ: Z=2JZ-纵Prog"ML": Prog"ZM"ZB-坐标A=360Frac((A+360)÷360):"ZB":X:Pause 1:Y:Pause 1ZH-直缓X=M:Y=N:A=S:I=K-L:S=ORZM-纵面Lbl 1:B=5:G=E-4F-1:Prog "D":T=Abs(O(P-L)÷2):K≤M+T =>Prog "SQ":Goto2△Dsz F:Goto 1:Lbl 2ZS-正算{U}:Lbl 1:{K}:U"1 B"=1=> Prog "P": Prog "PY" :Goto 1△Prog "PZ":Prog "ZB": Z▲Goto 1ZX-直线J=0:Z=0使用说明运行CL程序首先程序提示正反算1 ZS2 FS？输1正算、输2反算，其它重选。

CASIO4850计算器程序

缓和曲线段切线方位角1、A:B:C＂ZH＂:E＂L0＂:Q:R2、F=A+B（90（Q-C）2÷（πRE））▲字母说明：A—起始边方位角B—偏角，左偏为“-1”，右偏为“1”C＂ZH＂—ZH点里程E＂L0＂—缓和曲线长Q—所求点里程R—曲线半径F—切线方位角圆曲线段切线方位角1、A:B:C＂HY＂:E＂L0＂:Q:R2、F=A+B（90（E +2（Q-C））÷（πR））▲字母说明：A—起始边方位角B—偏角，左偏为“-1”，右偏为“1”C＂HY＂—HY点里程E＂L0＂—缓和曲线长Q—所求点里程R—曲线半径F—切线方位角1、D＂X0＂:F＂Y0＂: A:K2、Lbl 03、{W}: {Q}4、W=0=>Goto 1 △5、W≠0=>Goto 2 △6、Lbl 17、X = D + cos A（Q-K）▲8、Y = F + sin A（Q-K）▲9、Goto 010、Lbl 211、X = D + cos A（Q-K）+ W cos（A+90）▲12、Y = F + sin A（Q-K）+ W sin（A+90）▲13、Goto 0字母说明：D＂X0＂—起始点X坐标F＂Y0＂—起始点Y坐标A—起始边方位角K—起点里程W—边距，左边为“－”，右偏为“+”Q—所求点里程1、Defm 12、R : T : A : B: L : C : D : F3、E=C+L : M=（L－L 3÷（40R2）+ L 5÷（3456R4））÷cos （30L ÷（πR））4、B>0=>N=A+30L÷（πR）: ≠=>N= A－30L÷（πR）△5、Z[1]=D+McosN :G=F+MsinN:H=E+AbsBπR÷180－L: I=H+L6、M=2Tcos（AbsB÷2）: N=A+B÷2 : J=D+McosN: K=F+MsinN7、Lbl 38、{Q}: {W}: W : M =Q－C : M< 0=>Goto 3 △9、M=Q－I－0.005 : M>0=>Goto 3 △10、M=Q－E : M≤0=>Z=2 : Prog＂B＂△11、Q－E>0=>Goto 1 △12、Goto 213、Lbl 114、M=Q－H : M<Q=> Prog＂A＂△15、Q－H≥0=>Z=1 : Prog＂B＂△16、Lbl 217、W≠0=>Goto 5 △18、X=U + O cos P ▲19、Y=V + O sin P ▲20、W=0=> Goto 3 △21、Lbl 522、U=U+ O cos P: V= V+ O sin P: Prog＂D＂: X=U+Wcos（P+90）▲23、Y=V+Wsin（P+90）▲24、Goto 3子程序 A1、M=Q－E : O=2Rsin（90M÷（πR））: N=90÷（πR）×（L+M）2、B<0=>P=A－N △3、B>0=>P=A + N △4、U=Z[1] : V=G子程序 B1、Z = 2 =>M= Q－C △2、Z = 1 =>M= I－Q △3、X = M－M5÷（40R2L2）+ M9÷（3456R4L4）: Y= M3÷（6RL）－M7÷（336R3L3）+ M11÷（42240R5L5）4、O = √（X2+Y2）: N = 30M2÷（πRL）5、Z = 1 => Goto 1 △6、B<0=>P=A－N △7、B>0=>P=A + N △8、U=D: V=F9、Goto 210、Lbl 111、B<0=>P=A+B+N+180 △12、B>0=>P=A+B－N+180 △13、U=J : V=K : Lbl 2子程序 C1、Q－E<0=> Goto 1△2、Q－H<0=> Goto 2△3、B<0=>Z= 4 : Goto 3 △4、B>0=> Z= 6 : Goto 3 △5、Lbl 26、B<0=>Z= 1 : Goto 3 △7、B>0=> Z= 2 : Goto 3 △8、Lbl 19、B<0=>Z= 3 : Goto 3 △10、Z = 511、Lbl 3子程序 D1、N = 90L÷（πR）+180M÷（πR）: S=90M2÷（πRL）: Prog＂C＂2、Z = 1 => P = A－N △3、Z = 2 => P = A + N △4、Z = 3 => P = A－S △5、Z = 4 => P = A + S + B △6、Z = 5 => P = A + S △7、Z = 6 => P = A + B－S △备注：当曲线没有设缓和曲线时，需将“D”子程序中的“S=90M2÷（πRL）”更改为“S=0”，再进行计算。

卡西欧4850坐标计算程序说明

一．文件名:GSZFS（高斯-勒让德）Deg:I=0:I“JS.1,FS.2,SJK.3”：I=1=>GotoA⊿I=2=>GotoC⊿I=3=>GotoB⊿↙LbI A:{KZ}:QXYFLPROKZ“―0+”:Prog“A”:G=G-90:G<0=>G=G+360:G>36 0=>G=G-360⊿⊿“XYF=”:U:Pause0:V: Pause0:G→DMS▲(Q+J)=K=>Goto0: ≠> GotoA⊿↙LbI 0:Q=K:X=U:Y=V:F=G:P=R: “-----AC-----”: ▲ Prog“GS” ↙LbI B: {KZ}:K:Z“―0+”: Prog“SJK”: Prog“A”:C“C.X”:T“C.Y”:I=0:J=0:PoI(C-U.T-V: “FY=”:J=J+180→DMS I▲GotoB↙LbI C:{MN}:M“(X)”:N“(Y)”:K=0: C=0：I=0:J=0:PoI(M-X,N-Y:K=Q+Icos(F-J) ↙LbI 2:K=K-C: Prog“A”:PoI(U-M,V-N:T=G-90-J:C=IcosT:Int(1000C)=0=>Goto3:≠> Goto2⊿↙LbI 3:sinT<0=>A= -1:≠>A=1⊿“(K)=”:K+C: Pause0:“(―0+)=”:IA▲ GotoC↙二．文件名:A（内核）W=K-Q:D=OrW2(P-R)/2LPR:E=OrWP-1:A=.1184634425:B=.2393143352 :S=.046910077:H=.2307653449:Z[1]=F+SE+S2D:Z[2]=F+HE+H2D:Z[3]= F+.5E+.25D:Z[4]=F+(1-H)2D+E(1-H):Z[5]=F+(1-S)2D+E(1-S):G=F+E+D +90+(斜交) ↙U=X+ZcosG+W(AcosZ[1]+BcosZ[2]+.2844444444cosZ[3]+BcosZ[4]+Aco sZ[5]):V=Y+ZsinG+W(AsinZ[1]+BsinZ[2]+.2844444444sinZ[3]+BsinZ[4]+AsinZ[5])三.文件名:SJK（数据库）K≤下个起点桩号=>Q=起点桩号:X=起点（X）:Y=起点（Y）:F=前切线方位角:L=线型长度:P=起点（R）:R=终点（R）:O=转向:≠>K≤。

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。

CASIO4850全线坐标正反算程序

CASIO4850全线坐标正反算程序主程序：ZHANG LINGLbI 0：Deg：T“WJM”：C=20+40（T-1）：T=0：N“1.SR，2.JS”：N=2=>Goto A△M“JZ-XX”：O“QD-ZH”：U“X0”：V“Y0”：G“FWJ”：LbI B：{HRQP}：P“QD-R”：H“XX-CD”：R“ZD-R”：Q“）-1，——0，（+1”：W=H：Z=0：Prog“S1”：Prog“S4”：U=X：V=Y：O=O+H：G=F-90：P=R：“…SR…NEXT…”：Isz T：Dsz M：Goto B△LbI A：Fixm：{XYSZN}：T=0：N“1.JS-XY，2.XY-FS”：N=2=>Goto 2△S“JS-ZH=”：Z“D=”：LbI 1：S>Z[8T+C+1]= >Isz T：Goto 1△Prog“S3”：W=S-O：Prog“S1”：“JS-X=”：X：Pause 0：“JS-Y=”：Y▲Goto A：LbI 2：X“X…”：Y“Y…”：I=X：J=Y：LbI 3：Prog“S3”：W=Abs（（Y-V）Cos（G-90）-（X-U）Sin（G-90：W>H=>Isz T：Goto 3△Prog“S2”：“FS-ZH…”：S=O+W：Pause 0：“D…”：Z▲Goto A子程序S1A=0.1184634425：B=0.2393143352：Z[6]=0.2844444444：K=0.046910077：L=0.2307653449：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+0.5E+0.25D：Z[4]=G+(1-L)2D+E(1-L：Z[5]=G+(1-K)2D+E(1-K：F=G+E+D+90：X=U+ZCosF+W(ACosZ[1]+BCosZ[2]+Z[6]CosZ[3]+BCosZ[4]+ACosZ[5]：Y=V+ZSinF+W(ASinZ[1]+BSinZ[2]+Z[6]SinZ[3]+BSinZ[4]+ASinZ[5]子程序S2W=Abs（（Y-V）Cos（G-90）-（X-U）Sin（G-90：Z=0：LbI 4：Prog“S1”：L=G+E+D-90：Z=(J-Y) CosL-(I-X) SinL：AbsZ<E-6=>Goto 5：≠>W=W+Z：Goto4△LbI 5：Z=0：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=[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：[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编制的。

卡西欧4800和4850放样坐标计算程序

已知任意点的里程、半宽，计算该点平面坐标的程序前言1、本计算程序适用范围仅限于某ZH点至下一相邻ZH点间的线路范围。

2、使用本程序前需进行导线点坐标计算并求出各个曲线的曲线要素如第一切线方位角、转向角、曲线半径、缓和曲线长度、切线长度、曲线长度等。

3、使用过程中要注意长短链对里程传递的影响关系，否则将导致计算结果错误。

4、往计算器中输入程序符号时将⊿改为。

5、本程序仅供参考,实际使用时请先进行验证。

6、本程序适用的计算器型号为卡西欧4800/4850。

“wangDongping” (程序名称)Defm 13↙XYVWRLZTC：Fixm↙Z[4]=0.5×L－L3÷(240×R2)↙Z[2]=Z[4]+(R+L2÷24÷R)Abs tan(T÷2)：G=Z＋Ans：Fix3：G＂JDDK=＂◢Z[9]=R Abs Tπ÷180＋L：Z[9]＂QXL=＂◢Z[13]=Z＋Ans：Z[13]＂HZDK=＂◢D=C＋T：D＜0⇒D=D＋3600⊿D＞3600⇒D=D－3600⊿D◢Z[11]＂HZX＂=V＋Z[2]×(CosC＋CosD)◢Z[12]＂HZY＂=W＋Z[2]×(SinC＋SinD)◢Norm：Prog＂WEIXONG2＂(程序“WaangDongping”运行到此时，自动转到程序“WangDangping2”进行下一步运行)“WangDongping2”（程序名称）Lb1 0：{K,S,Q,E,F,N,U,H}：P=Q＋Abs(E－F)×N：K≤Z⇒Goto1⊿Z＜K⇒K≤Z＋L⇒Goto2⊿Z＋L＜K⇒K≤Z＋Z[9]÷2⇒Goto3⊿Z＋Z[9]÷2＜K⇒K＜Z[13]－L⇒Goto4⊿Z[13]－L≤K⇒K＜Z[13]⇒Goto5⊿ K≥Z[13]⇒Goto6⊿Goto0↙Lb1 1：Z[1]=V＋(Z－K)×cos(C＋1800)＋Pcos(C＋S)＋Ucos(C＋H)：Z[3]=W＋(Z －K)×sin(C＋1800)＋Psin(C＋S)+Usin(C+H)：Goto7⊿Goto0↙Lb1 2：Z[5]=K－Z－(K－Z)5÷(40R2L2)：Z[6]=（K－Z）3÷（6RL）：A=T÷AbsT：Z[7]=(K－Z)2×180÷(6πRL)：Z[1]=V＋√((Z[5])2＋(Z[6])2)×cos(C＋A×Z[7])＋Pcos(C＋3×A×Z[7]＋S)＋Ucos(C＋3×A×Z[7]＋H)：Z[3]=W＋√((Z[5])2＋(Z[6])2)×sin(C＋A×Z[7])＋Psin(C＋3×A×Z[7]＋S)＋Usin(C＋3×A×Z[7]＋H)：Goto7⊿Goto0↙Lb1 3：B=L×180÷(2πR)：Z[8]=B＋(K－Z－L)×180÷(πR)：Z[5]=Z[4]+R×sinZ[8]：Z[6]=R+L2÷(24R)－RcosZ[8]：Z[7]=tan-1(Z[6]÷Z[5])：Z[10]=(2×(K－Z)－L)×180÷(2πR)：A=T÷AbsT：Z[1]=V+ √((Z[5])2+(Z[6])2)×cos(C+A×Z[7])+Pcos(C+A×Z[10]+S)＋Ucos(C＋A×Z[10]＋H)：Z[3]=W+√((Z[5])2＋(Z[6])2)×sin(C+A×Z[7])+Psin(C+A×Z[10]+S)＋Usin(D＋A×Z[10]＋H)：Goto7⊿Goto0↙Lbl 4：B=L×180÷(2πR)：Z[8]=B+(Z[13]－K－L)×180(πR)：Z[5]=Z[4]+R×sinZ[8]：Z[6]=R+L2÷(24R)－RcosZ[8]：Z[7]=tan-1(Z[6]÷Z[5])：Z[10]=(2×(Z[13]－K)－L)×180÷(2πR)：A=－T÷AbsT：Z[1]=Z[11]+ √((Z[5])2+(Z[6])2)×cos(D+180+A×Z[7])+Pcos(D+A×Z[10]+S)＋Ucos(D＋A×Z[10]＋H)：Z[3]=Z[12]+√((Z[5])2+(Z[6])2)×sin(D+180+A×Z[7])+Psin(D+A ×Z[10]+S)＋Usin(D＋A×Z[10]＋H)：Goto7⊿Goto0↙Lbl 5：Z[5]=Z[13]－K－(Z[13]－K) 5÷(40R2L2)：Z[6]=(Z[13]－K) 3÷(6RL)：A=－T÷AbsT：Z[7]=(Z[13]－K)2×180÷(6πRL)：Z[1]=Z[11]+ √((Z[5])2+(Z[6])2)×cos(D+180+A×Z[7])+Pcos(D+3×A×Z[7]+S)+Ucos（D+3×A ×Z[7]+H）：Z[3]=Z[12]+√((Z[5])2+(Z[6])2)×sin(D+180+A×Z[7])+Psin(D+3×A×Z[7]+S)＋Usin(D+3×A×Z[7]+H)：Goto7⊿Goto0↙Lb1 6：Z[1]=Z[11]+(K－Z[13])×CosD+Pcos(D+S)＋Ucos(D+H)：Z[3]=Z[12]+(K －Z[13])×sinD+Psin(D+S)＋Usin(D+H)：Goto7⊿Goto0↙Lb1 7：M=Z[1]－X：O=Z[3]－Y：Pol(M,O)：J＜0 J=J+3600⊿Fix3：J◢I◢P◢E－F◢M+X◢O+Y◢Norm：Goto0↙坐标程序符号说明1、X，Y：测站坐标2、V，W：ZH点X,Y坐标3、R：曲线半径4、L：缓和曲线长度5、Z：ZH点里程6、T：曲线转向角（左负右正）7、C：第一切线方位角8、K：放样点里程9、Q：放样半宽（平行于S角非切线边在里程点左右侧的宽度，若已知实际半宽P，则Q=P；若放路基边桩则输入变坡点处的半宽P）10、S：放样点与其对应的中线点连线和该中线点切线的夹角（左负右正，即放样点在线路左侧S为负，同理在右侧为正）11、F：放样点地面标高12、N：路基边坡坡度13、E：变坡点设计标高14、P=Q+Abs（E－F）×N（Abs为绝对值），若P为已知，则Q的值直接输入P的值即可，E、F、N均为零。

CASIO4850程序

CASIO4850程序（一）坐标反算：（ZBFS）（1）A“XA”:B“YA”:C“XB”:D“YB”S=POL（（C-A）,（D-B））▲E=Frac（J÷360+1）×360“FWJ=”:E▲注：（1）计算XA？→第一点X坐标，YA? →第一点Y坐标，XB? →第二点X坐标，YB? →第一点Y坐标（2）结果r→两点间距离，J→象限角，FWJ→坐标方位角（二）直线坐标：（ZXZB）（1）Lbl 0：｛K,D｝:S=K-〖起点里程如:1188000〗：A=〖坐标方位角（十进制）如：163.3395411〗：X=〖起点X坐标〗+S×cos A:Y〖起点Y坐标〗+S×sin A:D“BZ”:D=1＝＞｛L,G｝:L“ZL”:A[1]=A-90:N=X+L×cosA[1]:E=Y+L×sinA[1]:G=“YL”:A[2]=A+90:X[2]=X+G×cosA[2]:Y[2]=Y+G×sinA[2]:“ZC”:“N=”N=N▲“E=”:E=E▲“YC”:“X[2]=”:X[2]=X[2] ▲“Y[2]=”:Y[2]=Y[2] ▲Goto 0:≠＝＞“X=”:X=X▲“Y=”:Y=Y▲Goto 0注：（1）计算K? →里程BZ? →任意数中桩坐标，1边桩坐标：ZL? →左边桩，RL?→右边桩（2）结果N=左X坐标，E＝左Y坐标，X[2]=右X坐标，Y[2]=右Y坐标注明：程序修改方位角163°20′22.35″=（十进制）163.3395411（三）曲线坐标｛缓和曲线，圆曲线｝（1）主程序（STAKEOUT）Prog “DA.SKIP”Lbl 1:｛H,Y｝H“（CE SHE）K”Y“（CE SHE）E”S=H-Q:Z=FYProg “P”Pol（V-K,W-N）T＜0＝＞T=T-360⊿U=T:K=IU“（CE SHE）X”=A+KcosT“（CE SHE）X=”:U=U▲K“（CE SHE）Y”=Ksin T+B“（CE SHE）Y=”:K=K▲Goto 1（2）子程序（3） 1.（DA.SKIP）DegC“DATA IN”C=0＝＞Fixm:≠＝＞Prog“DA.IN”△2.（DA.IN）R“（BAN JING）R”L“（HUAN HE）L”F“（QU XIAN）FX”Q“（ZH）K”X“（HZ）K”S=（X-Q）÷2:Z=0Prog “XY.BAS”Z[1]=V:Z[2]=W:Z[3]=KD“（CE ZHAN）K”O“（CE ZHAN）E”S=D-Q:Z=FOProg “P”M=V:N=WE“（HOU SHI）K”G“（HOU SHI）FW”S=E-Q:Z=OProg“P”Pol（V-M,W-N）Fixm:P=J3.（P）S≤X-Q-L＝＞Prog “XY.BAS”:≠＝＞S=X-Q-S:Prog“XY.BAS”:Pol（V-Z[1],W-Z[2]）:Fixm:V=Z[1]-Icos （2FZ[3]-J）:W=Z[2]-ISin （2FZ[3]-J）:｛E,U,G｝:K=2Z[3]-K△4.（XY.BAS）A=A:B=B:S≤0＝＞K=0:V=S:W=FZ:≠＝＞S≤L＝＞K=90S^2÷πRL：V=S-S^5÷40R^2 L^2-Zsin K:W=FS^3÷6RL-FS^7÷336R^3L^3＋FZcosK△△S＞L＝＞K=90（2S-L）÷ΠR:V=L-L^3÷40R^2-Rsin（90L÷πR）+（R-Z）sinK:W=FL^2÷6R-FL^4÷336R^3+FRcos（90L÷πR）+F（Z-R）cosK△注：计算：①DATA IN? →模式输0以外的任意数修改曲线要素要素，输入0继续上次的计算。

卡西欧4850计算器程序

Goto B
Lbl 3
Dsz F
Goto A
Q=G-Z[C+1]
Z=Z[C+2]+Z[C+4]Q
Goto B
Lbl B
子程序文件名:DL
J=0:I=0
Z[12]=0=>Goto Y⊿
Z[1]=-1=>Goto Y⊿
D=20+5Z[3]+5Z[4]+4Z[5]+4Z[6]:E=0:F=Z[12]:Goto A
A=Z[C+1]:{A}:A"JDN":Z[C+1]=A
E:Pause 0
A=Z[C+2]:{A}:A"JDE":Z[C+2]=A
E:Pause 0
A=Z[C+3]:{A}:A"JDR":Z[C+3]=A
E:Pause 0
A=Z[C+4]:{A}:A"JDA1":Z[C+4]=A
E:Pause 0
A=1=>Goto B⊿
J=J+I
Lbl B
G=G+J
Z[14]=J
G=-4=>D=20+5Z[3]+5Z[4]+4Z[5]+4Z[6]+2Z[12]:E=0:F=Z[13]: Prog "IN3": Goto B⊿
Z[9]=G
Prog "ZB"
D=20+5Z[3]+5Z[4]+4Z[5]:E=0:F=Z[6]-2
Prog "GC"

casio4850程序

CASIO4800、4850坐标正反计算程序[MC2008.01.29][/B][/fly][/color] 2008年1月27日更新，加入了坐标反算功能，坐标正算加入斜角计算功能。

程序特点：真正的全线贯通坐标正反计算！在曲线元要素输入时仅需要输入第一段全部曲线元要素，后面曲线元要素除起点半径、终点半径、曲线长、转向需输入外其他要素均从前一曲线按辛普森8等分计算得出，解决了主线坐标计算无法获得第二段及其以后曲线元起点参数的问题；辛普森公式任意等分，满足所有精度要求；全线曲线元数据一次性程序化输入，参数存储采用扩充变量数据库，无需修改程序内容；多功能采用单程序编程，避免频繁调用子程序，提高运算输度；傻瓜操作，特别适用初学者。

一、程序：MCW“1.ZS 2.FS 3.SZ”:W=1=>Z[2]=0:V=0：Goto 1 ΔW=2=> Goto 4ΔW=3=> O“KOU LING”:O≠123456=>O=0: “OUT” ◢Goto DΔO=0: V=0:Z[1]=0:Goto 0←┘Lbi 0←┘”N0.”：Z[1]+1 ◢←┘Z[1]=0=>{ABCREFGU}:A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G“LS”:U“G”:Z[Z[1]×8+3]=A:Z[Z[1]×8+4]=B:Z[Z[1]×8+5]=C:Z[Z[1]×8+6]=R-1:Z[Z[1]×8+7]=E-1:Z[Z[1]×8+8]=F: Z[Z[1]×8+9]=F+G: Z[Z[1]×8+10]=U: “NEXT”◢Isz Z[1]: Goto 0ΔZ[1]=1=>D=Z[9]:Z=0:Z[2]=0:GOTO 2ΔD=Z[(Z[1]-1)×8+9]:Z=0:Z[2]=Z[1]-1:GOTO 2←┘Lbi A←┘Z[Z[1]×8+3]=X:Z[Z[1]×8+4]=Y:Z[Z[1]×8+5]=J: Z[Z[1]×8+8]=D: {REGU}:R“Ro”：E“RN”: G“LS”:U“G”: Z[Z[1]×8+6]=R-1: Z[Z[1]×8+7]=E-1: Z[Z[1]×8+9]=D+G: Z[Z[1]×8+10]=U: “NEXT”◢Isz Z[1]: Goto 0←┘Lbi 1←┘{DZO }:D：Z：O“RJ”：Z[2]=0:Goto 2←┘Lbi 2←┘V≠1=>Z[2]>Z[1] =>GoToDΔΔD≤Z[Z[2]×8+9]=> A=Z[Z[2]×8+3]:B=Z[Z[2]×8+4]: C =Z[Z[2]×8+5]:R=Z[Z[2]×8+6]: E=Z[Z[2]×8+7]: F=Z[Z[2]×8+8]: G=Z[Z[2]×8+9]: U=Z[Z[2]×8+10]: Goto3ΔIsz Z[2]:Goto 2←┘Lbi 3←┘W=3 =>N=8：=≠>N=4ΔP＝U（E-R）÷Abs（G-F）：Q＝Abs（D-F）÷N：S=90Q÷π：J ＝C+(NPQ+2UR)NS：L=1←┘X＝A+Q÷6×(Cos C+Cos J +4∑（Cos (C+((L+0.5)PQ+2UR)×(L+.5)S),L,0,(N-1)）+2∑（Cos (C+((LPQ+2UR)LS,L,1,(N-1)))+ZCos（J+ O）←┘Y＝B+Q÷6×(Sin C+Sin J +4∑（Sin (C+((L+0.5)PQ+2UR)×(L+.5)S),L,0,(N-1)）+2∑（Sin (C+((LPQ+2UR)LS,L,1,(N-1)))+Z Sin（J+ O）:V=1=>Goto6ΔV=2=>Goto9ΔV=3=> GOTO CΔW=3=>GOTO AΔZ=0=>“X”:X:Pause 0: “Y” :Y◢ Goto 1ΔZ<0=>“XL”:X:Pause 0: “YL”:Y◢ Goto 1Δfx4850 ①Z>0=>“XR”:X:Pause 0: “YR”:Y ◢Goto 1 ←┘Z=0=> X “X” ◢Y “Y”◢ Goto 1ΔZ<0=> X “XL” ◢Y “YL”◢ Goto 1Δfx4800 ②Z>0=> X “XR” ◢Y “YR”◢Goto 1 ←┘Lbi 4←┘{MH} ：M“X” ：H“Y” ：Z[2]=0：GOTO 5←┘Lbi 5←┘V=1：D= Z[Z[2]×8+9]：Z=0 ：O=90：GOTO 2←┘Lbi 6←┘K=（（H -B）Cos（C-90）-（M-A）Sin（C-90））×（（H -Y）Cos（J-90）-（M-X）Sin （J-90））：K≤0=> Goto 7ΔIsz Z[2]：Goto5←┘Lbi 7←┘D=F+Abs((H -B)Cos(C-90)-(M-A)Sin（C-90）)：D>G=> Isz Z[2]: Goto5ΔGoto 8←┘Lbi 8←┘V=2 ：GOTO 3←┘Lbi 9 ←┘K=(H -Y)Cos(J-90)-(M-X)Sin(J-90)：AbsK<0.001=>Goto BΔD=D+K ：GOTO 8←┘Lbi B←┘V=3 ：Z=0 ：Goto 3←┘Lbi C←┘Z=(H-Y) ÷Sin(J+90)：“D”:D:Pause 0: “Z” : Z◢(Z=(H-Y) ÷Sin(J+90)：D“D” ◢“Z” Z◢)4800输出GOTO 4←┘Lbi D←┘二、说明a、编制说明本程序是运用复化辛普生公式根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距，对该曲线段范围内任意里程中边桩坐标进行计算，以及对卡西欧扩充变量的灵活应用，实现了真正意义上的的全线贯通及曲线要素输入程序化（在不修改程序内容的情况下可通过运行程序输入任意多段曲线元要素）。

casio4850全线坐标计算程序

该程序是本人自己编的适用于casio4850全线坐标计算，发布上来给大家参考，请大家试用指教，如有问题请修改，首先是方位角显示不出来，请高手修改。

一、LXZB（坐标计算主程序）LbI Z：｛P｝A=×××：V=××°××′××″：C=×××：D=×××：P≤直线的终点桩号＝＞Prog“ZX.SUB”：Goto Z（以上各字母含义详直线坐标计算子程序说明）A=×××：V=××°××′××″：C=×××：D=×××：R=×××：I=（1或-1）P≤圆曲线终点桩号＝＞Prog“YQX.SUB”：Goto Z（以上各字母含义详圆曲线坐标计算子程序说明）A=×××：V=××°××′××″：C=×××：D=×××：B=×××：W=××°××′××″：E=×××：F=×××：R=×××：L=×××：I=（1或-1）P≤HY点桩号＝＞Prog“HQ1.SUB”：Goto ZP≤YH点桩号＝＞Prog“HYQ.SUB”：Goto ZP≤HZ点桩号＝＞Prog“HQ2.SUB”：Goto Z（以上各字母含义详（HQ1.SUB、HYQ.SUB、HQ2.SUB）坐标计算子程序说明）该主程序是简单的示例说明，具体的应根据图纸的线型组合，随机应变的做顺序调整，并将其输入Casio4850计算器。

CASIO fx-4850计算器测量坐标程序(1)

直线坐标程序（prog“ZXY”）W=W+1μX=J+（K-D）×cosG+Abs W×cos（G+W×90÷Abs W）Y=O+（K-D）×sin G+Abs W×sin（G+W×90÷Abs W）A=GFix 3V=1→“X=”：X▲“Y=”：Y▲Fix 4：“FWA=”：A→DMS▲△说明：1、D、K分别为线元起点中桩桩号与所求点对应的中桩桩号。

2、J、O、G分别为线元起点X坐标，Y坐标与切线方位角。

3、W为偏距（左偏为“-”右偏为“+”）——计算点沿横断面方向到中线点的水平距离。

4、V坐标正算与反算的识别符。

圆曲线线元计算坐标程序（prog“YXY”）B=（（K-D）÷Abs R）r：W= W+1μX=J+Abs R(sinB×cosG+(1-cosB) ×cos (G+R×90÷Abs R))+ Abs W×cos (G+R ×B÷Abs R+W×90÷Abs W)Y=O+ Abs R(sinB×sin G+(1-cosB) ×sin (G+R×90÷Abs R))+ Abs W×sin (G+R ×B÷Abs R+W×90÷Abs W)A=G+R×B÷Abs RFix 3V=1→“X=”：X▲“Y=”：Y▲Fix 4：“FWA=”：A→DMS▲△说明：1、D、K分别为线元起点中桩桩号与所求点对应的中桩桩号。

2、R为线元半径值，左偏为“-”右偏为“+”。

3、J、O、G分别为线元起点X坐标，Y坐标与切线方位角。

4、W为偏距（左偏为“-”右偏为“+”）——计算点沿横断面方向到中线点的水平距离。

缓和曲线计算坐标程序（prog“HXY”）S=Int (B2÷Abs R×1000+0.5) ÷1000：W=W+1μC=S-Abs (E-D)：F=CG=G-R÷Abs R×（C2÷2÷B2）rP rog “HAB”：Q=M：P=NC=Abs (K-D)+F：Prog “HAB”X=J+（M-Q）×cos G+（N-P）×cos (G+R×90÷Abs R) + Abs W×cos (G+R÷Abs R×（C2÷2÷B2）r+W×90÷Abs W)Y=O+（M-Q）×sin G+（N-P）×sin (G+R×90÷Abs R) + Abs W×sin (G+R÷Abs R×（C2÷2÷B2）r+W×90÷Abs W)A=G+R÷Abs R×（C2÷2÷B2）rFix 3说明：1、D、E、K分别为线元起点、终点桩号与所求点对应的中桩桩号。

卡西欧4800 4850计算器程序

卡西欧4800\4850计算器计算程序及说明一：坐标反算：(ZBFS)1 程序步骤：C”X1”:D”Y1”:E”X2”:F”Y2”:Fixm:Pol(E-C,F-D:I”S1-2=”▲J≤Ｏ＝＞J=J+360△J”A1-2=”2 操作过程：ZBFS→EXE→输入X1值→EXE→输入Y1值→EXE→输入X2值→EXE→输入Y2值→EXE→EXE→EXE注：①：S1-2:计算得出的距离。

②：A1-2:计算得出的角度。

（按shift°′″转换为60进制的角度）③：此程序可循环计算。

3 案例：已知：D1坐标（7811.23，606.136），D2坐标（7805.915，431.910）求解： a D1-D2(A1-2)=268°15′09.56 ″S D1-D2(S1-2)=174.3071二：坐标正算：(ZBZS)1 程序步骤：C”X1”:D”Y1”:S”S1-2”:R”A1-2”:Fixm:X”X2”=C+Rec(S,R) ▲Y”Y2”=D+J2 操作过程：ZBZS→EXE→输入X1值→EXE→输入Y1值→EXE→输入S1-2距离值→EXE→输入A1-2角度（例268°15′09.56 ″）值→EXE→EXE注：①：X2计算得出坐标，Y2计算得出坐标。

②：此程序可循环计算。

3 案例：已知：D1坐标（7811.23，606.136），两点距离：S D1-D2(S1-2)=174.3071方位角：a D1-D2(A1-2)=268°15′09.56 ″ 求解：D2坐标（7805.915，431.910三：线路曲线计算程序:（XLQXJS）(万能公式)1 程序步骤：1.Lbl 0:{E G}：A“XA”：B“YA”：C“CA”：D“1/RA”：E“1/RB”：F“DKA”： G“DKA”2.Lbl 1:{H O R}:H“DKI”：O“DL”：R“DR”：H＞G=﹥Goto 3⊿3.P=(E-D)/Abs(G-F):Q=Abs(H-F):I=P×Q:T=D+I4.J=C+(I+2D)×Q×90/π ▲5.M=C+(I/4+2D)Q×45/(2π):N=C+(3I/4+2D)Q×135/(2π):6.K=C+(I/2+2D)Q×45/π7.X=A+Q(cosC+4(cosM+cosN)+2cosK+cosJ)/12 ▲8.Y=B+Q(sinC+4(sinM+sinN)+2sinK+sinJ)/12 ▲9.U“XL”=X+Ocos(J-90) ▲10.V“YL”=Y+Osin(J-90) ▲11.W“XR”=X+Rcos(J+90) ▲12.Z“YR”=Y+ Rsin(J+90) ▲13. A=X：B=Y：C=J：D=T：F=H：Goto 1Lbl 3: A=X：B=Y：D=E：F=G：C=J：Goto 0程序说明：A- 曲线元起点A的坐标；B- 曲线元起点B的坐标；C- 曲线元起点A的切线坐标方位角；F- 曲线元起点A的里程；G- 曲线元起点B的里程；H- 曲线上待求点i的里程；D- 曲线元起点A的曲率；E- 曲线元终点B的曲率；XL-左边线点位X坐标；YL-左边线点位Y坐标；XR-右边线点位X坐标；YR-右边线点位Y坐标；X- 中线点位纵坐标；Y- 中线点位横坐标；DL-左边线距中线平距；DR-右边线距中线平距；该程序需要输入的数据为：（1）.曲线元起点A的坐标及切线坐标方位角，计算器上用“XA”，“YA”，“CA”显示；（2）.曲线元起点A和B的曲率，计算器上用I÷RA，I÷RB 显示（曲线左偏时取“-”）；（3）.曲线元起点A和终点B的里程，计算器上用“DKA”，“DKB”显示；（4）.输入待求点里程和该点距左右的水平距离，计算器上用“DKI”，“DL”，“DR”显示；每算完一个待求点的中线及边线坐标，程序又让输入下一点的“DKI”，“DL”，“DR”当输入的“DKI”大于“DKB”时程序中显示“1÷RB”和“DKB”，此时输入下一个曲线元起点的曲率和里程，然后重复步骤（4），即可计算下一个曲线中线及边线点位坐标。

CASIO4850全线坐标正反算程序(辛普生公式)

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 BLbl A：Fixm:{XYSZN}:T=Ø:N"1.JS-XY,2.XY-FS": N=2=>Goto2S"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 ALbl 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=>Goto 5:≠=>W=W+Z:Goto 4Lbl 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. 本程序可根据多个线形数据计算平曲线任意点和边桩的坐标。

卡西欧4850程序计算简单绝对不是糊弄

主程序：<TYZDJS>“1。

（小数点）SZ=>XY”：“2。

XY=>SZ”：M：O“ZHK”：U“X0”（注意零和字母O的区别）：V“Y0”：G“QFWJ”：H“LS”：P“RA”：R“RB”：Q：D=(P-R)/(2HPR)：M=1=>Goto （符号）1：≠=> Goto2⊿Lb1 1：〔SZ〕（这个是大括号）:S“K”：Z:W=Abs（绝对值符号）（S-O）：Prog“SUB1”：“XS=”：X=X◣“YS=”：Y=Y◣“FS=”：F=F-90◣C“X1”:E“Y1”：I=0：J=0：Pol（X-C，Y-E）：“JULI=”:I=I◣J〈0=>J=J+360⊿“FWJ=”:J◣Goto1Lb12：〔XY〕（这个是大括号）：XY:I=X：J=Y：Prog“SUB2”：“K=”：S=O+W◣“B=”Z=Z◣Goto2 ⊿子程序<SUB1>A=011739274226:B=0.3260725774:K=0.0694318442:L=0.3300094782:F=1-L:M=1-KX=U+W(Acos(G+57.2958QKW(1/P+KWD))+Bcos(G+57.2958QLW(1/P+LWD))+Bcos(G+57.2958QFW (1/P+FWD))+Acos(G+57.2958QMW(1/P+MWD))):Y=V+W(Asin(57.2958QKW(1/P+KWD))+Bsin(G+ 57.2958QLW(1/P+LWD))+Bsin(G+57.2958QFW(1/P+FWD))+Asin(G+57.2958QMW(1/P+MWD))):F =G+57.2958QW(1/P+WD)+N“N=90”:X=X+ZcosF:Y=Y+ZsinF子程序<SUB2>T=G-90:W=Abs（绝对值符号）((Y-V)cosT-(X-U)sinT):Z=0: Lb10: Prog“SUB1”：L=T+57.2958QW(1/P+WD):Z=(J-Y)cosL-(I-X)sinL: Abs（绝对值符号）Z<1E-6=> Goto1：≠=>W=W+Z：Goto 0◣Lb1 1:Z=0: Prog“SUB1”：Z=(J-Y)/sinF子程序<QXZDJS>“1=>SJK1”：“2=>SJK2”:M=1=> Prog“SJK1”：Prog“TYQXJS”⊿M=2=> Prog“SJK2”◣M=n=> Prog“SJKn”◣注：⊿——按 Function 键 --3----3◣ ---- SHIFT +X2计算说明：运行主程序，屏幕显示M?输入1 正算 2 反算屏幕显示ZHK? 起点桩号屏幕显示X0? 起点X坐标屏幕显示Y0? 起点Y坐标屏幕显示QFWJ? 输入起点方位角屏幕显示LS？曲线长屏幕显示RA? 曲线起点的曲率（即起点半径）屏幕显示RB? 曲线终点的曲率（即终点半径）屏幕显示Q? 左偏-1 右偏1屏幕显示K? 带求点的里程桩号？屏幕显示Z ? 偏矩如果有不懂的+QQ273456237。