CASIO fx5800计算器 工程测量与公路测量程序.docx

CASIO fx5800计算器 工程测量与公路测量程序.docx

《CASIO fx5800计算器 工程测量与公路测量程序.docx》由会员分享，可在线阅读，更多相关《CASIO fx5800计算器 工程测量与公路测量程序.docx（21页珍藏版）》请在冰豆网上搜索。

CASIOfx5800计算器工程测量与公路测量程序

CASIOfx—5800计算器工程测量与公路测量程序

1、坐标正算

〝X0=〞?

X:

〝Y0=〞?

Y:

〝I=〞?

I:

〝J=〞?

J

X+ICos（J）→U:

Y+ISin（J）→V

〝X=〞:

U（待求点的X坐标）

〝Y=〞:

V（待求点的Y坐标）

说明：

X0Y0：

已知点坐标I：

两点的距离J：

方位角

２、坐标反算

Lbl0

〝X1=〞?

X:

〝Y1=〞?

Y:

〝X2=〞?

U:

〝Y2=〞?

V

Pol（U-X,V-Y）:

J<0J+360→J

〝I=〞:

I

〝J=〞:

JDMS

GOTO0

说明：

X1Y1：

第一点的坐标，X2Y2第二点的坐标，I：

两点的距离，J：

方位角

３、圆曲线

〝X0=〞?

X:

〝Y0=〞?

Y:

〝FWJ=〞?

D:

〝QDHAO=〞?

G:

〝ZDHAO=〞?

C:

〝R=〞?

R

Lbl1

〝DQHAO=〞?

L

L

〝PY=〞?

K

180（L-G）÷（πR）→E

2RSin（0.5E）→F

ReC（Abs（F）,D+0.5E）：

Cls

〝X=〞:

X+I+KCos（D+E+90）→A

〝Y=〞:

Y+J+KSin（D+E+90）→B

D+E→o:

o<0360+o→0o

〝FWJ=〞:

oDMS

Goto1

Lbl2

〝END〞

说明：

X0Y0：

起始点的坐标，FWJ：

起始点的方位角，QDHAO：

起点里程，ZDHAO终点里程，R：

半径，DQHAO：

待求点里程，PY：

偏移量

４、竖曲线

〝QZHAO=〞?

J:

〝H+B〞=?

B:

〝I1=〞?

C:

〝I2=〞?

D:

〝R=〞?

R:

〝T=〞?

T:

（（D-C）÷100）÷Abs（（D-C）÷100）→A:

Abs（0.5R（D-C）÷100）→S

〝T〞:

S

J-T→X

J＋T→Y

Lbl0

〝DQ=〞?

L

IfL

ThenGoto1:

ElseIf  L﹥Y:

ThenGoto1:

IfEnd:

IfEnd

〝GC=〞:

B+C（L-J）÷100+A（L-J+T）  ÷（2R）→H

Goto0

Lbl1

〝END〞

说明：

QZHAO：

曲中点里程，H+B：

曲中点高程，I1：

第一坡度，I2：

第二坡度，R：

半径，T：

切线长，DQ：

待求点里程。

5．正算偏移

〝X0=〞?

A:

〝Y0=〞?

B:

〝FWJ=〞?

D:

Lbl0

〝LS=〞?

L:

〝PY=〞?

C

ReC（L,D）

〝X=〞:

A+I+CCos（D+90）→X

〝Y=〞:

B+J+CSin（D+90）→Y

Goto0

说明：

X0Y0：

已知点坐标，LS：

两点距离，FWJ：

方位角

6.万能公式

〝X0=〞?

A:

〝Y0=〞?

B:

〝FWJ=〞?

C:

〝1÷R-QD〞?

D:

〝1÷R-ZD〞?

E:

〝QDHAO〞?

F:

〝ZDHAO〞?

G

Lbl1

〝DQHAO=〞?

H:

〝PY=〞?

o

IfH

ThenGoto2:

ElseIf  H﹥G:

ThenGoto2:

IfEnd:

IfEnd

（E-D）÷Abs（G-F）→P:

Abs（H-F）→Q:

PQ→I:

C+90Q（I+2D）÷π→J:

C+45Q（I÷8+D）÷π→M:

C+135Q（3I÷8+D）÷π→N:

C+45Q（I÷2+2D）÷π→K

