正算及反算——精选推荐
- 格式：pdf
- 大小：144.98 KB
- 文档页数：5
正算及反算
主程序文件名TYQXJS “1.SZ=>XY”：“2.XY=>SZ”:N:Porg”DAT1”:U”X0”:V”Y0
“:O”S0”:G”F0”:H”LS”:P”R0”:R”RN”:Q:C=1/P:D=(P-R)/(2HPR):E=180/ň:N=1=>Goto 1:≠
=>Goto 2⊿
Lbl 1:{SZ}:SZ:W=Abs(S-O):Prog”SUB1”:”XS=”:X=X◢ “YS=”:Y=Y◢“F=”；F=F◢
I=Pol(X-测站点X值，Y-测站点Y值：J<0=>J=J+360⊿“I=”：I=I◢
“J=”：J=◢J
Goto 1 Lbl 2:{XY}:XY:I=X:J=Y:Prog”SUB2”:”S=”:S=O+W◢
“Z=”:Z=Z◢ Goto 2
子程序SUB1(正算子程序)
A=0.1739274226:B=0.3260725774:K=0.0694318442:L=0.3300094782:F=1-L:M=1-K:X=U+W(ACOS(G+QEKW(C+KWD))+BC
OS(G+QELW(C+LWD))+BCOS(G+QEFW(C+FWD))+ACOS(G+QEMW(C+MWD)))：Y=V+W（Asin（G+QEKW（C+KWD））
+Bsin（G+QELW（C+LWD））+Bsin（G+QEFW（C+FWD））
+Asin（G+QEMW（C+MWD）））：F=G+QEW（C+WD）+90：X=X+ZCOSF：Y=Y+ZsinF
子程序SUB2(反算子程序) T=G-90:W=Abs((Y-V)COST-(X-U)sinT):Z=0:Lbl 0:Prog”SUB1”:
L=T+QEW(C+WD):Z=(J-Y)cos L-(I-X)sinL:Abs Z<1E-6=>Goto1：
≠=>W=W+Z:Goto 0⊿
Lbl 1:Z=0:Prog”SUB1”:Z=(J-Y)/sinF
DAT1(子程序数据库) J=1=>U=”X0”:V=”Y0:O=线点里程G=线元起点方位角:H=线
元长度:P=线元起点曲率:R=线元终点曲率:Q=(左偏-1.右偏+1.直线0)
N?选择计算方式，输入1表示进行由里程、边距计算计算坐
标：输入2表示由坐标反算里程和边距。 X0？线元起点的X坐标
Y0？线元起点的Y坐标 S0？线元起点里程
LS？线元长度
R0？线元起点曲率半径 RN？线元止点曲率半径
Q？线元左右偏标志（左偏Q=-1，右偏Q=1，直线段Q=0） S？正算时所求点的里程
Z？正算时所求距中线的边距（左侧取负值，右侧取正值，
在中线上取零） X？反算时所求点的X坐标
Y？反算时所求点的Y坐标
显示部分：
XS=xxx 正算时，计算得出的X坐标 YS= xxx 正算时，计算得出的Y坐标
I\J测站到反算点的距离和方位角
S= xxx 反算时，计算得出的所求点的里程 Z= xxx 反算时，计算得出的所求点的边距
数据库数字参照例题 J=144=>U=2541283.799：V=509849.156:O=24851.857:G=278
º57′40.9″:H=347.092:P=10E45:R=10E45:Q=0⊿
J=-151=>U=2541330.087：V=509505.075:O=25198.949:G=278
º57′40.9″:H=266.45:P=10E45:R=2000:Q=-1⊿
J=-152=>U=2541365.731：V=509241.073:O=25465.399:G=275
º08′41.9″:H=755.115:P=2000:R=2000:Q=-1⊿
J=-153=>U=2541291.567：V=508494.167:O=26220.454:G=253
º30′50.51″:H=266.45:P=2000:R=10E45:Q=-1⊿ J=-154=>U=2541204.702：V=508242.329:O=26486.904:G=249
º41′50.7″:H=692.698:P=10E45:R=10E45:Q=0⊿ J=-161=>U=2540964.350：V=507592.666:O=27179.602:G=249
º41′50.7″:H=284.462:P=10E45:R=2600:Q=1⊿
J=-162=>U=2540870.542：V=507324.157:O=27464.064:G=252 º49′54.26″:H=824.053:P=2600:R=2600:Q=1⊿
J=-163=>U=2540755.081：V=506511.712:O=28288.117:G=270
º59′28.54″:H=284.462:P=2600: R=10E45:Q=1⊿
J=-171=>U=2540770.366：V=506227.700:O=28572.579:G=274
º07′32.1″:H=278.261: P=10E45: R=2300:Q=-1⊿
J=-172=>U=2540784.784：V=505949.857:O=28850.84:G=270
º39′34.83″:H=378.438: P=2300: R=2300:Q=-1⊿ J=150=>U=2541345.644：V=509507.528:O=25198.949:G=278
º57′40.9″:H=17.423: P=10E45: R=10E45:Q=0⊿ J=151=>U=2541348.359：V=509490.317:O=25216.372:G=278
º57′40.9″:H=266.45: P=10E45: R=2000:Q=-1⊿
J=152=>U=2541384.003：V=509226.315:O=25482.822:G=275
º08′41.9″:H=615.177: P=2000: R=2000:Q=-1⊿
J=153=>U=2541344.812：V=508614.815:O=26097.999:G=257
º31′16.41″:H=266.45: P=2000: R=10E45:Q=-1⊿
J=154=>U=2541275.759：V=508357.522:O=26364.449:G=253
º42′16.6″:H=757.205: P=10E45: R=10E45:Q=0⊿ J=161=>U=2541063.295：V=507630.738:O=27121.654:G=253
º42′16.6″:H=1369.482: P=5500: R=5500:Q=1⊿ J=171=>U=2540845.797：V=506282.216:O=28491.136:G=267
º58′16″:H=1002.456: P=5500: R=5500:Q=-1⊿
说明： P=10E45:R=10E45表示直线段
“F=”；F=F◢ 求G时输入，正反算时删除 +90 求G时不输入，正反算时输入
