《矩阵与数值分析》课程数值实验报告.docxVIP

《矩阵与数值分析》课程数值实验报告 目录 TOC \o 1-5 \h \z 一、 Gauss消去法与Gauss列主元消去法求解方程组 2 二、 用Jacobi和Gauss-Seidel迭代法求解方程组 9 三、 用迭代法求方程的全部根 15 四、 三次样条插值函数的求解 17 五、 用最小二乘法对数据作三次多项式拟合 26 六、 复化求积公式求积分 30 七、 4阶Runge-Kutta法求解微分方程 33 一、Gauss消去法与Gauss列主元消去法求解方程组 1?1 ?问题 给定n阶方程组Ax = b,其中 6 1 、 厂7、 15 8 6 1 A二 ■ ? ? ? ? ? ? ? ? ,h = ? ■ ? 8 6 1 15 8 6丿 14, 则方程组有解x = (l,1)丁。对斤=10和〃 = 84,分别用Gauss消去法和列主元消去法解 方程组，并比较计算结果。 12算法描述 1.2.1算法基本思路 Gauss消去法基本思路：设有方程组Ax = b，设A是可逆矩阵。高斯消去法的基本思想 就是将矩阵的初等行变换作用于方程组的增广矩阵(A|b),将其化成一个上三角形矩阵，然 后通过冋代求这个上三角形方程组的解，就是原方程组的解。 Gauss列主元消去法基本思路：设有方程组Ax = hf设A是可逆矩阵。列主元高斯消去 法的基本思想就是在Gauss消去法屮增加选主元的过程，即在第k步消元时，首先在笫k 列主对角线以下(含主对角元)元素中挑选绝对值最大的数，并通过初等行变换，使得该数 位于主对角线上，然后再继续消元。 1.2.2算法实现步骤： 初始化矩阵A和b； 选择消元方法：1. Gauss消去法2. Gauss列主元消去法 当选择1时，进行Gauss消去：进行次消元，并实时输出每次消元结果，以便观测并 检查每次消元结果；消元后得到上三角形矩阵，回代求解； 当选择2时，进行Gauss列主元消去：进行n?l次消元，在每次消元之前查找列主元，并使 列主元位于对角线上，然后再继续消元并实时输出每次消元结果，以便观测并检查每次消元 结果； 输出最后的求解结果。 1.3程序代码 木程序是用C++语言在Microsoft visual C++编译环境下编译和运行的。由于代码太长，程序 代码见附录1 1.4计算结果： 1)Gauss消去法求解结果： 当矩阵阶数1尸10,即附录1程序代码屮的N设为10时，运行程序得到如下结果： 消去之前的a和b： 6 1 0 0 0 0 0 0 0 0 1 7 8 6 1 0 0 0 0 0 0 0 1 15 0 8 6 1 0 0 0 0 0 0 1 15 0 0 8 6 1 0 0 0 0 0 1 15 0 0 0 8 6 1 0 0 0 0 1 15 0 0 0 0 8 6 1 0 0 0 1 15 0 0 0 0 0 8 6 1 0 0 1 15  0 0 0 0 0 0 8 6 1 0 1 15 0 0 0 0 0 0 0 8 6 1 1 15 0 0 0 0 0 0 0 0 8 6 I 14 请选择消元方法： 1. Gauss消去法 2. Gauss列主元消去法 1\ （由键盘输入数字1,并冋车） 第1次消元结果: 6 1 0 0 0 0 0 0 0 0 1 7 0 4. 6667 1 0 0 0 0 0 0 0 I 5. 6667 0 8 6 1 0 0 0 0 0 0 I 15 0 0 8 6 1 0 0 0 0 0 1 15 0 0 0 8 6 1 0 0 0 0 1 15 0 0 0 0 8 6 1 0 0 0 1 15 0 0 0 0 0 8 6 1 0 0 1 15 0 0 0 0 0 0 8 6 1 0 1 15 0 0 0 0 0 0 0 8 6 1 15 0 0 0 0 0 0 0 0 8 6 I 14 第9次消元结果: 6 1 0 0 0 0 0 0 0 0 1 7 0 4. 6667 1 0 0 0 0 0 0 0 15. 6667 0 0 4. 2857 1 0 0 0 0 0 0 |5. 2857 0 0 0 4. 1333 1 0 0 0 0 0 |5. 1333 0 0 0 0 4. 0645 1 0 0 0 0 | 5. 0645 0 0 0 0 0 4. 0317 1 0 0 0 |5. 0317 0 0 0 0 0 0 4.0157 1 0 0 15.0157 0 0 0 0 0 0 0 4. 0078 1 0 15. 0078 0 0 0 0 0 0 0 0 4. 0039 1 15. 0039 0 0 0 0 0 0 0 0 0 4. 002 | 4. 002 求解结果： X[10]二{1 1 1 1 1 1 1 1 1 1 } 2?当矩阵阶数n=84,