卡西欧5800缓和曲线计算程序

卡西欧5800直线、缓和曲线程序1. 主程序：TYQXJS″1.SZ=>XY″：″2.XY=>SZ″：″N″？N：″XA″？U：″YA″？V：″DKA″？O：″CA″？G：″LS″？H：″RA″？P:″RB″？R：″Q″？Q： 1÷P→C:（P-R）÷（2HPR）→D： 180÷π→E：If N=1: Then Goto1:Else Goto 2:ifEndLBI 1: ″DKI″?S: ″D″?Z:Abs(S-O)→W:Prog″SUB1″: ″XS″:X⊿″YS″:Y⊿″FS″:F-90⊿Goto 1LbI 2：″X″?X: ″Y″?Y：X→I:Y→J：Prog″SUB2″：O+W→S:″S″:S⊿″Z″:Z⊿Goto 22. 子程序：SUB10.1739274226→A:0.3260725774→B:0.0694318442→K:0.3300094782→L:1-L→F:1-K→M:U+(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+ FWD))+Acos(G+QEMW(C+MWD))→X:Y+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD)))→Y:G+QEW(C+WD)+90→F:X+Zcos(F)→X:Y+Zsin(F)→Y3. 子程序：SUB2G-90→T： Abs((Y-V)cosT-(X-U)sin(T))→W:0→Z： Lbl 0：Prog″SUB1″:T +QEW(C+WD)→L :(J-Y)cos(L)-(I-X)sin(L)→Z:IfAbs(Z)<1*10-6: Then Goto1:ElseW+Z→W:Goto0:IfEndLbl 1: 0→Z: Prog″SUB1″: (J-Y) ÷sin(F)→Z4. 输入与显示说明N ：选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示由坐标反算里程和边距。

圆曲线缓和曲线线路中边桩放样卡西欧5800坐标计算程序程序名：FYXY"JD"?Q:"JDX"?W:"JDY"?K:R:"QDFWJ"?F:"A0"?N:"LS"?V:"-L,+R"?G:"DKQD"?C:"DK ZD" ?D:”ZJX"?A:”ZJY”?B8→DimzV÷2-V^(3)÷(240R2)→MV2÷(24R)－V^(4)÷(2688R^(3))→PpRN÷180+V→LM+(R+P)tan(N÷2)→TQ-T→Z[5]Z[5]+V→Z[6]Z[5]+L→Z[8]Z[8]-V→Z[7]Rec(T，F+180)W+I→Z[1]:K+J→Z[2]Rec(T，F+GN)W+I→Z[3]:K+J→Z[4]Lbl 0"DKI"?H:"BIAN JU"?S:"JIA JIAO"?EIf H<C:Then Goto 0:IfEndIf H>D:Then Goto 0:IfEndIf H<Z[5]:Then Goto 1:IfEnd:If H<Z[6]:Then Goto 2:IfEnd:If H<Z[7]:Then Goto 3:IfEnd:If H<Z[8]:Then Goto 4:Else Goto 5:IfEndLbl 1W+I→X: K+J→YRec(S，F+E)"X=":X+I→X◢"Y=":Y+J→Y◢Goto 7Lbl 2H-Z[5]→Z90Z2÷(pRV)→OZ-Z^(5)÷(40R2V2)+Z^(9)÷(3456R^(4)V^(4))→XZ^(3)÷(6RV)-Z^(7)÷(336R^(3)V^(3))+Z^(11)÷(42240R^(5)V^(5))→Y Rec(√(X2+Y2)，F+Gtan-1(Y÷X))Z[1]+I→X:Z[2]+J→YRec(S，F+GO+E)"X=":X+I→X◢"Y=":Y+J→Y◢Goto 7Lbl 3H-(Q-T)-V÷2→Z180Z÷(Rp)→ORsin (O)+M→XR(1-cos(O))+P→ZRec(X，F)Z[1]+I→X:Z[2]+J→YX+I→X:Y+J→YRec(S，F+OG＋E)"X=":X+I→X◢"Y=":Y+J→Y◢Goto 7Lbl 4Z[8]-H→Z90Z2÷(VRp)→OZ-Z^(5)÷(40R2V2)+Z^(9)÷(3456R^(4)V^(4))→XZ^(3)÷(6RV)-Z^(7)÷(336R^(3)V^(3))+Z^(11)÷(42240R^(5)V^(5))→Y Rec(√(X2+Y2)，F＋GN+180-Gtan-1(Y÷X))Z[3]+I→X:Z[4]+J→YRec(S，F+GN+180-OG-(180-E))"X=":X+I→X◢"Y=":Y+J→Y◢Goto 7Lbl 5Rec(H-Z[8]+T，F+GN)W+I→X:K+J→YRec(S，F+GN+E)"X=":X+I→X◢"Y=":Y+J→Y◢Goto 7Lbl 7Pol(X-A，Y-B)If J<0:Then J+360→J:IfEnd"FWJ=":J…DMS◢"CJ=":I◢Goto 0注：FYXY程序单代号Q"J D" 交点里程桩号W"JDX" 交点X坐标值K"JDY" 交点Y坐标值R 圆曲线半径F"QDFWJ" 第一直线段方位角N"A 0" 交点转向角V"LS" 缓和曲线长度G"-L,+R" 如果为左转角输入“－1”，如果为右转角输入“1”C"DKQD" 所求段起点里程桩号D"DKEZD" 所求段终点里程桩号ZJX 置镜点Ｘ坐标ZJY 置镜点Y坐标H"DKI" 所求点的里程桩号S"BIAN JU" 边桩至中桩的距离，求中桩时输入0 距离不为负E"JIA JIAO" 左边桩输入“－90”，右边桩输入“90”X Y 所求点的坐标值C J 置镜点至待求点的距离FWJ 置镜点至待求点的方位角。

