CASIO fx4850万能坐标计算程序

合集下载

高速MC公路测量CASIO4800&4850万能坐标计算程序(完整版)

高速公路测量CASIO4800&4850万能坐标计算程序(完整版）程序特点：真正的全线贯通坐标正反计算、任意斜角计算！！！程序中加入测站点，真正的实现了“坐标法”与“极坐标法”两种放样方法的同时显示的功能，使得放样操作方法选择时更加灵活！！！在曲线元要素输入时仅需要输入第一段全部曲线元要素，后面曲线元要素除起点半径、终点半径、曲线长、转向需输入外其他要素均从前一曲线按辛普森8等分计算得出，解决了主线坐标计算无法获得第二段及其以后曲线元起点参数的问题；辛普森公式任意等分，满足所有精度要求；全线曲线元数据一次性程序化输入，参数存储采用扩充变量数据库，无需修改程序内容；多功能采用单程序编程，避免频繁调用子程序，提高运算速度。

一、程序：ZBJSW“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 CΔO=0: V=0:Z[1]=0:Goto 0←┘Lbi 0←┘”N0.”：Z[1]+1 ◢Z[1]=0=>{ABCREFGUKO}:A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G “LS”:U“G”:K“X（00）”: O“Y（00）”: 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“R0”：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←┘{DZT }:D：Z：T“RJ”：Z[2]=0:Goto 2←┘Lbi 2←┘V≠1=>Z[2]>Z[1] =>GoToCΔΔ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：≠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+0.5)S),L,0,(N-1)）+2∑（Cos (C+((LPQ+2UR)LS,L,1,(N-1)))+ZCos（J+ T）←┘Y＝B+Q÷6×(Sin C+Sin J +4∑（Sin (C+((L+0.5)PQ+2UR)×(L+0.5)S),L,0,(N-1)）+2∑（Sin (C+((LPQ+2UR)LS,L,1,(N-1)))+Z Sin（J+T）:V=1=>Goto6ΔV=2=>Goto9ΔV=3=> GOTO CΔW=3=>GOTO AΔZ=0=>“X（Z）=”:X:Pause 0: “Y（Z）=”:Y◢Pol（（X－K），（Y－O））←┘“S（Z）=”:I ◢J<0=> J=J+360Δ“F（Z）=”: J→DMS◢Goto 1ΔZ<0=>“X（L）=”:X:Pause 0: “Y（L）=”:Y◢Pol（（X－K），（Y－O））←┘fx4850①“S（L）=”:I ◢J<0=> J=J+360Δ“F（L）=”: J→DMS◢Goto 1ΔZ>0=>“X（R）=”:X:Pause 0: “Y（R）=”:Y ◢Pol（（X－K），（Y－O））←┘“S（R）=”:J ◢J<0=> J=J+360Δ“F（R）=”: J→DMS◢Goto 1 ←┘Z=0=> X “X（Z）=”◢Y “Y（Z）=”◢Pol（（X－K），（Y－O））←┘I“S（Z）=”◢J<0=> J=J+360ΔJ“F（Z）=”◢Goto 1ΔZ<0=> X “X（L）=”◢Y “Y（L）”◢Pol（（X－K），（Y－O））←┘fx4800②I“S（L）=”◢J<0=> J=J+360ΔJ“F（L）=”◢Goto 1ΔZ>0=> X “X（R）=”◢Y “Y（R）=”◢Pol（（X－K），（Y－O））←┘I“S（R）=”◢J<0=> J=J+360ΔJ“F（R）=”◢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：T=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)：Abs K<（1÷E）^3=>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◢4850输出(Z=(H-Y) ÷Sin(J+90)：D“D”◢Z “Z”◢4800输出)GOTO 4←┘Lbi C←┘二、说明a、编制说明本程序是运用复化辛普生公式根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距，对该曲线段范围内任意里程中边桩坐标进行计算，以及对卡西欧扩充变量的灵活应用，实现了真正意义上的的全线贯通及曲线要素输入程序化（在不修改程序内容的情况下可通过运行程序输入任意多段曲线元要素）。

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”说明：计算时只需进入子程序（运算程序），无需进入主程序。

CASIOFX4850公路计算程序精品文档17页

