修改4800正反算2

合集下载

卡西欧4800P计算器程序

环运算语句J 子程序名称H=X+DcosG◣ 公式运算、数据显示语I=Y+DsinG◣ 公式运算、数据显示语T=X—EcosG◣ 公式运算、数据显示语U=Y—EsinG◣ 公式运算、数据显示语最后计算器状态设定语句是大家最容易忽视的。

如果将单位进行预设那么计算器就会默认其使用单位在进行下一单位换算时要一定要进行单位转换，否者会使计算结果错误。

在显示屏幕的左下角可以清楚地发现小提示符号：如D代表度为现在的缺省单位、R代表弧度为现在的缺省单位、G 代表梯度为现在的缺省单位。

首页┊ GPS测量┊工程测量┊遥感测量┊地籍测量┊控制测量┊数字测图┊地图制图┊数据处理┊测绘仪器┊软件应用┊房产测量┊地理信息┊测量管理┊程序开发┊理论文章┊其他┊ 文章更新列表┊游客发布论文所属分类：程序开发来源：作者：更新日期：2006-8-25 阅读次数：156 文章更新列表查看公告列表00-54700”Lbl 1:｛ZKB｝:Z:K:B:O=1=＞Goto 2△O=2=＞Goto 3△O=3=＞Goto 4△Lbl 2:Z＜47264=＞Goto 5△Z＞501▲Goto 0Lbl 6:Prog“DQXY”:“X=”:X:Pause 0:“Y=”:Y▲Pol（X-Q,Y-U）:J＜0=＞J=360+J△:I=Intg（1000I+.5）Z、待求点的桩号。

K、待求点在左幅输0右幅输1。

B、待求点距中距离。

Goto 5前的桩号为起止点或断链点桩号。

P T÷100T-Int T）÷.6X=Q+Rec（L,T）:Y=U+J:“X=”:X:Pause 0:“Y=”:Y▲Lbl 1:Prog“55”:Prog“YS”:Z＞48980 oto 4△Prog“FQ”: Goto ZLbl 4:Prog“57-1”:Prog“YS”:Z＞50354.96=＞Goto 5△Prog“FQ”Z＞50167.16 =＞Z og“FQ”: Goto ZLbl 7:Prog“60”:Prog“YS”:Z＞52546.8=＞Goto 8△Prog“FQ”: Goto ZLbl 8:Prog“61”:Pro B△Prog“FQ”: Goto ZLbl B:“CX ERROR”▲ Goto 0Lbl Z:Z=Intg（1000Z+.5）÷1000:B=Intg（1000B+.5）÷100号为第2个曲线组合的ZH点桩号，以后类推。

隧道坐标正反算

