数值计算方法实验报告_3.docVIP

PAGE PAGE 10 09-10（2）信计08级 数值计算方法课程设计实验报告 姓名： 王超 班级： 信计08-2班 学号： 200811011076 成绩： 实验一 实验名称：非线性方程（组）求解 实验目的：1.掌握二分法，迭代法，牛顿迭代法求方程近似根的基本思想与原理。 2.掌握常用迭代算法的程序实现。 实验题目：求方程f（x）=ｘ＾３－ｓｉｎｘ－１２ｘ＝１的全部实根， ｅｐｓ＝０.０００００１　　 实验过程： 设计思路：用简单迭代法求解。找到收敛的迭代公式，用迭代公式求解问题。 算法描述；由题目f（x）=ｘ＾３－ｓｉｎｘ－１２ｘ＝１得到收敛公式 ｘ（ｋ＋１）＝（ｓｉｎ（ｘｋ）＋１２ｘｋ－１）＾（１／３） 程序清单： function x=roots(x0,eps) k=1 xk=(sin(x0)+12*(x0)-1)^(1/3) while abs(xk-x0)eps k=k+1 p=xk xk=(sin(x0)+12*(x0)-1)^(1/3) x0=p end ４.调试： x=roots(3.1415926,0.000001) k = 1 xk = 3.3232 k = 2 p = 3.3232 xk = 3.3232 x0 = 3.3232 ５.实验小结： 基本掌握了简单迭代法求解非线性方程组的方法，对ｍａｔｌａｐ程序有了初步了解。 实验二 实验名称：解线性方程组的直接法 实验目的：1．掌握高斯（Gauss）消元法及选列主元素的技术 2．掌握追赶法 实验题目：用追赶法求求下述三对角线性方程组的解 实验过程： 设计思路：追赶法求解线性方程组。 算法描述：对三对角线矩阵的LU分解。 程序清单： function x = zgf（a，b，c，d） x = d; n = length(x); for j = 1:n-1 mu = a(j)/b(j); b(j+1) = b(j+1) - mu*c(j); x(j+1) = x(j+1) - mu*x(j); end x(n) = x(n)/b(n); for j = n-1:-1:1 x(j) = (x(j)-c(j)*x(j+1))/b(j); end 4.调试： -0.2667 -0.2000 -0.1333 -0.0667 5.实验小结：掌握了追赶法解决三对角线性方程的基本解法。 实验三 实验名称：解线性方程组的迭代法 实验目的：1.理解迭代法的基本原理及特点。 2.掌握jacobi迭代和Gauss-Seidel迭代格式的分量形式，矩阵形式及其各自特点。 3. 掌握jacobi迭代和Gauss-Seidel迭代算法的程序实现方法。 实验题目：用雅可比迭代法求下述线性方程组的解 要求：选取不同初值x0, 对于给定精度ε=10-6，用三种迭代法计算，观察迭代次数并分析计算结果给出结论。 实验过程： 设计思路：雅克比迭代。Ax=b，A=L+D+U，迭代公式：x=-D-1 （L+U）x+ D-1b 算法描述：Ax=b，A=L+D+U，迭代公式：x=-D-1 （L+U）x+ D-1 程序清单： function x=jacobiff(a,b,eps) n=length(a); d=diag(diag(a)); l=tril(a,-1); u=triu(a,1); x0=diag(zeros(n)); k=1 x=-inv(d)*(l+u)*x0+inv(d)*b while norm(x-x0)eps; k=k+1 p=x; x=-inv(d)*(l+u)*x+inv(d)*b x0=p; end 调试： a=[-3 -1 0 0;-1 2 -1 0;0 -1 2 -1;0 0 -1 2]; b=[1 0 0 0]; jacobiff(a,b,0.000001) k =1 x =-0.3333 0 0 0 k =2 x =-0.3333 -0.1667 0 0 k =3 x =-0.2778 -0.1667 -0.0833 0 k =4 x =-0.2778 -0.1806 -0.0833 -0.0417 k =5 x =-0.2731 -0.1806 -0.1111 -0.0417 k =6 x =-0.2731 -