CASIO_fx-5800P实用工程测量辛甫森公式放样程序

测量放样，适用于5800p计算器

一、QXFY 辛甫森公式放样程序

1. “X0”? U：“Y0”?V

2. “XA”? A：“YA”? B：“CA”? C：“1÷RA”?D：“1÷RB”?E：

“KA”?F：“KB”? G

3. Lb1 1：“KI”?H：“JJ”?L：“Y+Z－”?R

4. If H>G Or H

6. (E-D) ÷Abs(G-F) P: Abs(H –F)→Q：P×Q→Z：D+Z→T

7. C+（Z+2D）Q×90÷π→W :“QXFWJ”：

8. C+(Z÷4+2D)Q×22.5÷π→M：C+（3Z÷4+2D）Q×67.5÷π→

N：C+(Z÷2+2D)Q×45÷π→K

9. A+Q（cosC+4(cosM+cosN)+2cosK+cosW）÷12+Rcos(W+L) →X：

“X＝”：X

10. B+Q(sin(C)+4(sin(M)+sin(N))+2sin(K)+sin(W))÷12+Rsin (W+L)

→Y：“Y＝”：Y

11. Pol((X-U)，（Y-V）)

If J≤0 Then J+360→J: IfEnd : “FWJ”:J

12. “JU”: I

13. Goto 1

程序显示说明：

须输入参数：X0:置镜点X (对应变量U) Y0:置镜点Y (对应变量V)

XA:曲线起算点X (对应变量A)

YA:曲线起算点Y (对应变量B)

CA：曲线起算点切线方位角(对应变量C) 1÷RA：1÷半径，即起算点曲率，右偏为正左偏为负(对应变量D)

1÷RB：1÷半径，即曲终点曲率，右偏为正左偏为负(对应变量E)

KA:起算点里程(对应变量F)

KB:曲终点里程(对应变量G)

KI:待求点里程(对应变量H)

JJ:夹角(与前进方向切线向右的夹角)(对应变量L)

D“Y + Z－：偏距，右正左负(对应变量R)

计算结果: W: 待求点切线方位角(对应变量W) X:计算点x (对应变量X)

Y:计算点y (对应变量Y)

FWJ:放样方位角(对应变量J)

S:放样距离(对应变量I)

X0

Y0

KB

本程序依据复化辛甫森公式计算原理改进而成，特点是把曲线按曲率变化点分成若干计算单元单独计算，不论直线（曲率为0）、圆曲线（曲率为1/R）、卵形曲线，只要按曲率变化点分解弄清变化点曲率半径，右偏左偏曲线由曲率的正负号分别，既可求得该计算单元内任意里程中桩也可求斜交或正交的边桩。考虑到程序中调用变量较多与其它程序变量冲突，造成不方便来回交叉使用，便把反算程序嵌套其中，作为一个完整独立放样和坐标计算程序，不含子程序

当输入待求里程超出该曲线单元范围时程序会拒绝运算，并提示重新输入里程，避免计算错误结果；如要求另一曲线单元内里程须重新进入程序输入相应曲线元参数计算

由于应用中待求里程不一定按从小到大顺序便去掉了把计算里程替代为起算里程的语句，造成待求里程与起算里程距离较远时精度会略有影响，但对于一般工程中较大半径曲线，影响非常小仅在毫米上可以忽略不予考虑。

二、坐标反算ZBFS

1. “X0”？U：V“Y0”？V

2. Lbl 0：“X”？X：“Y”？Y：Pol((X－U),(Y－V))

3. If J≤0：Then J+360→J：IfEnd ：“FWJ”：J

4. “JU”：I

5. Goto 0

程序显示说明：

须输入参数：X0？: 起算点X

Y0？: 起算点Y

X ？: 计算点X

Y？: 计算点Y

计算结果: FWJ: 反算方位角

JU: 反算距离

X ？Y？: 重新输入计算点X Y

三、坐标正算ZBZS

1. “X0”?X：“Y0”? Y：“S”? S：“T”? T

2. “X1”:X+Scos(T)→M

“Y1”:Y+Ssin (T)→N

程序显示说明：

须输入参数：X0？: 起算点X

Y0？: 起算点Y

S？: 正算距离

T？: 方位角

计算结果: X1: 计算点X

Y1: 计算点Y

四、竖曲线SQX

1. “R＝”? A：“T＝”? B：“I1＝”? C：“I2＝”? D：“K0

＝”?E：“H0＝”? F

2. If C>D : Then －1→T : Else 1→T :IfEnd

3. Lbl 0：“K＝”?K：Abs(E－K）→S：

4. If S>B :Then Goto 0 : IfEnd

5. If E≥K : Then Goto 1: Else Goto 2 :IfEnd

6. Lbl 1：F－SC+T（B－S）2÷2÷A→G:“G”

: G

7. Goto 0

8. Lbl 2：F+SD+T（B－S）2÷2÷A→G:“G”

9. Goto 0

程序显示说明：

须输入参数：R＝？: 竖曲线半径

T＝？: 切线长

I1＝？: 变坡点前坡度值，上坡为正，下坡为负（如-2.5%则输-0.025）

I2＝？: 变坡点后坡度值，上坡为正，下坡为负

K0＝？：变坡点里程

H0＝？：变坡点坡线高程

K＝？：待求点里程

计算结果: G: 待求点曲线高程

K＝？：重新输入待求点里程

注：计算里程仅在竖曲线范围内

凹型竖曲线凸型竖曲线