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

缓和曲线各种迭代算法及比较 半只烟（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度制，取消小数位数