下载文档原格式
一、QXFY 辛甫森公式放样程序1. “X0”? U：“Y0”?V2. “XA”? A：“YA”? B：“CA”? C：“1÷RA”?D：“1÷RB”?E：“KA”?F：“KB”? G3. Lb1 1：“KI”?H：“JJ”?L：“Y+Z－”?R4. If H>G Or H<F : Then Goto 1: IfEnd6. (E-D) ÷Abs(G-F) P: Abs(H –F)→Q：P×Q→Z：D+Z→T7. 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÷π→K9. A+Q（cosC+4(cosM+cosN)+2cosK+cosW）÷12+Rcos(W+L) →X：“X＝”：X10. B+Q(sin(C)+4(sin(M)+sin(N))+2sin(K)+sin(W))÷12+Rsin (W+L)→Y：“Y＝”：Y11. Pol((X-U)，（Y-V）)If J≤0 Then J+360→J: IfEnd : “FWJ”:J12. “JU”: I13. 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)X0Y0KB本程序依据复化辛甫森公式计算原理改进而成，特点是把曲线按曲率变化点分成若干计算单元单独计算，不论直线（曲率为0）、圆曲线（曲率为1/R）、卵形曲线，只要按曲率变化点分解弄清变化点曲率半径，右偏左偏曲线由曲率的正负号分别，既可求得该计算单元内任意里程中桩也可求斜交或正交的边桩。
Fx-5800P万能程序核心计算程序摘自工测员网“yshf”,数学模型：辛普森公式一.程序清单:1. 正算子程序(SUB1)0.1739274226→Z[1]：0.3260725774→Z[2]：0.0694318442→Z[3]：0.3300094782→Z[4]：1- Z[4]→F：1- Z[3]→Z[5]：（将原来变量A、B、K、L、M,用Z[1] 、Z[2] 、Z[3] 、Z[4] 、Z[5]扩充变量来代替，以便嵌套其他程序）U+W(Z[1]cos(G+57.2958Q Z[3]W(1÷P+Z[3]WD))+Z[2]cos(G+57.2958QZ[4]W(1÷P+Z[4]WD))+Z[2]cos(G+57.2958QFW(1÷P+FWD))+Z[1]cos(G+57.2958Q Z[5]W(1÷P+Z[5]WD)))→X：V+W(Z[1] sin(G+57.2958Q Z[3]W(1÷P+Z[3]WD))+Z[2]sin(G+57.2958QZ[4]W(1÷P+Z[4]WD))+Z[2]sin(G+57.2958QFW(1÷P+FWD))+Z[1] sin(G+57.2958Q Z[5]W(1÷P+Z[5]WD)))→Y：G+57.2958QW(1÷P+WD)+M→F：（M为线路前进方向右交角）X+ZcosF→X：Y+ZsinF→Y2. 反算子程序(SUB2)G-90→T：Abs((Y-V)cosT-(X-U)sinT) →W：0→Z：Lbl 0：Prog "SUB1"：T+57.2958QW(1÷P+WD) →L：(J-Y)cosL-(I-X)sinL→Z：If AbsZ<10^(-6): Then Goto1：Else W+Z→W：Goto 0: IfEnd←┘Lbl 1：0→Z：Prog "SUB1"：(J-Y)÷sinF→Z←┘3.主程序（ZCX）(数据库程序)10→DimZ(扩充10个变量)xxxxxx.xxx→Cyyyyyy.yyy→E(坐标放样中测站点坐标)Lb1 4：”N”﹖N：”S”﹖SIf S≤下一线元起点里程:Then 本线元起点里程→O：本线元起点X坐标→U=：本线元起点Y坐标→V：本线元起点方位角→G：本线元长度→H：起点曲率半径→P：终点曲率半径→R：0（直线）或1（线路右转）、-1（线路左转）→Q： Goto0：IfEnd←┘(第一段线元数据要素)If S≤下一线元起点里程:Then 本线元起点里程→O：本线元起点X坐标→U=：本线元起点Y坐标→V：本线元起点方位角→G：本线元长度→H：起点曲率半径→P：终点曲率半径→R：0（直线）或1（线路右转）、-1（线路左转）→Q： Goto0：IfEnd←┘(第二段线元数据要素)If。
卡西欧fx5800P隧道测量程序…一、缓和曲线上的点坐标计算已知：①缓和曲线上任一点离ZH点的长度：l②圆曲线的半径：R③缓和曲线的长度：l0④转向角系数：K(1或－1)⑤过ZH点的切线方位角：α⑥点ZH的坐标：xZ，yZ计算过程：说明：当曲线为左转向时，K=1，为右转向时，K=-1，公式中n的取值如下：当计算第二缓和曲线上的点坐标时，则：l为到点HZ的长度α为过点HZ的切线方位角再加上180°K值与计算第一缓和曲线时相反xZ，yZ为点HZ的坐标切线角计算公式：二、圆曲线上的点坐标计算已知：①圆曲线上任一点离ZH点的长度：l②圆曲线的半径：R③缓和曲线的长度：l0④转向角系数：K(1或－1)⑤过ZH点的切线方位角：α⑥点ZH的坐标：xZ，yZ计算过程：说明：当曲线为左转向时，K=1，为右转向时，K=-1，公式中n的取值如下：当只知道HZ点的坐标时，则：l为到点HZ的长度α为过点HZ的切线方位角再加上180°K值与知道ZH点坐标时相反xZ，yZ为点HZ的坐标三、曲线要素计算公式公式中各符号说明：l——任意点到起点的曲线长度（或缓曲上任意点到缓曲起点的长度）l1——第一缓和曲线长度l2——第二缓和曲线长度l0——对应的缓和曲线长度R——圆曲线半径R1——曲线起点处的半径R2——曲线终点处的半径P1——曲线起点处的曲率P2——曲线终点处的曲率α——曲线转角值四、竖曲线上高程计算已知：①第一坡度：i1(上坡为“＋”，下坡为“－”)②第二坡度：i2(上坡为“＋”，下坡为“－”)③变坡点桩号：SZ④变坡点高程：HZ⑤竖曲线的切线长度：T⑥待求点桩号：S计算过程：五、超高缓和过渡段的横坡计算已知：如图，第一横坡：i1第二横坡：i2过渡段长度：L待求处离第二横坡点（过渡段终点）的距离：x求：待求处的横坡：i解：d=x/Li=(i2-i1)(1-3d2+2d3)+i1六、匝道坐标计算已知：①待求点桩号：K②曲线起点桩号：K0③曲线终点桩号：K1④曲线起点坐标：x0，y0⑤曲线起点切线方位角：α0⑥曲线起点处曲率：P0(左转为“－”，右转为“＋”)⑦曲线终点处曲率：P1(左转为“－”，右转为“＋”)求：①线路匝道上点的坐标：x,y②待求点的切线方位角：αT计算过程：注：sgn(x)函数是取符号函数，当x<0时sgn(x)=-1，当x>0时sgn(x)=1，当x=0时sgn(x)=0。
比积分公式快速CASIO fx-5800P放样程序（比积分公式快速）正算主程序[ZD-KB]“1,2,3,4,5”?Z:?X:?YLbl 1:Fix 4:?K:?B:?A:0→P:0→RZ=1＝＞Prog“YS-1”Z=2＝＞Prog“YS-2”Z=3＝＞Prog“YS-3”Z=4＝＞Prog“YS-4”Z=5＝＞Prog“YS-5”Prog“ZDYS”:M-X→I:N-Y→J:Prog“C”: Cls:Locate 1,1,“J=”:Locate 4,1, J°:Locate 1,2,“I=”:Locate 10,2,I:Locate 1,3,“M=”:Locate 7,3,M:Locate 1,4,“N=”:Locate 7,4,N◢Cls:Goto 1反算主程序[ZD-MN]“1,2,3,4,5”?Z:?KLbl 0:Fix 4:?M:?N:M→D:N→ELbl 1:90→A:0→B:0→P:0→RZ=1＝＞Prog“YS-1”Z=2＝＞Prog“YS-2”Z=3＝＞Prog“YS-3”Z=4＝＞Prog“YS-4”Z=5＝＞Prog“YS-5”Prog“ZDYS”:Pol(D-M+×10-9,E-N+×10-9:Isin(L-J→A:Icos(L-J→BIf Abs(A)＞0.001:Then K+A→K:Goto 1:Else Cls:Locate 1,1,“K=”:Locate 9,2, K:Loc ate 1,3,“B=”:Locate 10,4,B◢IfEnd:Goto 0子程序[C]I+×10-9→I:J+×10-9→J:Pol(I,J:J＜0＝＞J+360→J数据库子程序[YS-1]If K＜线元终点里程:Then线元起点切线方位角→C:线元起点X坐标→U:线元起点Y 坐标→V:线元起点里程→O:线元长度→F:±线元起点曲率半径(左转为负，右转为正，直线为零)→P:±线元止点曲率半径(左转为负，右转为正，直线为零)→R:Return:IfEn dIf K＜下一条线元终点里程:Then线元起点切线方位角→C:线元起点X坐标→U:线元起点Y坐标→V:线元起点里程→O:线元长度→F:±线元起点曲率半径(左转为负，右转为正，直线为零)→P:±线元止点曲率半径(左转为负，右转为正，直线为零)→R:Re turn:IfEnd子程序[ZDYS]P=0＝＞×1045→P:R=0＝＞×1045→R:If PR＞0:Then 1→Q:Else -1→Q:IfEnd:PR≥×1090＝＞0→Q:R＜0＝＞-1→Q:Abs(P→P:Abs(R→R:1÷P→H:180÷π→J:(P-R)÷(2FPR→I:Abs(K-O→W:0.1739274226→N:0.3260725774→M:0.0694318442→S:0.3300094782→L:1-L→G:C+QJSW(H+SWI→T:C+QJLW(H+LWI→L:C+QJGW(H+GWI→G:C+QJ W(1-S)(H+WI(1-S→F:A+C+QJW(H+WI→CU+W(Ncos(T)+Mcos(L)+Mcos(G)+Ncos(F→I:V+W(Nsin(T)+Msin(L)+Msin(G)+Nsin(F →J:C→L:I+Bcos(L→M:J+Bsin(L→N正算主程序[KB-MN]输入：1,2,3,4,5---选择待正算的线路X---测站X坐标Y---测站Y坐标 K---桩号B---距中（左负右正中零）结果：J=测站到放样点的方位角I=测站到放样点的平距M=放样点X坐标N=放样点Y坐标若果求斜交时边桩坐标，主程序[KB-MN]中90→A修改为斜交角度→A反算主程序[MN-KB]输入：1,2,3,4,5---选择待反算的线路K---待反算点附近的大概桩号M---待反算点X坐标 N---待反算点Y坐标结果：K=反算点的桩号 B=反算点的距中（左负右正中零）使用说明：(数据库[YS-1]中，当P或R的值为零时，可以省略输入0→P或0→R)输程序时，例“×1045→P”的小“×10”为计算器面板的白色键“×10X”子程序[ZDYS]中×1045→P相当于10∧(45)→PPR≥×1090＝＞0→Q即PR≥10∧(90)＝＞0→Q1、当线元为直线时，其起点、止点曲率半径为零。
fx－5800P计算器在工程测量中的应用(三)九、编程常用到数学公式及其注解：1、 sin(α)和cos(α)的定义：以坐标原点为圆心，1为半径，对应于α由0°～360°的变化，圆上各点的x值即为cos(α)值；圆上各点的y值即为sin(α)值；2、由sin(α)和cos(α)的定义可以引伸导出：①cos(α) =－cos(180－α) =－cos(180＋α) = cos(360－α) = cos(－α)②cos(α)= sin(90－α)=sin(90＋α)=－sin(270－α) =－sin(270＋α)③sin(α) =sin(180－α) =－sin(180＋α)=－sin(360－α) =－sin(－α)④sin(α)= cos(90－α) =－cos(90＋α) =－cos(270－α)= cos(270＋α)十、编程命令的输入操作及注解：1、输入操作①功能菜单内输入：（如，“Stop Stop)◢)，该命令的输入已经代替分隔符“：”键，也可代替分隔符“：”2、各命令的注解：（1）、“？”：输入提示。
语法①：？→变量，或“字符串”？→变量；该语法的功能是向变量赋值，当程序执行至此命令时，屏幕出现提示“？”或“字符串？”，且不显示变量的原值，如不输入任何数值而按EXE键时，计算器不会往下继续执行其他命令。
语法②：？变量，或“字符串”？变量；该语法的功能也是向变量赋值，当程序执行至此命令时，屏幕出现提示“变量？”或“字符串？”，且显示变量的原值，如不输入任何数值而按EXE键时，计算器默认变量原值，往下继续执行其他命令。
（2）、“：”：分隔符。
语法：语句：语句：…：语句。
注：语句可为某一数值、变量、表达式、命令或注释文本（“字符串”），如分隔符“：”前为数值、变量、表达式时，每当执行至分隔符“：”，会随之产生一个答案存储器的值Ans。
（3）、“→”：变量赋值。
语法：表达式→变量。
fx-5800p匝道(线元法)测量程序1.ZDZS(匝道正算程序)Do:〝1.AZD 2.BZD 3.CZD 4.DZD 5.EZ D↙ZD=〝？→Z[1]:Lpwhile Z[1]>6 Or Z[1]≤0(零):〝X0(零)=〞？J:〝Y0(零)=〞？N↙J→Z[2]:N→Z[3]↙LbI 0(零):〝K=〞？S:90→T(如将90→T改为〝Q°=〞？T可计算斜交)：〝D=〞？Z:Z[1]=1=＞Prog〝AZD〞: Z[1]=2=＞Prog〝BZD〞: Z[1]=3=＞Prog〝CZD〞: Z[1]=4=＞Prog〝DZD〞: Z[1]=5=＞Prog 〝EZD〞↙〝X=〞:X+Zcos(F+T)→X◢〝Y=〞:Y+Zsin(F+T)→Y◢Prog〝A〞:Goto 0(零)↙本程序中输入及显示说明：输入部分：ZD=? 输1为运行A匝道、输2为运行B匝道、输3为运行C匝道、输4为运行D匝道、输5为运行E匝道；X0=? 输入测站点X坐标；Y0=? 输入测站点Y坐标；K=? 输入待求点桩号；D=? 输入待求点距中桩距离。
(以线路前进方向区分左右、当待求点位于中桩时D=0(零)；当待求点位于中桩左侧时D取负值；当待求点位于中桩右侧时D取正值。
）显示部分：X= 显示待求点X坐标；Y= 显示待求点Y坐标；FWJ= 显示测站点至待求点方位角(即水平角)；HD= 显示测站点至待求点水平距离。
若测站X0坐标输入0(零)时、不论测站Y0坐标输入多少，都不显示FWJ=、HD=两项。
2.ZDFS(匝道反算程序)Do:〝1.AZD 2.BZD 3.CZD 4.DZD 5.EZ D↙ZD=〝？→Z[1]:Lpwhile Z[1]>6 Or Z[1]≤0(零)↙LbI 0(零):〝K=〞？S:〝X=〞？I:〝Y=〞？J: Z[1]=1=＞Prog 〝AZD〞: Z[1]=2=＞Prog〝BZD〞: Z[1]=3=＞Prog〝CZD〞: Z[1]=4=＞Prog〝DZD〞: Z[1]=5=＞Prog〝EZD〞:Porg〝E〞:O(欧)+W→S:〝K=〞:S◢〝D=〞:Z◢Goto 0(零)↙本程序中输入及显示说明：输入部分：ZD=? 输1为运行A匝道、输2为运行B匝道、输3为运行C匝道、输4为运行D匝道、输5为运行E匝道；K=? 输入待求点大概桩号；X=? 输入待求点X坐标；Y=? 输入待求点Y坐标。
(以下程序是专业人士编写，本店铺不对程序负责，仅供您参考使用。
)紫金龙逸zijinlongyi0820@正反算辛普森5800程序2010年4月22日修改（可回头）仲项飞1.主程序(ZXC)2→DimZ(增加变量)Lbl 0：Z[9] →E:“1-ZS,2-FS” ?E：E→Z[9]↲Lbl1: ”XL LC” ?Z:Prog“DA TE1”:E≥2=＞Goto2 : 5.375(设计线路在线路中线右边为“-”，设计线路在线路中线左边为“+”,数值为设计线路与线路中线的距离，所有偏距以隧道中线起算;若算斜交，则改为0，则所有偏距以设计线路起算。
)→N:“DZX”?W:“PJ”?S:W+N→W:Abs(Z-F)→Q:Prog“XZS”:W-N→W:J-S→J:Cls:Locate1,1,”X=”: Locate 3, 2,X:Locate 1, 3,”Y=”:Locate 3, 4, Y◢"FWJ=”:J:J►DMS◢Goto 3Lbl 2:?X:?Y:Prog“XFS”: F+Q→Z:-5.375→N:W-N→W:Cls:Locate 1, 1,”K=”: Locate 3, 2,Z: Locate 1, 3, ”W=”: Locate 3, 4, W◢Goto 3Lbl 3：Prog “SQX1”DA TE1(数据库)7704.149→A:5024.7762→B:264°42'42.62"→C:1÷10^(45)→D: 1÷10^(45)→R:122020.554→F:124140.538→G:IfZ≤G:Then GotoA:IfEnd↲7508.761→A:2913.8154→B:264°42'42.62"→C:1÷10^(45)→D: 1÷(-750)→R: 124140.538→F:124360.538→G:IfZ≤G:Then GotoA:IfEnd↲………………(有几个线圆依次添加几个判断语句，顺序从小里程到大里程，次序不得颠倒)Lbl A………………(有几条线路就添加几个数据库。
CASIO fx[/B]—5800计算器工程测量与公路测量程序[/B](1)[/B]小写o为大写字母“O”，与零区分。
[/B][/B](2)[/B]程序中没有空格符号。
[/B]1[/B]、坐标正算[/B]〝X0=〞?X：〝Y0=〞?Y：〝I=〞?I：〝J=〞?J X+ICos(J)→U：Y+ISin(J)→V〝X=〞：U◢〝Y=〞：V◢说明：X0 Y0：已知点坐标，I：两点的距离，J：方位角[/B]２、坐标反算[/B]Lbl 0-〝X1=〞?X：〝Y1=〞?Y：〝X2=〞?U：〝Y2=〞?VPol(U-X，V-Y)：J<0J+360→J〝I=〞：I◢〝J=〞：JDMS◢Goto 0-说明：X1Y1：第一点的坐标，X2Y2第二点的坐标，I：两点的距离，J：方位角[/B]３、圆曲线[/B]〝X0=〞?X：〝Y0=〞?Y：〝FWJ=〞?D：〝QDHAO=〞?G：〝ZDHAO=〞?C：〝R=〞?RLbl 1-〝DQHAO=〞?LL<G Or L﹥CGoto 2〝PY=〞?K180(L-G)÷(πR)→E2Rsin(0.5E)→FRec(Abs(F)，D+0.5E)〝X=〞：X+I+Kcos(D+E+90)→A◢〝Y=〞：Y+J+Ksin(D+E+90)→B◢D+E→o：o<0360+o→o〝FWJ=〞：oDMS◢Goto 1-Lbl 2-〝END〞说明：X0Y0：起始点的坐标，FWJ：起始点的方位角，QDHAO：起点里程，ZDHAO终点里程[/B]，R：半径，DQHAO：待求点里程，PY：偏移量[/B]Abs：绝对值function按键在按1 在按“下”第一个就是abs Cls 编程状态下按FUNCTION，按3，按两次下翻键，按5就可以了小写子母：编程状态下按FUNCTION，下翻按1，就可以输入小写子母和下标字符４、竖曲线[/B]〝QZHAO =〞?J：〝HB〞=?B：〝I1=〞?C：〝I2=〞?D：〝R=〞?R：〝T=〞?T：((D-C)÷100)÷Abs((D-C)÷100)→A：Abs(0.5R(D-C)÷100)→S〝T〞：S◢J-T→XJ＋T→YLbl 0-〝DQ=〞?LIf L<X：Then Goto1：Else If L﹥Y：Then Goto1：If End：If End〝GC=〞：B+C(L-J)÷100+（A(L-J+T)）²平方÷(2R)→H◢）Goto 0-Lbl 1-〝END〞说明：QZHAO：竖曲线中点里程，HB：交点高程（不是交点路面高），I1：第一坡度，I2：第二坡度，R：半径，T：切线长，DQ：待求点里程。
卡西欧5800P 测量程序汇总CASIO fx5800P公路施工放线测量公式及程序一、已知座标，求平距和方位角(座标反算)：公式：D=√（Xp-Xo）2+（Yp-Yo）2 α=arctg(Yp-Yo)/(Xp-Xo) 程序：“A”？→A：“B”？→B：Lbl 0:“X”？→X：“Y”？→Y：（X-A）→M：（Y-B）→N：“D=”：√（M2+N2）⊿ tan-1(N/M) →C:If M＜0:Then “Q=”：180+C →Q ⊿ Else If N＞0: Then “Q=”：C→Q ⊿ Else “Q=”：360+C→Q ⊿ If End : If End : Goto 0 说明：（A，B）为测站点坐标，（X，Y）为所求点坐标。
输出：D为平距，Q为方位角。
二、已知直线的坐标方位角Q和直线起点坐标（Xo,Yo），求直线上任一点的中桩坐标（X， Y），左右边桩坐标（XL，YL）、（XR，YR）：公式：X =Xo+LcosQ Y=Yo+LsinQ程序：“A”？→A：“B”？→B：“Q”？→Q：“Z”？→Z：Lbl 1: “P”?→P：“J”？→J：“M”？→M：Abs（M-Z）→W：“X=”：A+W*cos(Q)+P cos(Q+J)→X ⊿“Y=”：B+W*sin(Q)+P sin(Q+J) →Y ⊿Goto 1 说明：（C，D）为直线起点坐标，Q为直线方位角，Z为起点桩号，L为所求坐标点桩号。
“U”为左边距，“T”为右边距，“V”为偏角；U=0时不算边桩坐标。
输出：（X，Y）为中桩坐标，（XL，YL）为左边桩坐标，（XR，YR）为右边桩坐标。
三、已知圆曲线起点坐标（U，V），切线方位角Q，桩号Z和圆半径R，求圆曲线上桩号为 L的点中桩坐标（X，Y），左右边桩坐标（XL，YL）、（XR，YR）：公式：ψ=90L/(лR) (偏角公式） C=2Rsin ψ (对应弧的弦长公式）弦的方位角:Q=Qo±ψ（曲线左转时为“-”）程序：“U”？→U：“V”？→V：“Q”？→Q：“R”？→R：“Z”？→Z：“W=-1，1”：？→W：“ZJ=”：？→Z[1]：“YJ=”：？→Z[2]：“PIAN JIAO”：？→T：Lbl 2: “M”？→M：180*（M-Z）/（2π*R）→J：R*2sin(J)→K:If W=-1:Then“X=”:U+K*cos (Q-J）→X⊿ “Y=”：V+K*sin(Q-J)→Y ⊿“Q=”：Q-180*（M-Z）/（πR）→O⊿ Else If W=1:Then“X=”:U+K*cos(Q+J）→X⊿“Y=”：V+K*sin(Q+J)→Y ⊿ “Q=”：Q+180*（M-Z）/（πR）→O⊿ IfEnd:IfEnd:T=0=＞Goto 2：“XL=”:X+Z[1]*cos(O-T)→F⊿“YL=”：Y+Z[1]*sin(O-T)→P⊿“XR=”:X+Z[2]*cos(O+T)⊿“YR=”：Y+Z[2]*sin(O+T) ⊿ Goto 2 说明：W=-1时曲线左转, W=1时曲线右转。
“ZJ=”为左边距，“YJ=”为右边距，“PIAN JIAO”为偏角，偏角输0时不算边桩坐标。
输出：（X，Y）为中桩坐标，“Q=”为所求点方位角，（XL，YL）为左边桩坐标，（XR，YR）为右边桩坐标。
四、已知直缓点坐标（M，N）、方位角Q、桩号Z，缓和曲线全长S和连接圆半径R，求缓和曲线上任一点（桩号为L）的中桩坐标（X，Y），左右边桩坐标（XL，YL）、（XR，YR）：公式：θ=30L2/(πRLs) C=L-L5/(90R2Ls2) α=α±90Ls2/(πA2)程序：“M”？→M：“N”？→N：“R”→R：“A”？→A：“S”?→S：“Z”？→Z：“Q”？→Q：“W=-1，1”：？→W：“ZJ=”？→Z[1]：“YJ=”？→Z[2]：“PIAN JIAO”：？→ T：Lbl 3:“L”？→L：Abs(L-Z)→H:30*H2/(πRS)→D:H-H5/(90R2S2)→C:If W=-1:Then “X=”：M+C*cos(Q-D)→X⊿“Y=”：N+C*sin(Q-D)→Y ⊿ A =0 =＞Goto3:“Q=”:Q-90H/(πA）→E ⊿ Else If W=1:Then “X=”：M+C*cos(Q+D)→X⊿“Y=”：N+C*sin(Q+D)→Y ⊿ A=0 =＞Goto3:“Q=”:Q+90*H/(πA）→E ⊿IfEnd:IfEnd:T=0=＞Goto 3：“XL=”:X+Z[1]*cos(E-T)⊿“YL=”：Y+Z[1]*sin(E-T)⊿“XR=”:X+Z[2]*cos(E+T)⊿“YR=”：Y+Z[2]*sin(E+T)⊿ Goto 3说明：“A”为缓和曲线参数，W=-1时曲线左转, W=1时曲线右转。
“ZJ=”为左边距，“YJ=”为右边距，“PIAN JIAO”为偏角，偏角输0时不算边桩坐标。
输出：（X，Y）为中桩坐标，“Q=”为所求点方位角，（XL，YL）为左边桩坐标，（XR，YR）为右边桩坐标。
五、已知卵形曲线（即非完整缓和曲线）大圆半径R，小圆半径r，缓和曲线参数A，以大圆的圆缓点为起点，其坐标为（C，D），其方位角为Q，其桩号为Z，求桩号为L的点的中桩坐标（X，Y），左、右边桩坐标（XL，YL）、（XR，YR）。
公式：偏角：θ=arctg[(Yp-Yo)/(Xp-Xo)]-90LR/(πR)Xj=Lj-Lj5/(40A4)+Lj9/3456A8-…… ，Yj=Lj3/(6A2)-Lj7/(336A6)+Lj11/(42240A10)-…… Lj=A/Rj(Rj为曲率半径）S=√(Xp-Xo)2+(Yp-Yo)2（弦长公式）α=Q±θ（弦线方位角公式） Xp=Xo+S *cosα, Yp=Yo+sinα程序：“A”？→A：“C”？→C：“D”？→D：“R”？→R：“Q”？→Q：“Z”？→Z：“W=-1，1”？→W：A/R→J:J-J5/(40*A4)+J9/(3456*A8)→E:J3/(6*A2)-J7/(336*A6)+J11/(42240*A10)→F :90*J/(πR)→B:“ZJ=”？→Z[1]：“YJ=”？→Z[2]：“PIAN JIAO”：？→T：Lbl 4:“L”？→L: J+Abs(L-Z)→S:S-S5/(40*A4)+S9/(3456*A8)→M:S3/(6*A2)-S7 /(336*A6)+S11/(42240*A10)→N:Abs(M-E)→G:Abs(N-F)→H:√(G2+H2)→K:arctg(H/G)→I:If I＜0:Then I+360→U:Else I→U: IfEnd: I f W=-1:ThenQ-(U-B)→V:“X=”：C+K*cos(V)→X⊿“Y=”：D+K*sin(V)→Y⊿“Q=”:Q-90*S2/(πA2)-90*J2/(πA2) →O ⊿ Else If W=1:Then Q+(U-B)→V:“X=”：C+K*cos(V)→X⊿“Y=”：D+K*sin(V)→Y⊿“Q=”:Q+90*S2/(πA2)-90*J2/(πA2) →O ⊿IfEnd:IfEnd:T=0=＞Goto4:“XL=”:X+Z[1]*cos (O-T)⊿“YL=”：Y+Z[1]*sin(O-T)⊿“XR=”:X+Z[2]*cos(O+T)⊿“YR=”：Y+Z[2]* sin(O+T)⊿ Goto 4说明：“A”为缓和曲线参数，W=-1时曲线左转, W=1时曲线右转。
“ZJ=”为左边距，“YJ=”为右边距，“PIAN JIAO”为偏角，偏角输0时不计算边桩坐标。
输出：（X，Y）为中桩坐标，“Q=”为所求点方位角，（XL，YL）为左边桩坐标，（XR，YR）为右边桩坐标。
六、竖曲线计算公式：H=X2/(2R)程序：“A”？→A：“H”？→H：“D”？→D：“T”？→T：“R”？→R：“I1”？→U：“I2”？→V：“K(ZHONG-DIAN)”：D+T→B ⊿“LY=”：B-A→E ⊿ Lbl 6:“M”？→M：“G”？→G：If M≤A: Then D-M→L:H-LU→P:IfEnd:If M＞A And M＜D: ThenM-A→L:H-U(D-M)+GL2/(2R)→P: IfEnd:If M＞D And M＜B:ThenB-M→L:H+V(M-D)+GL2/(2*R)→P:IfEnd:If M≥B: Then M-D→L:H+LV→P: IfEnd: “H(SHEJI)=”：P ⊿ Goto 6说明：“A”为竖曲线起点桩号，“H”为交点高程，“D”为交点桩号，“T”为切线长，“R”为竖曲线半径，“I1”、“I2”为第一、二坡度。
输出：“H(SHEJI)=”为设计高程。
注意的是：计算范围不能超出到下一个竖曲线范围内。
G取+1，或者-1，当“I1”？→U《0，“I2”？→V 必然》0，即G=1反之G=-1卡西欧5800计算器正反算及缓和曲线计算程序时间:2011-06-15 21:19:47 来源:自编作者:杨连双我要投稿我要收藏投稿指南5800实用程序好用实用FW J -0（反算坐标）Lb1 0：“X1=”?A:”Y1=”?B↓“X2=”?C↓“Y2=”？D↓Po1(C-A,D-B):C1s↓J <0=>J+360->J↓“S=”:I▲“F(DMS)=”:J◆DMS▲Goto0↓备注;标A=起算点X1坐标 B=起算点Y1坐 C=起算点X2坐标 D=起算点Y2坐标ZX-1（直线程序）“QSX=”?A: “QSY=”?B: “QSL=”?Q: ?F↓Lb1 1: “JSL=”?J:? L:?P↓A+(J-Q)COS(F)+LC OS(F+P)→C↓“X=”:C▲B+(J-Q)Sin（F）+LSin(F+P)→D↓“Y=”:D▲Goto1备注：标A=起算点X1坐标 B=起算点Y1坐 Q=起算点里程.F=起算点方位角 J=计算点方位角YQX-2(圆曲线)“QSX=” ?A: “QSY=”?B ：“QSL=”?Q:?F:?R↓Lb1 2: “JSL=”?J:?L:?P↓J-Q→S:90S÷(兀R) →O: √(2RSin(0))2→I:A+ICOS(F+O)→E:B+ISin(F+O) →G↓E+LCOS(F+2O+P) →C↓“X=”：C▲G+LSin(F +2O+P ) →D↓“Y=”:D▲F+180×(J-Q) ÷（兀R）→G↓“T(DMS)=” :G◆DMS▲Goto2备注：A=起算点X1坐标 B=起算点Y1坐 Q=起算点里程 F：起算点里程 R：设计半径 J：计算点理程 L；左右偏距离 P:左右偏方位角HQX-4（缓和曲线）“QSX=”?A: “QSY=”？B: “QSL=”?Q: “HXC=”?H:?F:?R↓Lb1 4:“ JSL=”?J:?L:?P↓J-Q→N:N×√((1-(N2)2÷(40R 2H2)) 2+(N2÷(6RH)) 2)→M↓tan-1 （（20RHN2÷(120R2H2-( N2)2 )）→k:90N2 ÷(兀RH)→G↓A+Mcos(F+K)+Lcos(F+G+P)→C↓“X=”:C▲B+Msin(F+K)+Lsin(F+G+P) →D↓“Y=”:D▲Goto4A=起算点X1坐标 B=起算点Y1坐 Q=起算点里程H:缓曲线曲缓线的长度 F:计算点方位角 R：设计半径J：计算点理程 L；左右偏距离 P:左右偏方位角CASIO Fx-5800P计算器曲线万能程序时间:2011-06-09 09:42:20 来源:本站作者:李立军我要投稿我要收藏投稿指南CASIO Fx-5800P计算器万能程序QXJS---主程序〝L=〞?L:〝LI CHENG=〞?W <┘Lbl A <┘For 1→N To 20 <┘L=1＝＞Pros 〝1XY〞<┘L=2＝＞Pros 〝A〞<┘If W≥ReP(Z) And W≤ImP(Z):Then Pros 〝MP〞:Goto 0:If End:Next <┘〝--LI CHENG----?〞:Stop <┘Lbl 0 <┘If K<0:Then Q→X:0→Q:0→J:S+C→I:Abs(K) →K:Goto 3:If End <┘Goto 2 <┘Lbl 1 <┘Q<0＝＞360+Q→× <┘〝LI CHENG=〞?W <┘If W≥ReP(Z) And W≤ImP(Z):Then Goto 2:E lse〝>>>>〞:Goto A:If End <┘Lbl 2 <┘0→J<┘uo+Abs(W-Abs(K)) →I <┘W:K<0＝＞-W <┘If Ans < k : Then 0→U : -I : Goto 4 : If End <┘I>S+C＝＞S+C→I <┘Lbl 3 <┘If I<S:Then I²r÷(2RS) →U:I-I^(5)÷(40R²S²)+i(I^(3)÷(6RS)-I^(7)÷(336R^(3)S^(3)):Else (2I-S)ｒ÷(2R) →U:Z［2］+R(sin(U)-icos(U)+i): If End <┘Lbl 4 <┘J+Abs cos(Q)+Ans i sin(Q) →J <┘If I=S+C:Then Q+A→Q:-K-I-T→K:T→I:S→T:I→S:-R→R:Porg〝MP〞:Goto 3: If End <┘If Q+A=0:ThenX-A-180→Q:K+ReP(J)+ImP(J)÷tan(A) →K:O+ImP(J)(cos(Q)+isin(Q))÷sin(A) →O:Goto 2: If End <┘If I=S:Then Q-180→Q:O+J→O:Goto 2: If End <┘K<0⇒U-180→U <┘Q+U→U <┘〝BIAN JU=〞?E <┘ReP(O)+ReP(J)-Esin(U) →X:ImP(O)+ImP(J)+Ecos(U) →Y <┘Cls:Fix 3:Locate 1,1,〝X=〞:Locate 6,1,X:Locate 1,2, 〝Y=〞: Locate 6,2,Y◄X→Z［3］:Y→Z［4］<┘If E=×10-6:Then Prog〝BYXY〞:Goto 1: If End <┘Goto 1 <┘Prog〝ZBFS〞:Norm 2:Goto 1 <┘说明：1．L=12．输里程桩号3．输边距：左“负”，右“正”，中心桩“0”，*10-6时为反算，“M0=”为X；“P=”为Y,W OK为计算里程结果，I为左右距4.本程序为万能程序,只用编缉数据库即可运算.1XY—数据库If N=1:Then0→S:400→T:-4500→R:4523.084→C:300°24′42″→Q:69287.5839+506028.4 703i→O:15701.589→K:15701.589+21331.515i: If End <┘If N=2:Then400→S:400→T:3800→R:2914.665→C:245°22′6.4″→Q:68953.9506+500686 .4792i→0:21331.515→K:21331.515+26411.098i: If End <┘If N=3:Then…If N=4:Then… N=Ans⇒0→Z＜┘Ans⇒Z＜┘ Return＜┘说明：N=?—本段曲线；1.S—起点缓和曲线长度；2.T—止点缓和曲线长度；3.R—曲线半径；4.C—曲线长度；5.Q—计算方位角；6.O—起点坐标；7.K—起点桩号；8.本段曲线起止点桩号。
