2024年4月8日发(作者：巫马濡霈)

fx-5800P

易学易用程序集

卡西欧(中国)贸易有限公司

目录

序 内容

1

2

3

4

5

目录

页码

第1页

易学易用程序集

第3页

公路高铁程序

编程基础

第8页

第25页

实用用户公式集

第26页

第28页

第31页

第34页

第37页

6 基本操作

7 常见问题

8

9

编程技巧

数学功能

1

fx-5800P G2011-9-30

版程序列表

子程序

功能说明

极坐标放样程序（已知坐标计算边长、方位角计算程序）

已知边长方位角计算坐标程

建筑轴线偏移程序(或扩大基础程序)

测角前方交会坐标计算

测角后方交会坐标计算

坐标转换计算程序

直线中边桩坐标计算程序

圆曲线中边桩坐标计算程序

完整缓和曲线中边桩坐标计算程序

线元法匝道连续计算程序（万能程序）

交点法曲线连续计算程序

竖曲线高程连续计算程序

序 主程序

1

2

3

4

5

6

7

8

9

10

11

12

G1

G2

G3

G4

G5

G6

G7

G8

G9

G1

0

G1

1

G1

2

G1

0-1

；

G1

0-2

G1

1-1

；

G1

2-1

备注

序 参考书目作者

1 王中伟教授

参考书目名称

《道路施工放样综合程序》，王中伟教授在QQ（595077）空间把自己的程序电

子版贡献给广大用户，在此特别感谢！

《公路与铁路施工测量程序》。

《

fx-5800P

说明书》。

其他可能参考的书籍。

特别感谢广西路桥的李虎（东华理工学院97届毕业），他无私奉献了他的公路

高铁测量程序。

本书中程序仅供用户免费参考，卡西欧（中国）贸易有限公司不对用户使用本

书中程序发生的任何问题负责。

用户对用本书疑问，请发邮件到：gongchengji@。欢迎交流！

卡西欧的客服热线：400-700-6655。

2 覃辉教授

3 卡西欧

4

...

5 特别感谢

6

7

8

9

免责申明

联系方式

客服热线

程序传输本书的程序用户可以到全国省会城市的测绘店免费传输G版程序。

如果您有非常好的实用程序，希望广大用户能免费得到分享，请和本书编者

联系:gongchengji@

。

10 程序征集

本书欢迎翻印；希望更多的用户能得到，欢迎流通。本书电子稿及其他的测量

11 欢迎翻印

程序源代码请到卡西欧网站下载: /

欧易学易用程序集》下载。

客户支持

﹥下载中心﹥计算器﹥工程资料下载，或者到百度文库中搜索《fx-5800P卡西

易学易用程序集

 常用测量程序

1极坐标放样程序（已知坐标计算边长、方位角计算程序），程序名

1

Deg：Fix 3 

2

"XO="?A："YO="?B 

3

"XN="?C："YN="?D 

4

Pol((C–A),(D–B)) 

5

J<0J+360→J 

设置角度单位为十进制，3 位固定小数显示

提示输入起点的坐标（A,B）

提示输入端点的坐标（C,D）

调用极坐标函数计算距离和角度,自动保存到I和J

把小于0的角度换算成方位角

显示水平距离I

显示方位角J

6

"JL="：I

◢

7

"FWJ="：JDMS

◢

2已知边长方位角计算坐标程序，程序名

1

Deg：Fix 3 

2

"XO="?A："YO="?B 

3

"L="?L："C="?C 

4

A+Lcos(C) → X 

5

B+Lsin(C) → Y 

6

"X="：X

◢

7

"Y="：Y

◢

设置角度单位为十进制，3 位固定小数显示

提示输入起点的坐标（A,B）

提示输入所求点的距离L和方位角C

计算所求点的X坐标

计算所求点的Y坐标

显示所求点的X坐标

显示所求点的Y坐标

3建筑轴线偏移程序(或扩大基础程序)。程序名

1

Deg：Fix 3 

2

?A:?B:?C:?D: ?E 

3

A+Dcos(C) + Ecos(C-90) → X 

4

B+D sin(C) + Esin (C-90) → Y 

"X="：X

◢

设置角度单位为十进制，3 位固定小数显示

提示输入起点的坐标（A,B）、起点方位角C、轴线偏移距D,E.

计算所求点的X坐标

计算所求点的Y坐标

显示所求点的X坐标

显示所求点的Y坐标

（）

6

"Y="：Y

◢

7

（）

3

4测角前方交会坐标计算，程序名

1

Deg：Fix 3 

2

?A:?B:?C:?D: ?E : ?F

设置角度单位为十进制，3 位固定小数显示

提示输入两个已知点的坐标（A,B）, （C,D），观测的两

个角度E和F

计算所求点的X坐标

计算所求点的Y坐标

显示所求点的X坐标

显示所求点的Y坐标

（）

-1-1-1-1

3

(A(tanF))+C(tanE)+(D-B))÷((tanE)+(tanF)) → X 

-1-1-1-1

4

(B(tanF))+D(tanE)+(A-C))÷((tanE)+(tanF)) → Y 

5

"XP="：X

◢

6

"YP="：Y

◢

7

（）

（）

5测角后方交会坐标计算，程序名

1

2

3

4

5

6

7

8

9

Deg：Fix 3：11→DimZ 

?A:?B:?C:?D: ?E : ?F:?G: ?H : ?K 

Pol((C–A),(D–B)) : J<0

Pol((E–A),(F–B)) : J<0

Z[4]- Z[5]→Z[1] 

Z[1]<0 Z[1]+360→Z[1] 

J+360→J : J→Z[4] 

J+360→J : J→Z[5] 

J+360→J : J→Z[4] 

J+360→J : J→Z[5] 

设置角度单位为十进制，3位固定小数显示

提示输入两个已知点的坐标（A,B）, （C,D），（E,F）

观测的角度G,H,K

计算点（

计算点（

）的方位角

）的方位角

方位角计算顶角

计算点（

计算点（

计算点（

）的顶角并存贮到变量

）的方位角

）的方位角

Pol((E–C),(F–D)): J<0

Pol((A–C),(B–D)): J<0

Z[4]-Z[5]→Z[2]  方位角计算顶角

计算点（

计算点（

计算点（

）的顶角并存贮到变量

）的方位角

）的方位角

10

Z[2]<0Z[2]+360→Z[2] 

J+360→J :J→Z[4] 

J+360→J : J→Z[5] 

11

Pol((A–E),(B–F)) : J<0

12

Pol((C–E),(D–F)) : J<0

13

Z[4]-Z[5]→Z[3] 

14

Z[3]<0

方位角计算顶角

计算点（

计算点（

计算点（

计算点（

）的顶角并存贮到变量

）的参数

）的参数

）的参数

Z[3]+360→Z[3]

15

tan(Z[1]) tan(G) ÷(tan(Z[1])- tan(G))→Z[7]

16

tan(Z[2]) tan(H) ÷(tan(Z[2])- tan(H))→Z[8]

17

tan(Z[3]) tan(K) ÷(tan(Z[3])- tan(K))→Z[9]

18

(Z[7]A + Z[8]C+ Z[9]E) ÷(Z[7]+ Z[8]+ Z[9])→Z[10] 

20

"XP="：Z[10]

◢

21

"YP="：Z[11]

计算所求点的X坐标

19

(Z[7]B + Z[8]D+ Z[9]F) ÷(Z[7]+ Z[8]+ Z[9])→Z[11]  计算所求点的Y坐标

显示所求点的X坐标

显示所求点的Y坐标

4

22

（）

（）

（）

6坐标转换计算程序，程序名

1

Deg: Fix 3 

2

"X0"? A: "Y0"? B 

3

"ANGLE"?E 

4

Cls: "LEFT(1), RIGHT(ELSE) "?F 

5

"XPP"? C: "YPP"? D 

6

[[A] [B]] →Mat A 

7

[[cos(E), -sin(E)] [sin(E), cos(E)]] →Mat B 

8

F≠1-D→D 

9

[[C] [D]] →Mat C 

10

Mat A+Mat B×Mat C→Mat D

◢

两平面坐标系转换计算示意图：

11

设置角度单位为十进制，3位固定小数显示

输入施工坐标系原点在统一坐标系中的X、Y坐标（A,B）

输入统一坐标系的X轴顺时针旋转至施工坐标系的X轴的角度

值

判断施工坐标系是左手系还是右手系

待换算平面点在施工坐标系中的X、Y坐标（C,D）

向矩阵A赋值

向矩阵B赋值

根据施工坐标系的类型决定变量D的取值

向矩阵C赋值

计算转换后的坐标值，并将结果显示在矩阵D中

5

7直线中边桩坐标计算程序,程序名

1 Deg：Fix 3 

2

"X0="?X: "Y0="?Y: "K0="?E: "FWJ="?J 

3

Lbl 1

4

?K:?L:?A 

5

X+(K-E)cos( J ) +Lcos(J+A) → C 

6

"XP=": C

◢

7

Y+(K-E)Sin( J ) +LSin(J+A) → D 

8

"YP=": D

◢

9

Goto 1 

10 备注：该程序请从起点往终点计算。

设置角度单位为十进制，3 位固定小数显示

输入直线起点坐标(X,Y)，起点桩号E，起点方位角J

设置程序标签

提示输入所求点桩号K，边距L，和中线夹角A

计算所求边桩点的X坐标

显示所求边桩点的X坐标

计算所求边桩点的Y坐标

显示所求边桩点的Y坐标

转向程序标签1继续执行程序

8圆曲线中边桩坐标计算程序,程序名

1

Deg：Fix 3 

设置角度单位为十进制，3 位固定小数显示

输入圆曲线起点坐标（A,B），起点桩号K，起点切线方

位角C,半径R（左负右正）

设置程序标签

提示输入所求点桩号L，

提示输入所求点边距E和以及夹角N（左负右正）

计算圆心角的一半，并把弧度转换成角度值。的输入

[FUNCTION] [5](ANGLE) [2]()

计算所求点的弦长，并取绝对值。

计算所求边桩点的X坐标

计算所求边桩点的Y坐标

显示所求边桩点的X坐标

显示 所求边桩点的Y坐标

转向程序标签1循环执行程序

r

r

2

"X="?A:"Y="?B:"FWJ="?C:"KO="?K:"R="?R 

3

Lbl 1

4

"K="?L 

5

"B="?E:"N="?N 

6

(((L-K)÷R)÷2)→M 

7

Abs(2Rsin (M) )→F 

8

A+Fcos (C+M)+Ecos (C+2M+N)→G 

9

B+Fsin (C+M)+Esin (C+2M+N)→H 

10

"XP=":G

◢

11

"YP=":H

◢

12

Goto 1

13 备注：该程序请从起点往终点计算。

r

6

9完整缓和曲线中边桩坐标计算程序,程序名

1

Deg：Fix 3 

设置角度单位为十进制，3 位固定小数显示

输入缓和曲线起点坐标(X,Y)，起点切线方位角C，起点

桩号K，缓和曲线参数A,小半径R（左负右正）

设置程序标签

提示输入所求点桩号L，

提示输入所求点边距E和以及夹角D（左负右正）

计算偏转系数并存储到S

计算所求点的曲线长并存储到P

2

"X0="?X:"Y0="?Y:"FWJ="?C:"KO="?K:"A="?A :"R="?R 