CASIO FX-4800P圆曲线及缓和曲线上任意点切线方位角计算程序程序内容I=0：J=0：Norm：deg：Lbl 0:｛ABCDGXYR｝：A“KO”D“KOX”G“KOY”B“KN”X“KNX”Y“KNY”C“LS”R“R”：T“XX=Y：H”：V“FX=Y：Z”：L=Abs（B-A）：T=1=》Z=（L/2R）×（180/π）：≠》Z=19.10×L2/（R×LS）：POL（X-D，Y-G）：J＜0=》J+360ΔFix6：V=1=》N=J+Z：N“FWJ”=IntN+0.01Int (60 Frac N)+0.006 Frac (60 Frac N)▲：≠》N=J-Z：N“FWJ”=IntN+0.01Int (60 Frac N)+0.006 Frac (60 Frac N)▲：Goto0说明：KO：线元起点桩号KOX：线元起点X坐标KOY：线元起点Y坐标KN：所求点桩号KNX：所求点X坐标KNY：所求点Y坐标LS：缓和曲线长，如为圆曲线输为0R：线元半径“XX=Y：H”：线形，圆曲线输入1，缓和线输入2“FX=Y：Z”：方向，右转为1，左转为2FWJ：所求点切线方位角，为度分秒形式显示。

例如100.1122=100°11′22″CASIO FX-4800P运行通过，希望对大家有所帮助，如有什么问题请高手指正。

方位角计算通用公式曲线上任一点切线方位角计算通用公式：P=1/R1+(1/R2-1/R1)*ABS(DKP-DKA)/ABS(DKB-DKA)....式中P为所切点曲率；1/R1为起点曲率、1/R2为终点曲率、DKP为计算点桩号、DKA为起点桩号、DKB 为终点桩号。

J=E+28.64789*(1/R1+P)*ABS(DKP-DKA)......式中J为计算点切线方位角；E 为起点切线方位角。

由于公式为通用，故对直线、圆曲线、缓和曲线上任一点的切线方位角计算都有效，在使用公式时规定：1/R1、1/R2 均为曲率即设计曲线半径的倒数，直线、完整缓和曲线起点设计半径均输入“1E50”，曲率故为1/1E50；圆曲线和缓和曲线接圆曲线端为设计圆曲线半径值，取其倒数。

[正确] CASIO fx-5800P QXZB曲线坐标计算CASIO fx-5800P QXZB程序可计算不等缓和曲线、圆曲线上的任意中、边桩坐标：该程序适用于计算器 CASIO fx-5800P，可计算与线路中心成任意夹角的缓和曲线、圆曲线中、边桩坐标及待测点方位角和距离。

