二 Fx4800辛卜生公式坐标计算程式
一、程式：
Defm 1：N=12：
Lbl 1：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-1
Lbl 2 ‘将半径转换为曲率
Z[1]=0：Q=0：U=0:V=0：P=N+1：{L}：L“LC=”
R=1＝＞Goto3
L＞S＝＞Goto1 ‘当反算时，若计算点在线元里程范围外，则转至Lbl 1重新输入线元要素
L＜K＝＞Goto1
Goto4
Lbl 3 ‘当正算时，若计算点在线元里程范围外，则转至Lbl 1重新输入线元要素
L＜S＝＞Goto1
L＞K＝＞Goto1
Lbl 4
H=（L－S）÷N：D≠0＝＞Goto5 F= 0＝＞Goto 8
Lbl 5 判断线元是否为直线，是则转至直线计算程式Lab 8
Dsz P：Goto 6：Goto 7
Lbl 6 ‘计算各点曲率及方位角
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≠Int E＝＞Z[1]=Z[1]+CosM：Q=Q+SinM：≠＞U=U+CosM：V=V+SinM
Goto 5
X=A+Abs H÷3×（Cos C+4×Z[1]+2×U－CosO））▲
Y= B+Abs H÷3×（Sin C+4×Q+2×V－Sin O）▲
Z=－1＝＞Z=Z+180
O＞360＝＞O=O－360：≠＞O＜0＝＞O=O+360
O”A=”▲
Prog “LRZB”
Goto 2
Lbl 8 ‘直线上坐标计算程式
H=（L—S）×Z
X=A+H×Cos O▲ ‘线路中心点坐标
Y=B+H×Sin O▲ ‘线路中心点坐标
O”A=”▲ ‘线路中心点切线方位角
Prog “LRZB”
Goto 2
二、变量说明：
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＝＞Goto 2 R“D=”：O=O+W↓
X“[X]”=X+R×Cos O▲ ‘线路侧点坐标
Y“[Y]”=Y+R×Sin O▲ ‘线路侧点坐标
O＞360＝＞O=O－360：≠＞O＜0＝＞O=O+360
O“A=”▲ ‘侧点计算基点至该侧点连线方方位角
Goto 1
Lbl 2
说明：
本子程式计算曲线两侧点坐标,可以无限计算连续点坐标,前提是当提示”ANG”时,输入转向角度就是了.如果输入”-1”则回到计算中线坐标上来.重新计算下一点坐标
W----夹角（相当于曲线里面的转向角，为前一直线（或切线）的延长线至计算点的夹角
R----前一点至计算点的直线长度
前方交会是分别在已知坐标点A、B安置经纬仪向待定点P观测水平角 和检查角 ，以确定待定点P的坐标。为保证交会定点的精度，在选定P点时，应使交会角 位于30°~120°之间，最好近于90
变量对照表
数学模型变量 fx-4800P变量
P °′″ A点的水平角观测值
Q °′″ B点的水平角观测值
A
B
C °′″ 检查角
E,F m A点坐标
K,L m B点坐标
M,N m C点坐标
X,Y m P点坐标
W ″ 检查角差值
程序名：QFJH
E: F: K: L
{PQ}
A=1÷tan P: B=1÷tan Q
X=(EB+KA+L-F)÷(A+B)
Y=(FB+LA+E-K)÷(A+B)
G=X-E: H=Y-F:Pol(G,H)
J 0 Z=J+360: Z=J
M: N
G=M-E: H=N-F:Pol(G,H)
J 0 J=J+360
{C}
W=Z-J:W 0 W=W+360
W=3600(W-C)
案例：
前方交会数据
点名 坐标(m) 坐标(m) 、 水平角
A 781.227 1323.385 53°33′42″
B 769.862 1428.989 50°34′55″
C 866.294 1306.117 54°03′08″
程序步骤：
步骤 显示 按键 注释
1 E? ××× 781.227 输入A点坐标
2 F? ××× 1323.385
3 K? ××× 769.862 输入B点坐标
4 L？××× 1428.989
5 P? ××× 53 33 42 输入
6 Q？××× 50 34 55 输入
7 X=843.5338463 P点坐标
8 Y=1380.640577
9 J=42.58080857 中间计算结果
10 M? ××× 866.294 输入C点坐标
11 N？××× 1306.117
12 C? ××× 54 3 08 输入检查角
13 W=11.91808111 检查角差值
