2.1误差知识与算法知识要点.pptVIP

本节作业： 1、按照算法2.1编写列主元Gauss消去法的Matlab 程序。 2、分别用上述程序和求解线性方程组的Matlab命令（函数）： x=A/b或x=inv（A）*b 计算下列方程组： 小 结 ★基本思想： 消元与回代 ★顺序Gauss消去法能进行到底的条件 (1)主元 ★顺序Gauss消去法的缺点 Δ没有很好的数值稳定性。 Δ当A可逆时，AX=b有唯一解，但顺序Gauss消去 法不一定能进行到底。 1、顺序Gauss消去法 （2）矩阵A的前n-1个顺序主子式非零。 2、列主元Gauss消去法 ★列主元Gauss消去法的优点 Δ具有很好的数值稳定性； Δ当A可逆时，列主元Gauss消去法一定能进行到底。 ★基本思想：避免接近于零的数作分母。 ★列主元Gauss消去法能进行到底的条件 Δ具有与顺序Gauss消去法相同的计算量。 * 2002年NEC公司为日本地球模拟中心建造的一台“地球模拟器”，每秒能进行的浮点运算次数接近36万亿次，堪称目前的超级运算冠军。Blue Gene/L是一部由3.3万个处理器构成的超级计算机。将是一台拥有13万个处理器的系统，估计具有360万亿次的峰值速度。 * 一万万是一亿。一万亿是一兆。一万兆是一京。一万京是一恒河沙。 第一步： 第一步： 第一步： 第一步： 假设 第一步: 第二步: 第k步: 第一步: 第二步: 第k步: 第一步: 第二步: 第k步: 第一步: 第二步: 第k步: 第一步: 第二步: 第k步: ★顺序Gauss消去法的计算公式： 消元过程 对于 k=1,2,···,n-1 执行 （1）如果 则算法失效， （2）对于i=k+1,k+2,···n,计算 停止计算；否则转（2）; ★顺序Gauss消去法的计算公式： 回代过程 k=n-1,n-2,···,1 ★顺序Gauss消去法的计算公式： 消元过程 对于 k=1,2,···,n-1 执行 （1）如果 则算法失效， （2）对于i=k+1,k+2,···n,计算 停止计算；否则转（2）; 回代过程 k=n-1,n-2,···,1 ★顺序Gauss消去法的计算量： 消元过程 对于 k=1,2,···,n-1 执行 （1）如果 则算法失效， （2）对于i=k+1,k+2,···n,计算 停止计算；否则转（2）; ★顺序Gauss消去法 的计算量 消元过程的计算量 第k步(k=1,2,…,n-1) 需n-k次除法; 需 次乘法; 需 n-k 次乘法 ★顺序Gauss消去法的计算量 消元过程的计算量 第k步(k=1,2,…,n-1) 计算: 需n-k次除法; 需 次乘法; 需 n-k 次乘法 消元过程中乘除法总次数: ★顺序Gauss消去法的计算量 消元过程的计算量 第k步(k=1,2,…,n-1) 需n-k次除法; 需 次乘法; 需 n-k 次乘法 消元过程中乘除法总次数: 回代过程的计算量 回代过程的计算量 ★计算AX=b,顺序Gauss消去法的乘除次数 当n=10时,需430次乘除法; 当n=20时,需3060次乘除法; 消元过程中乘除法总次数 ★顺序Gauss消去法能进行到底的条件 (1)主元 证明关键： (2) ★顺序Gauss消去法的缺点 Δ没有很好的数值稳定性。 Δ当A可逆时，AX=b有唯一解，但顺序Gauss消去 法不一定能进行到底。 例：在四位十进制的限制下用 顺序Gauss消去法求解方程组 (2)列主元Gauss消去 法求解 消元: 回代: 精确解: (1)利用顺序Gauss消去法求解 交换两个方程: 消元 回代 解 (1)消元 (2)回代 精确解： 例1 用顺序Gauss消去法求解 用列主元素Gauss消去法解例1 消元过程 回代得 精确解： 3、列主元Gauss消去法 ★列主元Gauss消去法算法： 记 消元过程 ★列主元Gauss消去法算法： 记 消元过程 回代过程 算法2.1 （列主元Gauss消去法） （1）输入数据A和b，置det=1. （2）对于k=1，2，…，n-1，做 （a）按列选主元.令 如果 则停止计算（此时 ） （b）交换两行. 如果 ，则 （2）对于k=1，2，…，n-1，做 （a）按列选主元.令 （b）交换两行. （c）消元计算.置 （d）置 （3）如果 则停止计算（此时 ）； 否则进行回代， 对于 置 （4）置 （5）输出 为线性方程组 的