CASIOfx4800P万用坐标放样程序Word下载.docx

CASIOfx4800P万用坐标放样程序Word下载.docx

《CASIOfx4800P万用坐标放样程序Word下载.docx》由会员分享，可在线阅读，更多相关《CASIOfx4800P万用坐标放样程序Word下载.docx（31页珍藏版）》请在冰豆网上搜索。

⒌显示“FWJ＞”即所求点与测站点的方位角，如：

36.254121即为36°

25′41.21″。

“S=”即所求点与测站点的距离。

6.1/G，1/H可以输入G的-1次方，H的-1次方。

7.当所求点里程大于终点时，显示终点坐标。

之后要求输入下一线型终点里程，终点半径。

继续计算

casio4800曲线坐标计算（源程序）供测友们参考

曲线坐标计算（源程序）

File1:

QXZB

L1 

V“ZJZ（L-,R+）”:

R:

I“LS”:

C“ZH（ZY）”:

W“FWJ”:

A“XJD”:

B“YJD”

L2 

Z[1]=90I÷

R÷

Л：

Z[2]=I÷

2－I^3÷

240÷

R2:

Z[3]=I2÷

24÷

T=Z[2]+（R+Z[3]）tan（AbsV÷

2）:

H“L”＝ЛR（AbsV÷

2×

Z[1]÷

180+2I:

E=（R+Z[3]）（cos（AbsV÷

2））-1－R

L3 

Lbl0

L4 

｛S｝:

S“K”

L5 

P=W－V+180:

Z=R+E

L6 

L=S－C

L7 

L>

0=>

Prog“ZH”:

≠>

Prog“XZ” 

L8 

X=A+DcosU◢　Y=B+DsinU◢

L9 

Goto 

File2 

ZH

L≥I　=>

Prog“QZ”:

Prog“HY”

File3 

QZ

J=H－I

J 

=>

Prog“HZ”:

Prog“YH” 

File4 

HZ

L≥H=>

D=T+L－H:

U=W:

F=W◢　≠>

Prog“YZ” 

File5 

XZ

D=T+C－S:

U=P:

F=W－V◢

File6 

HY

M=L－L^5÷

40÷

R2÷

I2:

N=L^3÷

6÷

I－L^7÷

336÷

R^3÷

I^3:

K=T－M:

D=√（K2+N2）:

G=tan-1（N÷

K）

V<

0=>

U=P+G:

F=W－V－90（S－C）2÷

I÷

Л◢　 

U=P－G:

F=W－V+90（S－C）2÷

Л◢　

File7 

YH

Q=H÷

2－L

O=Abs180Q÷

Л:

D=√（R2+Z2－2RzcosO）

G=Abssin-1（RsinO÷

D）

0 

Prog“QX”:

Prog“QS” 

File8 

YZ

L=H－L

0 

U=W－G:

F=W+90（H+C－S）2÷

U=W+G:

F=W－90（H+C－S）2÷

Л◢

Casio 

fx-4800P平曲线坐标计算程序 

[推荐]

程序目的：

依平曲线要素计算直线、圆曲线、缓和曲线的任意中桩、左、右桩坐标。

程序说明：

K0：

起始桩号X0：

起始X坐标Y0：

起始Y坐标ALF：

起始方位角R：

半径LS：

缓和曲线长N：

曲线左转N=1，右转N=2K：

待求桩号LL、LR：

左、右桩距离Q：

左、右桩与中线斜交角求得XZ、YZ、XL、YL、XR、YR分别为中桩、左、右桩坐标。

一、直线段

文件名：

ZX（COMP）

程式：

L”K0”:

O”X0”:

P”Y0”:

W”ALF”:

Lbl0:

X”XZ”=O+（K-L）cosW◢

Y”YZ”=P+（K-L）sinW◢

{B}:

S”XL”=X-B”LL”cos（W+Q）◢

T”YL”=Y-Bsin（W+Q）◢

