计算是测量技能之笔，仪器应用是测量技能之眼

计算是测量技能之笔，仪器应用是测量技能之眼
我觉计算应放在测量的第一位，计算一天学不会，而用仪器最多一天可搞定，所以我们多专研测量计算，那么CISAO计算编程也是首当其冲的必要工作，现我把我的5800测量线路计算程序与大家分享一下
正算程序
01---“ZB”
"K="?K："B="?Q
6→N: Prog“SJK-PM”：
Lbl 1
Prog“YS-PM”
Z[6N]-Z[6] →C
IF K＞C+L+Z[8]+M Then N+1→N：Goto 1:IfEnd
Lbl 2
K≤C＝＞Goto 3⊿
K≤C+L＝＞Coto 4⊿
If K≤C+L+Z[8] Then Goto 5：IfEnd :Goto 6⊿
Lbl 3 进入直线段
K-C→S
Z[10]+SCosA→X： Z[11]+SSinA→Y
X+QCos（A+O）→X： Y+QSin（A+O）→Y
Goto Z
Lbl 4 进入一缓段
K-C→S
S-S5／（40(RL)2）+S9/(3456R4L4) →I:1→G:
S3／(6RL)-S7／(336(RL)3)+S11/(42240R5L5) →J
Z[10]+GICosA-JESinA→X
Z[11]+GISinA+JECosA→Y
X+QCos(A+90ES2／(RLπ)+O) →X
Y+QSin(A+90ES2／(RLπ)+O) →Y
Goto Z
Lbl 5 进入圆曲线段
K-C-L→S
90ES／（Rπ）→I： 2Rsin(AbsI) →J
Z[14]+JCos（A+90EL／（Rπ）+I）→X
Z[15]+JSin（A+90EL／（Rπ）+I）→Y
X+Qcos（A+90EL／（Rπ）+180ES／（Rπ）+O）→X
Y+QSin（A+90EL／（Rπ）+180ES／（Rπ）+O）→Y
Goto Z
Lbl 6 进入二缓段
C+L+Z[8]+M-K→S
S∧3／（6RM）-S∧7/(336(RM)∧3)+ S11/(42240R5M5) →J
S-S∧5/(40(RM)∧2) +S9/(3456R4M4):-1→G
Z[12]+GICosB-JESinB→X
Z[13]+GISinB+JECosB→Y
X+QCos(B-90ES∧2/(RMπ)+O) →X
Y+QSin(B-90ES∧2/(RMπ+O) →Y
Goto Z
Lbl Z:
Cls : "K="ocate 3，1，K
"B="ocate 3，2，Q
"X="ocate 3，3，X
"Y="ocate 3，4，Y

“YS-PM”
Lbl 1
Z[6N+3] → L: Z[6N+4] →R： Z[6N+5] →M
Pol(Z[6N+1]-Z[6(N-1)+1]，Z[6N+2]-Z[6(N-1)+2])
J→A:A＜0＝＞A+360→A⊿
Pol(Z[6(N+1)+1]-Z[6N+1]，Z[6(N+1)+2]-Z[6N+2])
J→B:B＜0＝＞B+360→B: B-A→D
If Sin D≥0 Then 1→E: Else -1→E :IfEnd⊿
If Abs D＞180 Then 360-Abs D →D: Else Abs D →D :IfEnd
L／2-L∧3／(240R∧2)+L∧5/（34560R∧4）-L∧7/(8386560R∧6) →Z[2]: 一缓切线增量
M／2-M∧3／(240R∧2)+M∧5/（34560R∧4）-M∧7/(8386560R∧6) →Z[3] 二缓切线增量
L∧2／(24R)-L∧4/（2688R∧3）-L∧6/(506880R∧5) →Z[4]:一缓内移值
M∧2／(24R)-M∧4/（2688R∧3）-M∧6/(506880R∧5) →Z[5] 二缓内移值
Z[2]+(R+Z[5]-(R+Z[4])COS(D)/SIN(D) →Z[6] 一缓切线长
Z[3]+(R+Z[4]-(R+Z[5])COS(D)/SIN(D) →Z[7]二缓切线长
( D-180L*(2Rπ)- 180M*(2Rπ))Rπ／180→Z[8]圆曲线长
Z[6N+1]-Z[6]CosA→Z[10]: Z[6N+2]-Z[6]SinA→Z[11] (ZH点XY坐标)
Z[6N+1]+Z[7]CosB→Z[12]: Z[6N+2]+Z[7]SinB→Z[13] (HZ点XY坐标)
Lbl 2
If L=0: Then Z[10] →Z[14]: Z[11] →Z[15]: Goto 4: IfEnd
L-L∧5／(40(RL)∧2)+ L∧9／(3456(RL)∧4) →I
L∧3／(6RL)-L∧7／(336(RL)∧3) +L∧11／(42240(RL)∧5) →J:
1→G
Z[10]+GICosA-JESinA→Z[14]: Z[11]+GISinA+JECosA→Z[15] (HY点XY坐标)
Lbl A
Z[14]+RCos(A+90EL／(Rπ)+90E) →Z[18] (圆心点X坐标)
Z[15]+RSin(A+90EL／(Rπ)+90E) →Z[19] (圆心点Y坐标)
Lbl 3
If M=0: Then Z[12] →Z[16]: Z[13] →Z[17]: Goto 4: IfEnd
M-M∧5／(40(RM)∧2)+ M∧9／(3456(RM)∧4) →I
M∧3／(6RM)-M∧7／(336(RM)∧3) +M∧11／(42240(RM)∧5) →J:
-1→G
Z[12]+GICosB-JESinB→Z[16]： Z[13]+GISinB+JECosB→Z[17] (YH点XY坐标)
Lbl 4
Return
★“SJK-PM” 平面曲线各交点DJ数据库
O；构筑物与线路交角
Z[6(N-1)]、Z[6(N-1)+1] 、Z[6(N-1)+2] 分别是起点里程、X、Y坐标
Z[6N]、Z[6N+1]、Z[6N+2]、Z[6N+3]、Z[6N+4]、Z[6N+5]分别是各DJ点里程、X、Y坐标、一缓长、半径、二缓长
……….
从N=6(Z[30]…..)开始…….

反算程序
“FQ”
"X="?→X："Y="?→Y:N=10: Prog“SJK-PM”：
Lbl 1
Prog“YS-PM”
X→I： Y→J： Z[12] →K： Z[13] →Q
Prog“FWJ”
K-B→K：If CosK＞0: Then N+1→N:Goto 1:IfEnd
Lbl 2
X→I： Y→J: Z[10] →K: Z[11] →Qrog“FWJ”： K-A→K：If CosK≤0:Then Goto 3：Else Goto 4:IfEnd
Lbl 3进入直线段
X→I： Y→J: Z[10] →K: Z[11] →Q：Prog“FWJ”
K-A→K： ISinK→Q：ICosK→H:Z[6N]-Z[6]+H→K:Goto S
Lbl 4进入HY点切线方位角判定
X→I： Y→J: Z[14] →K:Z[15] →Qrog“FWJ”
K-（A+90EL/（Rπ））→K [ HY点切线方位角A+Ls/（2R）*180/π]
If CosK≤0:Then Goto 5：Else Goto 6:IfEnd
Lbl 5进入一缓段
Prog“HQ”：Goto S
Lbl 6进入YH点切线方位角判定
X→I： Y→J： Z[16] →K：Z[17] →Q：Prog“FWJ”
K-（A+90EL/（Rπ）+180EZ[8]/（Rπ））→K [YH点切线方位角A+90EL/（Rπ）+180EZ[8]/（Rπ）]
If CosK≤0:Then Goto 7：Else Goto 8:IfEnd
Lbl 7进入圆曲线段
Z[14] →I： Z[15] →J: Z[18] →K:Z[19] →Qrog“FWJ”： K→O
X→I： Y→J:Z[18] →K:Z[19] →Qrog“FWJ”：K→P
E(R-I) →Q:Abs(P-O) →K:K＞180＝＞360-K→K
Z[6N]-Z[6]+L+RπAbsK/180→K：Goto S
Lbl 8
Prog“QH”：Goto S
Lbl S: Cls : "X="ocate 3，1，X
"Y="ocate 3，2，Y 显示ZH点X坐标
"K="ocate 3，3，K 显示ZH点Y坐标
"B="ocate 3，4，Q

“FWJ” Pol(I-K，J-Q)
K=J:K＜0＝＞K=K+360

“HQ”
I=X：J=Y：K=Z[10]：Q=Z[11]
Prog“FWJ”
P=K：K=P-A：H=ICosAbsK
I=（L-（L-L∧3/（40R∧2）））÷（L-L∧3/（40R∧2））∧5×H∧5
H=H+AbsI：V=30H∧2/（RLπ）：J=（H-H∧5/（40（RL）∧2））/CosV
Pol(Z[10]+Cos(A+VE)J -X， Z[11]+JSin(A+VE)-Y)
G=P-(A+VE):SinG＞0＝＞G=1:≠＞G=-1⊿
Z=H+(AbsI)Sin3V GE:Z=90Z∧2/(RLπ):
Z=H+(AbsI)SinZ GE:W=90Z∧2/（RLπ）：F=（AbsI）SinW
S=Ftan(Abs(W-3V)÷3×2)=(AbsI)CosW
U=F÷Cos(Abs(W-3V)÷3×2)=Q+SGE
K=Z[6N]-Z[6]+(H+UGE)=QG
Lbl S

“QH”I=X:J=Y:K=Z[12]=Z[13]
Prog“FWJ”
P=K：K=P-（B+180）：H=ICosAbsK
I=（M-（M-M∧3/（40R∧2）））/（M-M∧3/（40R∧2））∧5×H∧5
H=H+AbsI：V=30H∧2/(RMπ):J=(H-H∧5/(40(RM)∧2))/CosV
Pol（X-( Z[12]+Jcos（B+180-VE）)，Y-( Z[13]+JSin（B+180-VE）)）:
G=P-(B+180-VE)：SinG＞0＝＞G=-1：≠＞G=1⊿
Z=H+(AbsI)Sin3V GE：Z=90Z∧2/（RMπ）：Z=H+（AbsI）SinZ GE：W=90Z∧2/（RMπ）：F=（AbsI）SinW
S=Ftan(Abs(W-3V)÷3×2)=(AbsI)CosW
U=F÷Cos（Abs（W-3V）÷3×2）：Q=Q+SGE
K=Z[6N]-Z[6]+L+Z[8]+M-（H+UGE）：Q=QG：
Lbl S

 要素计算
02-“YSJS”
50→DimZ "DUANG="?N:5+N→N rog“SJK-PM”：区段选择后调用数据库文件计算可直接修改输入
Z[6(N-1)+1]→A　起点N坐标
Z[6(N-1)+2]→B　起点E坐标
Z[6N]→K　交点DJ点里程
Z[6N+1]→X　交点DJ点X坐标
Z[6N+2]→Y　交点DJ点Y坐标
Z[6N+3]→L　 一缓长
Z[6N+4]→R　 半径
Z[6N+5]→M
Z[6(N+1)+1]→W　下交点DJ点(终点)X坐标
Z[6(N+1)+2]→V　下交点DJ点(终点)Y坐标
Lbl 1
"ZD1-X="?A：A→Z[6(N-1)+1]:"ZD1-Y="?B:B→ Z[6(N-1)+2]
"JD-K="?K : K→Z[6N]:"JD-X="?X：X→Z[6N+1]:"JD-Y="?Y: Y→Z[6N+2]
"LS1="?L : L→Z[6N+3]:"R="?R：R→Z[6N+4]："LS2="?M： M→Z[6N+5]
"ZD2-X="?W：W→Z[6(N+1)+1]:"ZD2-Y="?V:V→ Z[6(N+1)+2]
Lbl 2
Prog“YS-PM”
Lbl 3
Cls : "a="ocate 3，1，ED 显示转角
"T="ocate 3，2，Z[6] 显示切线长
"LC=":Locate 4，3，Z[8] 显示圆曲线长
"E=":Locate 3，4，(R+Z[4])/ Cos(D/2)-R⊿ 显示外距长
Lbl 4
Cls : "ZH-K=":Locate 6，1，Z[6N]-Z[6] 显示ZH点里程
"ZH-X=":Locate 6，2，Z[10] 显示ZH点X坐标
"ZH-Y=":Locate 6，3，Z[11] 显示ZH点Y坐标
"ZH-Fa=":Locate 7，4，A⊿ 显示ZH点切线方位角
Lbl 5
Cls : "HY-K=":Locate 6，1，Z[6N]-Z[6]+L 显示HY点里程
"HY-X=":Locate 6，2，Z[14] 显示HY点X坐标
"HY-Y=":Locate 6，3，Z[15] 显示HY点Y坐标
"HY-Fa=":Locate 7，4，A+E180L/(2Rπ)⊿ 显示HY点切线方位角
Lbl 6
Cls : "YH-K=":Locate 6，1，Z[6N]-Z[6]+L+Z[8] 显示YH点里程
"YH-X=":Locate 6，2，Z[16] 显示YH点X坐标
"YH-Y=":Locate 6，3，Z[17] 显示YH点Y坐标
"YH-Fa=":Locate 7，4，A+E180L/(2Rπ)+ E180Z[8]/(Rπ)⊿ 显示HY点切线方位角
Lbl 7
Cls : "HZ-K=":Locate 6，1，Z[6N]-Z[6]+2L+Z[8] 显示HZ点里程
"HZ-X=":Locate 6，2，Z[12] 显示HZ点X坐标
"HZ-Y=":Locate 6，3，Z[13] 显示HZ点Y坐标
"HZ-Fa=":Locate 7，4，B 显示HZ点切线方位角

[ 本帖最后由 y01h02f03 于 2009-5-14 11:08 编辑 ]