谁知道卡西欧fx4800程序输入fx4850之间的转换代码啊？或怎么转变

有缓和曲线或无缓和曲线任意桩号中、边桩坐标4800计算程序

有缓和曲线或无缓和曲线任意桩号中、边桩坐标4800计算程序

四、QXZB(曲线坐标)

A”ZH(ZH)=〃:B”X(ZH)=〃:C”Y(ZH)=〃 ”X(HZ)=〃:E”Y(HZ)=〃:F”FWJ-1=〃:G”FWJ-2=〃:H”R=〃:I”L1=〃:J”L=〃:K”L2=〃:L”ZJ(Z=1,Y=-1)=〃

O=180/pi (pi为圆周率，电脑中找不到此符号以此代替)

LbI A

{MN}

M”ZHUANG HAO=〃

P=M-A

MM>(A+I+J+K)?GOTO A⊿

N”ZHUANG JU(Z+,Y-)=〃

M>(A+I+J)?GOTO D⊿

M>(A+I)?GOTO C⊿

LbI B(第一缓和曲线)

X=B+(P-(P5÷40÷H2÷I2))COSF+(P3÷6÷H÷I-P7÷336÷H3÷I3)LSINF

Y=C+(P-(P5÷40÷H2÷I2))SINF-(P3÷6÷H÷I-P7÷336÷H3÷I3)LCOSF

X=X+NCOS(F-(P2÷2÷H÷I)OL-90)◢

Y=Y+NSIN(F-(P2÷2÷H÷I)OL-90)◢

GOTO A

LbI C(圆曲线)

X=B+(HSIN(((P-I)÷H+I÷2÷H)O)+(I÷2-I3÷240÷H2))COSF-(H(1-COS(((P-I)÷H+I÷2÷H)O))+ I2÷24÷H)(-1)LSINF

Y=C+(HSIN(((P-I)÷H+I÷2÷H)O)+(I÷2-I3÷240÷H2))SINF+(H(1-COS(((P-I)÷H+I÷2÷H)O))+I2÷24÷H)(-1)LCOSF

X=X+NCOS(F-((P-I)÷H+I÷2÷H)OL-90)◢

Y=Y+NSIN(F-((P-I)÷H+I÷2÷H)OL-90)◢

GOTO A

LbI D(第二缓和曲线)

Q=A+I+J+K-M

X=D-(Q-Q5÷40÷H2÷K2)COSG+(Q3÷6÷H÷K-Q7÷336÷H3÷K3)LSING

Y=E-(Q-Q5÷40÷H2÷K2)SING-(Q3÷6÷H÷K-Q7÷336÷H3÷K3)LCOSG

X=X+NCOS(F-(I÷2÷H+J÷H+K÷2÷H-Q2÷2÷H÷K)OL-90)◢

Y=Y+NSIN(F-(I÷2÷H+J÷H+K÷2÷H-Q2÷2÷H÷K)OL-90)◢

GOTO A

变量说明：

以上变量依次为：

ZH(ZH)=直缓点桩号，无缓和曲线则为直圆点；

X(ZH)=直缓点X坐标；无缓和曲线则为直圆点；

Y(ZH)=直缓点Y坐标，无缓和曲线则为直圆点；

X(HZ)=缓直点X坐标，无缓和曲线则为圆直点；

Y(HZ)=缓直点Y坐标，无缓和曲线则为圆直点；

FWJ-1=第一切线方位角；

FWJ-2=第二切线方位角；

R=半径；

L1=第一缓和曲线长；

L=圆曲线长；

L2=第二缓和曲线长；

ZJ(Z=1,Y=-1)=转角形式(左=1，右=-1)

ZHUANG HAO=所求点桩号

ZHUANG JU(Z+,Y-)=边桩桩距，左边桩为正，右为负，中桩0

我已经用过了，没发现什么问题，请大家提意见。

附闭合导线计算

1、 源程序

F1 A1

L1 Defm 4N-2

L2 N:A:B:Pol(C-A,D-B):W<0=>W=W+360?T=W

L3 K=0=>M=T+180:E=C:F=D:GOTO 0: ≠>E:F:Pol(G-E,H-F):W<0=>W=W+360?M=W

L4 Lbl 0:L=0:U=0:I=0:R=2:Z[1]=T

L5 Lbl 1:{J}:Z[R]+360: ?R=N+1=>GOTO 2: ≠>R=R+1:GOTO 1

L6 Lbl 2:P”JB”=(Z[N+1]-M) ?Q”JL”=40√N?R=2