A+Q（Cos（C）+4（Cos（M）+Cos（N））+2Cos（K）+Cos（J））÷12→U:

B+Q（Sin（C）+4（Sin（M）+Sin（N））+2Sin（K）+Sin（J））÷12→V

〝X=〞:

U+oCos（J+90）→X

〝Y=〞:

V+oSin（J+90）→Y

〝FWJ=〞:

JDMS

Goto1

Lbl2

〝END〞

说明：

X0Y0：

起开点坐标，FWJ：

起开点方位角，1÷R-QD：

起开点曲率，1÷R-ZD终点曲率，QDHAO：

起点里程，ZDHAO终点里程，DQHAO：

待求点里程，XY：

待求点坐标

6、标高计算

〝QZHAO=〞?

J:

〝QDHAO=〞?

G:

〝ZDHAO=〞?

K:

〝B—HB=〞?

B:

〝SX=〞?

Q:

〝I1=〞?

C:

〝I2=〞?

D:

〝R=〞?

R:

〝T=〞?

T

（D-C）÷100→W:

W÷Abs（W）→A:

〝T〞:

0.5RAbs（W）→S

〝o-YIGAO=〞?

o

J-T→X:

J+T→Y

Lbl1

〝L—DQHAO=〞?

L:

〝PY=〞?

U:

〝XIELV=〞?

V

IfL

ThenGoto2:

ElseIf  L﹥Y:

ThenGoto3:

IfEnd:

IEnd

〝HL=〞:

B+0.01C（L-J）+A（L-J+T）  ÷（2R）+0.01UV+Q→Q

〝HO=〞:

o-H→Z

Goto1

Lbl2

L

〝HL=〞:

B+0.01（L-J）C+0.01UV+Q→H

〝H0=〞:

o-H→H

Goto1

Lab3

L>KGoto4

〝HL=〞:

B+0.01（L-J）D+0.01UV+Q→H

〝HO=:

o-H→Z

Goto1

Lbl4:

〝END〞

说明：

QZHAO：

曲中点里程（竖曲线），QDHAO：

起点里程，ZDHAO：

终点里程，B—HB曲中点高程，SX：

上下（路面标高下多少为水稳层面），I1：

第一坡度，I2：

第二坡度，R：

半径，T切线长，L—DQHAO：

待求点里程，PY：

偏移（中桩到边桩的距离），XIELV：

斜率（横坡的坡度，左负右正）o-YIGAO：

仪器高

卡西欧Fx5800P公路逐桩坐标计算程序

本程序先把整条线路要素先一次输入到计算器额外变量中，就能计算任意中、边桩坐标，适用高速公路等一般线型，能计算两边不等长缓和曲线。

ZXZB（文件名）

Lb1  1?

6→Z?

  （此处6只表示举例的数字，实际应输入要计算线路的交点个数）

2→W?

？

K?

  （程序运营时出现“K？

”表示要输入计算的桩号，如计算K5+334.23，则输入5334.23）

WhileK>Z[W+5Z+4]:

W+1→W：

WhileEnd?

K-Z[W+5Z+3]→S?

Z[W+2Z+2]→R：

Z[W+3Z+2]→L：

Z[W+4Z+2]→E?

Pol（Z[W]-Z[W-1],Z[W+Z+2]-Z[W+Z+1]）?

IfJ<0:

Then  J+360→J:

IfEnd?

J→A  ?

IfW=2:

Then  Goto  2:

IfEnd?

Pol（Z[W-1]-Z[W-2],Z[W+Z+1]-Z[W+Z]）?

IfJ<0:

Then  J+360→J:

IfEnd?

A-J→C：

J→A  ?

L÷2-L^3÷（240R^2）→M：

E÷2-E^3÷（240R^2）→N?

L^2÷（6R）-L^4÷（336R^3）-R（1-cos（90L÷（πR）））→P?

E^2÷（6R）-E^4÷（336R^3）-R（1-cos（90E÷（πR）））→Q?

（R+Q-（R+P）cos（Abs（C）））÷（sin（Abs（C）））+M→M?

（R+P-（R+Q）cos（Abs（C）））÷（sin（Abs（C）））+N→Q?

πRAbs（C）÷180+（L+E）÷2→N?

Z[W-1]-Mcos（A）→X?

Z[W+Z+1]-Msin（A）→Y?

Z[W-1]+Qcos（A+C）→M?

