数值分析实验报告一

1、精品资料实验报告课程名称：数值分析实验项目：解方程/解方程组专业班级：机械工程姓名：李奇学号： 201820065实验室号：实验组号：实验时间：批阅时间：指导教师：曲绍波成绩：可编辑修改精品资料沈阳工业大学实验报告（适用计算机程序设计类）专业班级：机械工程学号：201820065姓名：李奇实验名称：解方程可编辑修改精品资料1. 实验目的：掌握牛顿法和割线法求解方程能够使用 MATLAB编程求解2. 实验内容：(1) 学习掌握牛顿法，割线法原理。(2) 学习 MATLAB 语言的语法格式。(3) 编写 MATLAB程序 ,求解方程 cos x xex0 的最小正根，取初值x0 0，当xk 1xk

2、106时，迭代终止。3. 实验方案（程序设计说明）牛顿法是一种线性化的近似方法，其基本思想是将非线性方程转化为线性方程来进行迭代求解，当计算结果满足精度时，计算中止。割线法具有超线性的收敛速度，求解速度快。根据所掌握知识将牛顿法和割线法用计算机语言表达出来，并运行求解，比较两种求解方法的优劣。4. 实验步骤或程序（经调试后正确的源程序）见附件 A5 程序运行结果见附件 A6 出现的问题及解决方法出现的问题：在编辑器内编写程序时出现“字符数无效”的警告。解决方法：将输入法切换至英文模式，重新输入。可编辑修改精品资料可编辑修改精品资料附件 A沈阳工业大学实验报告（适用计算机程序设计类）专业班级：机

3、械工程学号：201820065姓名：李奇实验步骤或程序 ：1、 编写牛顿法 M 文件2 、调用 newton.m文件，输入求解可编辑修改精品资料3 、运行结果如下4 、编写割线法 M 文件（ 1）子程序（2 ）主程序可编辑修改精品资料5 、调用 secant.m文件，输入方程6 、运行结果如下可编辑修改精品资料沈阳工业大学实验报告（适用计算机程序设计类）专业班级：机械工程学号：201820065姓名：李奇可编辑修改1. 实验目的：精品资料掌握高斯消去法，列主元高斯消去法和LU 分解法求解方程组能够使用 MATLAB编程2. 实验内容：(1) 学习掌握高斯消元，列主元， LU 分解原理。(2)

4、学习 MATLAB 语言的语法格式。(3) 编写 MATLAB 程序求解下列方程组4 x5.3y5.6z3m3.4n100.165x2.1y3.2z4m8n75.722x4 y7.2z5m2.4n98.25x3y8z2.3m3n57.14.2x3y2n3.723. 实验方案（程序设计说明）高斯消去法的基本思想是通过消元把线性方程组化为等价的上三角方程组，再进行求解。高斯消去一般由“消元过程”和“回代过程“两部分组成。消元过程就是按确定的计算过程对方程组的增广矩阵进行初等变换，将原方程化为与之等价的上三角方程组；回代过程就是对得到的上三角方程组求解的过程。列主元高斯消去法的算法步骤：(1)对 k

5、 =1,2, ,n-1, 选主元素，选出 |a (k) kk | ，|a (k) k+1 ， k | , ,|a (k) nk |中绝对值最大的元素 |a (k) mk | ，将第 k 行和第 m 行交换。 (2) 消元过程。 (3) 回代过程。 (4) 输出方程组的解。LU 分解法的算法步骤： (1)首先进行 LU 分解，计算 U 的第一行和 L 的第一列。(2)对 k=2 ，3, ，n, 计算 U 的第 k 行和 L 的第 k 列。 (4) 计算 y, 求解 Ly=b 。(4)求解 x, 求解 Ux=y 。4. 实验步骤或程序（经调试后正确的源程序）见附件 B5 程序运行结果可编辑修改见附件 B6 出现的问题及解决方法精品资料实验名称：解方程组可编辑修改精品资料附件 B沈阳工业大学实验报告（适用计算机程序设计类）专业班级：机械工程学号：201820065姓名：李奇实验步骤或程序 ：1、 编写高斯消元法M 文件2、 调用 Gaussxiaoqu.m文件，输入方程组可编辑修改精品资料3 、运行结果如下可编辑修改精品资料4、 编写高斯列主元法M 文件可编辑修改精品资料5、 调用 Gauss.m文件，输入