隧道测量---快速坐标正反算程序(4800-4850计算器）2008-11-23 20:47:00点击：212ZB （坐标正算）1. Lb1 12. ｛L｝3.SMNARCL4. K=L-S5. R=0=>X“X”=M+KcosA▲6. Y“Y”=N+KsinA▲7. Z“FW”=A▲8.Goto 4△R≠0=>Abs K≤C=>Goto 2△R≠0=>K＞C=>Goto 3△R≠0=>K<0=>Abs K＞C=>Goto 1△9. Lb1 210. O=90K²÷（πRC）11. G=K³÷（6RC）-Kˆ7÷（336（RC）³＋Kˆ11÷（42240（RC）ˆ5-Kˆ15÷9676800（RC）ˆ7）＋Kˆ19÷（3530096640（RC）ˆ9）12. F=K－Kˆ5÷（（40（RC）²）－Kˆ9÷（3456（RC）ˆ4）＋Kˆ13÷（599040（RC）ˆ6）－Kˆ17÷（40320×4532（RC）ˆ8）13. J=0=>X“X”=M＋GsinA＋cosA▲14. Y“Y”=N＋FsinA- cosA▲15. Z“FW”=A-O▲Goto 4△16．J≠0=>X“X”=M＋F cosA－GsinA▲17．Y“Y”=N＋G cosA＋F sinA▲18．Z“FW”=A＋O▲Goto 4△19．LbI 320．B=90（2K-C）÷（лR）21.U=C÷2－C³÷（240R²）＋R sinB＋Cˆ5÷（34560Rˆ4）-Cˆ7÷（599040 Rˆ6）22．V=C²÷（24R）＋R（1－cosB）－Cˆ4÷（2688R³）＋Cˆ6÷（42240 Rˆ5）23．J=0=>X“X”=M＋UcosA＋VsinA▲Y“Y”=N＋U sinA－V cosA▲Z“FW”=A-B▲Goto 4△24．J≠0=>X“X”=M＋U cosA－VsinA▲Y“Y”=N＋U sinA＋V cosA▲Z“FW”=A＋B▲Goto 4△25．Lb1 426．｛D｝：｛T｝27．D=0 => Goto 1△28．D≠0=>X“LX”=X＋Dcos（T＋Z）▲Y“LY”=Y＋Dsin（T＋Z）▲Goto 1△SUB（坐反算）“1.SZ=>UV”: “2.UV =>SZ”: N（=2): X“X0”: Y“Y0”: O“S0”=L: G“F0”(=Z): Q：Q=0=>R=1E-45：≠> R△J=0：M=AbsQ：K=1-M：N=1=>Goto 1：≠>Goto 2△←┙Lbl 1：{SZ}：SZ：L=S-O：B=90QL÷π÷R：F=G+B：E=F+B+90：C=LK+2MRsinQB：U"US"= X+Rec(C,F)+ZcosE◢V "VS" =Y+J+ZsinE◢Goto 1←┙Lbl 2：W=G+90Q：E=X+Rec(R,W)：F=Y+J：W=W+180M：W>360=>W=W-360△←┙Lbl 3：{UV}：UV：D=Pol(U-E,V-F):T=J：J<0=>T=J+360△B=T-W：P=KB+MAbs B：S"S"=O+ DKcosP+ MπRP÷180◢H"H"=DKsinB+MQ(R-D)◢Goto 3←┙正算由下面输入：S： ZH或HZ点里程M： ZH或HZ点X坐标N： ZH或HZ点Y坐标A：切线方位角R：曲线半径C：缓和曲线长度L：求点里程X：求点X坐标Y：求点Y坐标FW：求点切线方位角D：左、右边桩距离T：边桩与线路夹角LX：边桩X坐标LY：边桩Y坐标J：右偏输“1”，左偏输“0”反算由下面输入：X0: 线路中心X坐标Y0: 线路中心Y坐标L: 相对应里程Z: 相对应方位角Q: ＝０U: 求点X坐标V: 求点Y坐标S= 求点里程H= 求点边距说明：１.可用于缓和曲线坐标正反算，反算时X0,Y0,L,Z利用正算储存的数据，速度快．２.反算时先输入近似里程正算．３.隧道断面炮位放样２０多个点位在操作熟练时只需15分钟左右，提高工程进度．。

4800计算器使用方法

4800L1 Defm 10L2 Lb1 0:{FG}:C"QX":D"QY":A"Q >":E"QDZH":F"ZDZH":H"QR":G"ZR":O"ZDX":U"ZDY" L2 Lb1 1:{K}:K>F＝＞Goto 3△L3 Lb1 2:P=A+(1/G+1/H)/Abs(F-E):L=Abs(K-E):J=P×LL4 I=A+90(J+2/H)L/π▲L5 Z[1]=A+45(J/4+2/H)L/(2π)L6 Z[2]=A+135(3J/4+2/H)L/(2π)L7 Z[3]=A+45(J/2+2/H)L/πL8 Z[4]=A+(J/8+2/H)L(90/8π)L9 Z[5]=A+(3J/8+2/H)L(3*90/8π)L10Z[6]=A+(5J/8+2/H)L(5*90/8π)L11Z[7]=A+(7J/8+2/H)L(7*90/8π)L12 {B}:X“X”=C+L(CosA+4(CosZ[4]+CosZ[5]+CosZ[6]+CosZ[7])+2(CosZ[1]+CosZ[2]+Cos Z[3])+CosI)/24+BCos(I+ 90)▲L13 Y“Y” =D+L(SinA+4(SinZ[4]+SinZ[5]+SinZ[6]+SinZ[7])+2(SinZ[1]+SinZ[2]+SinZ[3])+SinI)/24+BSin (I+90)▲L14 K=F＝＞Goto 4: ≠＞Goto 5△△L15 Lb1 3:K=F:Goto 2L16 Lb1 4:C=X:D=Y:E=F:H=G:A=I:Goto 0L17 Lb1 5:Pol(X-O,Y-U):W<0 ＝＞W=W+360 :Ｗ“FWJ＞”＝Int W+Int（100Frac W）/ 60+Frac （100W）/36▲Fix 3:V“S ”=V▲Norm: Goto 1△L18 ≠＞W =W:Ｗ“FWJ＞”＝Int W+Int（100Frac W）/ 60+Frac （100W）/36▲Fix 3:V“S ”=V▲Norm:Goto 1说明：本程序是根据复化辛普生公式（n=4）编制的，适合各种线型的正反算。

