数学建模数值算法案例.ppt

算法描述 输入平面方程一般式，螺旋线方程一般参数式，用3个欧拉角和映像z轴到螺旋线中心轴的变换向量 把输入平面方程一般式，螺旋线方程参数式 判断方程是否奇异 求极值点 有解否？ 输出无解标志 用二分法求唯一解 求周期解 解在极值附近否？ 用二分法求解 用牛顿法求解 是唯一解? y y y N N N N 二分法有根区间的确定 牛顿法初值的确定：(牛顿法局部收敛，而且会出现加收敛)以连接相邻两个极值点的直线与t轴交点横坐标为初始点。 数值积分和数值微分 1．黎曼和： 2．牛顿—柯特斯(Newton—Cotes)求积公式 将积分区间[a,b]等分n份,记,分点为,k=0,1,...,n, 用拉格朗日插值多项式代替被积函数f(x)，得到n阶Newton—Cotes公式 3.常用的Newton—Cotes公式 n=1时，得梯形公式: I? n=2时，得辛普森公式: I? 4 复化求积公式 随着n的增加可减少积分误差，但高阶Newton—Cotes公式又会造成数值不稳定，因而采用复化公式。 将积分区间[a,b]等分n等份,在每个小区间[xk,xk+1]上使用低阶Newton—Cotes公式, (k=0,1,...,n). 1).复化梯形公式 在每个子区间上用梯形公式,将它们累加得复化梯形公式： 2). 复化辛浦生公式 记xk+1/2是区间的中点, 在每个子区间上用辛