CASIOfx-4800P公路测量程序

1、主程序（XLZB）
Defm 10:Norm:O“X0”：G“Y0”:N≥1=＞Goto 1：≠＞N≥2=＞Goto 2：≠＞N≥3=＞Goto 3：≠＞N≥4=＞Goto 4△△△△
Lb1 1：｛KPM｝：KPM：Prog“C”：I=K-J：Prog“A”：X“XS”=X▲
Y“YS”=Y▲
Prog“W”
Goto 1
Lb1 2：｛KXY｝：KXY：A=X：B=Y：Prog“C”：Prog“B”：K=J+I▲
P=P▲
Goto 2
Lb1 3：｛KSQAB｝：KSQAB：Prog“J”: Prog“C”: I=K-J： Prog“A”: X“XS”=X▲
Y“YS”=Y▲
Prog“W”
Goto 3
Lb1 4：｛KE｝：KE： Prog“W”: Prog“C”: I=K-J：P=T: M=90:Prog“A”: X“XS”=X▲
Y“YS”=Y▲
Prog“W”
Goto 4
a、当N=1时，求一般坐标放样。K：待求点桩号  P：偏移距离，当求中桩时P=0（左负右正） M：与右侧夹角，当求中桩时M=0
b、当N=2时，求桩号及距中桩距离。K：待求点桩号曲线上随意的一个桩号   
X：反算时所求点的X坐标      Y：反算时所求点的Y坐标
c、当N=3时，求构造物防样几个角点坐标。K：待求点桩号  S：构造物的角点编号
Q：与路线右夹角   A：横向距离    B：纵向距离 
d、当N=4时，所测标高点进行边线放样。K：待求点桩号   E：待求点边桩位置标高



2、正算子程序（A）
Z[5]=90*I/п/R：Z=W+Z[3]*2*Z[5]+Z[4]*I*L*Z[5]：X=U+Z[2]*I*cosZ+Z[3]*（2RsinZ[5]*cos（W+Z[5]））+Z[4]（（I-I5/（40R2L2））cosW-（I3/（6RL）-I7/（336R3L3））sinW）：Y= V+Z[2]*I*csinZ+Z[3]*（2RsinZ[5]*sin（W+Z[5]））+Z[4]（（I-I5/（40R2L2））sinW+（I3/（6RL）-I7/（336R3L3））cosW）：Z＜0=＞Z=Z+360△Z＞360=＞Z=Z-360△Z“QA”=Z：X=X+Pcos（Z+M）：Y=Y+Psin（Z+M）
  
3、反算子程序（B）
T=W-90：I=Abs（（B-V）cosT-(A-U)sinT：P=0：M=90：Lb1 0：Prog“A”：L=T+Z[3]*180/п*I/R：P=（B-Y）cosL-（A-X）sinL：Abs P〈1E-6＝＞Goto 1：≠＞I=I+P：Goto 0△
Lb1 1：P=0：Prog“A”：P=（B-Y）/sin（Z-90）
式中的E为计算器中的“EXP”符号

4、平曲线要素（C）
K≥***=＞J=***：U=***：V=***：W=***：R=***：L=***：Z[2]=***：Z[3]=***：Z[4]=***：K≥***=＞J=***：U=***：V=***：W=***：R=***：L=***：Z[2]=***：Z[3]=***：Z[4]=***：K≥***=＞J=***：U=***：V=***：W=***：R=***：L=***：Z[2]=***：Z[3]=***：Z[4]=***：K≥***=＞J=***：U=***：V=***：W=***：R=***：L=***：Z[2]=***：Z[3]=***：Z[4]=***：K≥***=＞J=***：U=***：V=***：W=***：R=***：L=***：Z[2]=***：Z[3]=***：Z[4]=***：K≥***=＞J=***：U=***：V=***：W=***：R=***：L=***：Z[2]=***：Z[3]=***：Z[4]=***：
……
K：待测点桩号    J：线元起点桩号   U：线元起点的X坐标    V：线元起点的Y坐标
W：线元起点切线方位角    R：曲线半径   L=缓和曲线长度   
当为直线段时Z[2]=1、Z[3]=0、Z[4]=0；
当为圆曲线段时Z[2]=0、Z[3]=1、Z[4]=0；
当为缓和曲线段时Z[2]=0、Z[3]=0、Z[4]=1；



6、高程计算程序(W)
Lb1 1：Prog“X”：H=A+D*（K-B）2/2/R+（K-B）*C：H=H-（14-1.5）*0.02：S=H-E：Prog“Y”：T=（25.4-0.75（S+8））Z[6]+（15.9-S）Z[7]+（14+1.5S）Z[8]+（26+1.75（S-8））Z[9]
本程序计算得出的宽度

7、宽度计算子程序（Y）
S＞-16=＞Z[6]=1：Z[7]=0：Z[8]=0：Z[9]=0：S＞-8=＞Z[6]=0：Z[7]=1：Z[8]=0：Z[9]=0：S＞0=＞Z[6]=0：Z[7]=0：Z[8]=1：Z[9]=0：S＞8=＞Z[6]=0：Z[7]=0：Z[8]=0：Z[9]=1


7、竖曲线要素（X）
K≥***=＞A=***：B=***：C=***：D=***：         （竖曲线中直线段）
K≥***=＞A=***：B=***：C=***：D=***：R=***：  （竖曲线中曲线段）
K≥***=＞A=***：B=***：C=***：D=***：         （竖曲线中直线段）
K≥***=＞A=***：B=***：C=***：D=***：R=***：  （竖曲线中曲线段）
……
K：待测点桩号     A：起点桩号     B：起点桩号标高     C：坡度     R：曲线半径
当计算段为直线时D=0，为曲线时D=-1；当为凸形竖曲线时R取正值，为凹形竖曲线时R取负值。


8、构造物角点子程序（J）
S≥1=＞P=√（A2+B2+2ABcosQ）：M=cos-1（（B2+P2-A2）/2/B/P）：
S≥2=＞P=-√（A2+B2-2ABcosQ）：M=-cos-1（-（B2+P2-A2）/2/B/P）：
S≥3=＞P=-√（A2+B2+2ABcosQ）：M=cos-1（-（B2+P2-A2）/2/B/P）：
S≥4=＞P=√（A2+B2-2ABcosQ）：M=-cos-1（（B2+P2-A2）/2/B/P）：
S≥5=＞P=A：M=Q

9、算测点到待测点距离及转角的子程序（W）
I=0:J=0:Norm:Pol（O-X，G-Y）:Fix 3:A“S”=I▲
Fix 6：B=J+180：B“A”=Int B+Int(Frac B*60)/E2+Frac(60B)*6E-3▲
式中的E为计算器中的“EXP”符号

说明：X0：测站点X坐标   Y0：测站点Y坐标
XS=******    计算得出所求点的X坐标
YS=******    计算得出所求点的Y坐标
K=*******    计算得出所求点的桩号
P=*******    计算得出所求点的边距（当P为正值为左边，当P为负值为右边）
S=*******    计算得出测站点到待测点的距离
A=*******    计算得出测站点到待测点转角(小数点前为度，小数点后第一二位表示分，三四五六位表示为秒