曲线(含直线)任意里程中边桩坐标正反算

曲线（含直线）任意里程中边桩坐标正反算（CASIO fx-4800p&fx-4850）J-SQX(竖曲线数据输入)“J-SQX”：{NHUVMQP}：Z[9]=N“SJD”：Z[10]=H“JDZ”：Z[11]=U“I1”：Z[12]=V“I2”：Z[13]=M“R”：Z[14]=Q“QD”：Z[15]=P“ZD”：“TO J-JSMS”J-PQX(平曲线数据输入,自动切换到J-JSMS)A“JD”B“JDX”C“JDY”F“FJ”O“ZJ：Z-，Y+”RE“LS1”K“LS2”：E＜1=＞E=1E-9⊿K＜1=＞K=1E-9⊿Z[1]=EE÷(24R)-E^4÷(2688RRR)：Z[2]=E÷2-EEE÷(240RR)：X=(EE-KK)÷(24R)÷sinAbsO：Z[3]“T1”=(R+Z[1])tan(AbsO÷2)+Z[2]-X▲Z[4]“T2”=(R+KK÷(24R)-K^4÷(2688RRR))tan(AbsO÷2)+K÷2-KKK÷(240RR)+X▲L=AbsOπR÷180+(E+K)÷2▲J=tan-1((R+Z[1])÷(Z[3]-Z[2]))：X“E”=(R+Z[1])÷sinJ-R▲X=A-Z[3]：Y=X+E：E＜1=＞X“ZY”▲≠=＞X“ZH”▲Y“HY”▲⊿Y“QZ”=X+(L-K-E)÷2+E▲Y=X+L-K：X=X+L：K＜1=＞X“YZ”▲≠=＞Y“YH”▲X“HZ”▲⊿Prog“J-JSMS”J-JSMS(放样模式主程序)“1-ZS,2-FS,4-DMFY” Lb1 0：{Z}：Z“MS”≤1=＞Goto 1：≠=＞Z“MS”≤4=＞Goto 2⊿⊿Lb1 1：{PDW}：PD“BZ”W“BJ”：Prog“JP”：X“X=”▲Y“Y=”▲Goto 0⊿Lb1 2：{XYW}：XYW“BJ”： Prog“JF”：P“P=”▲D“BZ=”▲Z=3=＞Prog“JS”：Prog“DMFY”⊿(运行竖曲线高程计算程序)Z=4=＞Z[26] “H”=12.417+(P-75360)*5.2/1000: Prog“DMFY”⊿（运行单面坡比高程计算,语句中12.417为起点桩号DK75+360的高程,可以根据实际情况进行调整；5.2/1000为单面上坡率,上坡输正值，下坡输负值）Goto 0JS(竖曲线计算主程序)P＜Z[14] =＞Prog“J-SQX”⊿P＞Z[15] =＞ Prog“J-SQX”⊿N=Z[9]：U=Z[11]：V=Z[12]：H=Z[13]：G=H Abs(U-V)÷200：P＜N=＞I=U：M=N-G：M＞P=＞M=P⊿≠=＞I=V：M=N+G：M＜P=＞M=P⊿⊿J=(P-M)2÷(2H)：U-V＞0=＞J=-J⊿H=Z[10]+(P-N)×I÷100+JZ[26] “H”=HJP(平曲线正算子程序)FixmLb1 1：J=F：X=B-Z[3]cosF：Y=C-Z[3]sinF：G=EP≤A-Z[3] =＞I=A-Z[3]-P：M=-I：N=0：H=F+W：Goto5≠=＞P≤A-Z[3]+E=＞I=P-A+Z[3]：H=90II÷(REπ)：O＜0=＞H=-H⊿H=H+W+F： Goto 3≠=＞P≤A-Z[3]+L-K=＞I=P-A+Z[3]-E： Goto 4：≠=＞Goto 2⊿⊿Lb1 2：X=B+Z[4]cos（F+O）：Y=C+Z[4]sin（F+O）：J=F+O+180：G=KP≤A-Z[3]+L=＞I= A-Z[3]+L-P：H=90II÷(REπ)：O＞0=＞H=-H⊿H=H+J+W+180： Goto 3≠=＞I=P-A+Z[3]-L：M=-I：N=0：H=J+W+180： Goto 5Lb1 3：M=I-I^5÷(40RRGG)：N=III÷(6GR)-I^7÷(336RRRGGG)： Goto 5Lb1 4：H=(E+2I)×90÷(πR)：M=RsinH+Z[2]：N=R(1-cosH)+Z[1]：O＜0=＞H=-H⊿H=J+W+H： Goto 5Lb1 5：P≤A-Z[3]+L-K=＞O＜0=＞N=-N⊿≠＞O＞0=＞N=-N⊿⊿Goto 6Lb1 6：X=X+Mcos J-Nsin J+Dcos H：Y=Y+Msin J+N cosJ+Dsin HJF(平曲线反算子程序)FixmU=X：V=Y：D=0：J=F-W：P=A+(Y-C)cos J-(X-B)sin JLb1 1：Prog“JP”：J=H-180：I=(V-Y)cos J-(U-X)sin J：Abs I＜1E-4=＞Goto 2：≠=＞P=P+I： Goto 1⊿Lb1 2：D=(V-Y)÷sin HDMFY隧道开挖断面(超欠挖情况)计算程序{HIJMNG}：Z[20]= H“C”：Z[21]=I“A”：Z[22]= J“H1”：Z[23]= G“GC”：Z[24]= M“R1”：Z[25]= N“R2”D≤2.3=＞D“ZD”=Abs(D-2.3)▲≠=＞D＞2.3=＞D“YD”= D-2.3▲⊿⊿(C为圆心O1与O2的高差值,A为圆心O1与O2的宽度差值,H1为圆心O1至设计高程的高度,GC为实测高程,R1为上部第1个半径,R2为下部第2个半径)(显示值ZD为至隧道中心线左边的宽度,YD为至隧道中心线右边的宽度,2.3为设计线路距隧道中心线的宽度,可以实际情况进行改变调整)Z[26]“H”▲(显示所求桩号点设计高程值)Z[27]“O1” =Z[26]+ Z[22](第1个圆心高程)Z[28]“O2”=Z[26]+ Z[22]+ Z[20] (第2个圆心高程)Z[29]= Z[24]×Sin 60▲(此数据为上半弧60度范围的宽度值，可根据实际情况调整；数值可不显示) D>Z[29] =＞Goto 1⊿Z[30]=√（（Z[24]）2-D2）Z[31]= √（Abs(（Z[24]）2-(Z[23]-Z[27])2)）Z[32] “HGD”=Z[27]+ Z[30]- Z[23] ▲Z[33] “SKD”=Z[31]-D ▲Prog“J-JSMS”Lb1 1Z[23]<Z[27] =＞Goto 2⊿Z[31]= √(（Z[25]）2-(Z[23]-Z[28])2)Z[33] “SKD”=Z[21]+ Z[31]-D ▲Prog“J-JSMS”Lb1 2Z[31]=Z[25]-（Z[27]-Z[23]）*0.35/(0.78+Z[20]+Z[22])(0.35为下边墙往边墙底内缩的距离，0.78为设计高程至边墙底部的高度)Z[33] “SKD”= Z[31]-D ▲Prog“J-JSMS”一、程序简介1、本套程序共有2个主程序，5个子程序。

