数值分析实验报告一

实验报告实验报告 课程名称 课程名称 数数 值值 分分 析析 实验项目 实验项目 解解 方方 程程 解解 方方 程程 组组 专业班级 专业班级 机械工程机械工程 姓姓 名 名 李奇李奇 学学 号 号 实验室号 实验室号 实验组号 实验组号 实验时间 实验时间 批阅时间 批阅时间 指导教师 指导教师 曲绍波曲绍波 成成 绩 绩 1 沈阳工业大学实验报告沈阳工业大学实验报告 适用计算机程序设计类 专业班级 机械工程 学号 姓名 李奇 实验名称 解方程 1 实验目的 掌握牛顿法和割线法求解方程 能够使用 MATLAB 编程求解 2 实验内容 1 学习掌握牛顿法 割线法原理 2 学习 MATLAB 语言的语法格式 3 编写 MATLAB 程序 求解方程的最小正根 cos0 x xxe 取初值 当时 迭代终止 0 0 x 6 1 10 kk xx 3 实验方案 程序设计说明 牛顿法是一种线性化的近似方法 其基本思想是将非线性方程转化 为线性方程来进行迭代求解 当计算结果满足精度时 计算中止 割线法具有超线性的收敛速度 求解速度快 根据所掌握知识将牛顿 法和割线法用计算机语言表达出来 并运行求解 比较两种求解方法的优 劣 4 实验步骤或程序 经调试后正确的源程序 见附件 A 5 程序运行结果 见附件 A 6 出现的问题及解决方法 出现的问题 在编辑器内编写程序时出现 字符数无效 的警告 解决方法 将输入法切换至英文模式 重新输入 附件附件 A 沈阳工业大学实验报告沈阳工业大学实验报告 适用计算机程序设计类 专业班级 机械工程 学号 姓名 李奇 实验步骤或程序实验步骤或程序 1 编写牛顿法 M 文件 2 调用 newton m 文件 输入求解 1 3 运行结果如下 4 编写割线法 M 文件 1 子程序 2 主程序 5 调用 secant m 文件 输入方程 2 6 运行结果如下 沈阳工业大学实验报告沈阳工业大学实验报告 适用计算机程序设计类 专业班级 机械工程 学号 姓名 李奇 1 实验名称 解方程组 1 实验目的 掌握高斯消去法 列主元高斯消去法和 LU 分解法求解方程组 能够使用 MATLAB 编程 2 实验内容 1 学习掌握高斯消元 列主元 LU 分解原理 2 学习 MATLAB 语言的语法格式 3 编写 MATLAB 程序求解下列方程组 45 35 633 4100 16 52 13 24875 72 247 252 498 2 5382 3357 1 4 2323 72 xyzmn xyzmn xyzmn xyzmn xyn 3 实验方案 程序设计说明 高斯消去法的基本思想是通过消元把线性方程组化为等价的上三角方程组 再 进行求解 高斯消去一般由 消元过程 和 回代过程 两部分组成 消元过程 就是按确定的计算过程对方程组的增广矩阵进行初等变换 将原方程化为与之等 价的上三角方程组 回代过程就是对得到的上三角方程组求解的过程 列主元高斯消去法的算法步骤 1 对 k 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 实验步骤或程序 经调试后正确的源程序 见附件 B 5 程序运行结果 见附件 B 6 出现的问题及解决方法 附件附件 B 沈阳工业大学实验报告沈阳工业大学实验报告 适用计算机程序设计类 专业班级 机械工程 学号 姓名 李奇 实验步骤或程序实验步骤或程序 1 编写高斯消元法 M 文件 2 调用 Gaussxiaoqu m 文件 输入方程组 1 3 运行结果如下 2 4 编写高斯列主元法