《辛卜生曲线程式4850.doc》由会员分享，可在线阅读，更多相关《辛卜生曲线程式4850.doc（3页珍藏版）》请在金锄头文库上搜索。
1、CASIO fx-4850P计算器程序一、文件名：LXQX（主程序）DZG：A“XA”：B “YA”: C“FW”：D“1/RA”：Z “1/RB”: F“DKA”：G“DKB”LbI 1：H：H“DKI”：P=（Z-D）/AbS（G-F）：K=AbS（H-F）：S=PK：Q=C+（S+2D）K*90/:V=C+(S/4+2D)K*45/2: N=C+(3S/4+2D)K*135/2: J=C+(S/2+2D)K*45/: R=C+(S/8+2D)K*45/4:L=C+(3S/8+2D)K*135/4: T=C+(5S/8+2D)K*225/4: W=C+(7S/8+2D) K*315/4:X
2、=A+K/24*(cosC+4(cosR+cosL+cosT+cosW)+2(cosV+cosJ+cosN)+cosQ): Y=B+K/24*(sinC+4(sinR+sinL+ sinT+sinW)+2(sinV+sinJ+sinN)+sinQ)T=0: T: T0=Prog “HXY-JI” “X=”: X “Y=”: Y（T默认为零，输入不为零的数则计算线路中心外移点的坐标）M=X: U=Y: Prog “XY-SN”: Goto 1二、文件名:XY-SN(坐标反算)E “XO”:O “YO”: M “XI”: U “YI”: pol (M-E, U-O):J0=J=J+360 “QA=
3、”: J “QS=”: I三、文件名: HXY-JI(线路中线外移点坐标计算)M=X: U=Y: Lbl 3: V, D, L: V “QH”: D “ZY”: L “QA”:Pol(V, D):J0=J=J+360J=J+Q+L: J360=J=J-360X=M+IcosJ:Y=U+IsinJ: “X=”: X “Y=”: Ypol(X-E“XO”, Y-O“YO”):J0=J=J+360 “QA=”: J “QS=”: IGoto 3说明:1、该程序对任何曲线都适合,可以说是一万能公式，但需对线路进行曲线元的划分，然后分段计算； 2、对该程序要求输入数据的说明: 输入:XA? YA? :要
4、求输入曲线元计算起点X,Y坐标； FW? :要求输入曲线元计算起点的切线方位角(沿计算方向:) 1/RA? :要求输入曲线元起点处的曲率（半径的倒数，左偏为负，右偏为正）；1/RB? :要求输入曲线元终点处的曲率（半径的倒数，左偏为负，右偏为正）； DKA? :要求输入曲线计算段起点里程(以米为单位); DKB? :要求输入曲线计算段终点里程(以米为单位);DKI? :要求输入曲线计算段内待求点里程； QA？ ：结构物轴线与线路前进方向（曲线为该点切线方向）夹角，详见图示；QH？ ：点偏离前后方向的距离（前正后负，详见图示）；ZY？ ：点偏离左右方向的距离（右正左负，详见图示）；输出：Q= ：输出求得的曲线（某一里程的）切线方位角；X=，Y= ：输出求得的某点坐标；J ：为置镜点到待测点方位角；I ：为置镜点到待测点距离；3、利用此程式计算也可将曲线倒过来算，但一般无此必要，经实例计算，半径53米的圆，圆心角217度, 坐标差值最大18mm，因此，对于圆心角大于或等于180度的圆，可分两段（即：分两个曲线元）连续计算即可达到要求，差值可控制在2mm左右。4、程式三的输入图示：