3

Lbl 1

4

"K="?L 

5

"B="?E:"D="?D 

6

Abs (R)÷R→S 

7

L-K→P 

9

P÷(6R)- P

2

∧（

4

）

∧（

3

）

2

∧（

5

）∧（

4

）∧（

7

）∧（

6

）∧（

9

）∧（

8

）

÷ (40 R)+ P÷(3456 R)- P÷(599040R)+ P(÷175472640R)→M 

8

P- P

计算切线支距X坐标并存储到M

÷(336 R

∧（

3

）

)+ P

∧（

6

）

÷(42240 R

∧（

5

）

)- P

∧（

8

）

÷(9676800R

∧（

7

）

)+ P

∧（

10

）

÷(3530096640R

∧（

9

）

)→N 

计算切线支距Y坐标并存储到N

计算弦长并存储到Q

计算弦长对应的偏角值。并存储到T

计算切线夹角并存储到U

计算所求边桩点的X坐标

计算所求边桩点的Y坐标

显示所求边桩点的X坐标

显示 所求边桩点的Y坐标

转向程序标签1循环执行程序

10

M

2

+N

2

→Q

11

cos

-1

（M÷Q)→T



12

(P

2

÷A

2

÷2)

r

→U



13

X+Qcos (C+ST)+Ecos (C+SU+D)→G 

14

Y+Qsin (C+ST)+Esin (C+SU+D)→H 

15

"XP=":G

◢

16

"YP=":H

◢

17

Goto 1

18 备注：该程序只适合完整缓和曲线。

7

公路高铁程序

 fx-5800P程序使用流程：

本程序数据和主程序是分开的，编程时将不同的工程数据存放到不同的数据文件里，如A匝道，文件名

为A，将匝道A所有的曲线线元参数输入A文件里。运行时只要运行文件名A的程序就可以了，具体运行

流程见下图：

运行数据程序文件如ZX、B、C等

输入测站点坐标X0,Y0

第四位小数点不能是0.0001， 如1100.0第四位小数点后为0.0001，如1100.0001

输入计算区域附近桩号

输入计算桩号L

输入所求P点的坐标XP,YP

计算下一点

输入中桩切向右侧夹角M

输入距离所求点到所求桩号中线距

离K，左－右+

输出所计算位置坐标X,Y

输出所计算位置与测站坐标X,Y距

离和方位角DIS,FW

输出P点所在线路的桩号L和距离

K,左－右+

输入桩号L计算下一点

8

数据文件的编写

（一）交点法数据文件编辑

交点法编写数据文件必须是对称型的，即直线段→缓和曲线段→圆曲线段→缓和曲线段→直线段，（如

果任意一端没有直线段，则把直线段长度看做是0），另外圆曲线两侧缓和曲线的旋转常数必须相等，并且

和直线段连接处的半径必须是无穷大。

交点法数据文件编写一般是根据设计图纸提供的平面曲线参数一览表提供的参数来编写，每个弯道包

括：弯道起点方位角（C），交点X坐标（D），交点Y坐标（E），缓和曲线长度（F，当没有设缓和曲线

时，F=0），交点转交（G，向左转弯，G为负值，向右转弯，G取正值），交点桩号（H），弯道圆曲线半

径（R）。 下图是一段市政道路设计参数数据。

根据上图提供的数据，可以编辑成如下的数据文件：文件名（O表示变量字母O，其余的0均为零）

3→DimZ：“X0”?A：”Y0”?B：“Ln”?L：Abs(L)Int(Abs(1000L))÷1000→Z*3+ ：Lbl 0：If Z*3+≠0.0001：Then ?L：E

lse “Xp”?X：”Yp”?Y：X→Z*1+：Y→Z*2+：IfEnd：Lbl 1:If L>0 ：Then 98°39°35.12°→C：4774.384→D：2415

.861→E：140→F：31°17°23°→G：410.007→H：600→R：IfEnd：If L>1060 Then 129°56°58.19°→C：4206.

