《辛卜生曲线程式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、程式三的输入图示: