缓与曲线各种迭代算法及比较.docVIP

缓和曲线各种迭代算法及比较半只烟（850570455）关于缓和曲线的直接计算式都是采用近似计算，因而其计算精度和参数有关，不同的参数得到的计算精度是不一样的，那么很自然的会想到，有没有一种计算方法使计算结果达到一给定的精度后才结束过程而和参数无关，答案是肯定的。通常利用程序使用数值计算的迭代方法。下面给出常用的使用迭代原理进行计算的变步长辛普森积分法和高斯-勒让德求积法，至于别的方法，此处不再详述，有兴趣的可参阅数值计算方法方面的资料。变步长辛普森迭代求积法变步长辛普森积分法是计算定积分的经典方法，其计算步骤如下：用梯形公式计算，其中n=1,h=b-a,且令Sn=Tn。用变步长梯形法则计算用辛普森求积公式计算若，则令转到步骤（2）继续计算；否则结束，S2n即为所求的积分近似值。其中为事先给定的求积精度。由于需要对被积函数求值，先给出求解回旋线的函数值的子程序FX，用于求解回旋线上距起点x处的X坐标，求解Y坐标只需把cos(余弦函数)改成Sin(正弦函数),此处不在给出。FileName:FX计算f(x)行数代码说明3Cos(J+A×X+(B-A)×X2÷(2L)→X计算Y时，需把cos改成sin，变量见后FileSIMPSON复化Simpson迭代求回旋线坐标增量行数代码说明Deg:Norm设置360度制，取消小数位数“QiDianFangWei”?→J起点方位角J×180÷π→J:Rad转换成弧度，设置成弧度计算“QiDianBanJin”?→A起点半径，直线为0“ZhongDianBanJin”?→B终点半径，直线为0“HHQX Chang”?→L缓和曲线全长“JiSuanChang”?→D计算点到起点的长度If A0:Then 1÷A→A:IfEnd计算起点曲率If B0:Then 1÷B→B:IfEnd计算终点曲率1→N:D→H初始节点数N和积分上限HH→X:Prog “FX”计算f(b)的值,b=HX→Y暂存0→X:Prog “FX”计算f(a), f(a)= f(0)H×(X+Y)÷2→T计算Tn,T→S:1→E:0→wSn=Tn初始计算精度，W为迭代次数While E0.0001大于给定的精度就进行迭代0→PFor 0→K To N-1(K+0.5)×H→X:Prog “FX”P+X→P:W+1→W在For循环中计算累计和，并计数Next(T+H×P)÷2→U计算T2n(4×U-T)÷3→V计算结果Abs(V-S)→E修改迭代值U→T:V→S准备迭代参数N+N→N:H÷2→H准备迭代参数WhileEndDeg：V◢设回成角度模式，计算结果W◢迭代次数高斯-勒让德迭代求积法对定积分的积分变量x作变换，将原积分转化为区间[-1，1]上的积分，即由差值求积公式有其中，（k=0,1,2,……,n-1）,为区间[-1，1]上的n个求积结点，且，如果n个结点（k=0,1,2,……,n-1）取勒让德多项式在区间[-1，1]上的n个零点，则上述差值求值公式称为高斯-勒让德求积公式，其代数精度为2n+1。下表给出了n从2到10间的求积结点tk和求积系数λkn求积结点tk求积系数λk2±0.57735026921.0000000000 30.0000000000 0.8888888889 ±0.77459666920.5555555556 4±0.33998104360.6521451549 ±0.86113631160.3478548451 50.0000000000 0.5688888889 ±0.53846931010.4786286705 ±0.90617984590.2369268851 6±0.23861918610.4679139346 ±0.66120938650.3607615731 ±0.93246951420.1713244924 70.0000000000 0.4179591837 ±0.40584515140.3818300505 ±0.74153118560.2797053915 ±0.94910791230.1294849662 8±03626837834 ±0.52553240990.3137066459 ±0.79666647740.2223810345 ±0.96028985650.1012285363 90.0000000000 0.3302393550 ±0.32425342340.3123470070 ±0.61337143270.2606106964 ±0.83603110730.1806481607 ±0.96816023950.0812743884 10±0.14887433900.2955242247 ±0.43339539410.2692667193 ±0