计算方法实习3-4matlab部分2012东南大学.doc

实习题 3 用列主元消去法解方程组： 在GEpiv.m文件中程序段如下： function GEpiv(A,b) [m,n]=size(A); nb=n+1; Ab=[A b]; %........消元a for (i=1:m-1) [pivot,p]=max(abs(Ab(i:n,i))); ip=p+i-1; if(ip~=i) Ab([i ip],:)=Ab([ip i],:); end pivot=Ab(i,i); for (k=i+1:m) Ab(k,i:nb)=Ab(k,i:nb)-(Ab(k,i)/pivot)*Ab(i,i:nb); end end %........回代 x= zeros(n,1); x(n)=Ab(n,nb)/Ab(n,n); for i=n-1:-1:1 x(i)=(Ab(i,nb)-Ab(i,i+1:n)*x(i+1:n,1))/Ab(i,i); end for k=1:n fprintf(x%d=%f\n,k,x(k)); end end 在命令栏中输入程序如下： A=[ 1,1,0,3; 2,1,-1,1; 3,-1,-1,3; -1,2,3,-1 ] A = 1 1 0 3 2 1 -1 1 3 -1 -1 3 -1 2 3 -1 b=([4,1,-3,4]) b = 4 1 -3 4 GEpiv(A,b) x1=-1.333333 x2=2.333333 x3=-0.333333 x4=1.000000 这就是该方程组的解 用LU分解法解方程组Ax=b,其中： A=, b= 在LUpiv.m文件中程序段如下： function LUpiv(A,b) [m,n]=size(A); nb=n+1; Ab=[A b]; Ab(2:m,1)=Ab(2:m,1)/Ab(1,1); for k=2:m for j=k:nb Ab(k,j)=Ab(k,j)-Ab(k,1:k-1)*Ab(1:k-1,j); end for i=k+1:m Ab(i,k)=(Ab(i,k)-Ab(i,1:k-1)*Ab(1:k-1,k))/Ab(k,k); end end x=zeros(n,1); x(n)=Ab(n,nb)/Ab(n,n); for k=n-1:-1:1 x(k)=(Ab(k,nb)-Ab(k,k+1:n)*x(k+1:n,1))/Ab(k,k); end for k=1:n fprintf(x%d=%f\n,k,x(k)); end end 在命令栏中输入程序如下: A=[ 1,2,-12,8; 5,4,7,-2; -3,7,9,5; 6,-12,-8,3] A = 1 2 -12 8 5 4 7 -2 -3 7 9 5 6 -12 -8 3 b=([27,4,11,49]) b = 27 4 11 49 LUpiv(A,b) x1=3.000000 x2=-2.000000 x3=1.000000 x4=5.000000 4编写用追赶法解三对角线性方程组的程序，并解下列方程组： 2）Ax=b,其中： =, 5.分别用雅克比迭代法与高斯-赛德尔迭代法解下列方程组： 1） 雅可比迭代法： 在Jacobi.m文件中程序如下： function Jacobi(A,b,max,eps) %max为迭代最大容许次数 ，eps为容许误差 n=length(A);x=zeros(n,1);