L7 Lb1 3:{S}:Z[N+R]=S:L=L+S?

L8 Z[2N-1+R]=Rec(S,(Z[R]-P(R-1)/N)):U=U+V

L9 Z[3N-2+R]=W:I=I+W:N=R=>GOT 4: ≠>R=R+1:GOTO 3

L10 Lbl 4:P=U+C-E?Q=I+D-F?

L11 G”1:M”=L/Pol(P,Q) ?R=2

L12 Lbl 5:X”XI”=C+Z[2N-1+R]-PZ[N+R]/L?Y”YI”=D+Z[3N-2+R]-QZ[N+R]/L?

L13 R=N=>GOTO 6: ≠>R=R+1:C=X:D=Y:GOTO 5

L14 Lbl 6:”END”

2、 说明

（1）、本程序可计算附和导线和闭合导线的坐标，计算的坐标系经过角度闭合差及坐标增量闭合差分配后的结果，能显示角度闭合差、增量闭合差及导线全长的相对精度；

（2）、输入的观测角为导线的左角。

3、程序代号注释

N?导线观测角的折角数；

A、B?导线起始点所后视的已知点的坐标x,y；

C、D?导线起始点（即设站点）的坐标x,y；

E、F?导线终点（已知点）的坐标x,y；

G、H?在导线终点设站观测前视已知点的坐标x,y；

T?起始站后视至起始点的方位角；

M?终点站至前视已知点的方位角；

J?观测的左角值；

JB?角度闭合差；

JL?允许的角度闭合差，程序中是以40√n计算的，如和要求的不一致，可改一下L6语句中的有关部分。

S?所测导线的边长；

L?边长的累计数；

U?△x的累计数；

I?△y的累计数；

P?x坐标的闭合差；

Q?y坐标的闭合差；

K?转换符，当K=0时为计算闭合导线，当K≠0（任意数）时为计算附和导线。

铁路曲线坐标计算程序

( for CASIO fx-4500P )

一、主程序：M（M是文件名，下同)

L1 A”CX”B”CY”C”HX”D”HY”E”ZHX”F”ZHY”Q”HZX”L”HZY”G”ZHDK”H”ZHJD°”

I”A°”J”L0”K”L1”R

L2 Pol(C-A,D-B):M=W

L3 Lbl 2:W=90:{NOWZ}:Z”DK”N”Z/!/Y(1/2/3)”:N=2=>Goto 1⊿O”M”W”A°”

L4 Lbl 1:P=Z-G:S=E:T=F:U=H:V=H:Fixm

L5 P≤0=>Prog 3:≠=>P≤J=>Prog 1:≠=>P≤J+K=>Prog 2:V=90J/π/R+180π-

1R-1

(P-J)⊿⊿P≤J+K=>Prog 4:V=H+V⊿⊿

L6 P>J+K=>P=2J+K-P:S=Q:T=L:P>0=>Prog 1:X=-X:U=H+I:Prog 4: V=U-V:

≠=>

U=H+I:V=U:P=-P:Prog 3⊿⊿

L7 N≠2=>P=O:U=V-W:N=3=>P=-P⊿S=X:T=Y:Prog 3⊿

L8 Pol(X-A,Y-B):W=W-M:W<0=>W=W+360⊿

L9 W:”°°°=”◢V:”S=”◢X:”X=”◢Y:”Y=”◢Goto 2

二、子程序：1

L1 V=90P2(πRJ)-1

L2 X=P-Pxy5(40R2J2)-1

L3 Y=PVπ/540

三、子程序：2

L1 Y=180π-1R-1(P-0.5J)

L2 X=RsinY+J/2- Jxy3/240/R2

L3 Y=J/24/R+R-RcosY

四、子程序：3

L1 X=S+PcosU

L2 Y=T+PsinU

五、子程序：4

L1 I<0=>Y=-Y:V=-V⊿

L2 S=S+XcosU-YsinU

L3 Y=T+XsinU+YcosU

L4 X=S

说明：1、该程序适用于计算器 CASIO fx-4500PA。

2、程序符号定义说明：测站点坐标CX,CY；后视点坐标HX,HY；直缓点坐标

ZHX,ZHY；缓直

点坐标HZX,HZY；直缓点里程ZHDK；经过直缓点和交点的直线的方位角ZHJD°；

曲线偏

角（曲线左偏为负，右偏为正）A°；缓和曲线长L0；圆曲线长L1圆曲线半径R；

计算

