标准缓和曲线４８００程序
J＂JDZ＂：A＂JDX＂：B＂JDY＂：C＂FWJ1＂：D＂FWJ2＂：Ｒ：Ｌ＂ＬＳ＂：Ｎ：Ｔ：Ｅ＂ＺＨ＂：Ｆ＂ＨＺ＂
ＬＢ0：{ＫＵＺ}：Ｋ≤Ｅ=>ＧＯＴＯ１△Ｋ＞Ｅ=>Ｋ≤Ｅ+Ｌ：=>GOTO2△K＞E+L=>K＜F－L=>GOTO3△K≥F－L=>K＜F=>GOTO4△K≥F=>GOTO5 LBI1：P=C：X=A－(J-K)COSC▲Y=B-(J－K)SINC▲GOTO6△
LBI2 ：M=K－E：H==M－Mˇ5÷(40Rˇ2Lˇ2) ：O=Mˇ3÷(6RL) ：N=0=>S=C：P=C－28.7Mˇ2÷(RL) ：≠>S=C+180：P=C+28.7Mˇ2÷(RL) △X=A－TCOSC+HCOSC+OSINS▲Y=B-TSINC+HSIND－OCOSS▲GOTO6△
LBI3：M=K－(E+L) ：H=57.29579(M+2L)÷R：O=R(1－COSH)+Lˇ2÷(24R) ：N=0=>S=C：P=C-28.7(L+2M)÷R ：≠>S=C+180：P=C+28.7(L+2M)÷R△X=A－TCOSC+(RSINH+L÷2－Lˇ3÷240÷Rˇ2)COSC+OSINS▲Y=B-TSINC+(RSINH+ L÷2－Lˇ3÷240÷Rˇ2)SINC－OCOSS▲GOTO6△
LBI4：M=F－K：H==M－Mˇ5÷(40Rˇ2Lˇ2) ：O=Mˇ3÷(6RL) ：N=0=>S=C：P=C+28.7Mˇ2÷(RL) ：≠>S=C+180：P=C－28.7Mˇ2÷(RL) △X=A+TCOSD－HCOSD+OSINS▲Y=B+TSINC－HSIND－OCOSS▲GOTO6△
LBI5：P=D：X=A+(T+K－F)COSD▲Y=B+(T+K－F)SIND▲GOTO6△ LBI6：W＂BZX＂=X+ZCOS(P+U) ▲V＂ＢＺＹ＂=Y+ZSIN(P+U) ▲GOTO0△ 注:J为交点桩,AB分别为交点XY坐标,C为前一切线方位角,C为后一切张方位角,R半径,L缓和长,N(转角,左0右随意),T切线长,E为ZH桩号,F为HZ桩号,K为放样桩号,U为法线与切线夹角,Z为边桩与中桩距离,W为边桩X坐标,V为边桩Y坐标. ˇ为指数,下划线处为三角函数
(自己测试还得,要是大家计算正确我再加极坐标放样程序
667建筑网 www.667e.com 万份建筑资料免费下载！
CASIO fx—4800P型计算器对施工线路测量数据的连续计算
CASIO fx—4800P型计算器对施工线路测量数据的连续计算
如何适应施工现场快速、灵活的放样要求而适时快捷地计算施工线路上的任何测量数据？用《CASIOfx—4800P型计算器对施工线路测量数据的连续计算》能帮你实现这一设想。该程序对线路的三维坐标的计算不受单元、线路里程桩号、断面距离和断面斜交角度的影响，能方便快捷地得到施工线路的任何设计坐标值、高程值及其相应的放样数据。
该程序对于施工线路（回旋曲线相连的匝道线路除外）所需要的中桩、边桩坐标及放样参数以及它们的设计高程和高程的放样数据，只要输入其里程桩号、边桩距离、测站坐标或是视线高，这些数据即可获得。若是斜交构造物，只需输入其轴线与主线斜交的角度便也可获取相应的数据。
本程序稍加修改（即显示方式的修改）同样适合4850计算器。经过两个公路工程的运用，证明计算准确，使用方便，根据计算器的容量，可以输入数公里乃至十数公里的线路要素，真正实现整条施工线路的数据的连续计算。 本程序有计算例子！不但可以连续不间断的进行坐标或高程计算，还可以随意输入施工线路全线任意桩号里程而不受线路元素单元的影响即可得到相应的坐标或标高数据以及防样数据；不但可以计算施工线路数据，还可以计算施工线路以外的曲线坐标或高程，并且还能计算缓和曲线的超高值。 如果有不想进行线路外坐标计算的可将原程序稍做这样改动一下就会更方便了。
ZBFYJS（文件名） F：D：N\"XC\"：Q\"YC\"：B↙ Lbl0：{K，G}
K≤E2=>E=E1：U=X1：M=Y1：O=O1：R=R1：S=S1：T=T1：P=P1：F=±1：≠>↙
K≤E3=>E=E2：U=X2：M=Y2：O=O2：R=R2：S=S2：T=T2：P=P2：F=±1：≠>↙
667建筑网 www.667e.com 万份建筑资料免费下载！
K≤E4=>E=E3：U=X3：M=Y3：O=O3：R=R3：S=S3：T=T3：P=P3：F=±1：≠>↙
K≤En+1=>E=En：U=Xn:M=Yn：O=On：R=Rn：S=Sn：T=Tn：P=Pn：F=±1↙
Lbl1：K>E=>L=K-E：C=O：≠>L=E-K：C=O+180⊿↙ L>P=>Prog\"A\"↙
≠>R=0=>Rec（L，C）：X=U+I◢ Y=M+J◢ W=C+B：Goto5↙
≠>S=0=>A=180L/(πR)：X=R×SinA：Y=R（1-CosA）：Goto4↙ ≠>S≤L=>Goto2↙ ≠>S＞L=>Goto3↙ 接上“B”子程序第二行。
既将子程序“A”和“B”合并到主程序上。大家不妨试试！
如果你还想计算线路外坐标数据，可将所要计算的线元要素输入到数据库里，便也可计算该线元坐标数据。注意里程要与线路里程要有所区别！
667建筑网 www.667e.com 万份建筑资料免费下载！
常用CASIO4800计算程序集
FX4800P常用测量程序集
1、由一个已知点坐标和一个已知方位角计算支导线坐标（程序名：ZBJS） 程序清单：
X：Y：T：Lb1 0：Fixm：{AS}：T=T+A-180：T〈0 =〉T=T+360⊿ T〉360 =〉T=T-360⊿T◢ X=X+Scos T◢ Y=Y+Ssin T◢ Goto 0 程序运行说明：
运行“ZBJS”，每运行一步按EXE键一次，显示X？时输入已知点X坐标，显示Y？时输入已知点Y坐标，显示T？时输入已知点由后视到前视的已知坐 标方位角，运行一次本程序只需输入一次已知数据。显示A？时输入第一个左转角值，显示S？时输入测站到前视的水平距离。接下来计算器显示的T是测站到前视的坐标方位角，按 + 可将显示的十进制角度换算为度分秒。如此继续可计算出所有支导线未知坐标。
2、由两个已知点坐标计算支导线坐标（程序名：ZBJS2） 程序清单：
M“X0”：N“Y0”：D“XH”：E“YH” ↲
Pol(M-D，N-E)：Fixm：T=J：X=M：Y=N：T〈0=〉T=T+360⊿T“TH-J=” ◢ I“S=”◢
Lb1 0:{AS}: T=T+A-180：T〈0 =〉T=T+360⊿ T〉360 =〉T =T-360⊿T◢ X=X+Scos T◢ Y=Y+Ssin T◢ Goto 0 程序运行说明：
运行后，显示X0，Y0时输入已知置镜点坐标，显示XH，YH时输入已知后
667建筑网 www.667e.com 万份建筑资料免费下载！
视点坐标，其它同ZBJS的输入与显示。
3、由一个已知点坐标和一个已知方位角进行放射型导线坐标计算（程序名FSXZBJS） 程序清单：
C“X0”：D“Y0”：E“T0”：Fixm：Lb1 1:{SA}↲
T=E+A-180：T〈0=〉T=T+360⊿ T〉360 =〉T=T-360⊿T◢ X=C+Scos T◢ Y=D+Ssin T◢ Goto 1 程序运行说明：
变量输入参照ZBJS程序。
4、由两个已知点坐标进行放射型导线坐标计算（程序名FSXZBJS2） M“X0”：N“Y0”：D“XH”：E“YH” ↲
Pol(M-D，N-E)：Fixm：F=J：X=M：Y=N：F〈0=〉F=F+360⊿F“TH-J=” ◢ I“S=” ◢ Lb1 1:{SA}↲
T=F+A-180：T〈0=〉T=T+360⊿ T〉360 =〉T=T-360⊿T◢ X=X+Scos T◢ Y=Y+Ssin T◢ Goto 1
5、坐标计算优化程序（可进行连续式（LX）、放射式（FS）已知一点（1-D）或者两点（2-D）进行坐标计算）：
V“1=>LX，≠>FS”：W“1 =>1-D≠>2-D”：W=1=> Goto 1:≠> Goto2↲ Lb1 1：F“X0=”：G“Y0=”：H“T H-0=”：Goto 3↲
Lb1 2:F“X0=”：G“Y0=”：M“XH=”：N“YH=”:Pol(F-M，G-N)：Fixm: J〈0=〉J=J+360⊿H=J：H“TH-0=”◢ I“S-H=”◢ Goto 3↲
Lb1 3: X=F：Y=G：T=H ： V=1=> Goto 4:≠> Goto 5↲
667建筑网 www.667e.com 万份建筑资料免费下载！
Lb1 4: {AS}: A:S:T=T+A-180：T〈0 =〉T=T+360⊿ T〉360 =〉T =T-360⊿T◢ X=X+Scos T◢ Y=Y+Ssin T◢ Goto 4
Lb1 5：{AS}:A:S
T=H+A-180：T〈0=〉T=T+360⊿ T〉360 =〉T=T-360⊿T◢ X=F+Scos T◢ Y=G+Ssin T◢ Goto 5
6、已知置镜点坐标和置镜点到后视点坐标方位角求放样角度和放样距离的放样程序（程序名FY1） 程序清单：
F“X0”：G“Y0”：H“T0-H”：Fixm：Lb1 1:{XY}↲ Pol(X-F，Y-G) ↲
J〈0 =〉J=J+360⊿A=J-H↲ A〈0 =〉A=A+360⊿A“A=” ◢ I“S=” ◢ Goto 1 程序运行说明：
变量输入参照坐标计算，A=为放样顺时针角度，S=为放样距离。 7、已知置镜点和后视点坐标进行放样（FY2） 程序清单：
F“X0”：G“Y0”：M“XH”：N“YH”↲ Pol(M-F，N-G)：Fixm：H=J↲ H〈0=〉H=H+360⊿H“T-H=”◢ I“S-H=” ◢
Lb1 1:{XY}: Pol(X-F，Y-G) ↲ J〈0 =〉J=J+360⊿↲
667建筑网 www.667e.com 万份建筑资料免费下载！
A=J-H↲
A〈0 =〉A=A+360⊿A“A=”◢ I“S=” ◢ Goto 1
8、放样优化程序（可分已知一点或者两点坐标进行放样角度距离的计算） 程序名“FY” 程序清单： I“I=1 =>1-D≠>2-D”↲ I=1=> Goto 1:≠> Goto2↲
Lb1 1：F“X0”：G“Y0”：H“T0-H”：Goto 3↲
Lb1 2：F“X0”：G“Y0”：M“XH”：N“YH”:Pol(M-F，N-G)：Fixm:J〈0=〉J=J+360⊿H=J:H“T-H=”◢ I“S-H=”◢ Goto 3↲ Lb1 3：{XY}↲
Pol(X-F，Y-G) ：Fixm ↲ J〈0 =〉J=J+360⊿A=J-H↲ A〈0 =〉A=A+360⊿A“A=” ◢ I“S=” ◢ Goto 3 程序运行说明：
变量输入参照坐标计算，A=为放样顺时针角度，S=为放样距离。 9、已知两点坐标，反算两点之间的坐标方位角和距离（程序名ZBFS） 程序清单：
A“X1=”：B“Y1=”：C“X2=”：D“Y2=”：P ol（C-A，D-B）：Fixm：J〈0 =〉J=J+360⊿J“T（1-2）=” ◢ I“S（1-2）=”
10、分别已知两条直线上一点的坐标和坐标方位角，求这两条直线的交点坐标（程序名JDZB）
667建筑网 www.667e.com 万份建筑资料免费下载！
程序清单：
X=（tan A“T1”×B“X1”-C“Y1”-tan D“T2” ×E“X2”+F“Y2”） V÷（tan A- tan D）◢
Y= tan A（X-B）+C
11、直线线路坐标计算（ZXXY） 程序清单：
V“X0=”：W“Y0=”：T“T0=”：U“K0=”：Lb1 0：Fixm：{KD}：X=V+（K-U）cos T+D×cos(T+90)◢
Y=W+（K-U）sin T+D×sin(T+90)◢ Goto 0 程序运行说明：
已知某一直线线路上一点的坐标为X0，Y0，该点里程为K0以及该直线沿前进方向的坐标方位角T0，按照提示输入以上数据，如果输入某一里程K和法线外移距离D（向右移为正值，向左移为负值）即可计算出该里程点的法线外移桩坐标X、Y，当D=0时计算结果为中桩坐标。 12、“缓+圆+缓”型曲线坐标计算 主程序（QXJS） 程序清单：
Prog “BQXYS”：Q“K—ZH”：F“ZHX=”：G“ZHY=”：N“TZH—JD”：I“（Y1，Z-1）” ↲
Lbl 1:Fixm:{KD}:K“LICHENG=”：D“WAIYI=（Z-，Y+）”：Goto2↲ Lbl 2:K Goto3:≠> Goto4↲ Lbl 3:“K—OVER”: Goto1↲ Lbl 4:K>Q+L => Goto3:≠> Goto5↲ Lbl 5:K≤Q+B => Goto6:≠> Goto7↲
Lbl 6:S=K-Q：W=B：Prog “HQXY”：O=N+I×90S2÷(Л RB)：O<0 => O=O+360⊿O>360 => O=O-360⊿O“FANGWEI=” ◢ X=UcosN-VsinN+F+D×cos(O+90) ◢ Y=VcosN+UsinN+G+D×sin(O+90) ◢
667建筑网 www.667e.com 万份建筑资料免费下载！
Goto1↲
Lbl 7：K≤Q+L-C => Goto8:≠> Goto9↲
Lbl 8:S=K-Q：O=N+I×90(2S-B)÷(Л R)：O<0 => O=O+360⊿O>360 => O=O-360⊿O“FANGWEI=” ◢
Prog “YQXY”：X=UcosN-VsinN+F+D×cos(O+90) ◢ Y=VcosN+UsinN+G+D×sin(O+90) ◢ Goto1↲
Lbl 9: S= Q+L-K：W=C：Prog “HQXY”：E=N+IA+180：O=E-I×90S2÷(Л RC)+180：O<0 => O=O+360⊿O>360 => O=O-360⊿O“FANGWEI=”◢ V=-V↲
X=(U-H)cosE-VsinE+F+T cosN+D×cos(O+90) ◢ Y=VcosE+(U-H)sinE+G+T sinN+D×sin(O+90) ◢ Goto1↲
子程序“BQXYS” A:R:B“L01”:C“L02” ↲
M“M1”=B÷2-B3÷（240 R2）◢ J“M2”=C÷2-C3÷（240 R2）◢
P“P1”= B2÷（24 R）- B4÷（2688 R3）◢ Z“P2”= C2÷（24 R）- C4÷（2688 R3）◢ T“T1”=M+（R+P）tan(A÷2)+(Z-P) ÷sinA◢ H“T2”=J+（R+Z）tan(A÷2)+( P-Z) ÷sinA◢ L= B÷2+C÷2+A RЛ÷180◢ E= ◢
子程序“HQXY”
U=S-S5÷(40R2 W2)+ S9÷(3456 R4 W4)-S13÷(599040R6 W6):V=I×(S3÷(6RW)- S7÷(336 R3 W3) +S11÷(42240 R5 W5)) 子程序“YQXY”
U=Rsin(180(S-B÷2)÷Л÷R)+M:V=I×(R×(1-cos(180(S-B÷2)÷Л÷R))+P)
程序说明：本程序可以计算等长和不等长“缓+圆+缓”型曲线或者单圆曲线的中
667建筑网 www.667e.com 万份建筑资料免费下载！
线坐标、外移桩坐标以及中线的切线方位角。
先将主、子程序存入计算器，然后运行主程序“QXJS”，输入A（曲线偏角），曲线半径R，第一缓和曲线长L01，第二缓和曲线长L02，（等长缓和曲线时L01= L02），然后显示曲线逐要素，接着输入ZH坐标（ZHX，ZHY），ZH到JD坐标方位角TZH—JD，当曲线为左偏时（Y1，Z-1）输-1，右偏输1，输入计算里程K，外移距离D，计算中桩时D输0，计算左边桩时D输负数，算右边桩时D输正数，当所输里程不在本曲线计算范围内时，屏幕显示K—OVER，继续输入有效里程可继续计算。 对于单圆曲线，缓和曲线长输入0。
本程序缓和曲线X计算取了4项，Y计算取了3项，一般的曲线计算精度都可达到mm级。
13、采用积分进行任意曲线线路坐标计算程序（直线、圆曲线、卵型曲线） 程序名“XLJS”
程序名“XLJS”(V1.1，2005年12月7日星期三升级于秦岭公司)
G“XJ=”：H“YJ=”：R“XH=”：S“YH=”：A“XA=”：B“YA=”：C“TA=”：E“1÷RA ,Z-，Y+”：P“KA=”↲
Lbl 1：{FQ}：F“1÷RB, Z-，Y+”：Q“KB=” ↲ Lbl 2：{V}：V“K=”：O=V：O Goto3：≠> Goto4↲
Lbl 3：“K-SMALL”：Goto2↲ Lbl 4：O>Q => Goto5：≠> Goto7↲
Lbl 5：{I}：I“I=1=>ESC≠>NEXT”：I=1=> Goto2：≠> Goto 6↲ Lbl 6：“JISUAN：A-TXY”：V=Q：Z=0：D=0：Goto8↲ Lbl 7: {DZ}：Z“JIA-JIAO=”：D“WAIYI=（Z-，Y+）”↲
Lbl 8：T=C+90÷Л×（2E+（F-E）（V-P）÷（Q-P））（V-P）：T〈0=〉T=T+360⊿ T〉360 =〉T=T-360⊿↲
Rad：W=A+∫（cos（C°+0.5×(2E+(F-E)(X-P)÷(Q-P))(X-P)）,P,V,4）+ D×cos(T°+Z°) ↲
Y=B+∫（sin（C°+0.5×(2E+(F-E)(X-P)÷(Q-P))(X-P)）,P,V,4）+ D×sin(T°+Z°)：
667建筑网 www.667e.com 万份建筑资料免费下载！
Deg：Goto9↲
Lbl 9：O≤Q => GotoB：≠> GotoA↲
Lbl A：A=W：B=Y：C=T：E=F：P=Q： Goto1↲ Lbl B： W“X=”◢ Y“Y=”◢ T“T=”◢
Pol(R-G，S-H)：Fixm：M=J：Pol(W-G，Y-H) ：J=J-M↲ J<0 =>J=J+360⊿J>360 =>J=J-360⊿J“A=”◢ I“S=”◢ Goto 2
程序积分式中“°”的编辑输入为：“FUNCTION”键+5（DRG）+4（°）。 程序功能
由于程序根据线路曲线坐标计算通用数学模型编制，所以，对于任意线型的一段线路曲线（直线、圆曲线可看作是通用计算模型的特例），只要输入曲线信息和测站信息，都可以迅速地计算出该曲线上任意里程点处的切线方位、中桩或与该里程切线成任意夹角的外移桩坐标、极坐标法放样的角度和距离。 程序运行
运行程序“XLJS”，按提示输入置镜点坐标（XJ，YJ），后视点坐标（XH，YH），曲线起点A的坐标（XA，YA）、沿线路前进方向的切线方位角TA、里程KA、曲率（直接在计算器上以1÷RA的形式输入，R表示半径，左偏曲线输负值右偏曲线输正值），继续输入终点B的曲率（1÷RB输入同A点）、里程KB。当屏幕显示“K=？”时，输入该段曲线要计算放样数据的里程值，按提示输入夹角“JIA-JIAO”（夹角的范围为0-180度）和外移量“WAIYI”，“JIA-JIAO”值为该里程自线路前进方向的切线方位角顺时针至该里程中桩与右边桩连线间的角度，当为法线方向时JIA-JIAO为90度，“WAIYI”值中桩左侧输负值右侧输正值，计算中桩时“WAIYI”值输0即可，按“EXE键”屏幕即可依次显示计算中桩或外移桩的坐标、中桩切线方位角、放样角度A，放样距离S。 程序的输入输出都按通俗易懂的提示或测量习惯显示，为确保计算结果无误，编制程序时专门设置了错误里程输入提示，当输入的计算里程K小于该曲线起点里程时，程序将提示里程太小“K-SMALL”，继续输入正确的里程，当第一条曲线计算完毕需计算紧接的第二条曲线时，只需要给K输入一个大于第
667建筑网 www.667e.com 万份建筑资料免费下载！
一曲线终点里程的里程值，这时屏幕显示“I=1=>ESC≠>NEXT”，如果不小心误输入一个大于终点里程的K时（ESC表示取消），可输入1“取消”操作，继续计算当前曲线，如果确认要计算下一曲线（NEXT表示要计算下一曲线），输入一个不等于1的数值，接下来屏幕将显示“JISUAN：A-TXY”，程序将自动计算第二曲线的起点信息，接着自动提示输入第二曲线终点的曲率和里程，当显示“K=？”时，即可进行第二曲线放样数据的计算。依次类推，直至计算完所有需计算的曲线中、边桩坐标和放样数据。
迁站时，直接按“AC”键终止程序的运行，到新的测站后重新运行“XLJS”，输入新的测站点、后视点坐标并保持原曲线参数不变即可。
当遇到线路曲线的曲率发生突变（如直线直接接圆曲线或者两个不同半径的圆曲线直接相接）的特殊情况时，可以将突变点两侧的曲线分开单独进行计算，也可以在突变点处增加输入一段微小长度的回旋线把突变点两侧的曲线连接起来，增加一段长度为0.0001米长度的回旋曲线，不会影响计算结果的精度。比如某直线直接连接半径为R的圆曲线，计算完直线段需要接着计算圆曲线时，先输入微小长度的回旋线（1/RB=1/R，KB=直线终点里程+0.0001），然后再输入圆曲线的终点曲率和里程即可继续计算。两个不同半径（半径分别为R1，R2）的圆曲线直接相接时，先输入微小长度的回旋线（1/RB=1/R2，KB=第一圆曲线终点里程+0.0001），然后输入第二圆曲线终点的曲率和里程即可。 有时会遇到逆向测设曲线情况，即从大里程方向向小里程方向测设曲线，采用本程序进行计算时，只需将曲线倒过来进行输入即可，这时，输入和计算的量都应反号（只有“JIA-JIAO”值符号保持不变），曲线偏角反号输入（左偏变为右偏，右偏变为左偏），曲线里程输入时按负值对待，计算右边桩时外移量为负，左边桩时外移量为正（与正向刚好相反）。
14、四等水准测量记录计算程序（红黑面常数为4687、4787） Lbl 0:“H4687=>I=1”: Fixm: {IABCD}：I=1=> Goto1: ≠> Goto2↲ Lbl 1:U=4687:W=4787: Goto3: Lbl 2:U=4787:W=4687: Goto3: Lbl 3: G“H1”=A“HOU1”-B“QIAN1”◢ E“V2”=B-C“QIAN2”+W◢ F“V1”=A-D“HOU2”+U◢ P“H2”=D-C◢ H=（G+P+W-U）÷2
667建筑网 www.667e.com 万份建筑资料免费下载！
程序按后-前-前-后的观测顺序编制，当后尺为4687时I赋值1，否则赋不等于1的值。
15、高斯抵偿高程投影面边长改正计算FX-4850程序 E“YI-CHANG：hm=”：F“DI-CHANG：Hm=”↲
Lbl 1:{HYS}：H“PING-JUN：Hp=”:Y“PING-JUN：Ym =”:S“GUAN-CE:S=”:V=-S×(H+E-F)÷(6371000+H+E)×1000:W=S×Y2÷(2×63712)×1000:D=S+( V+W)÷1000↲ “V1=”:V◢ “V2=”:W◢ “D=”:D◢ Goto 1
说明：该程序可计算高斯投影或抵偿高程面高斯投影的边长投影，hm为异常高程值，当设计资料中没有提供或者未知时取hm=0，即认为高程异常可忽略，Hm为抵偿高程，当没有设抵偿高程面时取Hm=0，即为一般的高斯投影。S为观测水平距离，Hp为导线边两端高程概略平均值（m），Ym为导线边两端点Y坐标平均值减去500公里后的平均值（Km），D为改正到高斯投影面后的边长（m）。
hm=0，Hm取隧道平均高程面高程，Ym=0时，为计算隧道平均高程面的投影改正
667建筑网 www.667e.com 万份建筑资料免费下载！
曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序
曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序 一、程序功能
本程序由一个主程序(TYQXJS)和两个子程——正算子程序(SUB1)、反算子程序(
SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线
元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲
率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另
外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲
线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运行。
特别申明：(1). 适用于弧长小于2倍半径的各种线元坐标正反算，精度优 于1mm；
(2). 在引用该核心计算部分时，请注明来源。 二、源程序 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：W=Abs(S-O)：Prog \"SUB1\"：X\"XS\"=X◢ Y\"YS\"=Y◢ F\"FS\"=F-90◢ Goto 1←┘
667建筑网 www.667e.com 万份建筑资料免费下载！
Lbl 2：{XY}：XY：I=X：J=Y：Prog \"SUB2\"：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) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆 弧的半径。
667建筑网 www.667e.com 万份建筑资料免费下载！
(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坐标 显示部分：
667建筑网 www.667e.com 万份建筑资料免费下载！
XS=××× 正算时，计算得出的所求点的X坐标 YS=××× 正算时，计算得出的所求点的Y坐标 FS=××× 正算时，所求点对应的中线点的切线方位角 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.3 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.688 YS=28652.91379
667建筑网 www.667e.com 万份建筑资料免费下载！
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.352 S=940 Z= 3.009 计算得 XS=19733.47298 YS=28722.53168 2、 反算
X=19831.418 Y=28509.726 计算得 S=699.9999974 Z= -5 .000181 X=19827.336 Y=28506.838 计算得 S=699.99993 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.35 计算得 S=939.9999862 Z=- 0.000027710 X=19733.4730 Y=28722.5317 计算得 S=940.0000238 Z= 3.008694
667建筑网 www.667e.com 万份建筑资料免费下载！
因篇幅问题不能全部显示，请点此查看更多更全内容