数分实验一.docVIP

石家庄经济学院 实验报告 (学 院)系: 信息工程学院 专 业: 计算机科学与技术 姓 名: 仇经旭 班 级: 4081090301 学 号: 408109070302 任课教师: 罗 毅 《数值分析》实验报告姓名仇经旭学号408109070302日期2010-10-29 1-2节实验室计算机软件技术实验室(152机房) 设备编号78实验题目实验一 用雅可比迭代和高斯赛德尔迭代解方程组 一 实验目的掌握雅可比迭代和高斯赛德尔迭代的基本思想，编程实现雅可比迭代和高斯赛德尔迭代法求解线性方程组的根。二 实验内容1、用雅可比迭代法求方程组8*x1-3*x2+2*x3=20，4*x1+11*x2-x3=33，6*x1+3*x2+12*x3=36的解，要求精确到小数点后三位。（matlab）1)雅可比迭代方法算法： 1.1）适当提供迭代初值{ （X1）} 1.2）按雅可比公式X（k+1）=D-1（L+U）Xk+D-1b(k=0,1,2…) 将老值{Xi(k)}加工成新值{Xi（k+1）} 1.3) 若迭代偏差ek=max|Xi(k+1)-Xi(k)|e 则输出结果终止计算机，否则执行下一步1.4）若迭代次数k未达到最大迭代次数N 则转为步骤（2） 继续迭代，否则输出 失败，标志终止计算机。2）源程序A=[8 -3 2;4 11 -1;6 3 12];B=[20;33;36];Err_user=0.001;N=10;[m,n]=size(A);X=zeros(n,1);k=1;while k=N Xk=X; for i=1:n for j=1:n if i~=j AX(j)=A(i,j)*Xk(j); end end Sum_AX=sum(AX); AX=0; X(i)=(B(i)-Sum_AX)/A(i,i); end E=max(abs(Xk-X)); if EErr_user break; end k=k+1;enddisp(X);disp(k); 3.0003 1.9999 0.9997 9实验结果讨论和分析 本题使用二分法得到的x=1.61803，满足基本要求，题目要求精确到小数点后一位，告诉了本题二分法得应达到得精确度；二分法收敛性很好，收敛速度不快。2、用高斯赛德尔迭代法求方程组8*x1-3*x2+2*x3=20，4*x1+11*x2-x3=33，6*x1+3*x2+12*x3=36的解，要求精确到小数点后三位。（matlab）1）高斯赛德尔迭代方法算法3.1）input n，A,b,e,N 3.2) x[i]-0(i=1,2,…,n) 3.3)k-1 3.4)while(k=N)do(5~11) 3.5)err-0 3.6)for i=1 to n, do(7~9) 3.7）XI-x[i] 3.8) x[i]-(b[i]-∑A[i,j]*x[j])/A[i,i] 3.9) if(err|XI-x[i]|)then err=|XI-x[i]| 4.0) if(erre) then output (x[i],i=1,2,…,n),stop 4.1) k-k+1 4.2)output(“Maximum number of iterations exceeded.”),stop 2）源程序A=[8 -3 2;4 11 -1;6 3 12];B=[20;33;36];Err_user=0.01;N=10;[m,n]=size(A);X=zeros(n,1);k=1;while k=N Xk=X; for i=1:n for j=1:n if i~=j AX(j)=A(i,j)*X(j); end end Sum_AX=sum(AX); AX=0; X(i)=(B(i)-Sum_AX)/A(i,i); end Er=max(abs(Xk-X)); if ErErr_user b