4800坐标计算程序.docx

Fx4800辛卜生公式坐标计算程式

一、程式：

Defm1：N=12：

Lbl1：S“L1=”：A“X1=”：B“Y1=”：C“A1=”：O=C：I“R1=”：D=I：K“L2=”：T“R2=”：F=T：I=0＝＞T=0

＝＞J=1 J“L-R=”

K－S＞0＝＞Z=1：≠＞Z=－1 D≠0＝＞D=D-1 F≠0＝＞F=F-1Lbl2 ?将半径转换为曲率

Z[1]=0：Q=0：U=0:V=0：P=N+1：{L}：L“LC=”

R=1＝＞Goto3

L＞S＝＞Goto1 ?当反算时，若计算点在线元里程范围外，则转至Lbl1重新输入线元要素

L＜K＝＞Goto1Goto4

Lbl3 ?当正算时，若计算点在线元里程范围外，则转至Lbl1重新输入线元要素

L＜S＝＞Goto1L＞K＝＞Goto1Lbl4

H=（L－S）÷N：D≠0＝＞Goto5F=0＝＞Goto8

Lbl5 判断线元是否为直线，是则转至直线计算程式Lab8DszP：Goto6：Goto7

Lbl6 ?计算各点曲率及方位角

E=P÷2：G=D+（F－D）÷（K－S）×H×P：M=C+（G+D）×H×P×90÷π×J×Z

P=N＝＞O=M E≠IntE＝＞Z[1]=Z[1]+CosM：Q=Q+SinM：≠＞U=U+CosM：V=V+SinMGoto5

X=A+AbsH÷3×（CosC+4×Z[1]+2×U－CosO）▲

Y=B+AbsH÷3×（SinC+4×Q+2×V－SinO）▲Z=－1＝＞Z=Z+180

O＞360＝＞O=O－360：≠＞O＜0＝＞O=O+360

O”A=”▲Prog“LRZB”Goto2

Lbl8 ?直线上坐标计算程式

H=（L—S）×Z

X=A+H×CosO▲ ?线路中心点坐标Y=B+H×SinO▲ ?线路中心点坐标O”A=”▲ ?线路中心点切线方位角Prog“LRZB”

Goto2

二、变量说明：

N----曲线元N值S---曲线元起点里程A----起点X坐标B----起点Y坐标C 起点切线方位角

I、D----起点半径K---曲线元终点里程 T、F----终点半径J 线元左右偏判别（1右-1左）

R----曲线正反算判别（1为正算-1为反算） 注：正算即从起点算至终点，反算即从终点算至起点

L----计算点里程Z 坐标反算判别

子程式“LRZB”

?计算线路中线左右两侧点坐标

Lbl1：W=-1：{WR}：W“ANG=”：W=1＝＞Goto2 R“D=”：O=O+W↓

X“[X]”=X+R×CosO▲ ?线路侧点坐标Y“[Y]”=Y+R×SinO▲ ?线路侧点坐标O＞360＝＞O=O－360：≠＞O＜0＝＞O=O+360O“A=”▲ ?侧点计算基点至该侧点连线方方位角Goto1

Lbl2

说明：

本子程式计算曲线两侧点坐标,可以无限计算连续点坐标,前提是当提示”ANG”时,输入转向角度就是了.如果输入”-1”则回到计算中线坐标上来.重新计算下一点坐标

W 夹角（相当于曲线里面的转向角，为前一直线（或切线）的延长线至计算点的夹角

R 前一点至计算点的直线长度

前方交会是分别在已知坐标点A、B安置经纬仪向待定点P观测水平角和检查角，以确定待定点P的坐标。为保证交会定点的精度，在选定P点时，应使交会角位于30°~120°之间，最好近于90

变量对照表

数学模型变量fx-4800P变量P°′″A点的水平角观测值Q°′″B点的水平角观测值A

B

C°′″检查角E,FmA点坐标K,LmB点坐标M,NmC点坐标X,YmP点坐标W″检查角差值

程序名：QFJHE:F:K:L

{PQ}

A=1÷tanP:B=1÷tanQX=(EB+KA+L-F)÷(A+B)Y=(FB+LA+E-K)÷(A+B)G=X-E:H=Y-F:Pol(G,H)J0Z=J+360:Z=J

M:N

G=M-E:H=N-F:Pol(G,H)J0J=J+360

{C}

W=Z-J:W0W=W+360

W=3600(W-C)

案例：

前方交会数据

点名坐标(m)坐标(m)、水平角

A781.2271323.38553°33′42″

B769.8621428.989 50°34′55″

C866.2941306.11754°03′08″

程序步骤：

步骤显示按键注释

1E?×××781.227输入A点坐标

2F?×××1323.385

3K?×××769.862输入B点坐标

4