数值分析实验报告表一.doc

. 实验报告 课程名称： 数 值 分 析 实验项目： 解 方 程 / 解 方 程 组 专业班级： 机械工程 姓 名： 李奇 学 号： 201820065 实验室号： 实验组号： 实验时间： 批阅时间： 指导教师： 曲绍波 成 绩： . . 沈阳工业大学实验报告 （适用计算机程序设计类） 专业班级： 机械工程 学号： 201820065 姓名： 李奇 实验名称：解方程 实验目的： 掌握牛顿法和割线法求解方程 能够使用 MATLAB编程求解 实验内容： 学习掌握牛顿法，割线法原理。 学习 MATLAB语言的语法格式。 (3) 编写 MATLAB程序 , 求解方程 cos x xex 0 的最小正根， 取初值 x0 0 ，当 xk 1xk 10 6 时，迭代终止。 实验方案（程序设计说明） 牛顿法是一种线性化的近似方法，其基本思想是将非线性方程转化为线性方程来进行迭代求解，当计算结果满足精度时，计算中止。 割线法具有超线性的收敛速度，求解速度快。根据所掌握知识将牛顿法和割线法用计算机语言表达出来，并运行求解，比较两种求解方法的优劣。 实验步骤或程序（经调试后正确的源程序） 见附件 A 5．程序运行结果 见附件 A 6．出现的问题及解决方法 出现的问题：在编辑器内编写程序时出现“字符数无效”的警告。 解决方法：将输入法切换至英文模式，重新输入。 . . 附件 A 沈阳工业大学实验报告 （适用计算机程序设计类） 专业班级： 机械工程 学号： 201820065 姓名： 李奇 实验步骤或程序 ： 1、 编写牛顿法 M文件 2、调用 newton.m 文件，输入求解 . . 3、运行结果如下 4、编写割线法 M文件 （ 1）子程序 （2）主程序 5、调用 secant.m 文件，输入方程 . . 6、运行结果如下 . . 沈阳工业大学实验报告 （适用计算机程序设计类） 专业班级： 机械工程 学号： 201820065 姓名： 李奇 实验名称：解方程组 实验目的： 掌握高斯消去法，列主元高斯消去法和 LU分解法求解方程 能 使用 MATLAB 程 实验内容： 学 掌握高斯消元，列主元， LU 分解原理。 学 MATLAB 言的 法格式。 写 MATLAB程序求解下列方程 4 x 5.3y 5.6z 3m 3.4n 100.16 5x 2.1y 3.2z 4m 8n 75.72 2x 4 y 7.2z 5m 2.4n 98.2 5x 3y 8z 2.3m 3n 57.1 4.2x 3y 2n 3.72 实验方案（程序设计说明） 高斯消去法的基本思想是通 消元把 性方程 化 等价的上三角方程 ，再 行求解。高斯消去一般由“消元 程”和“回代 程“两部分 成。消元 程就是按确定的 算 程 方程 的增广矩 行初等 ，将原方程化 与之等价的上三角方程 ；回代 程就是 得到的上三角方程 求解的 程。 列主元高斯消去法的算法步 ： (1)k=1,2, ? ,n-1, 主元素， 出 {|a (k) kk| ， |a (k) k+1，k | , ?, |a (k) nk (k) mk 行和第 m 行交换。 (2) 消元 | } 中 最大的元素 |a | ，将第 k 程。 (3) 回代 程。 (4) 出方程 的解。 LU 分解法的算法步 ： (1) 首先 行 LU 分解， 算 U 的第一行和 L 的第一列。 (2) k=2， 3, ?， n, 算 U 的第 k 行和 L 的第 k 列。 (4) 算 y, 求解 Ly=b。 (4) 求解 x, 求解 Ux=y。 实验步骤或程序（经调试后正确的源程序） 附件 B 5．程序运行结果 附件 B 6．出现的问题及解决方法 . . 附件 B 沈阳工业大学实验报告 （适用计算机程序设计类） 专业班级： 机械工程 学号： 201820065 姓名： 李奇 实验步骤或程序 ： 1、 编写高斯消元法 M文件 2、 调用 Gaussxiaoqu.m 文件，输入方程组 . . 3、运行结果如下 . . 4、 编写高斯列主元法 M文件 5、 调用 Gauss.m文件，输入方程组 . . 6、运行结果如下 . . 7、编写 LU 分解法 M文件 . . 8、调用 LU.m文件，输入方程组 9、运行结果如下 .