MATLAB算方法3解线性方程组计算解法.pptxVIP

解线性方程组 第三章线性方程组数值解法 §3.1 直接法 一、 Gauss　消去法 设 有 消 元： % 回代 x=zeros(n,1); x(n)=a(n,n+1)/a(n,n); for k=n-1:-1:1 x(k,:)=(a(k,n+1)-a(k,(k+1):n)*x((k+1):n))/a(k,k); end 二 列主元素Gauss消去法---计算结果可靠 到此原方程组化为 (上三角方程组) （3.2) 直到(n-1) 原方程组化为 以上为消元过程。 （3.3) (n) 回代求解公式 系数矩阵为对称正定阵或严格对角占优阵的方程组按高斯消去法计算是数值稳定的，因而不必选主元。 严格对角占优阵：至少有一个主对角线元素的绝对值严格大于此行或此列其他元素的绝对值之和。 Gauss 列主元消去法： 优点 ------ 计算结果更可靠； 缺点 ------ 挑主元花机时更多, 次序有变动，程序复杂。 说明： (1)也可采用无回代的列主元消去法(叫Gauss- --Jordan消去法)，但比有回代的列主元消 去法的乘除运算次数多。 (2)有回代的列主元消去法所进行的乘除运算 次数为 ，量很小。 if flag==0,a,end end %回代 x=zeros(n,1); x(n)=a(n,n+1)/a(n,n); for k=n-1:-1:1 x(k,:)=(a(k,n+1)-a(k,(k+1):n)*x((k+1):n))/a(k,k); end 记 ，则 （三角因子分解） Gauss消元，初等行变换，化原方程组为上三角型。 三．矩阵三角分解法 定义3.1 叫 的三角（因子）分解，其中 是 是上三角。 下三角, 为单位下三角阵（对角元全为1）， 为上三角阵，则称 为Doolittle分解; 若 是下三角， 是单位上三角，则称 定义3.2 若 为Crout分解。 定义3.2 若 为什么要讨论三角分解？若在消元法进行前能实 现三角分解 ， 则 从而容易回代求解。 1．直接三角分解法(Doolittle分解为例） 　　 　 由矩阵乘法 ……………………… ……. (k) 例3．1 用Matlab实现LU分解 在Matlab程序编辑器中输入： function [L,U]=nalu(a) % a为可逆方阵；L返回单位下三角矩阵；U返回上三角矩阵 n=length(a); U=zeros(n,n);L=eye(n,n); U(1,:)=a(1,:);L(2:n,1)=a(2:n,1)/U(1,1); for k=2:n U(k,k:n)=a(k,k:n)-L(k,1:k-1)*U(1:k-1,k:n); L(k+1:n,k)=(a(k+1:n,k)-L(k+1:n,1:k-1)*U(1:k-1,k))/U(k,k); end 2．平方根法 定理3.1 设A对称正定，则有非奇异下三角阵L，使 分解方法：设 3.列主元的三角分解法 例4 对下面的增广矩阵用列主元三角分解法分解。 四、解三对角方程组——追赶法 利用Gauss消元法得到同解的三角方程为 , 故有 (3.1) 解 解 (3.1) —(3.3) 叫追赶法，工作量小， 非常有效。 Ex:第三章 4 (3.2) （3.3） % 回代 x(n)=a(n,n+1)/a(n,n); for k=n-1:-1:1 x(k)=(a(k,n+1)-a(k,k+1)*x(k+1))/a(k,k) end end a end for i=1:(n-1) a(n,n)=a(n,n)-a(i,n)*a(n,i); a(n,n+1)=a(n,n+1)-a(i,n+1)*a(n,i); end a x=zeros(n,1); x(n)=a(n,n+1)/a(n,n); for k=n-1:-1:1 x(k)=(a(k,n+1)-a(k,(k+1):n)*x((k+1):n))/a(k,k); end 基本要求 1.熟悉Gauss消去法及列主元Gauss消去法。 2.熟悉矩阵的直接三角分解法的过程； 2.熟悉平方根法的计算过程及其优缺点。 作业： 第三章 9，11 §3.2 向量范数和矩阵范数 一.简单迭代法 1.迭代法建立. 考虑 (矩阵B