点里程DK；计算曲线的左边点或者右边点）Z/!/Y(1/2/3)，曲线的右边输入1，左

边输

入3，线路上输入2；左边的点与当前里程点的连线的距离和当前里程点切线的逆

时针

方向的角度A°（如果上一项不选择2的话）。

3、计算结果显示的数据是：极角、极距和坐标。

4、该程序适用于极坐标法测量曲线时的数据（极角和极距）计算，以及曲线坐标

计算。

面积计算（多边形法）

1、源程序

F1 A2

L1 N:P=A:Q=B:S=0:I=2

L2 Lbl 0:{C,D}:F=(A+C)(B-D):S=S+F

L3 A=C:B=D:I=I+1

L4 I≤N=>GOTO 0⊿

L5 F=(C+P)(D-Q):S=S+F:S”W”=S/2◢

3、 说明：

（1）、本程序适用于所测断面为多边形闭合图形的面积计算。

（2）、折点坐标按顺时针方向输入，得出的面积为正，否则为负，绝对值是一样

的。

4、 程序代号注释

A、B—计算面积起始点纵横坐标；

C、D—各转折点的纵横坐标；

S—代表计算过程中的有关面积；

S“W”—为图形最后需要的计算面积。

N—多边形的折点个数。

体积计算

1、源程序

F1 A3

L1 J=0:H=0:WG

L2 Lbl 0:{NAB}:NAB:P=A:Q=B:S=0:I=1

L3 Lbl 1:{CD}:S=S+(A+C)(B-D)/2:A=C:B=D:I=I+1

L4 I<N=>GOTO 1⊿S=S+(C+P)(D+Q)/2◢

L5 J≠1=>GOTO 2: ≠>L=G-H:V=(R+S+√（R*S))*L/3◢⊿W=W+V◢

L6 Lbl 2:R=S:H=G:J=1:{G}:G:GOTO 0

2、说明

（1）、程序可自动计算每一断面面积，当进行到第二个断面时就会显示出1~2断

面间的体积，而后再进行第三断面面积计算，并累计出1~3断面之间的体

积。。。。。。，直到最后得出需算断面间的总体积。

（2）、坐标输入时，应按顺时针方向逐个输入折点坐标，这样得出的面积为正

值，一个桩号折点输入完后，程序自动进入下一桩号的输入状态。

3、程序代号注释

G—断面桩号；

A、B—断面起算折点的坐标；

C、D—断面上其他折点坐标；

S—断面面积；

L—断面间距；

V—本断面与前一断面之间计算出的体积。

N—G桩号断面上的折点个数；

W—本断面之前所有体积之和。

在任意控制点上测定直线上的任意位置与高程数据计算

1、源程序

F1 A4

L1 A”X”:B”Y”:E”Z0”:C”X1”:D”Y1”:F”Z1”:G”H0”:I

L2 Lbl 0:{LJV}:LJV

L3 K= tanF:T=tan(E+L)

L4 X”XP”=(B-D-AT+KC)/(K-T) ◢Y”YP”=B+TX-AT◢

L5 S=√（（X-A）2+（Y-B）2）

L6 H=G+StanJ+I-V◢GOTO 0

2、说明

（1）、本程序功能：在已知断面上某一点的坐标及断面方向的方位角后，可直接

将仪器架设在邻近控制点上，为测设各个断面上的点提供数据，不需要一定要将

仪器架设在断面桩上测设断面，而所测的断面点均有坐标和高程。

（2）、基本原理：如图所示，A、B为已知控制点，其方位角为Z0，断1-断2为断

面线控制桩，其方位角Z1可以算出，仪器架设在A点，后视B点，转角L1、L2、

L3。。。。。。，这时A~1、A~2，。。。。。。的方位角也为已知，根据解析

几何，两方位直线相交，可解出1#，2#。。。。。。的坐标及其与A的距离。

3、程序代号注释

A、B—测站点坐标

C、D—断面桩点坐标

L—观测断面点的水平角

J—观测断面点的竖直角

V—觇标高

S—测站至断面测点的距离

H—断面测点的高程

I—仪器高

Z0—测站至后视点的方位角

Z1—断面线的方位角

竖曲线计算

1、源程序

F1 A5

L1 BADTRZ

L2 Lbl 0:{C}:C

L3 Z≥1=>H=A+(B-A)/T*(C-D)-(C-D)2/(2R) ◢≠>H=A+(B-A)/T*(C-D)+ (C-D)

2/(2R) ◢

L4 GOTO 0

2、说明

