缓与曲线上各种迭代算法与比较.doc

缓和曲线各种迭代算法及比较 半只烟（850570455） 关于缓和曲线的直接计算式都是采用近似计算，因而其计算精度和参数有关，不同的参数得到的计算精度是不一样的，那么很自然的会想到，有没有一种计算方法使计算结果达到一给定的精度后才结束过程而和参数无关，答案是肯定的。通常利用程序使用数值计算的迭代方法。下面给出常用的使用迭代原理进行计算的变步长辛普森积分法和高斯-勒让德求积法，至于别的方法，此处不再详述，有兴趣的可参阅数值计算方法方面的资料。 变步长辛普森迭代求积法 变步长辛普森积分法是计算定积分的经典方法，其计算步骤如下： 用梯形公式计算，其中n=1,h=b-a,且令Sn=Tn。 用变步长梯形法则计算 用辛普森求积公式计算 若，则令转到步骤（2）继续计算；否则结束，S2n即为所求的积分近似值。其中为事先给定的求积精度。 由于需要对被积函数求值，先给出求解回旋线的函数值的子程序FX，用于求解回旋线上距起点x处的X坐标，求解Y坐标只需把cos(余弦函数)改成Sin(正弦函数),此处不在给出。 FileName:FX 计算f(x) 行数 代码 说明 3 Cos(J+A×X+(B-A)×X2÷(2L)→X 计算Y时，需把cos改成sin，变量见后 File SIMPSON 复化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和积分上限H H→X:Prog “FX” 计算f(b)的值,b=H X→Y 暂存 0→X:Prog “FX” 计算f(a), f(a)= f(0) H×(X+Y)÷2→T 计算Tn, T→S:1→E:0→w Sn=Tn初始计算精度，W为迭代次数 While E0.0001 大于给定的精度就进行迭代 0→P For 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 准备迭代参数 WhileEnd Deg：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和求积系数λk n 求积结点tk 求积系数λk 2 ±0.5773502692 1.0000000000 3 0.0000000000 0.8888888889 ±0.7745966692 0.5555555556 4 ±0.3399810436 0.6521451549 ±0.8611363116 0.3478548451 5 0.0000000000 0.5688888889 ±0.5384693101 0.4786286705 ?.9061798459 0.2369268851 6 ?.2386191861 0.4679139346 ?.6612093865 0.3607615731 ?.9324695142 0.1713244924 7 0.0000000000 0.4179591837 ?.4058451514 0.3818300505 ?.7415311856 0.2797053915 ?.9491079123 0.1294849662 8 ?.1834346425 0.3626837834 ?.5255324099 0.3137066459 ?.7966664774