单个的螺旋线.pptVIP

* 单个的螺旋线 [问题]帮助一家生物技术公司就位于空间中一般位置的一个螺旋线和一个平面交点的“实时”定位设计、证明、编程，并测试检验一个数学算法。 类似的计算机辅助几何设计（CAGD）程序能使工程师看到他们设计的物体——例如飞机喷气发动机、汽车悬挂物或医疗仪器——的平面截面部分。此外，工程师还可以把诸如气流、压力或温度等物理量用颜色或等位线标码显示在该平面截面上。更进一步，工程师可以让这种平截面迅速扫过整个物体来得到物体的三维成像以及物体对运动、力和热的反应。为取得这样的效果，计算机程序必须能以足够的速度和精度来定位所设计物体的每一部分和所观察平面的全部交点。 通过方程求解器（eqution solvers），原则上计算出这种交点。但对于特定的问题，已经证明特定的方法比通用的方法算得更快、更准确。特别就完成实时计算来说，计算机辅助几何设计的通用软件已被证明太慢了，或对由该公司研制的医疗器来说该通用软件太大了（杀鸡用牛刀）。上面所说的情况导致该公司要考虑下述问题：设计、证明、编程并测试检验一种算法，它们计算处于一般位置（即处于任何地方、任何指向）的平面和螺旋线的全部交点。 螺旋线的一段可以表示，例如一段螺旋状的弹簧或化学仪器或医疗仪器中的一段管状物，见图6.22 。 需要对所建议的算法给以某种理论证明是由于必须要用几种观点来证实解法的必要性，这可以通过对部分算法的数学证明，并用已知的例子来检验最终的程序，政府的医用部门要求这个文件和检验报告。 本题是1994年MCM的赛题，如题目所讲，对于特定的问题，特定的方法比通用的方法算得更快、更准确。这实际上体现了数学模型中具体问题具体分析的灵魂，增加了创造性。 本题的难点在于： ①求出全部交点，这实际上隐含了准确地计算出交点的个数，也隐含了已知交点的性质，如重根是否存在，重根有几个，还包含找到算法，并求出这些交点来； ②实时求出来，要求算法快，特别在交点很多的情况下这似乎很难得到解决； ③平面螺旋线处于空间任意位置，特别当问题的坐标系是给定的情况下，给求解带来困难，虽然经过坐标变换可以避免这一问题，但也增加了结果（也就是交点）再用原坐标系表达的问题。 一、数学模型 首先还是简化。坐标系先简化，取螺旋线的轴为z轴，x轴取通过螺旋线上一点，若螺旋线长度有限，取其起点为x轴上一点。这样，螺旋线与平面方程为 （6.8-1） （6.8-2） 式中r为螺旋线半径， h 是与螺旋距有关的常数，h为螺距的 ， θ为向径与x轴夹角， A、B、C、D为任意实数。 理论上，（6.8-1）式与（6.8-2）式联立求解即可求出全部交点，但这里4个方程，4个未知数x,y,zθ、又是非线性方程，快速求解显然是困难的，所以还应根据具体问题作进一步简化。 可否进一步简化？应向哪一个方向简化？这是具有创造性的问题。 但是如果我们设想将螺旋线置于一个平行于其轴的平面M的前方，且让截面所在平面与平面M垂直，将螺旋线与截平面向平面M作投影，则螺旋线投影后为一曲线，截平面投影为直线，z轴也投影为一直线，直线与曲线间可能会有交点，且这种交点和螺旋线与平面之间的交点是一一对应的。 因为根据（6.8-1）式螺旋线上任意两点z的坐标不同，而向平面M作投影，z坐标的值未发生变化。而由平面M上曲线与直线交点可以很容易根据z坐标定出θ，从而求出x,y坐标。 这说明可以将上述求三维空间中螺旋线与平面交点问题等价转化为求平面上曲线与直线交点问题。这样既解决了可否简化，又解决了向什么方向简化的问题。 根据上述分析，将（6.8-1）式代入（6.8-2）式，当 时得 其中有两个三角函数还可再化简为 其中 所以 由上式变为 （6.8-3） 不妨假设 ， 则否令 这样将求4个未知数方程组的解简化为求非线性方程的解，且 的根与方程组中θ的解仅差一个已知的常数 。 当 即A=0,B=0时，由平面方程Cz+D=0,可解出 ，代入（6.8-1）式，得唯一交点 。 若a=0时，即Ch=0 ,则(6.8-3)式变为