Z[W+Z+1]+Qsin（A+C）→V?

Abs（C）÷C→Q?

IfS≤L:

Then0→P：

Goto  3:

IfEnd?

IfS≤N-E:

Then  S-L→S：

Goto  4:

IfEnd?

IfS≤N:

ThenN-S→S?

-Q→Q：

A+C-180→A  ：

M→X：

V→Y：

E→L：

180→P：

Goto  3:

IfEnd  ?

A+C→P：

S-N→S：

M+Scos（P）→D：

V+Ssin（P）→F：

Goto  6  ?

Lb12  ?

A+C→P  ?

Z[W-1]+Scos（P）→D  ?

Z[W+Z+1]+Ssin（P）→F：

Goto  6  ?

Lb13  ?

S-S^5÷（40R^2）÷L^2+S^9÷（3456R^4）÷L^4→I  ?

Q（S^3÷（6RL）-S^7÷（336R^3）÷L^3）→J?

P+A+90QS^2÷（πRL）→P：

Goto  5  ?

Lb14  ?

90（2S+L）÷（πR）→M?

Rsin（M）+L÷2-L^3÷（240R^2）→I?

Q（L^2÷（24R）+R（1-cos（M））→J  ?

A+QM→P?

Lb15  ?

X+Icos（A）-Jsin（A）→D：

Y+Jcos（A）+Isin（A）→F  ?

Lb16  ?

Fix3?

″X=″D◢?

                  （显示中桩横坐标）

″Y=″F◢?

                  （显示中桩纵坐标）

″AT=″P◢?

                （显示中桩切线方位角）

″B=″?

T  ?

  （程序运营时出现“B=？

”表示要输入计算边桩的距离,左-，右+）

″α=″?

U?

  （程序运营时出现“α=？

”表示输入边桩斜交角，前进方向的右角）

P+U→P?

″XB=″:

D+Tcos（P）◢  ?

          （显示边桩横坐标）

″YB=″:

F+Tsin（P）◢  ?

          （显示边桩纵坐标）

Norm  2

Goto  1

使用说明

1、上述程序代码中“（）”里为说明，不输入5800中。

2、代码中下划线的字母表示要从计算器中调出，整体输入，见5800说明书。

3、程序使用前应设置好额外变量，额外变量输入前，应按5800说明书，先将额外变量的数目设置至够用的数目，一般大于线路交点数的6倍+6个，如5个交点，则为大于36个。

额外变量的形式为：

Z[N]，N—表示具体的数字

如要把1000存入到第1个额外变量中，具体操作为：

输入1000→Z[1]后，按，就把1000存入到Z[1]中，

如要把999999存入到第45个额外变量中，就是：

输入999999→Z[45]后，按，就把999999存入到Z[45]中了。

FX5800公路逐桩坐标计算程序操作说明

计算器显示

操作步骤

1按键，打开文件菜单，选择“ZXZB”文件名，按键确认。

“K？

”提示输入要计算的桩号

“950.000”为缺省的桩号

2出现此窗口

“1020.888”为输入要计算的桩号

3输入计算的桩号

如：

K1+020.888

则输入1020.888

显示X坐标

4按键确认后计算中桩X坐标，出现此窗口。

显示Y坐标

5按键确认后计算中桩Y坐标，出现此窗口。

显示中桩切线方位角

6按键确认后计算中桩切线方位角，出现此窗口。

提示输入边桩距离

缺省的边桩距离，如不变直接按键

7按键确认后，出现此窗口，提示输入计算边桩的距离，左-，右+

输入的边桩距离

8如要计算右边桩距离为15米，则输入15。

提示输入斜交角

缺省的边桩斜交角，如不变直接按键

9按键确认后，出现此窗口，提示输入计算边桩的斜交角，为前进方向的右角，正交为90。

输入的边桩斜交角

10如要计算边桩斜交角为90度，则输入90。

显示边桩X坐标

11按键确认后计算边桩X坐标，出现此窗口。

显示边桩Y坐标

12按键确认后计算边桩Y坐标，出现此窗口。

上个计算输入的桩号

提示输入下个桩号

13按键确认后提示计算下个桩号坐标，出现此窗口。

输入的下个桩号

14按键确认后计算下个桩号坐标

卡西欧5800计算器道路高程计算程序

Z[16]→D

For1→FToZ[19]-2

4F+D→C

Z[C+3]→R            ;输入半径

