5800简单全线坐标计算程序

5800全线任意坐标计算程序

1. 正算主程序(ZHCX) （不运行）

8→DimZ

1÷P→Z[4 ]：(P-R)÷(2HPR)→D: 180÷π→E

“Z=”?Z:”YJJ=”?A:Abs(S-O)→W

0.26→Z[1 ]: 0.74→B: 0.02→K: 0.82→Z[3 ]: 1-Z[3 ]→F:1-K→Z[2 ]

U+W(Z[1 ]cos(G+QEKW(Z[4 ]+KWD))+Bcos(G+Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+Bcos(G+QEFW (Z[4 ]+FWD))+ Z[1 ]cos(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→X:

V+W(Z[1 ] sin (G+QEKW(Z[4 ]+KWD))+B sin(G+ Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+B sin(G+QEFW (Z[4 ]+FWD))+ Z[1 ] sin(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→Y：

G+QEW(Z[4 ]+WD)→F：X+Zcos(F+A)→X：Y+Zsin(F+A)→Y:If F≧360:Then F-360→F:IfEnd

”X=”:X→X◢

”Y=”:Y→Y◢

If F﹤0:Then F+360→F:IfEnd

”QX FWJ=”:F▼DMS◢

“C=1=>XX: C=2=>XZ”: ”C=”?C: ”QHJU=”?L: If C=1：Then Goto 1:Else Goto 2: IfEnd 可以计算斜交斜做或斜交正做的桥涵坐标

Lbi 1

X+L cos(F)→X：Y+Lsin(F)→Y: Goto 3

Lbi 2

X+L cos(F+A-90)→X：Y+Lsin(F+A-90)→Y: Goto 3

Lbi 3

“QH-X=”: X →X◢

“QH-Y=”: Y →Y◢

Prog “FY”

2 . 参数子程序（直接运行）

M(主线) 一条线路一个名称

“S=”?S

If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd

…

…

If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd

Lbi 1

Prog “ZBJS”

3. 放样程序（FY）（不运行）

“X0=”?M:“Y0=”？N

Pol((X-M, Y-N)

If J﹤0:Then J+360→J:IfEnd

“FWJ=”:J▲DMS◢

“JU-LI=”:I◢

说明：

一、程序功能及原理

1.功能说明：

本程序由一个主程序和两个子程——正算主程序(ZH)、参数子程序（M）、放样程序（FY），可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径、左偏或右偏）及里程边距或坐标，对该线元段范围内任意里程中边桩坐标进行正算。本程序可以在CASIO fx-4800P计算器及CASIO fx-4850P计算器上运行。2．计算原理：

利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据。

利用待求点至线元起点切线作垂线，逐次迭代趋近原理反算里程及边距。

二、使用说明

1、规定

(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1。

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

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

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

(5) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。

2、运行参数子程序（要做哪条线就直接运行哪条线的参数子程序）

输入

S （S=?）为任意里程桩号

Z (Z=?) 为斜距（当所求点位于中线时，Z=0；当位于中线左侧时，Z 取负值；当位于中线右侧时，Z取正值）

A (YJJ=?) 为右夹角（0≦W≦180）

输出

X=********* 计算的X坐标值

Y=********* 计算的Y坐标值

QXFWJ= *°*′*″计算的里程桩号K的切线方位角

输入

C=? 当C=1时，为计算斜交斜做的桥涵坐标；当C=2时，为计算斜交正做的桥涵坐标

QHJU=? 相对于刚计算完坐标的那个点的大桩号（+）或小桩号（-）的距离输出

QH-X=*******.*** 计算的斜交斜做或斜交正做的X坐标值

QH-Y=********* 计算的斜交斜做或斜交正做的Y坐标值输入

X0=? 测站的X值

Y0=? 测站的Y值

输出

FWJ= *°*′*″计算的X 、Y对于测站的方位角

JU-LI=******* 计算的X 、Y对于测站的距离