1、1.列主元法（178頁實習1題）函數(shù)liezhu(A,b).m的頭文件的程序如下：n=length(b); detA=det(A); d=0; x=zeros(n,1); c=zeros(1,n); if abs(det(A)=eps error(系數(shù)矩陣是奇異的); %所求方程組系數(shù)矩陣必須是是非奇異的！ return; end %按列選主元 for i=1:n-1 max=abs(A(i,i); m=i; for j=i+1:n if max A=10,-7,0,1;-3,2.,6,2;5,-1,5,-1;2,1,0,2; b=8;5.;5;1; x,detA=liezhu(A,b)x =

2、 0.0000 -1.0000 1.0000 1.0000detA = -762.00012.LU分解法：（178頁實習1題）LUFJ.m文件中函數(shù)LUFJ(A,b)的程序如下： function LUFJ(A,b) %A為系數(shù)矩陣，b為右端項矩陣%UNTITLED2 Summary of this function goes here% Detailed explanation goes here m,n=size(A); %初始化矩陣A，b，L和Un=length(b); L=eye(n,n);U=zeros(n,n);U(1,1:n)=A(1,1:n); %開始進行LU分解L(2:n,1

3、)=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 L %輸出L矩陣 U %輸出U矩陣 y=zeros(n,1); %開始解方程組Ux=y y(1)=b(1);for k=2:n y(k)=b(k)-L(k,1:k-1)*y(1:k-1);endx=zeros(n,1);x(n)=y(n)/U(n,n);for k=n-1:-1:1 x(k)=(y(k)-U(k,k+1:n)*x(k+

4、1:n)/U(k,k);endfor k=1:n fprintf(x%d=%fn,k,x(k);end end運行結(jié)果如下： A=10,-7,0,1;-3,2.,6,2;5,-1,5,-1;2,1,0,2; b=8;5.;5;1; LUFJ(A,b)L = 1.0e+006 * 0.0000 0 0 0 -0.0000 0.0000 0 0 0.0000 -2.5000 0.0000 0 0.0000 -2.4000 0.0000 0.0000U = 1.0e+007 * 0.0000 -0.0000 0 0.0000 0 -0.0000 0.0000 0.0000 0 0 1.5000 0.

5、5750 0 0 0 0.0000x1=-0.x2=-1.x3=1.x4=1.3.雅可比迭代法：（179頁實習3題）Jacobi.m文件中函數(shù)Jacobi(A,b,eps)的程序如下：function Jacobi(A,b,eps) %A為系數(shù)矩陣，b為后端項矩陣，epe為精度m,n=size(A);D=diag(diag(A); %求矩陣DL=D-tril(A); %求矩陣LU=D-triu(A); %求矩陣Utemp=1;x=zeros(m,1);k=0;while abs(max(x)-temp)eps temp=max(abs(x); k=k+1; %記錄循環(huán)次數(shù) x=-inv(D)*

6、(L+U)*x+inv(D)*b; %雅克比迭代公式end for k=1:n fprintf(x%d=%fn,k,x(k);end運行結(jié)果如下： A=10 -1 2 0;0 8 -1 3;2 -1 10 0;-1 3 -1 11; b=-11;-11;6;25; Jacobi(A,b,0.00005)x1=-1.x2=-2.x3=0.x4=2.4 .Gauss-Seidel迭代程序（179頁實習3題）Gauss-Seidel.m文件中函數(shù)Gauss-Seidel(A,b,eps)的程序如下：function Gauss_Seidel(A,b,eps) %A為系數(shù)矩陣，b為后端項矩陣，epe為精度m,n=size(A);D=diag(diag(A); %求矩陣DL=D-tril(A); %求矩陣LU=D-triu(A); %求矩陣Utemp=1;x=zeros(m,1);k=0;while abs(max(x)-temp)eps temp=max(abs(x); k=k+1; %記錄循環(huán)次數(shù) x=inv(D-L)*U*x+inv(D-L)*b; %Gauss_Seidel的迭代公式end for k=1:n fprintf(x%d=%fn,k,x(k);end其運行結(jié)果如