（1）、本程序的功能是根据道路施工纵断面图上的设计数据，算出竖曲线上各加

桩点的高程；

（2）、本程序适用于由小桩号向大桩号端方向计算，在键入全部已知数据后，当

C出现时，只要键入该点的桩号，高程立即会显示出来。

3、程序代号注释

A—起点（或终点）高程

C—需计算点的桩号

T—竖曲线切线长

D—起点（或终点）高程

R—竖曲线半径

H—C桩号处的高程

B—切线交点的高程

Z—曲线凹凸判断符Z≥1时为凸曲线，z<0时为凹曲线。

两点测角前方交会坐标计算

1、源程序

F1 A6

L1 ABCDEF

L2 X“XP”=（A/tanF+C/tanE-B+D）/(1/tanE+1/tanF) ◢

L3 Y“YP”=（B/tanF+D/tanE-C+A）/(1/tanE+1/tanF) ◢

说明：

E—1#点的观测角

F—2#点的观测角

1#、2#点的编号时应注意：面向交会点P的左侧定为1#点，右侧定为2#点。

竖曲线计算

1、源程序

F1 A5

L1 BADTRZ

L2 Lbl 0:{C}:C

L3 Z≥1=>H=A+(B-A)/T*(C-D)-(C-D)2/(2R) ◢≠>H=A+(B-A)/T*(C-D)+ (C-D)

2/(2R) ◢

L4 GOTO 0

2、说明

（1）、本程序的功能是根据道路施工纵断面图上的设计数据，算出竖曲线上各加

桩点的高程；

（2）、本程序适用于由小桩号向大桩号端方向计算，在键入全部已知数据后，当

C出现时，只要键入该点的桩号，高程立即会显示出来。

3、程序代号注释

A—起点（或终点）高程

C—需计算点的桩号

T—竖曲线切线长

D—起点（或终点）高程

R—竖曲线半径

H—C桩号处的高程

B—切线交点的高程

Z—曲线凹凸判断符Z≥1时为凸曲线，z<0时为凹曲线。

两点测角前方交会坐标计算

1、源程序

F1 A6

L1 ABCDEF

L2 X“XP”=（A/tanF+C/tanE-B+D）/(1/tanE+1/tanF) ◢

L3 Y“YP”=（B/tanF+D/tanE-C+A）/(1/tanE+1/tanF) ◢

说明：

E—1#点的观测角

F—2#点的观测角

1#、2#点的编号时应注意：面向交会点P的左侧定为1#点，右侧定为2#点。 线路中、边桩测量放样程序 (Ver 3.2)

F1 XLCS （主程序，步数385）

L1 Norm:Deg:U=O〃A0〃:Prog 1:Q=U:C〃X-JD〃:D〃Y-JD〃:U=A〃A:R+L-

〃:Prog 1:

B=Abs U:R:S〃L0〃:E〃K-ZH〃

L2 Fix 3:M=.5S-Sxy3/240R2:P=S2/24R:T〃T〃=(R+P)tan.5B+M◢F〃L〃=

πRB/180

+S◢F=F+E:Norm

L3 K〃ZJ:XY＝＞1〃:K≠1＝＞{L}:V=L〃K〃:U=0:Prog 4: G=X:H=Y:≠＞G〃

X〃:H〃Y〃◣

{KJ}:K〃HS:XY＝＞1〃:K=1＝＞{I}:X=I〃X〃:Y=J〃Y〃:≠＞V=J〃K〃:U=0:Prog 4◣

Fix 3

