数值分析实验报告一(共13页).doc

《数值分析实验报告一(共13页).doc》由会员分享，可在线阅读，更多相关《数值分析实验报告一(共13页).doc（13页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上实验报告课程名称： 数 值 分 析 实验项目： 解 方 程 / 解 方 程 组 专业班级： 机械工程 姓 名： 李奇 学 号： 实验室号： 实验组号： 实验时间： 批阅时间： 指导教师： 曲绍波 成 绩： 沈阳工业大学实验报告（适用计算机程序设计类）专业班级： 机械工程 学号： 姓名： 李奇 实验名称：解方程1.实验目的：掌握牛顿法和割线法求解方程能够使用MATLAB编程求解2.实验内容：(1) 学习掌握牛顿法，割线法原理。(2) 学习MATLAB语言的语法格式。(3) 编写MATLAB程序, 求解方程的最小正根，取初值，当时，迭代终止。3. 实验方案（程序设计说明）

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

3、newton.m文件，输入求解3、运行结果如下4、编写割线法M文件 （1）子程序（2）主程序5、调用secant.m文件，输入方程6、运行结果如下沈阳工业大学实验报告（适用计算机程序设计类）专业班级： 机械工程 学号： 姓名： 李奇 1.实验目的：掌握高斯消去法，列主元高斯消去法和LU分解法求解方程组能够使用MATLAB编程2.实验内容：(1) 学习掌握高斯消元，列主元，LU分解原理。(2) 学习MATLAB语言的语法格式。(3) 编写MATLAB程序求解下列方程组3. 实验方案（程序设计说明）高斯消去法的基本思想是通过消元把线性方程组化为等价的上三角方程组，再进行求解。高斯消去一般由“消元过

4、程”和“回代过程“两部分组成。消元过程就是按确定的计算过程对方程组的增广矩阵进行初等变换，将原方程化为与之等价的上三角方程组；回代过程就是对得到的上三角方程组求解的过程。列主元高斯消去法的算法步骤：(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. 实验步骤或程序（经调试后正确的源程序）见附件B5程序运行结果见附件B6出现的问题及解决方法实验名称：解方程组附件B 沈阳工业大学实验报告（适用计算机程序设计类）专业班级： 机械工程 学号： 姓名： 李奇 实验步骤或程序：1、 编写高斯消元法M文件2、 调用Gaussxiaoqu.m文件，输入方程组3、运行结果如下 4、 编写高斯列主元法M文件5、 调用Gauss.m文件，输入方程组6、运行结果如下7、编写LU分解法M文件8、调用LU.m文件，输入方程组9、运行结果如下