匝道线元法

任意里程中边桩坐标正反算(CASIOfx-4800P计算器)程序时间:2007-12-18 13:13:15 来源:工测员网作者:未知1. 加编数据库作为主程序 , 计算中不必逐项输入 " 线元要素 ", 提高运算速度，避免现场忙中出错2. 将原来的主程序并入数据库3. 计算直观 , 人性化4. 正算直接输入里程和边距 , 反算输入近似里程便可5. 增加了“ 计算点与测站点” 的距离和方位角计算语句，方便直接放样6. 愿收获与大家共享7. 核心计算程序摘自“yshf”1. 正算子程序 (SUB1)A=0.1739274226 ： B=0.3260725774 ： K=0.0694318442 ： L=0.3300094782 ：F=1-L ：M=1-K ：X=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(G+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)+90 ： X=X+ZcosF ： Y=Y+ZsinF2. 反算子程序 (SUB2)T=G-90 ： W=Abs((Y-V)cosT-(X-U)sinT) ： Z=0 ： Lbl 0 ： Prog "SUB1" ：L=T+57.2958QW(1/P+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)÷sinF3、 . 增设数据库程序（ＳＪＫ主程序）Lb1 4 ： "1.SZ => XY" ： "2.XY => SZ" ：｛ NS ｝：S ∠ 下一线元起点里程＝＞O = 本线元起点里程： U= 本线元起点 X ： V= 本线元起点 Y ： G= 本线元起算方位角： H= 本线元长度： P= 起点曲率半径： R= 终点曲率半径： Q=0 或 1 、 -1 ： Goto0Δ←┘( 第一线元数据要素)S ∠ 下一线元起点里程＝＞O = 本线元起点里程： U= 本线元起点 X ： V= 本线元起点 Y ： G= 本线元起算方位角： H= 本线元长度： P= 起点曲率半径： R= 终点曲率半径： Q=0 或 1 、 -1 ： Goto0Δ←┘( 第二线元数据要素)S ∠ 下一线元起点里程＝＞O = 本线元起点里程： U= 本线元起点 X ： V= 本线元起点 Y ： G= 本线元起算方位角： H= 本线元长度： P= 起点曲率半径： R= 终点曲率半径： Q=0 或 1 、 -1 ： Goto0Δ←┘( 第三线元数据要素)。

