第4节 迭代法.pptVIP

第三章 线性方程组求解的数值方法 迭代法（iterative method ） Iterative method: attempts to solve a problem by finding successive approximations to the solution starting from an initial guess；（comes from wikipedia） 迭代法是从初始猜测值开始，通过依次逼近获得问题解的方法。 直接法：是利用有限次操作得到问题解的方法。 应用： 线性方程组求解、非线性方程求解、非线性方程组求解、特征值求解、最优化方法、分形… … 典型方法： 数值计算：Jacobi法、Gauss-Seidel法、Newton迭代法、最快下降法（Steepest Descent） 数字信号处理：维纳滤波、卡拉曼滤波… … 迭代法解方程 迭代序列的基本公式： 迭代法解方程 迭代函数的构造 例：用迭代法求解方程 2 x = 1 迭代法程序 分形 fractal 分形 fractal 分形艺术 迭代法求解线性方程组 迭代公式的构建： 将方程Ax=b改写为：x=Mx+c，M称为迭代矩阵。 迭代法求解线性方程组 迭代法求解线性方程组 迭代法求解线性方程组 迭代法求解线性方程组 迭代法求解线性方程组 迭代法求解线性方程组 迭代法示例 迭代法示例 作业： P. 66~70 思考题：1，2，4，8 习题：4，6，7 实验题：4 * 课程回顾 向量范数与矩阵范数： 向量范数的定义 1、2、无穷向量范数 矩阵范数： Frobenius 范数， 算子范数：1、2、无穷算子范数 线性算子的矩阵表示 范数的性质：相容性，范数和谱半径的关系，… 范数的应用——分析病态问题 病态问题的判定——条件数。 第四节 迭代法 f（x）为：任意n维空间到n维空间的函数。 如何利用迭代法解方程？ 迭代法基本原理：如果迭代序列｛x（k+1）= f（ x（k） ）｝收敛，则其极限点为方程 f（x）= x 的解。 x*称为函数 f（x）的不动点（fixed point） f（x）可为线性方程（组）、非线性方程（组）。 迭代法可用于各类方程（组）求解问题。 证明： 步骤二、构造迭代函数： 下列方程（组）等价。 步骤一、写出等价方程： clear all close all clc format long % 迭代法求解2x=1 x0 = 0 % 迭代初值 N = 30 % 迭代次数 % 方法1：f(x) = 1 - x x1(1) = x0; for iii = 2 : N x1(iii) = 1 - x1(iii - 1); end figure; plot(x1, o-) title(f(x) = 1 - x) 振荡序列 收敛序列 发散序列 收敛序列 迭代法的核心问题： 迭代函数收不收敛 收敛速度快不快 A fractal is a rough or fragmented geometric shape that can be split into parts, each of which is (at least approximately) a reduced-size copy of the whole, a property called self-similarity 分形是一种粗糙的或分裂的几何形状，其每一部分都是整体的缩小比例的复制（自相似性） 分形是具有自相似性的几何形状。 NOVA纪录片： hunting the hidden dimension “寻找隐藏的维度” 4’30‘’， 26’30’‘ ，30’00’‘ 分形应用： 计算机虚拟现实技术；天线设计；医学； 科赫雪花（Koch snowflake）： 1、初始为正三角形； 2、将每条边等分三段； 3、以中段为底向外做等边三角形 4、移除被三等分的边的中段。 5、如此迭代。 曼德勃罗集（Mandelbrot） 迭代公式一：最简单的构造 迭代公式二：最有效的构造 迭代公式三（Jacobi 迭代）： 迭代公式三（Gauss-Seidel迭代）： 一般认为新近似解要比老近似解更接近真实解，将已计算出的x（k+1）分量替换Jacobi 迭代公式中x（k）相应分量即可得到Gauss-Seidel迭代。 迭代公式三（Gauss-Seidel迭代）(续)： 迭代公式四（超松弛迭代 / SOR）： 相当于在高斯迭代法的基础上再向前（ω1）或向后（ω1）移动一段距离。 迭代公式四（超松弛迭代 / SOR）： 不同迭代算法收敛性不同。 *