工程测量计算器程序

CASIO-48X0程序

本程序适用于CASIO-48X0,可计算一个交点的五个部分的中边桩（直－缓－圆－缓－直），缓和曲线可以是不等长的．

P=C2/(24R)-C^4/(2384R^3):Q=C/2-C^3/(240R2):M=I2/(24R)-I^4/(2384R^3): N=I/2-I^3/(240R2):T=(R+M-(P+R)cosA)/sinA+Q:E=(R+P-(P+M)cosA)/sinA+N:

L=∏RA/180-(C+I)/2

LBL 0

[J,S,F]:J:S:F

Z=J-G+T:Z≤0=>Goto 1:≠>Z≤C=>Goto 2:≠>Z=J-G+T-C:Z≤L=>Goto 3:≠>Z=J-(G-T+L+I+C):Z≤C=>Goto 4:≠>Goto 5⊿⊿⊿⊿

LBL 1

W=Z:V=B:o=B

Goto 8

LBL 2

D=Z-Z^5/(40R2C2):V=Z3/(6RC):o=90Z2/(∏RC)

Goto 6

LBL 3

D=Rsin((180Z+90C)/(∏R))+Q:V=R(1-cos((180Z+90C)/(∏

R)))+P:o=(180Z+90C)/(∏R)

Goto 6

LBL 4

D=Z-Z^5/(40R2I2):V=Z3/(6RI):W=√(D2+V2):V=tan-1(V/D):A≤0=>V=B+A+V+180: o=B+A+90Z2/(∏RI):≠>V=B+A-V+180:o=B+A-90Z2/(∏RI)⊿

Goto 7

LBL 5

W=Z:V=B+A:o=V

Goto 7

LBL 6

W=√(D2+V2):V=tan-1(V/D):A≤0=>V=B-V:o=B-o:≠>V=B+V:o=B+o⊿

Goto 8

LBL 7

H=X+Ecos(B+A)+WcosV

K=Y+Esin(B+A)+WsinV

Goto 9

LBL 8

H=X-TcosB+WcosV

K=Y-TsinB+WsinV

Goto 9

LBL 9

"X=":H=H+Fcos(o+S)◢

"Y=":K=K+Fsin(o+S)◢

Goto 0

A─偏角（左负右正，国际惯例）

B─方位角（该交点与前一交点之间的）

C－第一缓和曲线长

I－第二缓和曲线长

R－半径

G－交点桩号

X－交点X坐标

Y－交点Y坐标

J－所求点桩号

S－所求边桩与路线夹角（左负右正）

F－所求边桩与路线间距离

2、

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

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

prog "PQX"

{mn}:M"X0":N"Y0":Defm 1:Z[1]=0

Lbl 1:{ABGHVWK}:K"JL":G"JX":H"JY":V"I0":W"J0":A"R":B"LS":Fixm

P=B^2/24/A:Q=B/2-B^3/240/A^2:T=(A+P)tan(Abs W/2)+Q:C=K-T:D=C+B:F=D+π

*A*AbsW/180:E=F-B

w<0=>s=-1:≠=>s=1△

Lbl 2:{L}:L:L=0=>GOTO 1△

Lbl 6:L≤C=>O=K-L:R=G-OcosV:U=H-OsinV:Z=V:GOTO 3△

L≥F=>O=L-F+T:Z=V+W:R=G+OcosZ:U=H+OsinZ:GOTO 3△

LO=L-C:prog

"XY":Z=V+SR:R=G+(I-T)cosV-SjsinV:U=H+(I-T)sinV+jscosV:GOTO 3△

L>E=>O=F-L:prog

"XY":O=V+W:Z=O-SR:R=G+(T-I)cosO-SjsinO:U=H+(T-I)sinO+SJcosO:GOTO 3△

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

Lbl 3:Z[1]=0.002=>GOTO 5:≠=>O=0△R"XL"◢U"YL"◢Z"ZL"◢

Lbl 4:{O}:O"YC":O=0=>GOTO 2△

O=-1=>GOTO 5△

X=R-OsinZ:Y=U+OcosZ:pol(X-M,Y-N:J<0=>J=J+360△

X"XZ"◢Y"YZ"◢J"A0"◢I"DD"◢GOTO 4

Lbl 5:{XY}:X"XF":X=-1=>Z[1]=0:GOTO 4△

Y"YF":pol(X-R,Y-U:O=Icos(J-Z):AbsO≤0.001=>O"DL"◢L+O◢

O=Isin(J-Z):O"YC"◢GOTO 5:≠=>Z[1]=0.002:L=L+O:GOTO 6△

prog "XY"

U=AB:I=O-O^5/40/U^2:J=O^3/6/U-O^7/336/U^3:R=90O^2/π/U

变量名称说明：

M"X0":N"Y0"为测站点坐标

K"JL":G"JX":H"JY":V"I0":W"J0":A"R":B"LS"分别为交点桩号、坐标、起始方位

角、转角(左-右+)、半径、缓长

L、R"XL"、U"YL"、Z"ZL" 分别为输入桩号、求出中桩坐标、即时方位角O"YC"：X"XZ"◢Y"YZ"◢J"A0"◢I"DD" 分别为输入边桩距离(左-右+)、求出边桩

坐标、边桩点到测站点的方位角和平距

O"YC"=0返回计算中桩坐标

O"YC"=-1进入计算桩号过程(输入坐标X"XF" 、Y"YF"分别为求出桩号L+O,及至中桩边距O"YC"，O"DL"为桩号计算的误差距)

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

X"XF"=-1返回计算坐标过程

计算桩号过程时，需多次确认输入的坐标X"XF"、Y"YF"直到出现O"DL"桩号误差

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

3、

坐标反算程序

A：B A,B：测LbI 0

{C：D} C,D：测点或后视点坐标

E=C-A

F=D-B

K=√（E^2+F^2）◢

K:测站到测点的跳离

F=0=>F=e-20◢

≠>F=F⊿J=180-90F/AbsF-Tan^-1(E/F) ◢

J:坐标方位角

Goto 0

直线坐标计算程序

A:B:Q:J A,B:测站坐标

LbI 1 Q：测点法线方向的左右偏移值

（左“-”右“+”）

{L} L：测站到测点的距离U=A+LCosJ◢ U,V:计算点中线坐标

V=B+LSinJ◢

X=U+Qcos(J+90) ◢ X,Y:计算点边线坐标

Y=V+Qsin(J+90) ◢

Goto 1

缓和曲线坐标计算程序