0前言：传统公路测量中，使用的仪器设备和方法都很落后，需带着数学用表、曲线用表、计算盘、计算尺和算盘等一类的工具，完成外业测量工作。

计算器的出现，改变了这一局面。

高速公路建设中，长大曲线比比皆是，传统中对公路中线的测设方法，被极坐标法彻底的否定与取代，但大量的计算工作，只能带着提前计算好的线路逐桩坐标、高程资料，进行外业测量工作，机动性很差，现场查找也不方便。

这些问题都能在CASIO系列可编程计算器上得到很好的解决，对C ASIO系列可编程计算器如何使用，直接影响到测量成果的质量和工作效率，本文将对CASIO系列可编程计算器快捷的计算方法进行分析与介绍。

1：以知线外任意点坐标，求对应线路里程在缓和曲线上，要计算任意里程的法线方向及任意宽度的边线坐标，非常简单。

但要计算任意一个已知坐标点，是对应哪一个里程法线方向上的点，就有一些困难。

很难推导一个这样的计算公式。

唯一的方法“渐进”，如果手工计算这可不是一个好方法。

但在有CASIO系列可编程计算器，如：FX-4500的情况下就变的非常简单了。

亦可用于直线和圆曲线的计算。

首先在缓和曲线上任选一点A为起始点，计算该点的坐标和切线方位角，通过坐标反算求起始点A与计算点B的方位角和距离，B点肯定对应A点切线方向上有一个垂足C点，把三点看成一个直角三角形，通过解直角三角形计算AC的距离，当该距离大于某一数值，如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 中线法线距离。

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反算，其它重选。

卡西欧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-fx-4850P-4800p-5800P-9750GⅡ程序运行示例及说明