{C}:

U”XR”=X+C”LR”cos（W+Q）◢

V”YR”=Y+Csin（W+Q）◢

Goto0

注：

在程序执行过程中，赋给的要素变数的值被固定不变，可对变数（K、LL、LR）赋予不同值，迅速求得所需坐标。

二、圆曲线段

YQX（COMP）

Lbl1:

J=（-1）^N（K-L）÷

R×

180÷

π:

D=2Rsin（（-1）^NJ÷

X”XZ”=O+Dcos（W+J÷

2）◢

Y”YZ”=P+Dsin（W+J÷

S”XL”=X-B”LL”cos（W+J+Q）◢

T”YL”=Y-Bsin（W+J+Q）◢

U”XR”=X+C”LR”cos（W+J+Q）◢

V”YR”=Y+Csin（W+J+Q）◢

Goto1

若没有直接HY点方位角，则ALFHY=ALFZH±

Ls/2/R×

180/π,（左转-，右转+）。

三、缓和曲线

HHQX（COMP）

M”LS”:

Lbl3:

I=（-1）^N×

（K-L）^2÷

M÷

6×

D=（K-L）-（K-L）^5÷

90÷

（RM）^2:

X”XZ”=O+Dcos（W+I）◢

Y”YZ”=P+Dsin（W+I）◢

S”XL”=X-B”LL”cos（W+3I+Q）◢

T”YL”=Y-Bsin（W+3I+Q）◢

U”XR”=X+C”LR”cos（W+3I+Q）◢

V”YR”=Y+Csin（W+3I+Q）◢

Goto3

1、坐标计算方法是根据偏角法原理；

2、缓和曲线（ZH~HY或YH~HZ）以ZH（或HZ）为起始点；

3、平曲线左转（ZH~HY段N=1，YH~HZ段N=2），曲线右转（ZH~HY段N=2，YH~HZ段N=1）。

Casiofx-4800P计算机具有强大的储存记忆、程序设计功能，能为众多专业技术工作提供大量复杂程序计算。

笔者根据工程测量中的需要汇编以上程序，本着严密准确、计算程式简洁，避免冗长复杂化的计算过程。

以上三部分平曲线计算程序有待继续研究成更简便的三合一程序。

笔者另编写Casiofx-4800P方位角计算程序、竖曲线设计高程计算程序，及依坐标求算里程桩号及至中桩距离的程序，Excel坐标计算程序、Excel平曲线要素表计算程序、Excel竖曲线计算程序均具大批量计算并打印功能，程序准确实用、方法简捷，下次见面。

望大家不吝赐教，谢谢！

本程序用于用全站仪测设中桩和边桩，采用坐标变换的方法编制而成，在使用fx-4500计算器时，可直接输入以下程序便可使用；

而在使用fx-4800计算器时，见fx-4800程序清单。

注意：

因程序的0（零）和O易混淆,所以字母O用较大的全角符号，零仍用标准字。

程序输入时应连续输入，无空格。

而且输入时文件名需和给出的一样，否则会出现调用文件错误。

变量说明：

<

程序运算符定义>

X输入摆站点的x坐标C输入第一个曲线角桩x坐标I输入中间曲线角桩的x坐标

Y输入摆站点的y坐标D输入第一个曲线角桩y坐标J输入中间曲线角桩的y坐标

E最后一个曲线角桩x坐标R中间曲线半径L中间曲线缓和曲线长

F最后一个曲线角桩y坐标K第一个K为中间角桩桩号S中间角桩曲线总长，检验用

K第二个K为要测设的桩号B测中桩时输入0，测边桩时输入宽度（左输负值，右输正值）

程序清单

注释

文件:

F1文件名：

1

L1

L2

L3

Defm8

XYCDIJEFRLK:

Fixm:

Z[2]=0:

ProgB

Z[5]=0.5L-LXy3/240R2:

Z[3]=Z[5]+（R+L2/24R）AbstanT┛2:

Z=K-Ans:

S=RAbsTπ/180+L:

S◢H=Z+Ans:

G=RL:

Prog2

F2文件名：

2

L3

L4

L5

L6

L7

L8

L9

L10

L11

L12

{K}

Lbl6:

Norm:

K（-1）XyZ[2]>

（H-L）（-1）XyZ[2]＝>

Goto4⊿FracZ[2]┛2=0=>

Z[6]=Z:

Z[6]=H⊿A=90T/AbsT

B<

0＝>

N=0:

N=1⊿B=AbsB

０=Abs（K-Z[6]:

Z[6]（-1）XyZ[2]＝>

Goto2⊿Q=Z[3]+０:

Z[4]=B

B≠0＝>

ProgC⊿ProgD:

P=0

ProgE:

Lbl2:

０≥L＝>

Goto3⊿

Q=Z[3]-０+０Xy5/40G2:

Z[4]=０Xy3/6G-０Xy7/336GXy3:

ProgD:

P=B

ProgF⊿ProgE:

Q=Z[3]-Z[5]:

Z[4]=R+L2/24R:

ProgD

Z[8]=U+A（2（０-L）/πR-1+L/πR）:

ProgG:

P=R⊿ProgE:

Lbl4:

M=C:

０=D:

C=E:

D=F:

E=M:

F=０:

Z[2]=Z[2]+1:

ProgB:

Goto6

F3文件名：

0

L1

Pol（M,０）:

W=W+360

F4文件名：

B

L2

M=E-I:

０=F-J:

Prog0:

T=W:

M=I-C:

０=J-D:

U=W:

T=T-W

T<

-180＝>

T=T+360⊿T>

180＝>

T=T-360

F5文件名：

C

L1A=-90（-1）Xy（N+Z[2]）

F6文件名：

D

L1Z[7]=I-QcosU+Z[4]cos（U+A）:

Z[1]=J-QsinU+Z[4]sin（U+A）

F7文件名：

E

M=Z[7]+PcosZ[8]-X:

０=Z[1]+PsinZ[8]-Y:

Prog0

W◢Fix3:

V◢Norm

F8文件名：

F

Z[8]=U+A+A０2/πG+90+A（-1）Xy（N+Z[2]）

F9文件名：

G

P=R+AB（-1）Xy（N+Z[2]）/90

此处L2句加入：

M+X◢０+Y◢可显示测点坐标值

casio4800超高及高程（源程序）供参考

超高及高程（源程序）

L1:

L“LS”:

Z“ZH”:

H“HZ”:

V“RS”:

F“I1”:

G“I2”:

O“BPH”:

N“BP”:

T=Vabs（G－F）÷

2:

X“X（R+1,L-1）”:

P“P（R+,L-1）”

L2:

Lbl6:

｛S｝

L3:

S+T－N<

D=O+（S－N）F:

S－N<

D=O+（S－N）F－U（S－N+T）2÷

2÷

V:

S－N≤T 

D=O+（S－N）G－U（T－S+N）2÷

D=O+（S－N）G 

D“H=”◢

L4:

R<

E= 

M= 

Goto1:

R<

E= 

L5:

Lbl 

1 

M≥L 

C=L:

Goto2:

C=M

L6:

2 

PX>

Goto3:

Goto4

L7:

4 

S≤Z 

I=0.02:

Goto5:

S<

Z+C 

I=0.02－（0.02－EPX）÷

C×

（S－Z）:

S<

H－C 

I=EPX:

H 

I=EPX+（0.02－EPX）÷

（S－H+C）:

Goto5

L8:

Lbl3:

Z[1]=0.04÷

（EPX+0.02）×

C:

S≤Z+Z[1] 

I=0.02:

I=0.02+（EPX－0.02）÷

（C－Z[1]）×

（S－Z－Z[1]）:

H－Z[1] 

I=EPX－（EPX－0.02）÷