Z[C]-Z[C+4]→M          ;M为正时是凸型竖曲线,负时为凹竖曲线

Abs（MR）÷2→T          ;切线长

IfG

Then            ;点在切点前

G-Z[C+1]→Q

QZ[C]+Z[C+2]→Z            ;计算高程Z

GotoE          ;跳转标签E

ElseIfG

Then                ;点在本交点的切点1和变坡点之间

G-Z[C+1]+T→Q

Q2÷2÷R→H

M>0=>-H→H

Z[C+2]-Z[C]T+Z[C]Q+H→Z          ;计算高程Z

GotoE          ;跳转标签E

ElseIfG

Then                  ;点在本交点的变坡和切点2之间

G-Z[C+1]-T→Q

Q2÷2÷R→H

M>0=>-H→H

Z[C+2]-Z[C+4]T+Z[C+4]Q+H→Z        ;计算高程Z

GotoE

IfEnd:

IfEnd:

IfEnd

Next          ;如果计算点在切点2后,就转入下一个交点交点,再搜索,程序将返回到（For1→FToZ[19]-2）                        ;如果点在最后一个交点的切线2执行以下计算,

G-Z[C+1]→Q

Z[C+2]+Z[C+4]Q+H→Z      ;计算高程Z

LblE

FX--4800P常用测量程序集

1、由一个已知点坐标和一个已知方位角计算支导线坐标（程序名：

ZBJS）

程序清单：

X：

Y：

T：

Lb10：

Fixm：

{AS}：

T=T+A-180：

T〈0=〉T=T+360⊿T〉360=〉T=T-360⊿T◢

X=X+ScosT◢

Y=Y+SsinT◢

Goto  0

程序运行说明：

运行“ZBJS”，每运行一步按EXE键一次，显示X？

时输入已知点X坐标，显示Y？

时输入已知点Y坐标，显示T？

时输入已知点由后视到前视的已知坐  标方位角，运行一次本程序只需输入一次已知数据。

显示A？

时输入第一个左转角值，显示S？

时输入测站到前视的水平距离。

接下来计算器显示的T是测站到前视的坐标方位角，按    +      可将显示的十进制角度换算为度分秒。

如此继续可计算出所有支导线未知坐标。

2、由两个已知点坐标计算支导线坐标（程序名：

ZBJS2）

程序清单：

M“X0”：

N“Y0”：

D“XH”：

E“YH”?

br/>Pol（M-D，N-E）：

Fixm：

T=J：

X=M：

Y=N：

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

I“S=”◢

Lb10:

{AS}:

T=T+A-180：

T〈0=〉T=T+360⊿T〉360=〉T=T-360⊿T◢

X=X+ScosT◢

Y=Y+SsinT◢

Goto  0

程序运行说明：

运行后，显示X0，Y0时输入已知置镜点坐标，显示XH，YH时输入已知后视点坐标，其它同ZBJS的输入与显示。

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

      程序清单：

  C“X0”：

D“Y0”：

E“T0”：

Fixm：

Lb11:

{SA}?

br/>T=E+A-180：

T〈0=〉T=T+360⊿T〉360=〉T=T-360⊿T◢

X=C+ScosT◢

Y=D+SsinT◢

Goto  1

程序运行说明：

变量输入参照ZBJS程序。

4、由两个已知点坐标进行放射型导线坐标计算（程序名FSXZBJS2）

M“X0”：

N“Y0”：

D“XH”：

E“YH”?

br/>Pol（M-D，N-E）：

Fixm：

F=J：

X=M：

Y=N：

F〈0=〉F=F+360⊿F“TH-J=”◢

I“S=”◢

Lb11:

{SA}?

br/>T=F+A-180：

T〈0=〉T=T+360⊿T〉360=〉T=T-360⊿T◢

X=X+ScosT◢

Y=Y+SsinT◢

Goto1

5、坐标计算优化程序（可进行连续式（LX）、放射式（FS）已知一点（1-D）或者两点（2-D）进行坐标计算）：

V“1=>LX，≠>FS”：

W“1=>1-D≠>2-D”：

W=1=>Goto1:

≠>Goto2?

br/>Lb11：

F“X0=”：

G“Y0=”：

H“TH-0=”：

Goto3?

br/>Lb12:

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=”◢

Goto3?

br/>Lb13:

X=F：

Y=G：

T=H：

