1 / 24 首次试做实验记录实验中心（室）：年月日实验课程名称数值分析面向专业信息与计算科学总学时数实验项目名称线性方程组的直接解法实验学时一、实验目的、要求目的：掌握解线性方程组直接法（特别是顺序高斯消去法）的基本思想，熟悉其算法，加强编程能力和编程技巧，练习从数值分析的角度看问题。要求： 针对给定的实验题目，根据所学的算法，能够熟练地使用某种语言上机编程，给出实验结果，注意上机编程的正确性。二、实验原理用高斯消去法解线性方程组的基本思想是设法消去方程组的系数矩阵A的主对角线下的元素， 而将线性方程组化为等价的上三角形方程组，然后再通过回代过程便可获得原方程组的解。三、使用仪器、材料计算机一台， Matlab 、C、Mathematica等软件。四、实验内容：用高斯消去法求解下列方程组Ax=b, 其中A = 【2 2 2 3 2 4 1 3 9】b = 【1.0000 0.5000 2.5000】五、实验过程原始记录 ( 数据、图表、计算等 ) Matlab 源代码：function x =Gauss_solve(A, b)[n, n] = size(A);for k = 1 : n-1if A(k,k)==0 disp('A is singular,stop.')endfor i = k+1 : n mik = A(i, k) / A(k, k); b(i) = b(i) - mik * b(k);for j=k+1:n2 / 24 A(i, j) = A(i, j) - mik * A(k, j);endendend% back substitutionx=zeros(n,1);x(n)=b(n)/A(n,n);for i=n-1:-1:1 sum=0;for j=i+1:n sum=sum+A(i,j)*x(j);end x(i)=(b(i)-sum)/A(i,i);End 六、实验结果及分析在命令窗口输入相关数据，运行上面定义的函数Gauss_solve>>A=[2 2 2;3 2 4;1 3 9]A = 2 2 2 3 2 4 1 3 9>>b=[1;1/2;5/2]b = 1.0000 0.5000 2.5000>>x = Gauss_solve(A, b)x = -0.5000 1.0000 0结果分析：根据顺序高斯消去法的算法，使用Matlab 编写了程序，对所给的题目进行了求解。利用高等代数中的基本内容，我们可以看出该程序能够成功求解所给的线性方程组。实验中心（室）验收审查意见实验中心（室）主任签字：3 / 24 首次试做实验记录实验中心（室）：年月日实验课程名称数值分析面向专业信息与计算科学总学时数实验项目名称线性方程组的迭代解法实验学时一、实验目的、要求目的： 掌握解线性方程组几种常见的迭代解法的基本思想，熟悉其算法， 练习引入迭代矩阵的形式来解决方程组。要求： 针对给定的实验题目，根据所学的算法，能够熟练地使用某种语言上...
VIP