421→D：3093.946→E：70→F：-33°50°48°→G：1285.437→H:600→R：IfEnd：Prog”XLJS” ：If Z[3+≠0.0001：

Then Goto 0：Else (Z[1]－X)cos(O)＋(Z[2]－Y)sin(O)→N：L＋N→L：Int(Abs(1000L))÷1000+0.0001→L：－(

Z[1]－X)sin(O)+(Z[2]－Y)cos（O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“：L

◢

“K=” ：K

◢

IfEnd：Goto 0: IfEnd

在程式中，有两个条件转移语句即

If L>0：Then 98°39°35.12°→C:4774.384→D: 2415.861→E:140→F:31°17°23°→G:410.007→H:600→R:IfEnd

If L>1060：Then 129°56°58.19°→C:4206.421→D: 3093.946→E:70→F:-33°50°48°→G:1285.437→H:600→R:IfEnd

„„

9

如果还有其他弯道，可以继续后加。在这些存放设计参数的语句程序表达式是固定的。编辑曲线参数时，每

个曲线参数放在一个If L>***.***（两个弯道中间直线段上的任意桩号） :Then ***°**°**.**°→C（弯道

起点方位角）：****.***→D（交点X坐标）：****.****→E（交点Y坐标）：***.***→F（缓和曲线长度：**°

**°**°→G（转角，向左转为负值，向右转为正值）：***.***→H（交点里程桩号）:*** →R（圆曲线半径） ：

IfEnd 条件式语句里。如果有多个弯道，一直按上述形式编下去。

变量说明：

序 语句

1

2

3

说明

L>***.***

***°**°**.**°→C

****.***→D、****.****→E

**°**°**°→G

***.***→H

***→R

弯道参数的起点范围，可以是弯道起点，也可以是弯

道前直线段里的任意点桩号。

弯道起点方位角

交点X坐标、Y坐标

缓和曲线长度，如果没设臵缓和曲线，则缓和曲线长度看做0

转角，向左转为负值，向右转为正值

交点里程桩号

圆曲线半径

4

***.***→F

5

6

7

【注意】运行程序时，第一步是输入测站坐标X0，Y0，然后输入近似桩号，主要是为了加快收敛，

输入时注意当根据桩号和距中距离计算所求点坐标时，请不要把小数点第四位输入0.0001，如

K21+369.0921，请把桩号输入21369.092即可，相反，当根据坐标反算桩号和距中距离时，请在输入

时把第四位小数数位0.0001，如K21+200，输入时请输入21200.0001。

10

（二）线元法数据文件编辑

线元法数据文件编辑适合于各种各样的线型组合，特别是在大型立交区或者是不规则的线型组合地方，

如小区道路、重丘区山路等有特别明显的有点，本程序为了减少使用者设计参数的输入，特别把直线与缓和

曲线以及圆曲线连接的地方只要输入缓和曲线参数就可以了，这样可以减少输入很多设计参数，同时节约大

量的内存。看下图：

该图是一个半喇叭型的互通立交，一共是5条匝道，我们以D匝道为例，来看看线元法数据文件的建立。

11

下表是D匝道的线元表

：

线元

名称

缓和曲线

圆曲线

缓和曲线

起点桩号

DK0+000.00

偏向

右偏

半径或缓和曲

线旋转常数A

55

40

65

线元长度

(m)

39.289

39.216

105.625

线元起点坐标、方位角

X(m)

494341.902

494347.926

494380.677

Y(m)

478027.825

477989. 873

477971.271

方位角

260°38′55.7″

302°18′26.0″

358°28′45.9″

DK0+039.289 右偏

DK0+078.504 右偏

D匝道线元数据分析：从上表中可以看出，D匝道由三个线元组成，第一段是缓和曲线段，根据上表判

断，该缓和曲线属于第一段缓和曲线，其最小半径处桩号为J=K0+039.289，该处坐标X=494347.926,Y=477989.

873，半径为R=40米，方位角W=302°18′26.0″，缓和曲线旋转常数P=55，在缓和曲线终点K0+039.289处，

有一圆曲线与其相连接，长度I=39.216。第二段是圆曲线，第三段又是一段缓和曲线，根据上表判断，该曲线

属于第二缓和曲线，曲线最小半径位于与园曲线连接处，桩号J=K0+78.504，坐标为X=494380.677，Y=477971.271,

方位角W=358°28′45.9″，最小半径R=40，缓和曲线旋转常数P=-65（为了判断是第一段还是第二缓和曲线，

我们通过P的正负来判断，第一缓和段P统一用正数表示，第二缓和段统一用负号表示。起点处相连的圆曲线

长度I可以看做0，也可以看做是39.216.由此我们整理下。从整个线型全部是向右转向，所以R全部取正数。

具体如下：

第一段缓和曲线线元参数说明：

最小半径处坐标

X=494347.926 Y=477989. 873

最小半径R=40

最小半径处方位角：W=302°18′26.0″

缓和曲线旋转常数P=55

与之相连的圆曲线长度I=39.216

最小半径处桩号J=39.289

此段数据计算范围可以从K0+000~ K0+78.504

第二段缓和曲线线元参数

最小半径处坐标

X=494380.677 Y=477971.271

最小半径R=40

最小半径处方位角：W=358°28′45.9″

缓和曲线旋转常数P=65

与之相连的圆曲线长度I=39.216

最小半径处桩号J= W=358°28′45.9″

此段数据计算范围可以从K0+78.504到匝道结束。

当桩号L位于K0+000~ K0+78.504时，曲线参数按第一段编写，当L>78.504时，按第二段缓和曲线编

写，当然第二段缓和曲线计算范围也可以从圆曲线起点K0+039.289到K0+78.504中间任意桩号开始。下面

是D匝道的数据程序文件（文件名“”）：（O表示变量字母O，其余的0均为零）

3→DimZ：“X0”?A: “Y0”?B: “Ln”?L: Abs(L)-Int(Abs(1000L))÷→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L:Else

“Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L≥0 :Then 302°18°26°→W：494347.926→X：477989.

873→Y：55→P: 40→R：39.216→I: 39.289→J：IfEnd: If L≥78.504：Then 494380.677→X：477971.271

→Y：358°28°45.9°→W：40→R： -65→P:39.216→I: 78.504→J：IfEnd:Prog””: Prog””:

If Z[3]≠0.0001:Then Goto 0:Else (Z[1]－X)cos(O)＋(Z[2]－Y)sin(O)→N:L＋N→L:-(Z[1]－X)sin(O)+(Z[2]－Y)c

os（O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“:L

◢

“K=”:K

◢

IfEnd:Goto 0: IfEnd

12

在上面程式中，每个缓和曲线的线元参数放在一个If „ Then „„ IfEnd循环语句中，If 后面紧跟计

算的起点范围，终点范围从下段的起点范围结束。

变量定义：

序 变量

1

X、Y

2

R

3

W

4

P

5

I

6

J

说明

输入缓和曲线半径最小处的X、Y坐标

缓和曲线最小半径，当线路左转时，半径数据前面加“-”号，以判断线路转向。

输入缓和曲线半径最小处的方位角

缓和曲线旋转常数，第一缓和曲线长度为正数，第二缓和曲线常数为负数。当不

设缓和曲线时，P=0

与缓和曲线最小半径端连接的圆曲线长度。

缓和曲线小半径段的桩号。

再看看E匝道数据文件编辑

E匝道线元设计数据表

半径或

缓和曲

线旋转

常数A

120

200

95

75

86.75

线元起点坐标、方位角

线元长

度(m)

71.296

72

60.070

45.125

64.841

17.800

线元名称 起点桩号 偏向

X(m)

494360.619

494373.205

494381.631

494372.527

494355.878

494337.373

Y(m)

477706.293

477776.470

477847.870

477907.018

477948.932

478010.658

方位角

79°49′57.0″

79°49′57.0″

90°08′44.6″

107°21′16.3″

113°49′05.6″

92°24′19.0″

直线

缓和曲线

圆曲线

缓和曲线

缓和曲线

圆曲线

EK0+000.00

EK0+071.296

EK0+143.296

EK0+203.366

EK0+248.491

EK0+313.333

右偏

右偏

右偏

左偏

左偏

13

在匝道线元表中，我们可以看出，E匝道式由一段直线段→缓和曲线段→圆曲线段→缓和曲线段→缓

和曲线段→圆曲线段连接起来的。值得注意的是，在位于K0+248.491处，是两段反向的缓和曲线连接点，

在编写数据程序文件时，分三段编写，第一段从K0+000到K0+203.366，第二段从K0+203.366到K0+248.919，

第三段从K0+248.919到线路终点K0+331.13。具体文件如下：（O表示变量字母O，其余的0均为零）

3→DimZ：“X0”?A:”Y0”?B: “Ln”?L:Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L: Else “Xp

”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L≥0 :Then 494381.6316→X：477847.870→Y：90°08°44.6°

→W：200→R：120→P:60.07→I: 143.296→J：IfEnd: If L≥203.366：Then 494372.527→X：477907.018→Y：

107°21°16.3°→W：200→R：-95→P:60.07→I: 203.366→J：IfEnd: If L≥248.491：Then 494337.373→X：

478010.658→Y：92°24°19°→W：－86.75→R：75→P:17.8→I: 313.333→J：IfEnd: Prog””: Prog”

：If Z[3]≠0.0001:Then Goto 0:Else (Z[1]－X)cos(Ｏ)＋(Z[2]－Y)sin(Ｏ) →N:L＋N→L:-(Z[1]－X)sin(Ｏ”

)+(Z[2]－Y)cos（Ｏ）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“:L

◢

”K=”:K

◢

IfEnd:Goto 0: IfEnd

14

为了让使用者熟悉掌握数据程序文件的编写线元法的数据程序文件，现在将C、F、G等匝道的线元数据

列出来，供使用者练习。C匝道数据；F匝道数据；G匝道数据

缓和曲线 FK0+111.505

FK0+202.434

FK0+303.548

GK0+000.000

GK0+133.152

GK0+282.286

右偏

右偏

右偏

右偏

右偏

右偏

100

110

120

175

230

190

90.909

101.114

130.909

133.152

149.124

149.72

494557.835

494542.812

494595.310

494718.802

494587.600

494472.335

478313.650

478224.689

478142.420

478112.963

478093.368

478002.890

252°32′00.2″

276°12′33.6″

32852′35.7″

182°58′11.8″

199°33′17.5″

236°42′21.1″

直线 FK0+000.00 111.525 494591.311 478420.033 252°32′00.2″

圆曲线

直线

缓和曲线

圆曲线

缓和曲线

直线

缓和曲线

圆曲线

CK0+000.00

CK0+056.876

CK0+078.153

CK0+108.332

CK0+138.763

CK0+168.941

CK0+590.964

CK0+675.964

右偏

右偏

右偏

右偏

右偏

右偏

右偏

32

65

140

65

85

85

85

56.876

21.277

30.179

30.431

30.179

422.023

85

72.739

495005.991

494991. 262

494971. 511

494943.127

494913.167

494883.008

494461.552

494379.485

478056.677

478104. 123

478112.036

478122.242

478127.219

478126.740

478104.874

478086.681

56°19′42.3″

158°09′54.4″

158°09′54.4″

164°20′25.7″

176°47′40.6″

182.°58′11.8″

182.°58′11.8″

211°37′04.2″

半径或

线元名

称

起点桩号 偏向

缓和曲

线旋转

常数A

线元长

度(m)

X(m)

线元起点坐标、方位角

Y(m) 方位角

圆曲线

缓和曲线

缓和曲线

圆曲线

缓和曲线

15

道路计算程序

计算程序是核心程序，但是道路计算程序是固定的，使用者不需要做任何修改，也不直接运行道路计算

程序分为两部分，第一部分是常规的对成型线型组合计算的程序（文件名：

件直接调用；第二部分是适用于非对称的线型组合计算文件（文件名：

文件后，还要调用“

Deg:Fix 4

G÷Abs(G)→V: F÷2-F^(3)÷（240R

2

）+F^（5）÷(34560R^（4）)-F^（7）÷(8386560R^（6）)+F^（9）÷(3158507

2

520R^（8）)→Q:F÷(24R)-F^(4)÷(2688R^(3))+F^(6)÷(506880R^(5))-F^(8)÷(154828800R^(7))→P:Q+(R+P)Vtan(G÷2)

→T:πRGV÷180+F→I:

IF L≤H-T:Then L-H→J:Goto 1:IfEnd:

If L

），交点法编辑的数据文

），在数据文件调用完“”

”文件。程序代码：

（为了输入方便代码分行显示，使用者不必分行输入。）

If L

If L

If L≥H-T+I:Then L-H-I+2T→J:C+G→C:Goto 1:IfEnd:

Lbl 1:

D+Jcos (C)→X:

E+Jsin(C)→Y:

C→O: Goto C:Lbl 2:

(J-J^(5)÷(40(RF)

2

)+J^(9)÷(3456(RF)^(4))-J^(13)÷(599040(RF)^(6))+J^(17)÷(175472640(RF)^(8))-T)S→M:

(J^(3)÷(6RF)-J^(7)÷(336(RF)^(3))+J^(11)÷(42240(RF)^(5))-J^(15)÷(9676800(RF)^(7))+J^(19)÷(3530097000(RF)^(9)))

V→N:

C+90J

2

SV÷(πRF)→O:

Goto B

Lbl 3:

180(J-0.5F) ÷(πR)→O:

Q+Rsin (O)-T→M:(P+R(1-cos (O)))V→N:C+VO→O：

Goto B

Lbl B:

D+Mcos(C)-Nsin(C)→X:

E+Msin(C)+Ncos(C)→Y:

Goto C:

Lbl C:If Z*3+≠0.0001:Then ?K: "M"?U:K<0=>180-U →U:O+U×Abs(K)÷(K+10^(-10)) →Z:X+Abs(K)×cos(Z)→X:

:Y+Abs(K)×sin(Z)→Y:

"X":X

◢

"Y":Y

◢

tan-1((Y-B)÷(X-A+10^(-10))→W: If X-A<0:Then W+180→W: IfEnd: If W<0:Then W+360→W:IfEnd:√((X-A)

2

+(Y-B

2

)

)→S:"DIS": S

◢

"FW": W

DMS

◢

IfEnd

16

非对称线型组合计算程序（）（O表示变量字母O，其余的0均为零）

Abs(P+10^(-10)）÷(P+10^(-10)) →S：Abs (R)÷R→V:Abs(R) →R:P

2

÷R→F:

180(I+F)V÷(πR) →G:90FSV÷(πR) →O:

If S=1:Then W-O→C: IfEnd: If S=-1:Then W-O-G→C :IfEnd:C<0C+360→C:C≥360C-360→C:

F

2

÷(24R)-F^(4)÷(2688R^(3))+F^(6)÷(506880R^(5))-F^(8)÷(154828800R^(7))→P:F÷2-F^(3)÷(240R

2

)+F^5÷(

34560R^(4))-F^(7)÷(8386560R^6)+F^9÷(3158507520R^(8))→Q:

Q+(R+P)Vtan(G÷2)→T:If S=1:Then J-F+T→H:C→Z: Else J-I-F+T→H:C+G→Z:IfEnd:

(F-F^(3)÷(40R

2

)+F^(5)÷(3456R^(4))-F^(7)÷(599040R^(6))+F^(9)÷(175472640R^(8))-T)S→M:(F2÷(6R)-F^(4)

÷(336R^(3))+F^(6)÷(42240R^(5))-F^(8)÷(9676800R^(7))+F^(10)÷(3530097000R^(9)))V→N:

X-Mcos(Z)+Nsin(Z)→D:Y-Msin(Z)-Ncos(Z)→E

最小半径处坐标X=494347.926 Y=477989. 873 最小半径R=40 最小半径处方位角：W=302°18′

26.0″ 缓和曲线旋转常数P=55 与之相连的圆曲线长度I=39.216 最小半径处桩号J=39.289

17

程序运行

程序运行分为两种，第一种是根据根据使用者输入的桩号和距中距离计算出所求点的坐标和测站点到该

点的方位角和水平距离，可以直接利用计算器提供的距离和方位角进行极坐标放样。另一种是根据使用者输

入的实测坐标计算出该点的桩号和距中距离，可以根据桩号和距中距离进行高边坡、隧道断面点的设计高程、

以及进行路面高程控制计算方面很方便。

（一）据桩号和距中距离计算所求点的坐标、及方位角和水平距离。要进行程序的线路计算，直接运行数

程序文件就可以了，不需要理会计算程序。如果计算某桩号L距离中线K处的坐标X，Y和该点到测站点的

距离和方位角。运行提示如下：

序 变量

1

X0?

2

Y0？

3

Ln?

4

L？

5

说明

提示输入测站点X坐标

提示输入测站点Y坐标 （测站点的位数和数据程序文件里坐标取位要一致）

计算点近似桩号，近似桩号输入第四位小数后不能等于0.0001，

输入计算点准确桩号，近似桩号输入第四位小数后不能等于0.0001。

提示输入距离中桩的桩号如果计算点位于线路的左幅，则在数据前加入“-”号，如

果是计算中桩坐标，则输入0

K？

6

7

M？

X

FW

输入该点与线路的前进方向的按顺时针转过的角度，常于斜交的涵洞，桥梁等构造

物

显示所求点的X坐标

显示所求点的Y坐标

显示所求点至测站点的方位角

显示所求点至测站点的距离

8

Y

9

10

DIS

要继续计算则从L?开始，计算新点的坐标，具体流程见程序流程图的左边部分。例：要计算D匝道的

DK0+020中桩坐标、DK0+080左边5.3米的坐标、Dk0+240斜交75°58′11.8″、距中右侧5米处的坐标，

并计算出该点至坐标为X=495005.991，Y=478056.677控制点的水平距离和方位角。

18

运行文件“

序 变量

1 X0?

2 Y0？

3 Ln?

4 L？

5

6

7

K？

M？

X

说明

495005.991 测站点X坐标

478056.677 测站点Y坐标

”提示如下：

0 输入小输点第四位起不等于0.0001的桩号

20 输入所求点桩号

0 输入距中距离

90 输入斜交角度

494341.482 显示所求点X坐标

478007.903 显示所求点Y坐标

666.2959 显示所求点至测站平距

184°11′52.42″显示所求点至测站方位角

80 输入下点的桩号

-5.3

90

494382.2290

477965.9594

630.324

188°16′29.5″

240

5.0

75°58°11.8°

494455.3110

478101.5491

552.505

175°20′29.6″

8 Y

9 DIS

10 FW

11 L？

12 K？

13 M？

14 X=

15 Y=

16 DIS

17 FW

18 L？

19 K？

20 M？

21 X=

22 Y=

23 DIS

24 FW

19

（二）据坐标反算桩号和距中，根据某点的坐标，计算出该点在线路中的桩号和距离，程序运行时提示如下：

序 变量

1 X0?

2 Y0？

3 Ln?

4 Xp？

5

6

7

Yp？

L

K

说明

提示输入测站点X坐标

提示输入测站点Y坐标（测站点的位数和数据程序文件里坐标取位要一致）

计算点近似桩号，近似桩号输入第四位小数后不能等于0.0001，

所求点P的X坐标

所求点P的Y坐标

显示所求点的桩号

显示所求点的距中距离，如果位于线路的左边，则显示结果为负值，如果K为0时，则该点在线路

中线上。

计算新点的桩号和距中距离，则从?Xp 开始。如此循环计算。具体流程见流程图的右边部分

例：求点X=494382.2290，Y=477965.9594在D匝道中的位置，运行文件“

序 变量

1 X0?

2 Y0？

3 Ln?

4 Xp？

5

6

7

Yp？

L

K

说明

495005.991 测站X坐标(可输入0)

478056.677 测站Y坐标(可输入0)

60.0001 输入小输点第四位起等于0.0001的桩号

494382.2290 输入所求点X坐标

477965.9594 输入所求点Y坐标

79.999999 显示所求点桩号

-5.3000 显示所求点距中距离, “-“ 表示该点位于中桩左侧

”如下：

20

程序应用

（一）本程序在高边坡放样中的应用

高边坡放样是目前公路道路测量工作普遍要做的一项工作，边坡坡顶线或者填方边坡坡脚线是边坡放样

最繁琐的工作，利用本程序进行边坡放样，工作步骤如下：

① 司镜员走到目的地后，先测出一个点的坐标和高程；

② 利用本程序求出该点在线路中的桩号和距中距离，然后根据计算出来的桩号和设计图纸，内插出道

路边桩高程，根据道路边桩高程和道路设计宽度、坡度、设计碎落平台的宽度和间隔高度、实测点

的高程计算出按实测高程，边坡距中距离。

③ 根据高程算出的距中距离和程序计算出立镜点至距中距离差值dK，指挥司镜员靠近道路中线或者

往外偏移多少米。

④ 重新做1~3工作，直到dK 小于边坡宽度控制精度要求为止。

dK

m

h

m

h

K1 K2

如上图，司镜员在P点立镜，测出P点坐标和高程后，用本程序计算出P点到的中桩桩号和距离

中桩距离K，根据中桩桩号和有关设计数据内插出边坡坡脚点的设计高程Hs和距离中桩的距离K1，

根据Hc、和Hs的高差计算出K2，按照P点的高程，边坡顶距中桩的理论宽度K

理论

，具体计算如下：

K

理论

=K1+h*n1+m+h*n2+m+（Hc-Hs-2h）*n3

则dK=K- K

理论

dK 求出来后，可以指挥司镜员靠近道路dK米，然后重复以上步骤，直到dK值能满足边坡控制精度。

21

（二）本程序在隧道断面欠挖、超挖放样中的应用。利用本程序进行隧道断面开挖的欠挖和超挖计算，看下图：

R2'

R2

用全站仪测出P点的坐标和标高后，根据坐标算出P点的中桩桩号和距中距离K，根据中桩桩号查看P

点的路面高程，算出P点至设计路面的高差h，根据P点所在圆弧的圆心O2在断面中的位置（K2，h2）算

出P点到O2点距离R2’，用R2’-R2的差距就是P点的超挖或者欠挖的数据。

（三）本程序在路面施工控制放样中的应用

路面施工时，路面的轴线偏位、路面宽度、路面高程是测量的重要控制指标，也是道路施工控制最严格

的一道测量控制。利用本程序输入桩号和距中距离，可以直接计算出所求点至测站点的水平距离和方位角，

从而直接通过极坐标法直接放出该点在实地的位置。

（四）处理线路设计中的“断链”问题

断链是线路设计中经常出现的问题，断链分为长链与短链，如果是短链，那么只要弯道的计算范围从短

链处分界就可以了，如果是长链，则需要把线路分成两个文件，从开始出现断链处分，这样可以避免桩号相

同，但位置不同的情况。

（五）本程序在道路竣工验收过程中的应用

道路竣工验收，其主要方面就是轴线偏位和宽度，可以通过现场测量道路中线的坐标，通过坐标反算出

该点和设计中线的偏移值K（轴线偏位值）。

22

公路纵断面高程计算

公路纵断面计算也是采用计算程序和数据程序分开的原则，数据文件是专门存放道路纵断面设计参数，

对于多个工程，可以建立多个数据文件，这样每个数据文件通过调用计算程序就可以计算出该点的高程。

① 数据文件：每个纵断面由纵坡和竖曲线组成，本程序每一个变坡点计算范围是从竖曲线前的纵坡任意桩

号开始到竖曲线后纵坡上的任一桩号结束。下表是某线路路基纵断面设计的参数一览表：

路基纵断面设计

桩号

K037+950.00

K038+650.00

K038+900.00

K039+245.00

K039+410.00

K039+590.00

变坡点高程

80.56

75.93

86.01

87.75

91.86

93.21

99.84

半径

3000.00

4500.00

6000.00

6000.00

10000.00

切线长

70.40

79.372

59.597

52.227

52.214

竖曲线起点桩号 竖曲线终点桩号

K038+579.599

K038+820.628

K039+185.403

K039+357.773

K039+532.786

K038+720.401

K038+979.372

K039+304.597

K039+462.227

K039+647.214

纵坡

-0.6614%

4.0320%

0.5043%

2.4909%

0.7500%

1.8943%

根据上表编程数据文件如下：文件名

"Dh"?E:Lbl 0:"L="?L:"I="?I:"K="?K 

If L≥37950:Then 75.93→C:38650→D:-0.006614→M:0.04032→N:3000→R:IfEnd 

If L≥38820.628:Then 86.01→C:38900→D:0.04032→M:0.005043→N:4500→R:IfEnd 

If L≥39185.403:Then 87.75→C:39245→D:0.005043→M:0.024909→N:6000→R:IfEnd 

If L≥39357.773:Then 91.86→C:39410→D:0.024909→M:0.0075→N:6000→R:IfEnd 

If L≥39532.786:Then 93.21→C:39590→D:0.0075→M:0.018943→N:10000→R:IfEnd 

Prog "":Goto 0

从上面程序看，"Dh"?E:Lbl 0:"L="?L:"I="?I:"K"?K 语句是固定的，“Dh”提醒使用者输入所求点距离路面

的厚度，比如路面底基层距离路面厚度为37cm，则输入0.37，表示所求点是路面底基层顶面高程，“L=”，

输入所求点桩号,“I”=表示输入路面横坡值，“K=”表示输入路面某点至中桩距离，均为正值。以上数据均

在程序运行时输入。在每一个变坡点处，要求输入的参数有：

序 变量

1 C；

2 D：

3 M：

4 N：

5 R：

说明

输入变坡点高程

输入边坡点桩号

输入变坡点小里程段纵坡坡度

输入变坡点大里程侧纵坡坡度

竖曲线半径

L表示竖曲线前纵坡任意里程的桩号，本程序每段竖曲线参数计算范围是从竖曲线前的纵坡上任意点到竖曲

线后纵坡上任意点的桩号。

23

② 高程计算程序（文件名：

③ 运行程序显示如下

说明

），高程计算文件代码如下：

N-M→P: Abs(P)÷P→S:(Abs(P)×R)÷2→T:If L＜D-T:Then C+(L-D)M→H:IfEnd 

If L＞D-T:Then C-MT+(M+(L-D+T)S÷(2R))×(L-D+T)→H:IfEnd 

If L≥D+T:Then C+(L-D)N→H:IfEnd 

H-E-KI→H:"H=":H

◢

序 变量

1 Dh？

2 L=?

3 i？

4 K？

5 H=

提示输入至路面厚度

提示输入桩号

提示输入横坡度

提示输入距中距离，均为正数

提示输出所求点高程

例：根据上表计算K38+660中桩高程、K38+900右侧4.5米高程、K39+600底基层面左侧4.5米高程（假设

路面沥青厚度为6cm，封釉层1cm，水泥稳定层基层厚度为30cm），运行

序 变量

1 Dh？

2 L=？

3 i？

4 K？

5

6

7

8

9

H=

L=？

i？

K？

H=

说明

0

38660

0

0

76.9412

38900

0.02

4.50

85.2200

文件显示如下：

退出程序，重新运行。

序 变量

1 Dh？

2 L=？

3 i？

4 K？

5 H=

说明

0.370

39600

0.02

4.50

93.0509

24

编程基础

 fx-5800P的编程

采用卡西欧三步编程法写出您自己需要的程序。步骤如下：

① 按【MODE】【5】(PROG)【1】(NEW)(输入文件名如P1)后【EXE】【1】(COMP)

② 按以下步骤输入程序源代码。

序 程序源代码

1

?A

:

?B

:

?C

:

?D: ?E 

2

A+Dcos(C) + Ecos(C-90) → X 

3

B+D sin(C) + Esin (C-90) → Y 

4

"X="：X

◢

5

"Y="：Y

◢

编程思路说明

第①步：

提示用户输入已知数据,

语法『?+变量+:（断开）』

 是换行符号。

第②步：

数据计算处理，

语法『直接输入公式，→+变量（存储结果数据）』

第③步：成果数据显示在屏幕上。

语法『"变量="+：（断开）+变量+

◢

』

③ 按【FILE】键再按

调用程序P1。

在案例中，输入公式的部分可以是用户的任意公式，用户只要直接输入公式（公式是什么样就

输入什么样，没有任何语法）。表明卡西欧计算器编程操作简单，功能强大。

 fx-5800P的公式自定义

没有语法,变量可带脚标，α β θ等符号直接作为变量使用；

请直接写出您自己需要的公式并调用。步骤如下：

① 按【MODE】【5】(PROG)【1】(NEW) (输入文件名如F1)后【EXE】【3】(Formula)

α

◢

α

② 公式编写好后，按【EXIT】退出到程序菜单。再按[MODE] [1]回到普通计算界面，按[FMLA]

[1] (Original)移动光标选择程序名如F1，调用程序。

【注意】此操作过程常见问题：

1：脚标的输入是[FUNCTION] [4](ALPHA) [5]

[7]。

2：在MATH状态下按[FMLA] [1] (Original)直接调用自定义公式。

3：以上公式是已知直线起点坐标,里程求算直线上中桩坐标。

4：从以上公式可以看出，公式自定义可以在一个文件中连续输入多个公式，只要用

◢

断开即可。

25

实用用户公式集

 实用公式集

① 新建公式按【MODE】【5】(PROG)【1】(NEW) (输入文件名如F1)后【EXE】【3】(Formula)

② 公式运行在普通计算界面，按[FMLA] [1] (Original)移动光标选择程序名如F1运行。

1

测角前方交会坐标计算公式

β

β

2

α

α

α

α

β

β

◢

测角后方交会坐标计算公式

◢

其中

、、

需要另外建个公式文件

α

◢

θ

θ

α

◢

αθ

◢

β

θβθ

γγ

α

θ

3 直线中边桩计算坐标公式（也是坐标反算公式）

4 坐标正算公式

2024年4月8日发(作者：巫马濡霈)

fx-5800P

易学易用程序集

卡西欧(中国)贸易有限公司

目录

序 内容

1

2

3

4

5

目录

页码

第1页

易学易用程序集

第3页

公路高铁程序

编程基础

第8页

第25页

实用用户公式集

第26页

第28页

第31页

第34页

第37页

6 基本操作

7 常见问题

8

9

编程技巧

数学功能

1

fx-5800P G2011-9-30

版程序列表

子程序

功能说明

极坐标放样程序（已知坐标计算边长、方位角计算程序）

已知边长方位角计算坐标程

建筑轴线偏移程序(或扩大基础程序)

测角前方交会坐标计算

测角后方交会坐标计算

坐标转换计算程序

直线中边桩坐标计算程序

圆曲线中边桩坐标计算程序

完整缓和曲线中边桩坐标计算程序

线元法匝道连续计算程序（万能程序）

交点法曲线连续计算程序

竖曲线高程连续计算程序

序 主程序

1

2

3

4

5

6

7

8

9

10

11

12

G1

G2

G3

G4

G5

G6

G7

G8

G9

G1

0

G1

1

G1

2

G1

0-1

；

G1

0-2

G1

1-1

；

G1

2-1

备注

序 参考书目作者

1 王中伟教授

参考书目名称

《道路施工放样综合程序》，王中伟教授在QQ（595077）空间把自己的程序电

子版贡献给广大用户，在此特别感谢！

《公路与铁路施工测量程序》。

《

fx-5800P

说明书》。

其他可能参考的书籍。

特别感谢广西路桥的李虎（东华理工学院97届毕业），他无私奉献了他的公路

高铁测量程序。

本书中程序仅供用户免费参考，卡西欧（中国）贸易有限公司不对用户使用本

书中程序发生的任何问题负责。

用户对用本书疑问，请发邮件到：gongchengji@。欢迎交流！

卡西欧的客服热线：400-700-6655。

2 覃辉教授

3 卡西欧

4

...

5 特别感谢

6

7

8

9

免责申明

联系方式

客服热线

程序传输本书的程序用户可以到全国省会城市的测绘店免费传输G版程序。

如果您有非常好的实用程序，希望广大用户能免费得到分享，请和本书编者

联系:gongchengji@

。

10 程序征集

本书欢迎翻印；希望更多的用户能得到，欢迎流通。本书电子稿及其他的测量

11 欢迎翻印

程序源代码请到卡西欧网站下载: /

欧易学易用程序集》下载。

客户支持

﹥下载中心﹥计算器﹥工程资料下载，或者到百度文库中搜索《fx-5800P卡西

易学易用程序集

 常用测量程序

1极坐标放样程序（已知坐标计算边长、方位角计算程序），程序名

1

Deg：Fix 3 

2

"XO="?A："YO="?B 

3

"XN="?C："YN="?D 

4

Pol((C–A),(D–B)) 

5

J<0J+360→J 

设置角度单位为十进制，3 位固定小数显示

提示输入起点的坐标（A,B）

提示输入端点的坐标（C,D）

调用极坐标函数计算距离和角度,自动保存到I和J

把小于0的角度换算成方位角

显示水平距离I

显示方位角J

6

"JL="：I

◢

7

"FWJ="：JDMS

◢

2已知边长方位角计算坐标程序，程序名

1

Deg：Fix 3 

2

"XO="?A："YO="?B 

3

"L="?L："C="?C 

4

A+Lcos(C) → X 

5

B+Lsin(C) → Y 

6

"X="：X

◢

7

"Y="：Y

◢

设置角度单位为十进制，3 位固定小数显示

提示输入起点的坐标（A,B）

提示输入所求点的距离L和方位角C

计算所求点的X坐标

计算所求点的Y坐标

显示所求点的X坐标

显示所求点的Y坐标

3建筑轴线偏移程序(或扩大基础程序)。程序名

1

Deg：Fix 3 

2

?A:?B:?C:?D: ?E 

3

A+Dcos(C) + Ecos(C-90) → X 

4

B+D sin(C) + Esin (C-90) → Y 

"X="：X

◢

设置角度单位为十进制，3 位固定小数显示

提示输入起点的坐标（A,B）、起点方位角C、轴线偏移距D,E.

计算所求点的X坐标

计算所求点的Y坐标

显示所求点的X坐标

显示所求点的Y坐标

（）

6

"Y="：Y

◢

7

（）

3

4测角前方交会坐标计算，程序名

1

Deg：Fix 3 

2

?A:?B:?C:?D: ?E : ?F

设置角度单位为十进制，3 位固定小数显示

提示输入两个已知点的坐标（A,B）, （C,D），观测的两

个角度E和F

计算所求点的X坐标

计算所求点的Y坐标

显示所求点的X坐标

显示所求点的Y坐标

（）

-1-1-1-1

3

(A(tanF))+C(tanE)+(D-B))÷((tanE)+(tanF)) → X 

-1-1-1-1

4

(B(tanF))+D(tanE)+(A-C))÷((tanE)+(tanF)) → Y 

5

"XP="：X

◢

6

"YP="：Y

◢

7

（）

（）

5测角后方交会坐标计算，程序名

1

2

3

4

5

6

7

8

9

Deg：Fix 3：11→DimZ 

?A:?B:?C:?D: ?E : ?F:?G: ?H : ?K 

Pol((C–A),(D–B)) : J<0

Pol((E–A),(F–B)) : J<0

Z[4]- Z[5]→Z[1] 

Z[1]<0 Z[1]+360→Z[1] 

J+360→J : J→Z[4] 

J+360→J : J→Z[5] 

J+360→J : J→Z[4] 

J+360→J : J→Z[5] 

设置角度单位为十进制，3位固定小数显示

提示输入两个已知点的坐标（A,B）, （C,D），（E,F）

观测的角度G,H,K

计算点（

计算点（

）的方位角

）的方位角

方位角计算顶角

计算点（

计算点（

计算点（

）的顶角并存贮到变量

）的方位角

）的方位角

Pol((E–C),(F–D)): J<0

Pol((A–C),(B–D)): J<0

Z[4]-Z[5]→Z[2]  方位角计算顶角

计算点（

计算点（

计算点（

）的顶角并存贮到变量

）的方位角

）的方位角

10

Z[2]<0Z[2]+360→Z[2] 

J+360→J :J→Z[4] 

J+360→J : J→Z[5] 

11

Pol((A–E),(B–F)) : J<0

12

Pol((C–E),(D–F)) : J<0

13

Z[4]-Z[5]→Z[3] 

14

Z[3]<0

方位角计算顶角

计算点（

计算点（

计算点（

计算点（

）的顶角并存贮到变量

）的参数

）的参数

）的参数

Z[3]+360→Z[3]

15

tan(Z[1]) tan(G) ÷(tan(Z[1])- tan(G))→Z[7]

16

tan(Z[2]) tan(H) ÷(tan(Z[2])- tan(H))→Z[8]

17

tan(Z[3]) tan(K) ÷(tan(Z[3])- tan(K))→Z[9]

18

(Z[7]A + Z[8]C+ Z[9]E) ÷(Z[7]+ Z[8]+ Z[9])→Z[10] 

20

"XP="：Z[10]

◢

21

"YP="：Z[11]

计算所求点的X坐标

19

(Z[7]B + Z[8]D+ Z[9]F) ÷(Z[7]+ Z[8]+ Z[9])→Z[11]  计算所求点的Y坐标

显示所求点的X坐标

显示所求点的Y坐标

4

22

（）

（）

（）

6坐标转换计算程序，程序名

1

Deg: Fix 3 

2

"X0"? A: "Y0"? B 

3

"ANGLE"?E 

4

Cls: "LEFT(1), RIGHT(ELSE) "?F 

5

"XPP"? C: "YPP"? D 

6

[[A] [B]] →Mat A 

7

[[cos(E), -sin(E)] [sin(E), cos(E)]] →Mat B 

8

F≠1-D→D 

9

[[C] [D]] →Mat C 

10

Mat A+Mat B×Mat C→Mat D

◢

两平面坐标系转换计算示意图：

11

设置角度单位为十进制，3位固定小数显示

输入施工坐标系原点在统一坐标系中的X、Y坐标（A,B）

输入统一坐标系的X轴顺时针旋转至施工坐标系的X轴的角度

值

判断施工坐标系是左手系还是右手系

待换算平面点在施工坐标系中的X、Y坐标（C,D）

向矩阵A赋值

向矩阵B赋值

根据施工坐标系的类型决定变量D的取值

向矩阵C赋值

计算转换后的坐标值，并将结果显示在矩阵D中

5

7直线中边桩坐标计算程序,程序名

1 Deg：Fix 3 

2

"X0="?X: "Y0="?Y: "K0="?E: "FWJ="?J 

3

Lbl 1

4

?K:?L:?A 

5

X+(K-E)cos( J ) +Lcos(J+A) → C 

6

"XP=": C

◢

7

Y+(K-E)Sin( J ) +LSin(J+A) → D 

8

"YP=": D

◢

9

Goto 1 

10 备注：该程序请从起点往终点计算。

设置角度单位为十进制，3 位固定小数显示

输入直线起点坐标(X,Y)，起点桩号E，起点方位角J

设置程序标签

提示输入所求点桩号K，边距L，和中线夹角A

计算所求边桩点的X坐标

显示所求边桩点的X坐标

计算所求边桩点的Y坐标

显示所求边桩点的Y坐标

转向程序标签1继续执行程序

8圆曲线中边桩坐标计算程序,程序名

1

Deg：Fix 3 

设置角度单位为十进制，3 位固定小数显示

输入圆曲线起点坐标（A,B），起点桩号K，起点切线方

位角C,半径R（左负右正）

设置程序标签

提示输入所求点桩号L，

提示输入所求点边距E和以及夹角N（左负右正）

计算圆心角的一半，并把弧度转换成角度值。的输入

[FUNCTION] [5](ANGLE) [2]()

计算所求点的弦长，并取绝对值。

计算所求边桩点的X坐标

计算所求边桩点的Y坐标

显示所求边桩点的X坐标

显示 所求边桩点的Y坐标

转向程序标签1循环执行程序

r

r

2

"X="?A:"Y="?B:"FWJ="?C:"KO="?K:"R="?R 

3

Lbl 1

4

"K="?L 

5

"B="?E:"N="?N 

6

(((L-K)÷R)÷2)→M 

7

Abs(2Rsin (M) )→F 

8

A+Fcos (C+M)+Ecos (C+2M+N)→G 

9

B+Fsin (C+M)+Esin (C+2M+N)→H 

10

"XP=":G

◢

11

"YP=":H

◢

12

Goto 1

13 备注：该程序请从起点往终点计算。

r

6

9完整缓和曲线中边桩坐标计算程序,程序名

1

Deg：Fix 3 

设置角度单位为十进制，3 位固定小数显示

输入缓和曲线起点坐标(X,Y)，起点切线方位角C，起点

桩号K，缓和曲线参数A,小半径R（左负右正）

设置程序标签

提示输入所求点桩号L，

提示输入所求点边距E和以及夹角D（左负右正）

计算偏转系数并存储到S

计算所求点的曲线长并存储到P

2

"X0="?X:"Y0="?Y:"FWJ="?C:"KO="?K:"A="?A :"R="?R 

3

Lbl 1

4

"K="?L 

5

"B="?E:"D="?D 

6

Abs (R)÷R→S 

7

L-K→P 

9

P÷(6R)- P

2

∧（

4

）

∧（

3

）

2

∧（

5

）∧（

4

）∧（

7

）∧（

6

）∧（

9

）∧（

8

）

÷ (40 R)+ P÷(3456 R)- P÷(599040R)+ P(÷175472640R)→M 

8

P- P

计算切线支距X坐标并存储到M

÷(336 R

∧（

3

）

)+ P

∧（

6

）

÷(42240 R

∧（

5

）

)- P

∧（

8

）

÷(9676800R

∧（

7

）

)+ P

∧（

10

）

÷(3530096640R

∧（

9

）

)→N 

计算切线支距Y坐标并存储到N

计算弦长并存储到Q

计算弦长对应的偏角值。并存储到T

计算切线夹角并存储到U

计算所求边桩点的X坐标

计算所求边桩点的Y坐标

显示所求边桩点的X坐标

显示 所求边桩点的Y坐标

转向程序标签1循环执行程序

10

M

2

+N

2

→Q

11

cos

-1

（M÷Q)→T



