数值分析编程作业.docxVIP

《数值分析》实验报告第二章：解线性方程组的直接方法2、试用MATLAB软件编程实现追赶法求解三对角方程组的算法，并考虑梯形电阻电路问题，电路如下：其中电路中的各个电流{，，…，}须满足下列线性方程组：设，，运用求各段电路的电流量。解：上述方程组可用矩阵表示为：MatLab程序： %赋初值；a=[0 -2 -2 -2 -2 -2 -2 -2];b=[2 5 5 5 5 5 5 5];c=[-2 -2 -2 -2 -2 -2 -2];d=[8.1481 0 0 0 0 0 0 0];%三对角方程的追赶法for i=2:8%“追”的过程；a(i)=a(i)/b(i-1);b(i)=b(i)-c(i-1)*a(i);d(i)=d(i)-a(i)*d(i-1);end;d(8)=d(8)/b(8);%“赶”的过程；for i=7:-1:1d(i)=(d(i)-c(i)*d(i+1))/b(i);end; x=d; x程序运行结果：x = 8.1477 4.0737 2.0365 1.0175 0.5073 0.2506 0.1194 0.0477即。1、试分别用（1）Jacobi迭代法；（2）Gauss-Seidel迭代法；（3）共轭梯度法解线性方程组迭代初始向量取=0，0，0，0，0。解：实验步骤及程序、结果取要求达到的精度。以下程序中的均表示迭代次数。（1）Jacobi迭代法MatLab源程序。format longA=[10,1,2,3,4;1,9,-1,2,-3;2,-1,7,3,-5;3,2,3,12,-1;4,-3,-5,-1,15];b=[12,-27,14,-17,12]; x0=[0,0,0,0,0];x1=x0;Nmax=1000;k=0;for i=1:5 sum=0;for j=1:5 if j~=isum=sum+A(i,j)*x0(j); end; end; x1(i)=(b(i)-sum)/A(i,i); end; while abs(norm(x1-x0,inf))1e-8 kNmaxx0=x1;for i=1:5 sum=0; for j=1:5 if j~=isum=sum+A(i,j)*x0(j); end; end; x1(i)=(b(i)-sum)/A(i,i); end; k=k+1; end;x1k输出结果：x1 = 1.000000024731627 -2.000000022663167 2.999999958446139 -1.999999993335823 0.999999970325168k =88迭代了88次。（2）Gauss-Seidel迭代法MatLab源程序。 A=[10,1,2,3,4;1,9,-1,2,-3;2,-1,7,3,-5;3,2,3,12,-1;4,-3,-5,-1,15];b=[12,-27,14,-17,12]; x0=[0,0,0,0,0];x1=x0; Nmax=1000;k=1;sum=0; for j=2:5 sum=sum+A(1,j)*x0(j); end; x1(1)=(b(1)-sum)/A(1,1);for i=2:4 sum=0;for j=1:(i-1) sum=sum+A(i,j)*x0(j); end;x1(i)=b(i)-sum;sum=0;for j=(i+1):5 sum=sum+A(i,j)*x0(j);end;x1(i)=(x1(i)-sum)/A(i,i); end; sum=0;for j=1:4 sum=sum+A(5,j)*x0(j); end; x1(5)=(b(5)-sum)/A(5,5); while abs(norm(x1-x0,inf))1e-8 kNmaxx0=x1;sum=0; for j=2:5 sum=sum+A(1,j)*x0(j)