1、DK(JD)？输入交点桩号2、X(JD)？输入交点坐标X3、Y(JD)？输入交点坐标Y4、T1？输入第一切线长(如果只有一条切线两者都输入一致)5、T2？输入第二切线长(如果只有一条切线两者都输入一致)6、FWJ？输入直线方位角(ZH→JD)7、A？输入转角：左转为负，右转为正8、R？输入圆曲线半径9、LS1？输入第一缓和曲线长(如果只有一条缓和曲线两者都输入一致)10、LY？输入圆曲线长(L-LS1-LS2)11、LS2？输入第二缓和曲线长(如果只有一条缓和曲线两者都输入一致)12、X(ZJD)？输入置镜点坐标X13、Y(ZJD)？输入置镜点坐标Y14、JSDK？输入前视点里程15、PL？输入偏距16、PA？输入偏角程序名称：QXZB＂DK(JD)＂?R：＂X(JD)＂?S：＂Y(JD)＂?T：＂T1＂?U：＂T2＂?V：＂FWJ＂?F：＂A＂?W：＂R＂?H：＂LS1＂?I：＂LY＂?J：＂LS2＂?K：＂X(ZJD)＂?D：＂Y(ZJD)＂?EF+W→GIf W>0： Then -1→L ：IfEnd(式中为数字0…括号内不输入)If W<0： Then 1→L ：IfEnd (式中为数字0…括号内不输入)If W=0： Then 0→L ：IfEnd (式中为数字0…括号内不输入)R-U→AS-Ucos( F) →BT-Usin( F) →CS+Vcos( G) →Z[1]T+Vsin( G) →Z[2]180÷π→O(式中为字母0…括号内不输入)Lb1 A＂JSDK＂?MM-A→PIf M<A： Then Goto A ：IfEndIf M>(A+I+J+K)： Then Goto A ：IfEnd＂PL＂?NIf M>(A+I+J) ：Then Goto D：IfEndIf M>(A+I) ：Then Goto C：IfEndLb1 BB+(P-(P^(5)÷40÷H^2÷I^2))cos( F)+(P^(3)÷6÷H÷I-P^(7)÷336÷H^(3)÷I^(3))Lsin( F)→XC+(P-(P^(5)÷40÷H^2÷I^2))sin( F)-(P^(3)÷6÷H÷I-P^(7)÷336÷H^(3)÷I^(3))Lcos( F)→Y“X=”：X+Ncos(F-(P^2÷2÷H÷I)OL+Z) →X ◢(式中OL为字母O…括号内不输入)“Y=”：Y+Nsin(F-(P^2÷2÷H÷I)OL+Z) →Y ◢(式中OL为字母O…括号内不输入)Prog＂ZBFS＂(此处为调用直线程序…括号内不输入)Goto ALb1 CB+(Hsin(((P-I)÷H+I÷2÷H)O)+(I÷2-I^(3)÷240÷H^2))cos(F)-(H(1-cos(((P-I)÷H+I÷2÷H)O))+I^2÷24÷H)(-1)Lsin( F)→X(中间O为字母O)C+(Hsin(((P-I)÷H+I÷2÷H)O)+(I÷2-I^(3)÷240÷H^2))sin(F)+(H(1-cos(((P-I)÷H+I÷2÷H)O))+I^2÷24÷H)(-1)Lcos( F)→Y(式中为字母0…括号内不输入)“X=”：X+Ncos(F-((P-I)÷H+I÷2÷H)OL+Z) →X ◢(式中OL为字母O…括号内不输入)“Y=”：Y+Nsin(F-((P-I)÷H+I÷2÷H)OL+Z) →Y ◢(式中OL为字母O…括号内不输入)Prog ＂ZBFS＂(此处为调用直线程序…括号内不输入)Goto ALb1 D(A+I+J+K-M) →QZ[1]-(Q-Q^(5)÷40÷H^2÷K^2)cos( G)+(Q^(3)÷6÷H÷K-Q^(7)÷336÷H^(3)÷K^(3))Lsin( G)→XZ[2]-(Q-Q^(5)÷40÷H^2÷K^2)sin( G)+(Q^(3)÷6÷H÷K-Q^(7)÷336÷H^(3)÷K^(3))Lcos( G)→Y “X=”：X+Ncos(F-(I÷2÷H+J÷H+K÷2÷H-Q^2÷2÷H÷K)OL+Z) →X ◢(OL为字母O…括号内不输入) “Y=”：Y+Nsin(F-(I÷2÷H+J÷H+K÷2÷H-Q^2÷2÷H÷K)OL+Z) →Y ◢(OL为字母O…括号内不输入) Prog＂ZBFS＂(调用直线程序…括号内不输入)Goto A--------------------------------------------------------------------------------说明：本程序需要扩展变量，方法：20→DimZ附件：本程序附曲线要素坐标计算核对表。