4850计算器万能程序(公路专用)

4850计算器万能程序(公路专用)

主程序 TYQXJS

Defm10：Lbl 0: “1.LC=>XY”: “2.XY=>LC”: {V}:V=1=>GOTO 1: ≠> GOTO 2 △

Lbl 1: {DZG}:D:Z”JL=”:G”J=”:Z=0=>G=0：△Prog “QXYS”;Prog “ZSZB” : ”X=”：X ▲ ”Y=”：Y▲"FWJ="：J→DMS▲

{HU}：H“XC”：U“YC”：I=0：J=0：Pol（X-H，Y-U）：J<0＝＞J=J+360：Δ“FW=”：J→DMS◢

“S=”：S◢

Prog “SQX”：

“H0=”：Z[8] ◢

GOTO 0

Lbl 2: D=***:Z=0:G=0:{MI}:M”XO”:I”YO”: Prog “FSZB”:”D=”：D ▲”JL=”：Z▲

Prog “SQX”：“H0=” ：Z[8]◢

GOTO 0

D=有效里程最好是起点这样计算需要3秒

想计算的更快就在{MI} 添加一个D 就是{DMI}这时候的D＝计算点的大概里程越近越快不要忘记把前面的D=***:去掉

2.正算坐标ZSZB

N=5:P＝（1\E-1\R）÷Abs（K-F）：Q＝Abs（D-F）÷N：S=90Q÷π： J＝

C+(NPQ+2\R)NS：L=1←┘

X＝A+Q÷6×(Cos C+Cos J +4∑（Cos

(C+((L+0.5)PQ+2\R)×(L+.5)S),L,0,(N-1)）+2∑（Cos

(C+((LPQ+2\R)LS,L,1,(N-1))）+ZCos（J+G)：←┘

Y＝B+Q÷6×(Sin C+Sin J +4∑（Sin

(C+((L+0.5)PQ+2\R)×(L+.5)S),L,0,(N-1)）+2∑（Sin

(C+((LPQ+2\R)LS,L,1,(N-1))）+Z Sin（J+G)：←┘

3.反算坐标：FSZB

Lbl 0:Prog “QXYS” :Prog “ZSZB”:O=J-90:

P= (I-Y)cosO-(M-X) sinO :D=D+P:

AbsP≥0.001=> GOTO 0 :≠> GOTO 1 △

Lbl 1: Z= (I-Y)cosJ-(M-X) sinJ

4. 《SQX》高程计算子程序

Prog“SB”：W=Z[1]－Z[2]:W＞0=>Z[6]＝－1：≠>Z[6]＝1△Z[3]：T“T”＝Abs(Z[3]W/2) ：E“E”=T＾2÷2Z[3]：C=D-Z[4]：D≦Z[4]=>Z[7]=Z[1]：

≠>Z[7]=Z[2]△Z[5]：AbsC≦T=>Z[8]“H0”=Z[5]+CZ[7]+Z[6](T- Abs C)＾

2÷2Z[3]：≠>Z[8]“H0”=Z[5]+CZ[7]：

5、曲线元要素数据库：QXYS

D≥F=>D＜K=> F=**:A=**:B=**:C=**:R=**:K=**E=**⊿⊿←┘

D≥F=>D＜K=> F=**:A=**:B=**:C=**:R=**:K=**E=**⊿⊿←┘

D≥F=>D＜K=> F=**:A=**:B=**:C=**:R=**:K=**E=**⊿⊿←┘

D≥F=>D＜K=> F=**:A=**:B=**:C=**:R=**:K=**E=**⊿⊿←┘

D≥F=>D＜K=> F=**:A=**:B=**:C=**:R=**:K=**E=**⊿⊿←┘

……………………………

D≥F=>D＜K=>

F=**:A=**:B=**:C=**:R=**:K=**E=**⊿⊿←┘

(注：如有多个曲线元要素继续添加入数据库QXYS中)

6、《SB》竖曲线数据库

D≥竖曲线起点里程=>Z[1]=？：Z[2]=？：Z[3]=？：Z[4]=？：Z[5]=？△

D≥？=>Z[1]=？：Z[2]=？：Z[3]=？：…………△

附注：

1、D所求的桩号：Z[1]=坡度1:Z[2]=坡度2：Z[3]=半径：Z[4]=交点桩号：Z[5]=变坡点控制高程：

2、注意Z[1]、Z[2]的正负。

3、本程序可以计算上百公里的路线中桩高程。

变量说明：

F……..起点里程A…..起点X坐标B….起点Y坐标C….起点方位角R….….起点半径K…..终点里程 E.…终点半径D….计算点里程G…输入左右（左-，零，右+）角度Z….中桩到边桩的距离

注意：1:计算中桩坐标 J,JL 等于零。

2:如果建立数据库（主线坐标计算ZXZB），子程序中的

3：这里的方位角是指“主点作为起点” 的方位角

例如：（ZY , YZ , ZH , HY , YZ）每个点都是起点所以特别注意

1、规定

(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，半径＝负半径；当线元往右偏时,半径=半径；当线元为直线时，半径=以10的45次代替。

(2) 当所求点位于中线时，Z,J=0；当位于中线左侧时，J取负值；当位于中线右侧时，J取正值。

(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45

次代替。

(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。

(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。

(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明

V=1的时候

JL ？正算时所求点距中线的边距(在中线上取零)

J？边桩时左右边桩连线与线路前进中线的交角线元左右偏标志

(左偏J=-1*角度，右偏J=角度，直线段J=0) 角度

X=×××正算时，计算得出的所求点的X坐标

Y=×××正算时，计算得出的所求点的Y坐标

FWJ=*** 计算得出所求点的切线方位角°′″

FW=*** 计算得出置镜点到测点方位角

V=2的时候

XO=×××所求点的X坐标

YO=×××所求点的Y坐标

D=**** 求点所对应的里程

JL=**** 求点到所对应的里程的距离：负就是左边相反就是右边

本人有个人意见:

一般公路测量的时候可以将主程序中的

Prog “SQX”：

去掉可以提高计算器的运算速度一秒左右