4850 4500 4800正反算程序

1．中边桩计算程序<1>主程序名—XYLbl 1:{KM}：K：M：Prog”3”:X=C+Mcos(W-90:Y=D+Msin(W-90:”X=”:X:Pause 0:”Y=”:Y▲Goto1<2>子程序—3PrOg”QX”：Prog”20”<3>曲线要素—QX数据实例：某高速公路中线设计如下图所示K<2000=>B:E:F:O:P:Q:R: K≤4570=>B=25:E=1000:F=1000:O=100P=66:Q=4000:R=800: K≤5040=>B=-17:E=1086.4378:F=1555.1451:O=90:P=91:Q=4570:R=800: K≤5670=>B=35:E=1168.2556:F=2016.0909:O=120:P=74:Q=5040:R=1000: B=-25:E=1180.9927:F=2636.2334:O=150:P=109:Q=5670:R=1500注：无缓和曲线时，O=0.00001当输入里程小于2000时，为新曲线要素输入，可修改值域范围<4> 计算过程—20Z=O：Prog”8”：A=B/AbsB:L=ABπR/180+O:T=(D+RcosN)tan(AB/2)+C-RsinN：K≤Q-T=>C=(K-Q)cosP+E：D=(K-Q)sinP+F:W=P： K≤Q-T+O=>Z=K-Q+T：Prog”9”：K≤Q-T+L-O=>Z=O: Prog”9”:V=90(K-Q+T-O)/ πR:S=2R sinV:N=W+AV:C= C+S cosN:D=D+S sinN:W=W+2AV: Z=Q-T+L-K: Prog”8”:V=P+B-180-AV:W=P+B-AN :C=S cosV+E+T cos(P+B:D=S sinV+F+T sin(P+B: Z=K-Q+2T-L:W=P+B:C=E+ZcosW:D=F+Z sinW△△△△N=Z[3]<5> 缓和曲线方程式—8V=RO:C=Z-Z^5/40V2+Z^9/3456V^4:D=Z^3/6V-Z^7/336V^3+Z^11/42240V^5:N= 90Z^2/πV:V=tan-1(D/C:S=(D^2+C^2)^0.5<6> 子程序—9Prog”8”:C=Scos(P+AV)+E-TcosP:D=Ssin(P+AV)+F-TsinP:W=P+AN2.桩号、边桩距返算程序—KM方式一：试算法（速度慢，程序简单）Fixm:K=Q-T+L-O:Prog”20”Lbl 1:{XY}Lbl 2:POL(X-C,Y-D:K=Isin(J-W+90)+K:Abs sin(J-W+90)<0.0001＝>”K=”:K:Pause0:”M=”:M=Icos(J-W+90▲≠>Prog”20”:Goto2△方式二：解算法（程序复杂，计算速度快）Defm5: Lbl 1:Fixm:{XY}:Pol(X-E,Y-F:Z[1]=I:Z[2]=J:cos(J-P<0＝>K=Q+Icos(J-P: ≠>K=Q+Icos(J-P-B)+2T+L△K≤Q-T＝>M=Icos(J-P+90:Goto5≠>Goto2△△Lbl 2:Z=O:Porg”8”:U=(D+RcosN)+cos(B/2):Z=P+B/2+90A:C=UcosZ+E:D=UsinZ+F:Pol(X-C,Y-D:Asin(J-P-AN+90 A<0＝>K=K-AZ[1]cos(Z[2]-P+90)tan(90 (K-Q+T)^2/πRO):Goto3:≠>Asin(J-P-B+AN+90A>0 ＝>K=K-AZ[1]cos(Z[2]-P-B-90 )tan(90(K-Q+T-L)^2/πRO):Goto4:≠>K=Q-T+O+Asin-1(sin(J-P-AN+90 A))×πR/180 :M=A(I-R:Goto5△△Lbl 3:Z=K-Q+T:Prog”9”Pol(X-C,Y-D:K=K+Isin(J-W+90 :Abs sin(J-W+90<.001＝>M=Icos(J-W+90 :Goto5: ≠>Goto3△Lbl 4:Z=Q-T+L-K:Prog”8”:V=P+B-180-AV:W=P+B-AN:C=cosV+E+Tcos(P+B:Pol(X-C,Y-D:K=K+Isin(J-W+90 :Abs sin(J-W+90 <.001＝>M=Icos(J-W+90:Goto5: ≠>Goto4△Lbl 5:”K=”:K:Pause0:”M=”:M▲注：程序“20”为中桩计算程序Q—直缓点里程；T—切线长，L—曲线长；O—缓和曲线长；W—切线方位角Pol—极坐标转换，为涵数特别注意数字0和字母O的区别一、程序功能本程序由一个主程序(TYQXjs)和两个子程——正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。

卡西欧4800计算器程序

LX、RX：左侧或右侧点X坐标
LY、RY：左侧或右侧点Y坐标
I：置镜点至左或右点距离
J：置镜点至左或右点方位角
CONT：输入0计算下一里程；输入1则从本程序第一步开始
QX1计算范围：
DK214+832.674~DK216+053.229
QX2计算范围：
DK215+707.196~DK217+253.771(反算)
LbI 5 ↙
H=1=〉Y= -1×Y△ U=B+X cosT-YsinT：V=C+XsinT+YcosT：↙
X=U▲
Y=V▲
PoL（（X-M，Y-N）↙
I ▲
J＜0=〉J=J+360△↙
J ▲
{W}：W“NEXT”：W=1=〉Goto 1△↙
LHale Waihona Puke I 6 ↙ {O}：O“O=”：O=1=〉Goto 7△↙
QX3计算范围：
DK216+623.229~DK219+480.262
QX4计算范围：
DK217+923.771~DK222+114.707(反算)
J ▲
F=T ▲
Goto 6 △↙
LbI 3 ↙
L＞A=〉Goto 4 △ X=L-L5÷40÷（RA）2：Y=L3÷（6RA）-L7÷336÷（RA）3：Z=60×L2÷（πRA）+tan-1（Y÷X）：H=1=〉Z= -1×Z △F=T+Z▲
Goto 5 △↙
LbI 4 ↙
Z=180（Abs（E-D）-A）÷（πR）+90A÷（πR）：X=R sinZ+A÷2-A3÷（240R2）：Y=R（1-cos Z）+A2÷（24R）-A4÷（2688R3）：H=1=〉Z= -1×Z △ F=T+Z▲

线元法线路坐标正反算程序

经苦心钻研，奋战多日，终于编写出了代码短，速度快，精度高, 功能全的线路坐标正反算程序，欢迎试用并提出宝贵意见。

功能简介及特点：1、选用高斯-勒让德公式作计算内核，保证精度，模块化设计，便于扩充功能。

2、线元数据可自动从数据库调用，也可手工输入。

3、可管理多条线路，如里程不在线路或线元范围，将警告里程偏大、偏小。

4、边桩计算设计为导线式递推方式，可用于由一个中桩推出结构物所有角点坐标。

5、反算实现了智能化操作，只需输入线路号（或手工输线元资料）坐标，不需近似里程，即可自动从起点向后开始试算出里程、位置，如对算出里程、位置表示怀疑，还可以让计算器从终点起再向前试算下一个可能的位置（匝道、回头曲线同一坐标可能会有一个以上结果）。

第三次及以后试算才要求输入近似里程。

6、程序代码规范简洁，便于阅读、理解。

完整程序清单：ZFS %正反算主程序B=.1739274226:C=.5-B:Lbl 1:U"0 ZS 1 FS"=0=>Prog "ZS":工>U=1=>Prog"FS":工>Goto 1ZS %正算子程序{K}:Prog"ZZ":l=0:{l}:l"L" 丰 0=>"Prog"WY":丰 >Prog"ZB"FS %反算子程序{KVW}:V"XC"W" YC":Lbl 2:Prog "ZZ":I=V-S:J=W-T:Pol(I,J: J=J~F:K二K+Rec(l,J:Absl<1m=>Prog"WZ": 半 >Gota\2M=O:{M}:M"O NEXT"=O二>U=U+1:Goto 2: 半 >U=1ZZ %高斯法中桩子程序(4节点)Prog"XL":M=K-L: O=(P-R) 2PQR:D=.0694318442:E=.3300094782:F=1:G=1-E:H=1-D:l=5:Lbl 1:C[l]=A+MrC[l](1 P+OMC[l]:Dsz l:Goto 1:S=X+M(BcosD+CcosE+CcosG+BcosH:T=Y+M(Bsi nD+Csi nE+Csi nG+Bsi nHWY %外移点计算子程序Lbl 1:J=90:{J}:J=F+J"<":F=J:S=S+Rec(l,J:T=T+J: Prog"ZB":l=0:{l}:l"L" 工0=>Goto 1WZ %位置显示子程序"KJ":K:Pause 1:J丄ZB %坐标显示子程序"XY ":S:Pause 1:T AYC %异常处理子程序U=1=>K=L:U=2 △U=3=>K=M:U=4AU=5=>{K}:U=4 △K<L=>"V<!":Z=1△K>M=>">>!":Z=1DL %断链处理子程序"DL":K=L:I>0=>K=L+Q------------------ 以上为程序运算部分,以下为数据库部分XL %线路数据库选择子程序Lbl 1:Z=0:N"0 SD"=0=>Prog"0" △N=1=>Prog"1" △N=2=>Prog"2" △有几条线路仿上行格式输几行Z=1=>{NLX YOPQRK}:Goto 10 %手工输入子程序L"K0"X YAQ"LS"P"R0"R"RN":M二L+Q:Prog" YC"1 %线路一数据库子程序①Lbl B:L=线路起点里程:M=线路终点里程:Prog"YC":Z=1=>Goto E △②Q=线元长:P=起点半径:R=终点半径:K< L+Q二>X=起点X坐标:Y二起点Y坐标:A=起点方位角:Goto E △ L=L+Q:③…… @Q=短链长:K<L+Q=>Prog "DL":Goto B △ L=L+Q:⑤Q=线元长:P=起点半径:R=终点半径:K＜L+Q二＞X=起点X坐标:Y二起点Y坐标:A=起点方位角:Goto E △ L=L+Q:⑥……炉0=线元长:P=起点半径:R=终点半径:X=起点X坐标:丫二起点Y坐标:A=起点方位角丄bl E2 %线路二数据库子程序输入要求和线路一相同。

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

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

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

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

二、源程序
1.主程序（TYQXjs）
"1.SZ => XY"："2.XY => SZ"：N：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：{SZ}：SZ：S<O=> Goto 1：≠>Goto A:Δ
Lbl A：S>(O+H) =>Goto 1：≠>Goto 3Δ←┘
Lbl 3 :W=Abs（S-O）：Prog "SUB1"：X"XS"=X◢
Y"YS"=Y◢
Goto 1←┘
Lbl 2：{XY}：XY：I=X：J=Y：Prog "SUB2"：(O+W)<O=>Goto 2：≠>Goto BΔ←┘
Lbl B：(O+W)>(O+H) =>Goto 2：≠>Goto 4Δ←┘
Lbl 4：：S"S"=O+W◢
Z"Z"=Z◢
Goto 2
2. 正算子程序（SUB1）
A=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：F=1-L：
M=1-K：X=U+W（Acos（G+QEKW（C+KWD））+Bcos（G+QELW（C+LWD））+Bcos（G+QEFW
（C+FWD））+Acos（G+QEMW（C+MWD）））：Y=V+W（Asin（G+QEKW（C+KWD））+Bsin（G+ QELW（C+LWD））+Bsin（G+QEFW（C+FWD））+Asin（G+QEMW（C+MWD）））：F=G+QEW（C+ WD）+90：X=X+ZcosF：Y=Y+ZsinF
3. 反算子程序（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）÷sinF
三、使用说明
1、规定
（1）以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

（2）当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。

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

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

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

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

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

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

2、输入与显示说明
输入部分：
1. SZ => XY
2. XY = > SZ
N ？选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示由坐标反算里程和边距。

X0 ？线元起点的X坐标
Y0 ？线元起点的Y坐标
S0 ？线元起点里程
F0 ？线元起点切线方位角
LS ？线元长度
R0 ？线元起点曲率半径
RN ？线元止点曲率半径
Q ？线元左右偏标志（左偏Q=-1，右偏Q=1，直线段Q=0）
S ？正算时所求点的里程
Z ？正算时所求点距中线的边距（左侧取负，值右侧取正值，在中线上取零）
X ？反算时所求点的X坐标
Y ？反算时所求点的Y坐标
显示部分：
XS=×××正算时，计算得出的所求点的X坐标
YS=×××正算时，计算得出的所求点的Y坐标
S=×××反算时，计算得出的所求点的里程
Z=×××反算时，计算得出的所求点的边距
四、算例
某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程S0、起点坐标X0 Y0、起点切线方位角F0、线元长度
LS、起点曲率半径R0、止点曲率半径RN、线元左右偏标志Q）如下：
S0 X0 Y0 F0 LS R0 RN Q
500.000 19942.837 28343.561 125 16 31.00 269.256 1E45 1E45 0 769.256 19787.340 28563.378 125 16 31.00 37.492 1E45 221.75 -1 806.748 19766.566 28594.574 120 25 54.07 112.779 221.75 221.75
-1
919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9579.228
-1
999.812 19744.038 28781.659 80 40 50.00 100.000 1E45 1E45 0
1、正算
（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）
S=700 Z=-5 计算得XS=19831.41785 YS=28509.72590
S=700 Z=0 计算得XS=19827.33592 YS=28506.83837
S=700 Z= 5 计算得XS=19823.25398 YS=28503.95084
S=780 Z=-5 计算得XS=19785.25749 YS=28575.02270
S=780 Z=0 计算得XS=19781.15561 YS=28572.16358
S=780 Z= 5 计算得XS=19777.05373 YS=28569.30446
S=870 Z=-5 计算得XS=19747.53609 YS=28654.13091
S=870 Z=0 计算得XS=19742.68648 YS=28652.91379
S=870 Z= 5 计算得XS=19737.83688 YS=28651.69668
S=940 Z=-5.123 计算得XS=19741. 59118 YS=28722.05802
S=940 Z=0 计算得XS=19736.47687 YS=28722.35642
S=940 Z= 3.009 计算得XS=19733.47298 YS=28722.53168
2、反算
X=19831.418 Y=28509.726 计算得S=699.9999974 Z= -5 .00018164 X=19827.336 Y=28506.838 计算得S=699.9996493 Z= 0.000145136
X=19823.25398 Y=28503.95084 计算得S=699.9999985 Z= 5.000003137 X=19785.25749 Y=28575.02270 计算得S=780.0000035 Z= -5 .000001663 X=19781.15561 Y=28572.16358 计算得S=780.0000025 Z=- 0.000002979 X=19777.05373 Y=28569.30446 计算得S=780.0000016 Z= 4.99999578
X=19747.536 Y=28654.131 计算得S=870.0001137 Z= -4.99941049
X=19742.686 Y=28652.914 计算得S=870.0003175 Z=- 0.00041814
X=19737.837 Y=28651.697 计算得S=870.0002748 Z= 4.999808656
X=19741.5912 Y=28722.0580 计算得S=939.9999786 Z= -5.123024937 X=19736.4769 Y=28722.3564 计算得S=939.9999862 Z=- 0.000027710 X=19733.4730 Y=28722.5317 计算得S=940.0000238 Z= 3.00898694。