12

(P

2

÷A

2

÷2)

r

→U



13

X+Qcos (C+ST)+Ecos (C+SU+D)→G 

14

Y+Qsin (C+ST)+Esin (C+SU+D)→H 

15

"XP=":G

◢

16

"YP=":H

◢

17

Goto 1

18 备注：该程序只适合完整缓和曲线。

7

公路高铁程序

 fx-5800P程序使用流程：

本程序数据和主程序是分开的，编程时将不同的工程数据存放到不同的数据文件里，如A匝道，文件名

为A，将匝道A所有的曲线线元参数输入A文件里。运行时只要运行文件名A的程序就可以了，具体运行

流程见下图：

运行数据程序文件如ZX、B、C等

输入测站点坐标X0,Y0

第四位小数点不能是0.0001， 如1100.0第四位小数点后为0.0001，如1100.0001

输入计算区域附近桩号

输入计算桩号L

输入所求P点的坐标XP,YP

计算下一点

输入中桩切向右侧夹角M

输入距离所求点到所求桩号中线距

离K，左－右+

输出所计算位置坐标X,Y

输出所计算位置与测站坐标X,Y距

离和方位角DIS,FW

输出P点所在线路的桩号L和距离

K,左－右+

输入桩号L计算下一点

8

数据文件的编写

