《高斯消元法》课件2.pptVIP

对比得： 计算流程： Ax=f LUx=f 记Ux=y,Ly=f 先由Ly=f求出y,再由Ux=y求出x 追赶法公式： 1.计算的递推公式： 2.解Ly=f： 3.解Ux=y： 由 例P230.8 由公式 解:设 由 得 function x=threedia(a,b,c,f) n=length(f); L=eye(n); x=zeros(1,n); y=zeros(1,n); d=zeros(1,n); u=zeros(1,n); d(1)=b(1); for i=1:n-1 u(i)=c(i)/d(i); d(i+1)=b(i+1)-a(i+1)*u(i); end y(1)=f(1)/d(1); for i=2:n y(i)=(f(i)-a(i)*y(i-1))/d(i); end x(n)=y(n); for i=n-1:-1:1 x(i)=y(i)-u(i)*x(i+1); end 解方程组 a=[0 -1 -1 -3]; b=[2 3 2 5]; c=[-1 -2 -1 0]; f=[6 1 0 1] ‘ threedia(a,b,c,f) X= 5 4 3 2 §4 误差分析 一、矩阵范数 设 ，给