V=1=>Goto4:

≠>Goto5?

br/>Lb14:

{AS}:

A:

S:

T=T+A-180：

T〈0=〉T=T+360⊿T〉360=〉T=T-360⊿T◢

X=X+ScosT◢

Y=Y+SsinT◢

Goto4

Lb15：

{AS}:

A:

S

T=H+A-180：

T〈0=〉T=T+360⊿T〉360=〉T=T-360⊿T◢

X=F+ScosT◢

Y=G+SsinT◢

Goto  5

6、已知置镜点坐标和置镜点到后视点坐标方位角求放样角度和放样距离的放样程序（程序名FY1）

程序清单：

F“X0”：

G“Y0”：

H“T0-H”：

Fixm：

Lb11:

{XY}?

br/>Pol（X-F，Y-G）?

br/>J〈0=〉J=J+360⊿A=J-H?

br/>A〈0=〉A=A+360⊿A“A=”◢

I“S=”◢

Goto  1

程序运行说明：

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

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

程序清单：

F“X0”：

G“Y0”：

M“XH”：

N“YH”?

br/>Pol（M-F，N-G）：

Fixm：

H=J?

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

I“S-H=”◢

Lb11:

{XY}:

Pol（X-F，Y-G）?

br/>J〈0=〉J=J+360⊿?

br/>A=J-H?

br/>A〈0=〉A=A+360⊿A“A=”◢

I“S=”◢

Goto  1

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

程序名“FY”

程序清单：

I“I=1=>1-D≠>2-D”?

br/>I=1=>Goto1:

≠>Goto2?

br/>Lb11：

F“X0”：

G“Y0”：

H“T0-H”：

Goto3?

br/>Lb12：

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=”◢

Goto3?

br/>Lb13：

{XY}?

br/>Pol（X-F，Y-G）：

Fixm?

br/>J〈0=〉J=J+360⊿A=J-H?

br/>A〈0=〉A=A+360⊿A“A=”◢

I“S=”◢

Goto  3

程序运行说明：

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

9、已知两点坐标，反算两点之间的坐标方位角和距离（程序名ZBFS）

程序清单：

A“X1=”：

B“Y1=”：

C“X2=”：

D“Y2=”：

Pol（C-A，D-B）：

Fixm：

J〈0=〉J=J+360⊿J“T（1-2）=”◢

I“S（1-2）=”

10、分别已知两条直线上一点的坐标和坐标方位角，求这两条直线的交点坐标（程序名JDZB）

程序清单：

X=（tanA“T1”×B“X1”-C“Y1”-tanD“T2”×E“X2”+F“Y2”）V÷（tanA-tanD）◢

Y=tanA（X-B）+C

11、直线线路坐标计算（ZXXY）

程序清单：

V“X0=”：

W“Y0=”：

T“T0=”：

U“K0=”：

Lb10：

Fixm：

{KD}：

X=V+（K-U）cosT+D×cos（T+90）◢

    Y=W+（K-U）sinT+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）”?

br/>Lbl1:

Fixm:

{KD}:

K“LICHENG=”：

D“WAIYI=（Z-，Y+）”：

Goto2?

br/>Lbl2:

KGoto3:

≠>Goto4?

br/>Lbl3:

“K—OVER”:

Goto1?

br/>Lbl4:

K>Q+L=>Goto3:

≠>Goto5?

br/>Lbl5:

K≤Q+B=>Goto6:

≠>Goto7?

br/>Lbl6:

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?

br/>Lbl7：

K≤Q+L-C=>Goto8:

≠>Goto9?

br/>Lbl8:

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?

br/>Lbl9:

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?

br/>X=（U-H）cosE-VsinE+F+TcosN+D×cos（O+90）◢

Y=VcosE+（U-H）sinE+G+TsinN+D×sin（O+90）◢

Goto1?

br/>子程序“BQXYS”

A:

R:

B“L01”:

C“L02”?

br/>M“M1”=B÷2-B3÷（240R2）◢

J“M2”=C÷2-C3÷（240R2）◢

P“P1”=B2÷（24R）-B4÷（2688R3）◢

Z“P2”=C2÷（24R）-C4÷（2688R3）◢

T“T1”=M+（R+P）tan（A÷2）+（Z-P）÷sinA◢

H“T2”=J+（R+Z）tan（A÷2）+（P-Z）÷sinA◢