程序运行示例及说明运行主程序“LQXJS”一、任意点坐标计算（道路、桥、涵）第一步：Z? 1(进入道路、桥、涵任意点坐标计算).第二步：A? 742589.425(输入起点的“X”坐标).第三步：B? 463404.387(输入起点的“Y”坐标).第四步：C? 15°23′31.7″(输入起点的方位角).第五步：D?1÷340 (输入起点曲率).第六步：E? 1÷2286.5(输入终点曲率).第七步：F? 714.188(输入起点桩号).第八步：G? 890.019 (输入终点桩号)第九步：H? 800 (输入待求点桩号)第十步：S? 0((输入待求点至中桩的距离(中心桩为“0”、左为“-”右为“+”) 第十一步：U=26°50′56.17″ (待求点所计算的切线方位角).第十二进步：X=742669.0657 (显示待求点所计算的“X”坐标).Y=463435.9536(显示待求点所计算的“Y”坐标).第十三步：Q? 1(输入“1”时返回第九步，再重复步骤第十步~十三步计算各点的道路坐标；当输入“0”时进入桥、涵坐标计算,再继续下面步骤).第十四步：P? 135(输入桥、涵中心点处纵向轴线(路线前进方向)与横向轴线(从左到右方向)之间的夹角.如图).第十五步：K? 0(输入桥、涵中心点处纵向轴线(路线前进方向)偏离切线方向的偏角(左“-”右“+”)如图). 第 107頁第十六步：J? 8(输入桥、涵的垂直跨距(桥、涵各点至中心处横向轴线的垂直距离的2倍)).第十七步：W? 20 (输入桥、涵上各点至纵向轴线(路线前进方向)的垂直距离(左“-”右“+”)如图). 第十八步： M? 1(计算以桥涵中心处横轴线为界在前面各点的坐标;(M=1为前面、M=-1为后面、M=0 为横轴线)如图).第十九步：L=742647.2359(显示待求点所计算的“X”坐标).N=463447.3197(显示待求点所计算的“Y”坐标).第二十步：返回第十六步,再根据需要计算各点坐标.第二十一步：J、W同时输入“0”则可以退出桥、涵坐标计算，并返回到程序第九步状态.二、坐标反算桩号、距离第一步：Z? 2(进入坐标反算桩号、距离).第二步：A? 742589.425(输入起点的“X”坐标).第三步：B? 463404.387(输入起点的“Y”坐标).第四步：C? 15°23′31.7″(输入起点的方位角).第五步：D?1÷340 (输入起点曲率).第六步：E? 1÷2286.5(输入终点曲率).第七步：F? 714.188(输入起点桩号).第八步：G? 890.019 (输入终点桩号)第九步：J? 742669.0657(输入已知点X坐标).第十步：P? 463435.9536(输入已知点Y坐标).第十一步：H=800(显示待求点所计算的里程桩号)..第十二步：S=0 (显示待求点所计算的距离左为“-”、右为“+”)第十三步：返回第九步.三．曲率判断方法1．直线段计算时： D=0：E=02．圆曲线段计算时：D=1÷R：E=1÷R（左“-”，右“+”）3．直缓段计算时： D=0：E=1÷R（左“-”，右“+”）4．缓直段计算时： D=1÷R：E=0（左“-”，右“+”）5．非完整型缓和曲线段计算时：D=1÷R1：E=1÷R2（左“-”，右“+”）H H（X、Y）N=1 时U ZCR纵向轴线横向轴线KPM=1M=-1M=OJS（J、P）N=2 时H UUw(L、N)起点（A，B）FD G E。

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编制的。

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编制的。

4850计算器坐标放样程序

一般直线，曲线任意桩号放样程序A7QX(文件名)LbI3:Fixm:{L}:L≤53459.53=>Prog”7”:GotoAΔLbIA:Prog”2”:T≤0=>Goto7ΔT>E=>Goto0ΔProg”3”:Prog”4”:LbI1:{SU}:S=0=>Goto3ΔRec(S,Z+QF+U):Prog”5”:Goto1:LbI0:T>H-P-V=>Goto6Δprog”6”:Prog”4”:LbI2:{SU}:S=0=> Goto3ΔRec(S,Z+OF+U):Prog”5”:Goto2:LbI7:Rec(-T,Z+180):Prog”4”: LbI4:{SU}:S=0=> Goto3ΔRec(S,Z+U:Prog”5”: Goto4: LbI6:I=W:G=Z+180(H-P-(V+E)÷2)÷兀RF:Z=G+180: Rec(I,G:K=I+C:G=J+D:T ≥H-P=> Goto5ΔT=H-L:F=-F:E=V:Prog”3”:F=-F:N=K+I▲M=G+J▲LbI8:{SU}:S=0=>Goto3ΔRec(S,Z+180-QF+U:Prog”5”:Goto8: LbI5 :Rec(T-H+P,Z+180):N=K+I▲M=G+J▲LbI9: {SU}:S=0=>Goto3ΔRec(S,Z+180+U):Prog”5”:Goto92(文件名) T=L-P:Pol(C-A,D-B:Z=J3(文件名) W=ER:X=T-T^5÷40W^2+T^9÷3456W^4:Y=T^3÷6W-T^7÷336W^3: Pol(X,Y: Rec(I,Z+JF:Q=90T^2÷R÷E÷兀4(文件名) N=A+I▲M=B+J▲5(文件名) X=N+I▲Y=M+J▲6(文件名) V=E÷2-E^3÷240R^2:G=E^2÷24R-E^4÷2688R^3:O=180(T-0.5E)÷R兀:X=V+RsinO:Y=G+(1-CosO)R:Pol(X,Y:Rec(I,Z+JF7（文件名） (数据库)L=任意桩号 S=居中距离 U＝左-90，右+90 A＝ZH(ZY)点X坐标B= ZH(ZY)点Y坐标 C=交点X坐标 D＝交点Y坐标 E＝LS1 V＝LS2 R=半径 H=YH(YZ)点桩号 P=ZH(ZY)点桩号 F=左转-1，右转+1 W=切线长如果两弯道间没有直线段：L≤ZY(YZ)点桩号 =>Prog”数据库文件名”如果有直线段：L≤下一弯道ZY点桩号 =>Prog”数据库文件名”一般直线，曲线任意桩号放样程序B HYQPXG（文件名）Prog”1”:LbI3:Fixm:{L}:Prog”2”:T≤0=>Goto7ΔT>E=>Goto0ΔProg”3”:Prog”4”:LbI1:{SU}:S=0=>Goto3ΔRec(S,Z+QF+U):Prog”5”:Goto1:LbI0:T>H-P-E=>Goto6Δprog”6”:Prog”4”:LbI2:{SU}:S=0=> Goto3ΔRec(S,Z+OF+U:Prog”5”:Goto2:LbI7:Rec(-T,Z+180:Prog”4”: LbI4:{SU}:S=0=> Goto3ΔRec(S,Z+U:Prog”5”: Goto4:LbI6: G=Z+180(H-P-E)÷兀RF:Z=G+180: Rec(I,G:K=I+C:G=J+D:T≥H-P=> Goto5ΔT=H-L:F=-F: Prog”3”:F=-F:N=K+I▲M=G+J▲LbI8:{SU}:S=0=>Goto3ΔRec(S,Z+180-QF+U:Prog”5”:Goto8: LbI5: Rec(T-H+P,Z+180:N=K+I▲M=G+J▲LbI9: {SU}:S=0=>Goto3ΔRec(S,Z+180+U:Prog”5”:Goto91（文件名) A”X1”B”Y1”C”X2”D”Y2”E”LS”H”HZ”RFP:2(文件名) T=L-P:Pol(C-A,D-B:Z=J3(文件名) W=ER:X=T-T^5÷40W^2+T^9÷3456W^4:Y=T^3÷6W-T^7÷336W^3: Pol(X,Y: Rec(I,Z+JF:Q=90T^2÷R÷E÷兀4(文件名) N=A+I▲M=B+J▲5(文件名) X=N+I▲Y=M+J▲6(文件名) V=E÷2-E^3÷240R^2:G=E^2÷24R-E^4÷2688R^3:O=180(T-0.5E)÷R兀:X=V+RsinO:Y=G+(1-CosO)R:Pol(X,Y:Rec(I,Z+JFL=任意桩号 S=居中距离 U＝左-90，右+90 X1＝ZH(ZY)点X坐标Y1= ZH(ZY)点Y坐标 X2=交点X坐标 Y2＝交点Y坐标 E＝LS1 V＝LS2 R=半径 H=YH(YZ)点桩号 P=ZH(ZY)点桩号 F=左转-1，右转+1任意曲线坐标正反算CASIO-4850计算程序含数据库、测站放样计算（第五次修改-正式版）1.加编数据库及计算总调度程序,计算中不必逐项输入"线元要素",提高运算速度，避免现场忙中出错2.可加入多条线路的数据库,内业输入数据,外业一目了然3.计算直观,人性化4.正算直接输入里程和边距,反算输入近似里程便可5.增加了“计算点与测站点”的距离和方位角计算语句，方便直接放样6.愿收获与大家共享,同时也希望大家提出心得和宝贵建议7.核心计算程序摘自“yshf”一.改动后的程序清单，增设数据库程序,可输入n条线路的数据库1.(QXZDJS 计算总调度程序)M=1=>Prog "SJK1":Prog "SUBSJK"△←┘...........M=n=>Prog "SJKn":Prog "SUBSJK"△←┘2.(SJK1 数据库程序)"1.SZ => XY"："2.XY => SZ"：｛N,S｝：N：S"DKI"←┘S≤本线元终点里程＝＞U=本线元起点X坐标：V=本线元起点Y坐标：G=本线元起点正切线方位角：P=本线元起点曲率半径：R=本线元终点曲率半径：O=本线元起点里程：H=本线元终点里程：Q=本线元左直右偏向(0或1、-1):Goto0Δ←┘S≤下一线元终点里程=＞O=。

实用程序Casio4850 坐标方位角距离

这是本人自编的Casio-4850测量程序,实用!请测友提出保贵意见! 愿和测友共享!符号说明-------------------A ？交点桩号B ？转角(左-右+)C ？缓和曲线长D = Li (距ZH点长度)E ？X [起算点]F ？Y [起算点]G ？X [交点]H ？Y [交点]I = POL（）：r极距J = POL（）：θ极角K ？DK [桩号］L = 曲线长（ZH～HZ）M ？X[后视点]N ？Y[后视点]O= 方位角[ZH～JD]P = 方位角[JD～HZ]Q = Li(距HZ点长度）R ？半经S = HD(水平距离)T = T (切线长)U ？X[测站]V ？Y[测站]W = 测站至后视方位角X = X坐标Y = Y坐标Z = 测站至测点方位角Z[1]= β0Z[2]= mZ[3]= pZ[4]= [ZH]点X坐标Z[5]= [ZH]点Y坐标Z[6]= [HZ]点X坐标Z[7]= [HZ]点Y坐标Z[8]= 圆曲线所夹圆心角Z[9]= 待测点X坐标Z[10]= 待测点Y坐标=============================================================== ================程序-------------------------------------------------------------------------------DEG: Defm10：A“DK[JD]”：B“JD°°°Z[-]Y[+]”→DMS： R：C“L0”：E“X[QSD]”：F“Y[QSD]”：G“X[JD]”：H“Y[JD]”：U“X[CE ZHAN]”：V“Y[CE ZHAN]”：M“X[HOU SHI]”：N“Y[HOU SHI]”： Z[1]=90C/(πR)：Z[2]=C/2-C^3/(240R^2): Z[3]=C^2/24R-C^4/(2688R^3):T=（R+Z[3]）tg（AbsB/2）+Z[2]：L=πR（AbsB-2Z[1]）/180+2C：I=0：J=0： Pol（M-U，N-V）：J<0=>J= J+360：△：“FWJ[CZ～HS]=”：J→DMS▲LbI A｛K｝：K“DK”：I=0：J=0：Pol（G-E，H-F）：J＜0=>J=J+360：△：P=J+B： P＜0=>P=P+360：△：Z[4]=G+Tcos（J+180）：Z[5]=H+Tsin（J+180）：Z[6]=G+TcosP：Z[7]=H+Tsi nP：D=K-A+T：D＜0=>X=﹣D: Y=0：O=J+180：Goto 1：△D＜C =>X=D- D^5/(40R^2C^2) + D^9/(3456R^4C^4)-D^13/(599040R^6C^6): Y=D^3/(6RC) -D^7/(336R^3C^3) +D^11/(42240R^5C^5):B＜0=>Y=-Y:△: O=J: Goto1:△:D≤L-C=>Z[8]=180(D-0.5C)/(πR): X=RsinZ[8]+Z[2]: Y=R-RcosZ[8]+Z[3]: B＜0=>Y=-Y:△: O=J: Goto1: △:Q=K-A+T-L: Q＜0 => Q=-Q：X=Q- Q^5/(40R^2C^2) + Q^9/(3456R^4C^4)-Q^13/(599040R^6C^6):Y=Q^3/(6RC) -Q^7/(336R^3C^3) +Q^11/(42240R^5C^5):B＞0=> Y=–Y：△：P=P+180：Goto2：△：Q≥0=>X=Q：Y=0：P=P： Goto2：△LbI 1 :Z[9]=Z[4]+XcosO-YsinO: Z[10]=Z[5]+XsinO+YcosO: goto 3LbI 2 :Z[9]=Z[6]+XcosP-YsinP: Z[10]=Z[7]+XsinP+YcosP: Goto 3LbI 3 :I=0: J=0: Pol(Z[9]-U,Z[10]-V): J<0=>J=J+360:△“X=”:Z[9]▲“Y=”:Z[10]▲“FWJ[CZ～DK]°°°=”: J→DMS“HD=”:I▲Goto A=============================================================== =================注:①本程序用于CASIO fx-4850p型计算器。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

CASIO fx4850万能坐标计算程序[Machao 2006.10.23]
一、程序内容
程序：Mc
W“1.js 2.SZ”:W=1=>Goto 1 ΔW=2=> O“KOU LING”:O≠123456=>O=0“OUT”◢Goto 4 ΔO=0:
Defm 83:Z[81]=1:Goto 0←┘
Lbi 0←┘
{ABCDEFGU}:A“X0”:B“Y0”:C“F0’:D“R0”:E“RN’:F“D0”:G“LS”:U“G”←┘
Z[Z[81]]=A:Z[Z[81]+10]=B:Z[Z[81]+20]=C:Z[Z[81]+30]=1÷D:Z[Z[81]+40]=1÷E:Z[Z[81]+50]=F: Z[Z[81]+60]=F+G: Z[Z[81]+70]=U←┘
A=0=> Dsz Z[81]:Goto 1ΔIsz Z[81]:Goto 0←┘
Lbi 1←┘
{HS}:H“D’:S“Z”:Z[82]=1:Goto 2←┘
Lbi 2←┘
Z[81]>Z[82] =>”OUT”◢:Goto5←┘
H≤Z[Z[82]+60]=>A=Z[Z[82]]:B=Z[Z[82]+10]:Z[Z[82]+20]=C:D=Z[Z[82]+30]: E=Z[Z[82]+40]: F=Z[Z[81]+50]: G=Z[Z[81]+60]: U=Z[Z[81]+70]: Goto3:ΔIsz Z[82]:Goto 2
Lbi 3←┘
P＝U（E-D）÷Abs（G-F）：Q＝Abs（H-F）：I＝PQ：J＝C+90 Q（I+2UD）/π：J<0=>J=J+360ΔM＝C+45 Q（I÷4+2UD）÷2π：N＝C+135Q（3I÷4+2UD）÷2π：K＝C+45Q（I÷2+ 2UD）÷π←┘
S=0=>Z=0：Goto4ΔS＜0=>Z＝－1：Goto4ΔZ＝1：Goto4←┘
Lbi 4
X＝A+Q(Cos C+4（Cos M+Cos N）+2Cos K+Cos J) ÷12 +ZSCos（J+90Z) ←┘
Y＝B+Q(Sin C+4（Sin M+Sin N）+2Sin K+Sin J) ÷12 +ZSSin（J+90Z) ←┘
S=0=>“X”:X:Pause 0: “Y”:Y◢Goto 1Δ
S<0=>“XL”:X:Pause 0: “YL”:Y◢Goto 1Δ
S>0=>“XR”:X:Pause 0: “YR”:Y ◢Goto 1 ←┘
Lbi5←┘
二、变量及说明
X0：Y0：F0――――曲线元起点X、Y坐标及起点正切线方位角
R0：RN――――――曲线元起点及终点半径
D0：LS：Q―――曲线元起点桩号、路线长度及线路左右偏标志（左=-1,直线=0,右=1）D――――――――曲线元中待求点桩号
Z―――――――――计算边桩距中线平距，左边输入负值，右边输入正值，中桩输入0 X：Y
XL:YL
XR:YR ―――――――待求点中边桩桩号的X，Y坐标
扩充变量：Z[Z[81]]: Z[Z[81]+10]: Z[Z[81]+20]: Z[Z[81]+30]: Z[Z[81]+40]: Z[Z[81]+50]: Z[Z[81]+60]: Z[Z[81]+70]: 分别为各线元X0:Y0:F0;R0:RN:D0:LS:Q
(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，
Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

(2) 当所求点位于中线时，Z=0，坐标显示X Y；当位于中线左侧时，Z取负值,坐标显示XL YL,；当位于中线右侧时，Z取正值,坐标显示XR YR。

(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。

(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。

(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45
次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。

止点与直线相接时，曲率半
径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。

(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。

止点与直线相接时，曲率半径等
于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。

本程序可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距，对该曲线段范围内任意里程中边桩坐标进行正算。

本程序是对扩充变量的应用，实现了真正意义上的的全线贯通。

程序分为两部分：1为计算2为设置。

首先对内存变量进行扩充，增加了83个扩充变量，为Z[1]~Z[83]，为防止误操作在输入正确口令123456时方可进行设置。

本程序设计为十段线元要素，如有需要可以根据内存容量进行增加。

在执行程序中将各个线元要素按规律输入到扩充内存变量中，在坐标计算时首先判断其在哪一线元内，并把其线元要素调出进行计算。

在进行计算时直接输入桩号、宽度就可以进行全线坐标计算。

本程序也可在4800上运行，只需将输出部分进行修改。

四、示例
郑石高速路面9标主线线元要素及相关数据如下：
X0 Y0 F0 R0 RN D0 LS Q
744450.244 510454.431 220º41 º31.8 º7000 7000 116075.889 3386.293 1
742505.351 507722.676 248º24 º33.6 º1045 1045 119462.182 1980.825 0
741776.461 505880.832 248º24 º33.6 º8500 8500 121443.007 1517.999 1
741346.546 504427.086 258º38 º30.1 º1045 1045 122961.007 1985.977 0
桩号宽度X Y
117000 0 743791.308 509807.481
118000 -20 743164.306 509027.467
119000 20 742707.550 508137.933
119462.182 0 742505.349 507722.673
120000 -15.253 742293.266 507228.206
121000 15.260 741953.665 506287.142
121443.007 0 741776.461 505880.833
122000 -10.20 741578.902 505359.701
122961.007 50.355 741395.915 504417.169
123650 0 741210.853 503751.587。