数值计算方法与算法.ppt

例题10雅可比迭代求解k01220x10?0.20.20.0240893x200?0.566667?0.578343x3022.11.6987401Jacobi迭代 收敛条件:02Gauss-Seidel迭代收敛条件:讨论方程组 的迭代收敛条件。例题11原理：若模最大特征值λ是实根且唯一，则

x(k+1)=Ax(k)/||Ax(k)||收敛到λ的特征向量。运算量：O(n2log(1/ε))幂法原理：对A-1应用幂法，计算A的模最小的特征值(向量)。运算量：O(n3+n2log(1/ε))反幂法第7章矩阵的特征值和特征向量数值计算方法与算法误差范数第0章绪论例题1计算矩阵 的范数 计算矩阵 的谱半径 Lagrange插值Newton插值多项式插值第1章插值插值函数误差Runge现象差商第1章插值三次样条函数01M关系式02m关系式03第1章插值第1章插值例题2给定f(-1)=3,f(2)=5,f(3)=7,f(4)=5，作出差商表，写出Newton插值公式。x-1234Δ03575Δ12/32-2Δ21/3-2Δ3-7/15例题3给定数据f(0),f(1),f(3),f(3),作出插值多项式，并写出插值余项。第2章数值微分和数值积分数值微分第2章数值微分和数值积分数值积分1复化数值积分2Romberg积分第2章数值微分和数值积分例题4构造积分 的数值积分公式I(f)=a0f(-h)+a1f(0)+a2f(2h)。例题5logo用Romberg公式计算积分：Ri,j01205/2119/87/3275/327/37/3最小二乘原理设A为m×n列满秩实矩阵，b为列向量，则 使函数 取最小值。A多项式拟合B第3章曲线拟合的最小二乘法例题6logo给定数据用最小二乘法求形如y=aebx的经验公式。x-0.70-0.500.250.75y0.991.212.574.23对分法求实根原理：若f(a)f(b)0则f在[a,b]中至少有一零点。收敛步数：O(log(1/ε))迭代法原理：若|f(x)|≤L1，则xk+1=f(xk)收敛。收敛步数：O(log(1/ε))第4章非线性方程求根原理：f(x)≈f(xk)+f(xk)(x-xk)公式：xk+1=xk－f(xk)／f(xk)收敛步数：O(loglog(1/ε))Newton迭代法原理：f(x)≈f(xk)+(x-xk)(f(xk)-f(xk-1))／(xk-xk-1)收敛步数：O(loglog(1/ε))弦截法第4章非线性方程求根02Weierstrass-Durand-Kerner方法：同时逼近一元多项式的所有根原理：F(X)≈F(Xk)+J(Xk)(X－Xk)公式：Xk+1=Xk－J(Xk)-1F(Xk)收敛步数：O(loglog(1/ε))方程组的Newton迭代法01第4章非线性方程求根例题7logo用Newton迭代法求解非线性方程组取(x0,y0)=(0.8,0.6)，误差控制量10-3。k0123xk0.80.8270490.8260320.826031yk0.60.5639340.5636240.563624原理：通过初等行变换，化方程组Ax=b为上三角。1由此得到Dolittle分解和Courant分解A=LU。2条件：A的顺序主子式非零。3运算量：O(n3)4Gauss消元法第5章解线性方程组的直接法A=PLU，P为置换方阵，L中各元素的模都≤1。条件：A行列式非零。运算量：O(n3)原理：在Gauss消元过程中，先选取一列中模最大元素， 将其换行至左上角位置，再作消元。由此得到分解列主元消元法第5章解线性方程组的直接法条件：A行列式非零。原理：在Gauss消元过程中，先选取所有元素模最大者， 将其换行至左上角位置，再作消元。由此得到分解A=PLUQ，P和Q为置换方阵，L中各元素的模