纯手打 求支持
本套程序包含的主程序及子程序：
1、主程序：0---ZXJS
2、坐标计算判断数据库：ZXDATA
3、坐标正算子程序：SUB1
4、坐标反算子程序：SUB2
5、竖曲线判断数据库程序：ZXGC
6、超高段判断数据库程序：ZXGC
7、竖曲线高程计算子程序：SUB3
8、竖曲线超高计算子程序：SUB4
要更换新路线数据，只需要更改ZXDATA、ZXGC、ZXGC这三个数据库即可
1、主程序0---ZXJS：
26→DimZ↙ 为Z数组变量定义范围Z[1]~Z[26]
180÷π→K↙ 弧度转换到角度系数
If N=1:Then Goto A:IfEnd
If N=11:Then Goto B:IfEnd
Lbl A 主线计算主程序起始
"ZH="?S:"E="?E:"JL="?Z:Prog "ZXDATA":Abs(S-O)→W:1÷P→C:(P-R)÷(2HPR)→D:Prog "SUB1"↙
"ZX=":X◢ ←计算出的X坐标
"ZY=":Y◢ ←计算出的Y坐标
Prog "ZXGC" ←判断所求高程点所在竖曲线
Prog "ZXCG" ←判断所求高程点所在超高段
Prog "SUB3" ←计算路线中心高程
Prog "SUB4" ←计算路线边桩高程
"H=":H◢←计算出的Z坐标
"TJ=":Z[5]?DMS◢ 曲线当前点的切线←方位角
Goto A
Lbl B
?X:?Y?:90→E:X→I:Y→J:Prog "SUB2":"Kxxx+xxx=":O+W→S◢:"D=":Z◢
Goto B
2、坐标计算判断数据库：ZXDATA
遵循 直线-缓和曲线-圆曲线-缓和曲线-直线 法则 无缓和曲线可以省略缓和曲线参数
参数含义：
U:平曲线X起点坐标
V:平曲线Y起点坐标
O:平曲线起点里程
G:平曲线起点方位角
H:曲线长度(直线就是直线长度，缓和曲线就是缓和曲线长度，圆曲线就是圆曲线长度)
P:曲线起点半径(直线为与缓和曲线为10^45，圆曲线为圆曲线半径)
R:曲线终点半径 当曲线为圆曲线时，P=R(直线为10^45，缓和曲线与圆曲线为圆曲线半径)
Q:路线转角方向(左转为负，右转为正，直线为0)
输入格式：
If S≥曲线起点里程 And S≤曲线终点里程:Then 曲线起点X坐标→U:曲线终点坐标→V:曲线起点里程→O:曲线起点方位角→G:曲线长度→H:曲线起点半径→P:曲线终点半径→R:+1、-1或0→Q:IfEnd
3、坐标正算子程序：SUB1
0.1739274226→Z[1]: 0.3260725774→Z[2]:0.0694318442→Z[3]:0.3300094782→Z[4]:1-Z[4] →F:1-Z[3] →M
U+W(Z[1]cos(G+QKZ[3]W(C+Z[3]WD))+Z[2]cos(G+QKZ[4]W(C+Z[4]WD))+Z[2]cos(G+QKFW(C+FWD))+Z[1]cos(G+QKMW(C+MWD))→X
V+W(Z[1]sin(G+QKZ[3]W(C+Z[3]WD))+Z[2]sin(G+QKZ[4]W(C+Z[4]WD))+Z[2]sin(G+QKFW(C+FWD))+Z[1]sin(G+QKMW(C+MWD))→Y
G+QKW(C+WD)→Z[5]
Z[5]+E→F
X+Zcos(F)→X
Y+Zsin(F)→Y
注：Z[1] Z[2] Z[3] Z[4] 为计算参数 不需要改变 因为我也不知道怎么定的 应该是因为本身计算公式是积分法 里面的常数相乘的得数
4、坐标反算子程序：SUB2
G-90→T:Abs((Y-V)cos(T)-(X-U)sin(T)) →W:0→Z
Lbl 0
Prog “SUB1”:T+QKW(C+WD) →L:(J-Y)cos(L)-(I-X)sin(L)→Z
If Abs(Z)<1e^(-4):Then Goto 1:Else W+Z→W
Goto 0
Lbl 1
0→Z:Prog “SUB1”:(J-Y)÷sin(F
)→Z
5：竖曲线判断数据库程序：ZXGC
参数含义：
B：变坡点桩号
G：变坡点高程
I：变坡点前坡度(区分正负，带符号)
J：变坡点后坡度(区分正负，带符号)
R：竖曲线半径
输入格式：
If S≥桩号1 And S≤桩号2:Then 桩号1与桩号2之间的变坡点桩号→B:变坡点高程→G:变坡点前坡度→I:变坡点后坡度→J:竖曲线半径→R:IfEnd
If S≥桩号2 And S≤桩号3:Then 桩号2与桩号3之间的变坡点桩号→B:变坡点高程→G:变坡点前坡度→I:变坡点后坡度→J:竖曲线半径→R:IfEnd
If S≥桩号3 And S≤桩号4:Then 桩号3与桩号4之间的变坡点桩号→B:变坡点高程→G:变坡点前坡度→I:变坡点后坡度→J:竖曲线半径→R:IfEnd
注：桩号1、桩号2、桩号3、桩号4 为两个变坡点之间的[直线]上任意一点，不可以是边坡半径曲线上的桩号
6、超高段判断数据库程序：ZXGC
参数含义：
I：超高渐变起点坡度(1.5%就输入1.5)
J：超高渐变终点坡度(-3%就输入-3)
O：超高渐变起点里程
A：超高渐变段长度
Z[10]：超高渐变参数：0为单向横坡无变化，1为起点、终点坡度有变化
输入格式：
①单向横坡无变化的(例如：起点横坡为2%，终点横坡为2%)：
If S≥超高渐变段起点里程 And S≤超高渐变段终点里程：Then 超高渐变起点坡度→I:0→Z[10]
②起点终点坡度有变化的(例如：起点坡度为-1.5%，终点横坡为1.5%)：
If S≥超高渐变段起点里程 And S≤超高渐变段终点里程：Then 超高渐变起点坡度→I:超高渐变终点坡度→J:超高渐变起点里程→O:超高渐变段长度→A:1→Z[10]
注：超高渐变起终点横坡正负的判断是相对变化的，如果路线左侧高程大于路线右侧高程，那么横坡就为正，反之为负
7、竖曲线高程计算子程序：SUB3
I÷100→I:J÷100→J:J-I→W:RAbs(W)÷2→T:B-T→P:B+T→Q:If W>0:Then F=1:Else F=-1:IfEnd
B-S→U
If S≤P:Then G-IU→H:IfEnd
If S≤B:Then G-IU+F(S-P)2÷2R→H:IfEnd
If S
8、竖曲线超高计算子程序：SUB4
If Z[10]=0:Then Goto 0:IfEnd
If Z[10]=1:Then Goto 1:IfEnd
Lbl0
H-ZI→H
Return
Lbl1
H-Z((J-I)(3((S-O)÷A)^2-2((S-O)÷A)^3)+I)→H
Return