二 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   检查角差值  |
|【路桥隧】
( 湘ICP备16018960号-1 ) QQ群:
GMT+8, 2026-4-4 19:41
Powered by Discuz! X3.2
© 2001-2013 Comsenz Inc.