L9:

Lbl5:

{B}:

W“HB”=D－IB◢

L10:

Goto6

平曲线半径 

缓和曲线长 

P：

偏转，Z=－1,Y=1

X“XL”：

线路（Z,Y）左线＝－1,右线＝1 

Z“ZH”:

直缓点里程

H“HZ”：

缓直点里程 

V“RS”：

竖曲线半径 

F“I1”：

前坡度值（带符号）

G“I2”：

后坡度值（带符号） 

O“BPH”：

变坡点高程 

N“BP”：

变坡点里程

T：

切线长 

E：

最大超高值 

C：

渐变段长度（LC） 

I：

横坡值 

M：

最大渐变长度（Lmax） 

B：

距中心距离 

W“HB”：

边部高程

竖曲线标高计算程序（4800）

（适用于4800） 

此程序为竖曲线标高计算程序，程序用变坡点高程和桩号来计算坡度，使高程计算比较精确，程序运行时，

先按顺序输入三点的桩号、高程和中间曲线的半径，然后可输入桩号求相应点的高程，当桩号到达曲线尾和下一曲线间的直线段时，

会出现“AFTERQX”的提示，这时可继续输入桩号求标高，也可输入一负值进行输入下一变坡点数据，输入负值后会出现提示要求输入R，

此R为第三个变坡点的半径，继续要求输入的EF则为第四个变坡点的桩号和高程。

程序自动将前次的后两个变坡点的数据交换到计算时所

用的第一第二个变坡点，而不用每次都输入。

但要注意的是，不能输入超过下一个曲线的桩号，否则会出现标高错误，最好在出现

“AFTERQX”提示后便进行下一变坡点数据的输入。

本程序在输出标高后可设计简单的横坡度、超高坡度计算程序，可一次输出断面中的标高。

变量说明：

<

A 

前一个变坡点桩号 

B 

前一个变坡点高程 

R 

中间变坡点半径

C 

中间变坡点桩号 

D 

中间变坡点高程 

G 

前段坡度

E 

后一个变坡点桩号 

F 

后一个变坡点高程 

后段坡度

L 

中间曲线的切线长 

I 

所求点的高程 

切线至竖曲线的竖向距离

K 

所求桩号 

程 

序 

清 

单

文件1 

程序名：

SQX

ABCD:

{ERF}:

REF:

G=（D-B）÷

（C-A）:

H=（F-D）÷

（E-C）:

L=R×

Abs（G-H）÷

2←

Lbl1:

K≤0＝>

Goto2⊿K≤C+L＝>

L=-AbsL:

Prog"

SHU"

M=D+（K-C）H+J:

M=D+（K-C）H:

AFTERQX"

⊿L=AbsL:

K≤C＝>

M=B+（K-A）G+J⊿K≤C-L＝>

M=B+（K-A）G⊿M=1000M:

SSWR"

I=M÷

1000◢

K=K+20:

A=C:

B=D:

文件2 

程序名：

SHU

J=（K-C+L）2÷

（2R）:

G-H>

J=-J⊿

文件3 

SSWR

M-IntM<

0.5＝>

M=IntM:

M=IntM+1⊿

CASIO中计算结果直接显示为度、分、秒

在4500和4800中将计算结果直接显示为度、分、秒

本人所写，请转帖时保持完整！

！

IntE+Int（FracE×

e2）/60+Frac（E×

e2）/36

显示结果是将E中的数据显示为度、分、秒的形式。

E是计算器中的存储器，e是计算器下的EXP键

它可以放在其它程序中，也可以单独做为一个子程序供其它程序调用！

平曲线坐标、反算桩号计算程序

prog"

PQX"

{mn}:

M"

X0"

N"

Y0"

Defm1:

Z[1]=0

{ABGHVWK}:

K"

JL"

JX"

JY"

V"

I0"

W"

J0"

R"

B"

LS"

Fixm

P=B^2/24/A:

Q=B/2-B^3/240/A^2:

T=（A+P）tan（AbsW/2）+Q:

C=K-T:

D=C+B:

F=D+π*A*AbsW/180:

E=F-B

w<

0=>

s=-1:

≠=>

s=1△

{L}:

L:

L=0=>

GOTO1△

L≤C=>

O=K-L:

R=G-OcosV:

U=H-OsinV:

Z=V:

GOTO3△

L≥F=>

O=L-F+T:

Z=V+W:

R=+OcosZ:

U=H+OsinZ:

L<

D=>

O=L-C:

XY"

Z=V+SR:

R=G+（I-T）cosV-SJsinV:

U=H+（I-T）sinV+JScosV:

L>

E=>

O=F-L:

O=V+W:

Z=O-SR:

R=G+（T-I）cosO-SJsinO:

U=H+（T-I）sinO+SJcosO:

O=180（L-D+B/2）/π/A:

I=AsinO+Q:

J=P+A（1-cosO）:

Z=V+SO:

R=G-（T-I）cosV-SJsinV:

U=H-（T-I）sinV+SJcosV

Z[1]=0.002=>

GOTO5≠=>

O=0.00000△:

XL"

?

YL"

Z"

ZL"

{O}:

YC"

O=0:

=>

GOTO2△

O=-1=>

GOTO5△

X=R-OsinZ:

Y=U+OcosZ:

pol（X-M,Y-N:

J<

J=J+360△

X"

XZ"

Y"

YZ"

J"

A0"

I"

DD"

GOTO4

Lbl5:

{XY}:

XF"

X=-1=>

Z[1]=0:

GOTO4△Y"

YF"

pol（X-R,Y-U:

O=Icos（J-Z）:

AbsO≤0.001=>

DL"

L+O?

O=Isin（J-Z）:

Z[1]=0.002:

L=L+o:

GOTO6△

U=AB:

I=O-O^5/40/U^2:

J=O^3/6/U-O^7/336/U^3:

R=90O^2/π/U

变量名称说明：

M、N为测站点坐标 

分别为交点桩号、坐标、起始方位角、转角（左-右+）、半径、缓长

L、R"

、U"

、Z"

分别为输入桩号、求出中桩坐标、即时方位角 

：

分别为输入边桩距离（左-右+）、求出边桩坐标、边桩点到测站点的方位角和平距

=0返回计算中桩坐标

=-1进入计算桩号过程（输入坐标X"

、Y"

分别为求出桩号L+O,及至中桩边距O"

，O"

为桩号计算的误差距）

=约0.00001可求中桩点到测站点方位角和平距 

=-1返回计算坐标过程

计算桩号过程时，需多次确认输入的坐标X"

、Y"

直到出现O"

桩号误差距较小时即可得出准确的桩号

［摘要］程序中采用的计算方法是基于不受线型限制的复化辛甫生公式，该程序适合于线路各种线型的中线和边线点位坐标计算。

［关键词］复化辛甫生公式坐标计算程序

［分类号］P204

1前言

在高速公路的路面施工中，为了保证路面的平整度和宽度，施工监理部门常要求施工单位5米加密一个中线点，并且中线点和边线点要用全站仪或光电测距仪放样。

设计文件中仅给了部分中线点位坐标，这些点位坐标远远不能满足施工放样的要求，由此，要进行大量的内业计算工作，为了充分发挥现有施工单位CASIOfx-4500P可编程计算器的作用，本文以复化辛甫生公式为基础，编制了适合各种线型的线路中线和边线点位坐标计算通用程序。

该程序在我院举办的高级测量工培训班上推广受到了良好的效果。

2线路中线和边线点位坐标计算公式

2.1计算线路中线坐标的复化辛甫生公式

式中

上式各符号的意义：

XA——曲线元起点A的x坐标；

YA——曲线元起点A的y坐标；

αA——曲线元起点A切线的坐标方位角；

αk——里程