本人急求算缓和曲线准的程序（5800），并且能放匝道。

我的程序在算缓和曲线时差很多，缓和曲线长时最多差9毫，所以很想求一套算的准的程序。全线的。

 

YANG6815475:互通式立交、匝道多线路正、反算综合程序(4800P优化程序)(四) 2007-11-25 21:33 (十六)ZDY-8程序(互通式立交、匝道多线路综合计算程序)   ⑴有关资料    采用复化辛普森公式：　　   X=XA+H / 6（cos aA+4∑cos aK+1/2 +2 ∑cos aK +cos ai）（K=0～n-1）　　   Y=YA+ H / 6（sin aA+4∑sin aK+1/2 +2 ∑sin aK +sin ai） （K=0～n-1）  　式中：aA为回旋曲线起点A的切线方位角；aK+1/2为里程DXK+1/2点切线方位角；aK为里程DXK点切线方位角；ai为里程DKi点切线方位角。　对于上式，虽然是由回旋曲线导出的，但该式也适用直线段和圆曲线段。    ⑵ ZDY-8程序(主程序)AC   MODE   5   1   ZDY   ALPHA — 7   EXE   1 Lbl 0 ∶ “1.KB=>XY”: “2.XY=>KB”: :J=1=> Goto   1∶≠> GOTO 2∶LbI 1∶｛H L O｝∶H“KP”L“B”O∶Prog ″ABCD1″∶Prog ″XY 9 ″∶Prog ″XY10 ″∶ Goto   1∶Lbl 2∶Prog ″XY 8 ″∶ GOTO   0   EXE   ⑵ 支程序    ①XY8程序(反算程序)AC   MODE   5   1   X Y   ALPHA   8   EXE   1 Lbl 0 ∶H“KP”∶ Prog ″ABCD1″∶Prog ″XY9 ″∶∶L”X0” ∶O”Y0” ∶U= Z - 90∶P= (O-N)cos U -(L-M) sin U   ∶AbsP≥0.001=>H=H+P∶GOTO 0 : ≠>H” KI =” ▲   V= (O -N)cos Z - (L -M) sin Z ▲      ②XY9程序(正算程序)AC   MODE   5   1   XY   ALPHA   9   EXE   1 P=（E－D）÷（G－F）∶ Q= H－F∶    I=PQ ∶R =8∶LbI 3∶R［R］=C＋RQ(RI÷8＋2D)r÷16∶DSZ R∶Goto   3∶R＝C∶M“X”=A＋Q∑（cosR［2K］＋4 cos R［2K＋1］＋cos R［2K＋2］ ，K,0,3）÷24 ▲ N“Y”=B＋Q∑（sinR［2K］＋4 sin R［2K＋1］＋sin R［2K＋2］ ，K,0,3）÷24   ▲ Z“QI”=360 Fra