Ch10[1].数值算法实现.ppt

Ch10. 数值算法实现 §1. 线性方程组解法 1、三角形线性方程组解法 以上三角形线性方程组 为例。 回代： %文件 uptri.m function u = uptri ( a, b ) n= size(a,1) ; x(n)= b(n) / a(n, n) ; for i = n-1:-1:1 s=0 ; for j = i+1: n s=s+a(i, j) * x( j ) ; end x(i) = ( b(i) –s) / a(i, i) ; end u = x ; 2、顺序Gauss消去法 （1）消去过程： 第 k 步，计算 （2）回代过程： %文件 gauss.m function u = gauss (a, b) n = length (b) ; for k=1: n –1 for i = k+1 : n mult = a ( i, k) / a (k, k) ; for j =k +1: n % if abs ( a( k, k) ) 1e–6 a (i, j) = a (i, j)