L4 V〃D0〃=Pol(X-G,Y-H◢Fix 4:N=W:W<0＝＞W=W+360◣Prog 2:W〃

A0〃=W◢

Norm:U=0:{U}:U〃AB〃:Prog 1:K=U

L5 Lb1 0:U=0:{U}:U〃CS:XY＝＞1〃:U=0＝＞{UZ}:V=Z〃K〃: U〃BZ:R+L-

〃:A<0＝＞U

=-U◣Prog 4: ≠＞{XY}:X:Y◣Prog 5: Prog 6:Goto 0

F2 1 （十进制→六十进制子程序，步数19）

L1 U=Int U+Frac U/.6+Frac 100U/90

F3 2 （六十进制→十进制子程序，步数22）

L1 60Frac W:W=Int W+.01Int Ans+.006Frac Ans

F4 3 （缓和曲线上任意点坐标计算子程序，步数57）

L1 Y=RS:X=V-Vxy5/40Y2:Y=Vxy3(1-Vxy4/56Y2)/6Y+URec(1,90V2/π

Y:X=X-UW

F5 4 （中桩、边桩坐标计算子程序，步数195）

L1 V＞F＝＞X=T-UsinB+Rec(V-F+T,B:Y=W+UcosB:Goto 0◣

L2 V＞F-S＝＞V=F-V:Prog 3: Rec(1,B:U=X:X=T+TV-XV-YW:Y=TW-

UW+YV:

Goto 0◣

L3 V＞E+S＝＞Y=P+R-Rec(R-U,180(V-E-.5S)/πR:X=M+W:Goto 0◣

L4 V＞E＝＞V=V-E:Prog 3: ≠＞X=V-E:Y=U◣

L5 Lbl 0: A<0＝＞Y=-Y◣Rec(1,Q:U=X:X=C-TV+XV-YW:Y=D-

TW+UW+YV

F6 5 （测设数据输出子程序，步数76）

L1 Pol(X-G,Y-H:W=W-N:W＜0＝＞W=W+360◣W=K+W:W≥360＝＞

W= W-360◣

Fix 4:Prog 2:W〃AC〃=W◢Fix 3:V〃DC〃=V◢Norm

F7 6 （测设时移桩数据计算子程序，步数92）

L1 Lbl 0:Norm:L=0:W=0

L2 Lbl 1:U=0:{U}:U〃SC〃:U≠0＝＞W=W+U:L=L+1:Goto 1◣W≠0＝＞

{U}:U〃V〃:

Prog 1:Fix 3:W〃SD〃=W/L×sin U◢W〃MOVE〃=V-W◢Goto 0◣

【使用说明】

1、 程序说明

(1) 本程序特点：可置镜任意点放样任意点（中桩点、边桩点、导线

点）；人性化设计，提示信息全面，并充分考虑工程术语和习惯；专门设计的角

度输入输出方式；绝对优秀的存储器分配方案。

(2) 本程序适用于放射螺旋线作为缓和曲线且前、后缓和曲线长度相同的

单曲线放样，可放样线路中桩和边桩，按里程和坐标放样均可。

(3) 角度输入（出）方式：度.分秒，如 56°3′19〃 输入（出）为

56.0319。注意：分、秒必须为两位数，当小于10时，应在前面加“0”。

(4) 放样前应准备的数据：建立统一的施工坐标系（正测量坐标系、右手

坐标系），计算控制点和曲线交点坐标、曲线第一切线方位角；整理曲线要素和

放样点桩号、相应的边桩距离（或放样点坐标）。

(5) 若按坐标放样，则当程序提示输入曲线要素时，可输入任意值，仅在

提示输入置镜点、后视点和放样点属性时输入“1”，然后输入相应坐标值。

(6) 如果要实时显示放样点的坐标，可在子程序F6行L1的“Fix 4:”后增加以

下语句(双引号内)：“X◢Y◢”或者“X〃X〃=X◢Y〃Y〃=Y◢”。

(7) 子程序F1行L4中“:U=0:{U}:U〃AB〃:Prog 1:K=U”（输入仪器后视读

数）和行L5中“Prog 6:”（调用移桩计算子程序）、子程序F6行L1中“W=K+W:W

≥360＝＞W= W-360◣”（计算仪器前视读数，必须与前述“输入仪器后视读数”

同时取舍）、子程序F8（移桩计算子程序）可根据需要取舍。

(8) 本程序按CASIO fx-4500P计算器设计。若使用CASIO fx-4800P计算

器，应作以下改动：①变量“V”与“I”调换，变量“W”与“J”调换；②调用子程序时，

子程序名应加双引号（〃），如：PROG 1应改为PROG 〃1〃；③“xy”应为“＾”；④

在子程序F1行L1或L2中任意位置添加“I=0:J=0”语句。

(9) 程序步数说明：每个子程序实际占用空间=程序步数+行数+程序名字

符数+2。

by 中铁二十四局福建公司工程部•洪淮斌

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” �6�8

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时输入已知后视点坐标，其它同ZBJS的输入与显示。

3、由一个已知点坐标和一个已知方位角进行放射型导线坐标计算（程序名FSXZBJS）

程序清单：

C“X0”：D“Y0”：E“T0”：Fixm：Lb1 1:{SA}�6�8

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” �6�8

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}�6�8

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�6�8

Lb1 1：F“X0=”：G“Y0=”：H“T H-0=”：Goto 3�6�8

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�6�8

Lb1 3: X=F：Y=G：T=H ： V=1=>Goto 4:≠>Goto 5�6�8

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}�6�8

Pol(X-F，Y-G) �6�8

J〈0 =〉J=J+360⊿A=J-H�6�8

A〈0 =〉A=A+360⊿A“A=” ◢

I“S=” ◢

Goto 1

程序运行说明：

变量输入参照坐标计算，A=为放样顺时针角度，S=为放样距离。

7、已知置镜点和后视点坐标进行放样（FY2）

程序清单：

F“X0”：G“Y0”：M“XH”：N“YH”�6�8

Pol(M-F，N-G)：Fixm：H=J�6�8

H〈0=〉H=H+360⊿H“T-H=”◢

I“S-H=” ◢

Lb1 1:{XY}: Pol(X-F，Y-G) �6�8

J〈0 =〉J=J+360⊿�6�8

A=J-H�6�8

A〈0 =〉A=A+360⊿A“A=”◢

I“S=” ◢

Goto 1

8、放样优化程序（可分已知一点或者两点坐标进行放样角度距离的计算）

程序名“FY”

程序清单：

I“I=1 =>1-D≠>2-D”�6�8

I=1=>Goto 1:≠>Goto2�6�8

Lb1 1：F“X0”：G“Y0”：H“T0-H”：Goto 3�6�8

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�6�8

Lb1 3：{XY}�6�8

Pol(X-F，Y-G) ：Fixm �6�8

J〈0 =〉J=J+360⊿A=J-H�6�8

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）

程序清单：

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）” �6�8

Lbl 1:Fixm:{KD}:K“LICHENG=”：D“WAIYI=（Z-，Y+）”：Goto2�6�8

Lbl 2:K Goto3:≠>Goto4�6�8

Lbl 3:“K—OVER”: Goto1�6�8

Lbl 4:K>Q+L =>Goto3:≠>Goto5�6�8

Lbl 5:K≤Q+B =>Goto6:≠>Goto7�6�8

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) ◢

Goto1�6�8

Lbl 7：K≤Q+L-C =>Goto8:≠>Goto9�6�8

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�6�8

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�6�8

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�6�8

子程序“BQXYS”

A:R:B“L01”:C“L02” �6�8

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)

程序说明：本程序可以计算等长和不等长“缓+圆+缓”型曲线或者单圆曲线的中线坐标、外移桩坐标以及中线的切线方位角。

先将主、子程序存入计算器，然后运行主程序“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=”�6�8

Lbl 1：{FQ}：F“1÷RB, Z-，Y+”：Q“KB=” �6�8

Lbl 2：{V}：V“K=”：O=V：O

Goto3：≠>Goto4�6�8

Lbl 3：“K-SMALL”：Goto2�6�8

Lbl 4：O>Q =>Goto5：≠>Goto7�6�8

Lbl 5：{I}：I“I=1=>ESC≠>NEXT”：I=1=>Goto2：≠>Goto 6�6�8

Lbl 6：“JISUAN：A-TXY”：V=Q：Z=0：D=0：Goto8�6�8

Lbl 7: {DZ}：Z“JIA-JIAO=”：D“WAIYI=（Z-，Y+）”�6�8

Lbl 8：T=C+90÷Л×（2E+（F-E）（V-P）÷（Q-P））（V-P）：T〈0=〉T=T+360⊿ T〉360 =〉T=T-360⊿�6�8

Rad：W=A+∫（cos（C°+0.5×(2E+(F-E)(X-P)÷(Q-P))(X-P)）,P,V,4）+ D×cos(T°+Z°) �6�8

Y=B+∫（sin（C°+0.5×(2E+(F-E)(X-P)÷(Q-P))(X-P)）,P,V,4）+ D×sin(T°+Z°)：Deg：Goto9�6�8

Lbl 9：O≤Q =>GotoB：≠>GotoA�6�8

Lbl A：A=W：B=Y：C=T：E=F：P=Q： Goto1�6�8

Lbl B： W“X=”◢

Y“Y=”◢

T“T=”◢

Pol(R-G，S-H)：Fixm：M=J：Pol(W-G，Y-H) ：J=J-M�6�8

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输入一个大于第一曲线终点里程的里程值，这时屏幕显示“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�6�8

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

程序按后-前-前-后的观测顺序编制，当后尺为4687时I赋值1，否则赋不等于1的值。

15、高斯抵偿高程投影面边长改正计算FX-4850程序

E“YI-CHANG：hm=”：F“DI-CHANG：Hm=”�6�8

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�6�8

“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时，为计算隧道平均高程面的投影改正。

---