最速降线优秀培训书.pptVIP

即 所以 可令 要使 最小，相当于求多元函数的极值，因此，令 科学出版社 南通大学理学院计算科学与统计学系 * * 又由 于是有 所以 同理得 …… 科学出版社 南通大学理学院计算科学与统计学系 * * +…… 于是 得到： 科学出版社 南通大学理学院计算科学与统计学系 * * ，（ 为给定的误差），则 为所求的结果。 再将所求出的 值代入上述方程中，求出相应的 ，由各点 由上式，可用二分法求 值，设 ，且 ，计算 为了求 ，可令： 所连成的曲线就是所求的“最速降线”。 科学出版社 南通大学理学院计算科学与统计学系 * * 用Matlab编程如下： H=10.0;n=16;h=H/n; g=9.8;f=1.0; a =0; b=2/(sqrt(2*g*(n-1)/n*H)+sqrt(2*g*H)); c=(a+b)/2; i=1; while abs(f)1e-10 s=0; for j=1:n v=(sqrt(2*g*j*h)+sqrt(2*g*(j-1)*h))/2; s=s+v/sqrt(1.0-c^2*v^2); end f=c-10.0/(h*s); if f0 科学出版社 南通大学理学院计算科学与统计学系 * * b=c; else a=c; end c=(a+b)/2; i=i+1; end x(1)=sqrt(9.8*h/2)*c*h/sqrt(1.0-c*c*v*v); for k=2:n v=(sqrt(2*g*k*h)+sqrt(2*g*(k-1)*h))/2; x(k)=x(k-1)+c*v*h/sqrt(1.0-c*c*v*v); end plot(x,-(0.625:.625:10),*r) n=16 （x,-(0.3125:.3125:10),’*r’）n=32 科学出版社 南通大学理学院计算科学与统计学系 * * 运行后所得图形如下： 科学出版社 南通大学理学院计算科学与统计学系 * * 用Mathematica编程得出的结果如下 n=16时，T=1.84534 n=32时，T=1.8446 n=100时，T=1.84429 上面我们用计算机模拟实验的方法得到了最速降线的近似形状， 曲线方程是什么呢？下面我们就来推导这一方程。 越大，步长越小，曲线越接近最速降线，那么最速降线的精确 科学出版社 南通大学理学院计算科学与统计学系 * * 将 代入上式得： 于是质点滑行时间应表为 的泛函 A B x y 任一点处的速度 满足（ 为弧长） 为了使问题简化，将 点取为坐标原点，在如图所示的坐 标系中 点坐标为 根据能量守恒定律，质点在曲线 上 四、理 论 推 导 用变分法求解最速降线的精确方程 科学出版社 南通大学理学院计算科学与统计学系 * * 端点条件为 最速降线满足欧拉方程 等价于 作一次积分得 只依赖于 ，而 不含 ，欧拉方程可展开为 科学出版社 南通大学理学院计算科学与统计学系 * * 即 化为参数方程可表为 这是旋轮线方程，由 可知 ， 由另一端 点条件 确定。 为任意常数，将 代入得 科学出版社 南通大学理学院计算科学与统计学系 * * 按上述参数方程作出的旋轮线经过原点A（0，0）， 求出适当的 B（a,h）时解方程组 用Mathematica语句 求出的 t=2.412, c1=5.72917 用Plot[]语句画出左边函数图象：发现它与横坐标的交点 在2.5左右 当a=10,h=10时，两式相除消去c1，整理得t的一元 方程 此时旋轮线方程为： 科学出版社 南通大学理学院计算科学与统计学系 * * 五、总 结 （1）时间的比较 n 16 18 32 100 极值法 1.85125 1.85681 1.9354 8 —— 模拟法 1.84534 —— 1.8446 1.84429 理论推导 1.84422 科学出版社 南通大学理学院计算科学与统计学系 * * （2）将前面的用求多元函数最小值的方法与模拟出的最速降线与 求极小值 旋轮线 旋轮线进行比较如图所示（红色的表示模拟的，绿色表示 用多元函数求最小值的方法，蓝色表示旋轮线） n=100时模拟出的曲线 科学出版