复化辛卜生公式可以算是公路线形坐标计算的万能公式.docVIP

复化辛卜生公式可以算是公路线形坐标计算的万能公式。它不仅对直线、圆曲线、缓和曲线，不完整的缓和曲线（卵形曲线）都能用一套公式进行计算。当然它的精度可以人为判断，不过它使用起来还是比较繁琐，一不小心更容易出错，对直线、圆曲线的坐标计算更无必要用此公式。但是如果在Casio上编一套程序，对任何线形只要用这一套程序，只要输入起算点坐标、起始方位角、桩号几个参数，使用起来就很方便了。笔者编程如下：File1FHXPS（复化辛卜生）L1???????Lb1??0????L2??????｛L｝L3??????W=（1∕B-1∕A）∕(S-K)L4???????O=1∕A＋W×(L-K)L5???????I=N＋(O＋1∕A)(L-K)×180∕π∕2L6???????V=1∕A＋W×(L-K)/2?L7???????U=N＋(V＋1∕A)(L-K)×180∕π∕4L8???????C=1∕A＋W×(L-K)∕4L9???????D=N＋(C＋1∕A)(L-K)×180∕π∕8L10??????E=1∕A＋W×(L-K)×3∕4L11??????F=N＋(E＋1∕A)×(L-K)×180×3∕π∕8L12??????R＂YPPQX:1,ZPPQX:-1＂L13??????R=1??=ProgYPPQX:=\??R=-1=??ProgZPPQX⊿⊿L14??????Goto?0File2??????YPPQX??L1??????X=G＋(L-K)∕6∕2×(cosN＋4(cosD＋cosF)＋2cosU＋cosI)???L2???????Y=H＋(L-K)∕6∕2×(sinN＋4(sinD＋sinF)＋2sinU＋sinI)?File3??????ZPPQX??L1???????X=G＋(K-L)∕6∕2×(cosN＋4(cosD＋cosF)＋2cosU＋cosI)???L2???????Y=?H＋(K-L)∕6∕2×(sinN＋4(sinD＋sinF)＋2sinU＋sinI)?A??----??起点半径（曲率）B??----??终点半径（曲率）S??----??终点桩号K??----??起点桩号L??----??待求点桩号I??----??待求点切线方位角N??----??起始切线方位角V??----??1∕2等分点曲率U??----??1∕2等分点切线方位角C??----??1∕4等分点曲率D??----??1∕4等分点切线方位角E??----??3∕4等分点曲率F??----??3∕4等分点的切线方位角G---起算点X坐标H---起算点Y坐标当然你在程序里可以加上一个计算边桩的子程序。演示一、?直线例1：一直线起点桩号为0（X=4000，Y=3000）起始方位角N=100°要求桩号为K0＋050的中桩坐标????????????????操作?????????????????????显示????????????说明??????????????AC?FILE???????????????????FHXPS????????????????EXE??????????????????????B？????????????10xy99?EXE???????????????????A?????????????(由于直线的曲率是+∞，而计算器没有+∞的表示方式，只用一个相当大的数值代替)????????????10Xy99?EXE???????????????????S？?????????????100?EXE?????????????????????K？??????????(在直线计算中终点桩号你可以输入一个真实桩号也可以输入只要大于计算点桩号就可以)???????????????0?EXE??????????????????????L？??????????????50?EXE??????????????????????N？??????????????100°?EXE??????????????????I=?100???????????????EXE???????????????????YPPQX:1,ZPP????????????????1?EXE??????????????????????G？????????????????????4000??EXE??????????????X=3991.317591????????(设计值3991.317591)???????????????????EXE????????????????????????H？???????????????3000?EXE???????????????Y=3049.2