[建筑]辛卜生曲线程式4850.docVIP

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=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}: T≠0=＞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): J＜0=＞J=J+360⊿ “QA=”: J◢ “QS=”: I◢ 三、文件名: HXY-JI(线路中线外移点坐标计算) M=X: U=Y: Lbl 3: {V, D, L}: V “QH”: D “ZY”: L “QA”: Pol(V, D):J＜0=＞J=J+360⊿J=J+Q+L: J≥360=＞J=J-360⊿X=M+IcosJ:Y=U+IsinJ: “X=”: X◢ “Y=”: Y◢ pol(X-E“XO”, Y-O“YO”): J＜0=＞J=J+360⊿ “QA=”: J◢ “QS=”: I◢Goto 3 说明:1、该程序对任何曲线都适合,可以说是一万能公式，但需对线路进行曲线元的划分，然后分段计算； 2、对该程序要求输入数据的说明: 输入:XA? YA? :要求输入曲线元计算起点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、程式三的输入图示：