（一）交点法数据文件编辑

交点法编写数据文件必须是对称型的，即直线段→缓和曲线段→圆曲线段→缓和曲线段→直线段，（如

果任意一端没有直线段，则把直线段长度看做是0），另外圆曲线两侧缓和曲线的旋转常数必须相等，并且

和直线段连接处的半径必须是无穷大。

交点法数据文件编写一般是根据设计图纸提供的平面曲线参数一览表提供的参数来编写，每个弯道包

括：弯道起点方位角（C），交点X坐标（D），交点Y坐标（E），缓和曲线长度（F，当没有设缓和曲线

时，F=0），交点转交（G，向左转弯，G为负值，向右转弯，G取正值），交点桩号（H），弯道圆曲线半

径（R）。 下图是一段市政道路设计参数数据。

根据上图提供的数据，可以编辑成如下的数据文件：文件名（O表示变量字母O，其余的0均为零）

3→DimZ：“X0”?A：”Y0”?B：“Ln”?L：Abs(L)Int(Abs(1000L))÷1000→Z*3+ ：Lbl 0：If Z*3+≠0.0001：Then ?L：E

lse “Xp”?X：”Yp”?Y：X→Z*1+：Y→Z*2+：IfEnd：Lbl 1:If L>0 ：Then 98°39°35.12°→C：4774.384→D：2415

.861→E：140→F：31°17°23°→G：410.007→H：600→R：IfEnd：If L>1060 Then 129°56°58.19°→C：4206.

421→D：3093.946→E：70→F：-33°50°48°→G：1285.437→H:600→R：IfEnd：Prog”XLJS” ：If Z[3+≠0.0001：

Then Goto 0：Else (Z[1]－X)cos(O)＋(Z[2]－Y)sin(O)→N：L＋N→L：Int(Abs(1000L))÷1000+0.0001→L：－(

Z[1]－X)sin(O)+(Z[2]－Y)cos（O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“：L

◢

“K=” ：K

◢

IfEnd：Goto 0: IfEnd

在程式中，有两个条件转移语句即

If L>0：Then 98°39°35.12°→C:4774.384→D: 2415.861→E:140→F:31°17°23°→G:410.007→H:600→R:IfEnd

If L>1060：Then 129°56°58.19°→C:4206.421→D: 3093.946→E:70→F:-33°50°48°→G:1285.437→H:600→R:IfEnd

„„

9

如果还有其他弯道，可以继续后加。在这些存放设计参数的语句程序表达式是固定的。编辑曲线参数时，每

个曲线参数放在一个If L>***.***（两个弯道中间直线段上的任意桩号） :Then ***°**°**.**°→C（弯道

起点方位角）：****.***→D（交点X坐标）：****.****→E（交点Y坐标）：***.***→F（缓和曲线长度：**°

**°**°→G（转角，向左转为负值，向右转为正值）：***.***→H（交点里程桩号）:*** →R（圆曲线半径） ：

IfEnd 条件式语句里。如果有多个弯道，一直按上述形式编下去。

变量说明：

序 语句

1

2

3

说明

L>***.***

***°**°**.**°→C

****.***→D、****.****→E

**°**°**°→G

***.***→H

***→R

弯道参数的起点范围，可以是弯道起点，也可以是弯

道前直线段里的任意点桩号。

弯道起点方位角

交点X坐标、Y坐标

缓和曲线长度，如果没设臵缓和曲线，则缓和曲线长度看做0

转角，向左转为负值，向右转为正值

交点里程桩号

圆曲线半径

4

***.***→F

5

6

7

【注意】运行程序时，第一步是输入测站坐标X0，Y0，然后输入近似桩号，主要是为了加快收敛，

输入时注意当根据桩号和距中距离计算所求点坐标时，请不要把小数点第四位输入0.0001，如

K21+369.0921，请把桩号输入21369.092即可，相反，当根据坐标反算桩号和距中距离时，请在输入

时把第四位小数数位0.0001，如K21+200，输入时请输入21200.0001。

10

（二）线元法数据文件编辑

线元法数据文件编辑适合于各种各样的线型组合，特别是在大型立交区或者是不规则的线型组合地方，

如小区道路、重丘区山路等有特别明显的有点，本程序为了减少使用者设计参数的输入，特别把直线与缓和

曲线以及圆曲线连接的地方只要输入缓和曲线参数就可以了，这样可以减少输入很多设计参数，同时节约大

量的内存。看下图：

该图是一个半喇叭型的互通立交，一共是5条匝道，我们以D匝道为例，来看看线元法数据文件的建立。

11

下表是D匝道的线元表

：

线元

名称

缓和曲线

圆曲线

缓和曲线

起点桩号

DK0+000.00

偏向

右偏

半径或缓和曲

线旋转常数A

55

40

65

线元长度

(m)

39.289

39.216

105.625

线元起点坐标、方位角

X(m)

494341.902

494347.926

494380.677

Y(m)

478027.825

477989. 873

477971.271

方位角

260°38′55.7″

302°18′26.0″

358°28′45.9″

DK0+039.289 右偏

DK0+078.504 右偏

D匝道线元数据分析：从上表中可以看出，D匝道由三个线元组成，第一段是缓和曲线段，根据上表判

断，该缓和曲线属于第一段缓和曲线，其最小半径处桩号为J=K0+039.289，该处坐标X=494347.926,Y=477989.

873，半径为R=40米，方位角W=302°18′26.0″，缓和曲线旋转常数P=55，在缓和曲线终点K0+039.289处，

有一圆曲线与其相连接，长度I=39.216。第二段是圆曲线，第三段又是一段缓和曲线，根据上表判断，该曲线

属于第二缓和曲线，曲线最小半径位于与园曲线连接处，桩号J=K0+78.504，坐标为X=494380.677，Y=477971.271,

方位角W=358°28′45.9″，最小半径R=40，缓和曲线旋转常数P=-65（为了判断是第一段还是第二缓和曲线，

我们通过P的正负来判断，第一缓和段P统一用正数表示，第二缓和段统一用负号表示。起点处相连的圆曲线

长度I可以看做0，也可以看做是39.216.由此我们整理下。从整个线型全部是向右转向，所以R全部取正数。

具体如下：

第一段缓和曲线线元参数说明：

最小半径处坐标

X=494347.926 Y=477989. 873

最小半径R=40

最小半径处方位角：W=302°18′26.0″

缓和曲线旋转常数P=55

与之相连的圆曲线长度I=39.216

最小半径处桩号J=39.289

此段数据计算范围可以从K0+000~ K0+78.504

第二段缓和曲线线元参数

最小半径处坐标

X=494380.677 Y=477971.271

最小半径R=40

最小半径处方位角：W=358°28′45.9″

缓和曲线旋转常数P=65

与之相连的圆曲线长度I=39.216

最小半径处桩号J= W=358°28′45.9″

此段数据计算范围可以从K0+78.504到匝道结束。

当桩号L位于K0+000~ K0+78.504时，曲线参数按第一段编写，当L>78.504时，按第二段缓和曲线编

写，当然第二段缓和曲线计算范围也可以从圆曲线起点K0+039.289到K0+78.504中间任意桩号开始。下面

是D匝道的数据程序文件（文件名“”）：（O表示变量字母O，其余的0均为零）

3→DimZ：“X0”?A: “Y0”?B: “Ln”?L: Abs(L)-Int(Abs(1000L))÷→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L:Else

“Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L≥0 :Then 302°18°26°→W：494347.926→X：477989.

873→Y：55→P: 40→R：39.216→I: 39.289→J：IfEnd: If L≥78.504：Then 494380.677→X：477971.271

→Y：358°28°45.9°→W：40→R： -65→P:39.216→I: 78.504→J：IfEnd:Prog””: Prog””:

If Z[3]≠0.0001:Then Goto 0:Else (Z[1]－X)cos(O)＋(Z[2]－Y)sin(O)→N:L＋N→L:-(Z[1]－X)sin(O)+(Z[2]－Y)c

os（O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“:L

◢

“K=”:K

◢

IfEnd:Goto 0: IfEnd

12

在上面程式中，每个缓和曲线的线元参数放在一个If „ Then „„ IfEnd循环语句中，If 后面紧跟计

算的起点范围，终点范围从下段的起点范围结束。

变量定义：

序 变量

1

X、Y

2

R

3

W

4

P

5

I

6

J

说明

输入缓和曲线半径最小处的X、Y坐标

缓和曲线最小半径，当线路左转时，半径数据前面加“-”号，以判断线路转向。

输入缓和曲线半径最小处的方位角

缓和曲线旋转常数，第一缓和曲线长度为正数，第二缓和曲线常数为负数。当不

设缓和曲线时，P=0

与缓和曲线最小半径端连接的圆曲线长度。

缓和曲线小半径段的桩号。

再看看E匝道数据文件编辑

E匝道线元设计数据表

半径或

缓和曲

线旋转

常数A

120

200

95

75

86.75

线元起点坐标、方位角

线元长

度(m)

71.296

72

60.070

45.125

64.841

17.800

线元名称 起点桩号 偏向

X(m)

494360.619

494373.205

494381.631

494372.527

494355.878

494337.373

Y(m)

477706.293

477776.470

477847.870

477907.018

477948.932

478010.658

方位角

79°49′57.0″

79°49′57.0″

90°08′44.6″

107°21′16.3″

113°49′05.6″

92°24′19.0″

直线

缓和曲线

圆曲线

缓和曲线

缓和曲线

圆曲线

EK0+000.00

EK0+071.296

EK0+143.296

EK0+203.366

EK0+248.491

EK0+313.333

右偏

右偏

右偏

左偏

左偏

13

在匝道线元表中，我们可以看出，E匝道式由一段直线段→缓和曲线段→圆曲线段→缓和曲线段→缓

和曲线段→圆曲线段连接起来的。值得注意的是，在位于K0+248.491处，是两段反向的缓和曲线连接点，

在编写数据程序文件时，分三段编写，第一段从K0+000到K0+203.366，第二段从K0+203.366到K0+248.919，

第三段从K0+248.919到线路终点K0+331.13。具体文件如下：（O表示变量字母O，其余的0均为零）

3→DimZ：“X0”?A:”Y0”?B: “Ln”?L:Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L: Else “Xp

”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L≥0 :Then 494381.6316→X：477847.870→Y：90°08°44.6°

→W：200→R：120→P:60.07→I: 143.296→J：IfEnd: If L≥203.366：Then 494372.527→X：477907.018→Y：

107°21°16.3°→W：200→R：-95→P:60.07→I: 203.366→J：IfEnd: If L≥248.491：Then 494337.373→X：

478010.658→Y：92°24°19°→W：－86.75→R：75→P:17.8→I: 313.333→J：IfEnd: Prog””: Prog”

：If Z[3]≠0.0001:Then Goto 0:Else (Z[1]－X)cos(Ｏ)＋(Z[2]－Y)sin(Ｏ) →N:L＋N→L:-(Z[1]－X)sin(Ｏ”

)+(Z[2]－Y)cos（Ｏ）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“:L

◢

”K=”:K

◢

IfEnd:Goto 0: IfEnd

14

为了让使用者熟悉掌握数据程序文件的编写线元法的数据程序文件，现在将C、F、G等匝道的线元数据

列出来，供使用者练习。C匝道数据；F匝道数据；G匝道数据

缓和曲线 FK0+111.505

FK0+202.434

FK0+303.548

GK0+000.000

GK0+133.152

GK0+282.286

右偏

右偏

右偏

右偏

右偏

右偏

100

110

120

175

230

190

90.909

101.114

130.909

133.152

149.124

149.72

494557.835

494542.812

494595.310

494718.802

494587.600

494472.335

478313.650

478224.689

478142.420

478112.963

478093.368

478002.890

252°32′00.2″

276°12′33.6″

32852′35.7″

182°58′11.8″

199°33′17.5″

236°42′21.1″

直线 FK0+000.00 111.525 494591.311 478420.033 252°32′00.2″

圆曲线

直线

缓和曲线

圆曲线

缓和曲线

直线

缓和曲线

圆曲线

CK0+000.00

CK0+056.876

CK0+078.153

CK0+108.332

CK0+138.763

CK0+168.941

CK0+590.964

CK0+675.964

右偏

右偏

右偏

右偏

右偏

右偏

右偏

32

65

140

65

85

85

85

56.876

21.277

30.179

30.431

30.179

422.023

85

72.739

495005.991

494991. 262

494971. 511

494943.127

494913.167

494883.008

494461.552

494379.485

478056.677

478104. 123

478112.036

478122.242

478127.219

478126.740

478104.874

478086.681

56°19′42.3″

158°09′54.4″

158°09′54.4″

164°20′25.7″

176°47′40.6″

182.°58′11.8″

182.°58′11.8″

211°37′04.2″

半径或

线元名

称

起点桩号 偏向

缓和曲

线旋转

常数A

线元长

度(m)

X(m)

线元起点坐标、方位角

Y(m) 方位角

圆曲线

缓和曲线

缓和曲线

圆曲线

缓和曲线

15

道路计算程序

计算程序是核心程序，但是道路计算程序是固定的，使用者不需要做任何修改，也不直接运行道路计算

程序分为两部分，第一部分是常规的对成型线型组合计算的程序（文件名：

件直接调用；第二部分是适用于非对称的线型组合计算文件（文件名：

文件后，还要调用“

Deg:Fix 4

G÷Abs(G)→V: F÷2-F^(3)÷（240R

2

）+F^（5）÷(34560R^（4）)-F^（7）÷(8386560R^（6）)+F^（9）÷(3158507

2

520R^（8）)→Q:F÷(24R)-F^(4)÷(2688R^(3))+F^(6)÷(506880R^(5))-F^(8)÷(154828800R^(7))→P:Q+(R+P)Vtan(G÷2)

→T:πRGV÷180+F→I:

IF L≤H-T:Then L-H→J:Goto 1:IfEnd:

If L

），交点法编辑的数据文

），在数据文件调用完“”

”文件。程序代码：

（为了输入方便代码分行显示，使用者不必分行输入。）

If L

If L

If L≥H-T+I:Then L-H-I+2T→J:C+G→C:Goto 1:IfEnd:

Lbl 1:

D+Jcos (C)→X:

E+Jsin(C)→Y:

C→O: Goto C:Lbl 2:

(J-J^(5)÷(40(RF)

2

)+J^(9)÷(3456(RF)^(4))-J^(13)÷(599040(RF)^(6))+J^(17)÷(175472640(RF)^(8))-T)S→M:

(J^(3)÷(6RF)-J^(7)÷(336(RF)^(3))+J^(11)÷(42240(RF)^(5))-J^(15)÷(9676800(RF)^(7))+J^(19)÷(3530097000(RF)^(9)))

V→N:

C+90J

2

SV÷(πRF)→O:

Goto B

Lbl 3:

180(J-0.5F) ÷(πR)→O:

Q+Rsin (O)-T→M:(P+R(1-cos (O)))V→N:C+VO→O：

Goto B

Lbl B:

D+Mcos(C)-Nsin(C)→X:

E+Msin(C)+Ncos(C)→Y:

Goto C:

Lbl C:If Z*3+≠0.0001:Then ?K: "M"?U:K<0=>180-U →U:O+U×Abs(K)÷(K+10^(-10)) →Z:X+Abs(K)×cos(Z)→X:

:Y+Abs(K)×sin(Z)→Y:

"X":X

◢

"Y":Y

◢

tan-1((Y-B)÷(X-A+10^(-10))→W: If X-A<0:Then W+180→W: IfEnd: If W<0:Then W+360→W:IfEnd:√((X-A)

2

+(Y-B

2

)

)→S:"DIS": S

◢

"FW": W

DMS

◢

IfEnd

16

非对称线型组合计算程序（）（O表示变量字母O，其余的0均为零）

Abs(P+10^(-10)）÷(P+10^(-10)) →S：Abs (R)÷R→V:Abs(R) →R:P

2

÷R→F:

180(I+F)V÷(πR) →G:90FSV÷(πR) →O:

If S=1:Then W-O→C: IfEnd: If S=-1:Then W-O-G→C :IfEnd:C<0C+360→C:C≥360C-360→C:

F

2

÷(24R)-F^(4)÷(2688R^(3))+F^(6)÷(506880R^(5))-F^(8)÷(154828800R^(7))→P:F÷2-F^(3)÷(240R

2

)+F^5÷(

34560R^(4))-F^(7)÷(8386560R^6)+F^9÷(3158507520R^(8))→Q:

Q+(R+P)Vtan(G÷2)→T:If S=1:Then J-F+T→H:C→Z: Else J-I-F+T→H:C+G→Z:IfEnd:

(F-F^(3)÷(40R

2

)+F^(5)÷(3456R^(4))-F^(7)÷(599040R^(6))+F^(9)÷(175472640R^(8))-T)S→M:(F2÷(6R)-F^(4)

÷(336R^(3))+F^(6)÷(42240R^(5))-F^(8)÷(9676800R^(7))+F^(10)÷(3530097000R^(9)))V→N:

X-Mcos(Z)+Nsin(Z)→D:Y-Msin(Z)-Ncos(Z)→E

最小半径处坐标X=494347.926 Y=477989. 873 最小半径R=40 最小半径处方位角：W=302°18′

26.0″ 缓和曲线旋转常数P=55 与之相连的圆曲线长度I=39.216 最小半径处桩号J=39.289

17

程序运行

程序运行分为两种，第一种是根据根据使用者输入的桩号和距中距离计算出所求点的坐标和测站点到该

点的方位角和水平距离，可以直接利用计算器提供的距离和方位角进行极坐标放样。另一种是根据使用者输

入的实测坐标计算出该点的桩号和距中距离，可以根据桩号和距中距离进行高边坡、隧道断面点的设计高程、

以及进行路面高程控制计算方面很方便。

（一）据桩号和距中距离计算所求点的坐标、及方位角和水平距离。要进行程序的线路计算，直接运行数

程序文件就可以了，不需要理会计算程序。如果计算某桩号L距离中线K处的坐标X，Y和该点到测站点的

距离和方位角。运行提示如下：

序 变量

1

X0?

2

Y0？

3

Ln?

4

L？

5

说明

提示输入测站点X坐标

提示输入测站点Y坐标 （测站点的位数和数据程序文件里坐标取位要一致）

计算点近似桩号，近似桩号输入第四位小数后不能等于0.0001，

输入计算点准确桩号，近似桩号输入第四位小数后不能等于0.0001。

提示输入距离中桩的桩号如果计算点位于线路的左幅，则在数据前加入“-”号，如

果是计算中桩坐标，则输入0

K？

6

7

M？

X

FW

输入该点与线路的前进方向的按顺时针转过的角度，常于斜交的涵洞，桥梁等构造

物

显示所求点的X坐标

显示所求点的Y坐标

显示所求点至测站点的方位角

显示所求点至测站点的距离

8

Y

9

10

DIS

要继续计算则从L?开始，计算新点的坐标，具体流程见程序流程图的左边部分。例：要计算D匝道的

DK0+020中桩坐标、DK0+080左边5.3米的坐标、Dk0+240斜交75°58′11.8″、距中右侧5米处的坐标，

并计算出该点至坐标为X=495005.991，Y=478056.677控制点的水平距离和方位角。

18

运行文件“

序 变量

1 X0?

2 Y0？

3 Ln?

4 L？

5

6

7

K？

M？

X

说明

495005.991 测站点X坐标

478056.677 测站点Y坐标

”提示如下：

0 输入小输点第四位起不等于0.0001的桩号

20 输入所求点桩号

0 输入距中距离

90 输入斜交角度

494341.482 显示所求点X坐标

478007.903 显示所求点Y坐标

666.2959 显示所求点至测站平距

184°11′52.42″显示所求点至测站方位角

80 输入下点的桩号

-5.3

90

494382.2290

477965.9594

630.324

188°16′29.5″

240

5.0

75°58°11.8°

494455.3110

478101.5491

552.505

175°20′29.6″

8 Y

9 DIS

10 FW

11 L？

12 K？

13 M？

14 X=

15 Y=

16 DIS

17 FW

18 L？

19 K？

20 M？

21 X=

22 Y=

23 DIS

24 FW

19

（二）据坐标反算桩号和距中，根据某点的坐标，计算出该点在线路中的桩号和距离，程序运行时提示如下：

序 变量

1 X0?

2 Y0？

3 Ln?

4 Xp？

5

6

7

Yp？

L

K

说明

提示输入测站点X坐标

提示输入测站点Y坐标（测站点的位数和数据程序文件里坐标取位要一致）

计算点近似桩号，近似桩号输入第四位小数后不能等于0.0001，

所求点P的X坐标

所求点P的Y坐标

显示所求点的桩号

显示所求点的距中距离，如果位于线路的左边，则显示结果为负值，如果K为0时，则该点在线路

中线上。

计算新点的桩号和距中距离，则从?Xp 开始。如此循环计算。具体流程见流程图的右边部分

例：求点X=494382.2290，Y=477965.9594在D匝道中的位置，运行文件“

序 变量

1 X0?

2 Y0？

3 Ln?

4 Xp？

5

6

7

Yp？

L

K

说明

495005.991 测站X坐标(可输入0)

478056.677 测站Y坐标(可输入0)

60.0001 输入小输点第四位起等于0.0001的桩号

494382.2290 输入所求点X坐标

477965.9594 输入所求点Y坐标

79.999999 显示所求点桩号

-5.3000 显示所求点距中距离, “-“ 表示该点位于中桩左侧

”如下：

20

程序应用

（一）本程序在高边坡放样中的应用

高边坡放样是目前公路道路测量工作普遍要做的一项工作，边坡坡顶线或者填方边坡坡脚线是边坡放样

最繁琐的工作，利用本程序进行边坡放样，工作步骤如下：

① 司镜员走到目的地后，先测出一个点的坐标和高程；

② 利用本程序求出该点在线路中的桩号和距中距离，然后根据计算出来的桩号和设计图纸，内插出道

路边桩高程，根据道路边桩高程和道路设计宽度、坡度、设计碎落平台的宽度和间隔高度、实测点

的高程计算出按实测高程，边坡距中距离。

③ 根据高程算出的距中距离和程序计算出立镜点至距中距离差值dK，指挥司镜员靠近道路中线或者

往外偏移多少米。

④ 重新做1~3工作，直到dK 小于边坡宽度控制精度要求为止。

dK

m

h

m

h

K1 K2

如上图，司镜员在P点立镜，测出P点坐标和高程后，用本程序计算出P点到的中桩桩号和距离

中桩距离K，根据中桩桩号和有关设计数据内插出边坡坡脚点的设计高程Hs和距离中桩的距离K1，

根据Hc、和Hs的高差计算出K2，按照P点的高程，边坡顶距中桩的理论宽度K

理论

，具体计算如下：

K

理论

=K1+h*n1+m+h*n2+m+（Hc-Hs-2h）*n3

则dK=K- K

理论

dK 求出来后，可以指挥司镜员靠近道路dK米，然后重复以上步骤，直到dK值能满足边坡控制精度。

21

（二）本程序在隧道断面欠挖、超挖放样中的应用。利用本程序进行隧道断面开挖的欠挖和超挖计算，看下图：

R2'

R2

用全站仪测出P点的坐标和标高后，根据坐标算出P点的中桩桩号和距中距离K，根据中桩桩号查看P

点的路面高程，算出P点至设计路面的高差h，根据P点所在圆弧的圆心O2在断面中的位置（K2，h2）算

出P点到O2点距离R2’，用R2’-R2的差距就是P点的超挖或者欠挖的数据。

（三）本程序在路面施工控制放样中的应用

路面施工时，路面的轴线偏位、路面宽度、路面高程是测量的重要控制指标，也是道路施工控制最严格

的一道测量控制。利用本程序输入桩号和距中距离，可以直接计算出所求点至测站点的水平距离和方位角，

从而直接通过极坐标法直接放出该点在实地的位置。

（四）处理线路设计中的“断链”问题

断链是线路设计中经常出现的问题，断链分为长链与短链，如果是短链，那么只要弯道的计算范围从短

链处分界就可以了，如果是长链，则需要把线路分成两个文件，从开始出现断链处分，这样可以避免桩号相

同，但位置不同的情况。

（五）本程序在道路竣工验收过程中的应用

道路竣工验收，其主要方面就是轴线偏位和宽度，可以通过现场测量道路中线的坐标，通过坐标反算出

该点和设计中线的偏移值K（轴线偏位值）。

22

公路纵断面高程计算

公路纵断面计算也是采用计算程序和数据程序分开的原则，数据文件是专门存放道路纵断面设计参数，

对于多个工程，可以建立多个数据文件，这样每个数据文件通过调用计算程序就可以计算出该点的高程。

① 数据文件：每个纵断面由纵坡和竖曲线组成，本程序每一个变坡点计算范围是从竖曲线前的纵坡任意桩

号开始到竖曲线后纵坡上的任一桩号结束。下表是某线路路基纵断面设计的参数一览表：

路基纵断面设计

桩号

K037+950.00

K038+650.00

K038+900.00

K039+245.00

K039+410.00

K039+590.00

变坡点高程

80.56

75.93

86.01

87.75

91.86

93.21

99.84

半径

3000.00

4500.00

6000.00

6000.00

10000.00

切线长

70.40

79.372

59.597

52.227

52.214

竖曲线起点桩号 竖曲线终点桩号

K038+579.599

K038+820.628

K039+185.403

K039+357.773

K039+532.786

K038+720.401

K038+979.372

K039+304.597

K039+462.227

K039+647.214

纵坡

-0.6614%

4.0320%

0.5043%

2.4909%

0.7500%

1.8943%

根据上表编程数据文件如下：文件名

"Dh"?E:Lbl 0:"L="?L:"I="?I:"K="?K 

If L≥37950:Then 75.93→C:38650→D:-0.006614→M:0.04032→N:3000→R:IfEnd 

If L≥38820.628:Then 86.01→C:38900→D:0.04032→M:0.005043→N:4500→R:IfEnd 

If L≥39185.403:Then 87.75→C:39245→D:0.005043→M:0.024909→N:6000→R:IfEnd 

If L≥39357.773:Then 91.86→C:39410→D:0.024909→M:0.0075→N:6000→R:IfEnd 

If L≥39532.786:Then 93.21→C:39590→D:0.0075→M:0.018943→N:10000→R:IfEnd 

Prog "":Goto 0

从上面程序看，"Dh"?E:Lbl 0:"L="?L:"I="?I:"K"?K 语句是固定的，“Dh”提醒使用者输入所求点距离路面

的厚度，比如路面底基层距离路面厚度为37cm，则输入0.37，表示所求点是路面底基层顶面高程，“L=”，

输入所求点桩号,“I”=表示输入路面横坡值，“K=”表示输入路面某点至中桩距离，均为正值。以上数据均

在程序运行时输入。在每一个变坡点处，要求输入的参数有：

序 变量

1 C；

2 D：

3 M：

4 N：

5 R：

说明

输入变坡点高程

输入边坡点桩号

输入变坡点小里程段纵坡坡度

输入变坡点大里程侧纵坡坡度

竖曲线半径

L表示竖曲线前纵坡任意里程的桩号，本程序每段竖曲线参数计算范围是从竖曲线前的纵坡上任意点到竖曲

线后纵坡上任意点的桩号。

23

② 高程计算程序（文件名：

③ 运行程序显示如下

说明

），高程计算文件代码如下：

N-M→P: Abs(P)÷P→S:(Abs(P)×R)÷2→T:If L＜D-T:Then C+(L-D)M→H:IfEnd 

If L＞D-T:Then C-MT+(M+(L-D+T)S÷(2R))×(L-D+T)→H:IfEnd 

If L≥D+T:Then C+(L-D)N→H:IfEnd 

H-E-KI→H:"H=":H

◢

序 变量

1 Dh？

2 L=?

3 i？

4 K？

5 H=

提示输入至路面厚度

提示输入桩号

提示输入横坡度

提示输入距中距离，均为正数

提示输出所求点高程

例：根据上表计算K38+660中桩高程、K38+900右侧4.5米高程、K39+600底基层面左侧4.5米高程（假设

路面沥青厚度为6cm，封釉层1cm，水泥稳定层基层厚度为30cm），运行

序 变量

1 Dh？

2 L=？

3 i？

4 K？

5

6

7

8

9

H=

L=？

i？

K？

H=

说明

0

38660

0

0

76.9412

38900

0.02

4.50

85.2200

文件显示如下：

退出程序，重新运行。

序 变量

1 Dh？

2 L=？

3 i？

4 K？

5 H=

说明

0.370

39600

0.02

4.50

93.0509

24

编程基础

 fx-5800P的编程

采用卡西欧三步编程法写出您自己需要的程序。步骤如下：

① 按【MODE】【5】(PROG)【1】(NEW)(输入文件名如P1)后【EXE】【1】(COMP)

② 按以下步骤输入程序源代码。

序 程序源代码

1

?A

:

?B

:

?C

:

?D: ?E 

2

A+Dcos(C) + Ecos(C-90) → X 

3

B+D sin(C) + Esin (C-90) → Y 

4

"X="：X

◢

5

"Y="：Y

◢

编程思路说明

第①步：

提示用户输入已知数据,

语法『?+变量+:（断开）』

 是换行符号。

第②步：

数据计算处理，

语法『直接输入公式，→+变量（存储结果数据）』

第③步：成果数据显示在屏幕上。

语法『"变量="+：（断开）+变量+

◢

』

③ 按【FILE】键再按

调用程序P1。

在案例中，输入公式的部分可以是用户的任意公式，用户只要直接输入公式（公式是什么样就

输入什么样，没有任何语法）。表明卡西欧计算器编程操作简单，功能强大。

 fx-5800P的公式自定义

没有语法,变量可带脚标，α β θ等符号直接作为变量使用；

请直接写出您自己需要的公式并调用。步骤如下：

① 按【MODE】【5】(PROG)【1】(NEW) (输入文件名如F1)后【EXE】【3】(Formula)

α

◢

α

② 公式编写好后，按【EXIT】退出到程序菜单。再按[MODE] [1]回到普通计算界面，按[FMLA]

[1] (Original)移动光标选择程序名如F1，调用程序。

【注意】此操作过程常见问题：

1：脚标的输入是[FUNCTION] [4](ALPHA) [5]

[7]。

2：在MATH状态下按[FMLA] [1] (Original)直接调用自定义公式。

3：以上公式是已知直线起点坐标,里程求算直线上中桩坐标。

4：从以上公式可以看出，公式自定义可以在一个文件中连续输入多个公式，只要用

◢

断开即可。

25

实用用户公式集

 实用公式集

① 新建公式按【MODE】【5】(PROG)【1】(NEW) (输入文件名如F1)后【EXE】【3】(Formula)

② 公式运行在普通计算界面，按[FMLA] [1] (Original)移动光标选择程序名如F1运行。

1

测角前方交会坐标计算公式

β

β

2

α

α

α

α

β

β

◢

测角后方交会坐标计算公式

◢

其中

、、

需要另外建个公式文件

α

◢

θ

θ

α

◢

αθ

◢

β

θβθ

γγ

α

θ

3 直线中边桩计算坐标公式（也是坐标反算公式）

4 坐标正算公式