数值分析实验报告一

实验报告课程名称：数值分析实验项目：解方程/解方程组专业班级：机械工程姓名：李奇学号：201820065实验室号：实验组号：实验时间：批阅时间：指导教师：曲绍波成绩：1沈阳工业大学实验报告（适用计算机程序设计类）专业班级：机械工程学号：201820065姓名：李奇实验名称：解方程1.实验目的：掌握牛顿法和割线法求解方程能够使用MATLAB编程求解2.实验内容：(1)学习掌握牛顿法，割线法原理。(2)学习MATLAB语言的语法格式。(3)编写MATLAB程序,求解方程cos0xxxe的最小正根，取初值00x，当6110kkxx时，迭代终止。3.实验方案（程序设计说明）牛顿法是一种线性化的近似方法，其基本思想是将非线性方程转化为线性方程来进行迭代求解，当计算结果满足精度时，计算中止。割线法具有超线性的收敛速度，求解速度快。根据所掌握知识将牛顿法和割线法用计算机语言表达出来，并运行求解，比较两种求解方法的优劣。4.实验步骤或程序（经调试后正确的源程序）见附件A5．程序运行结果见附件A6．出现的问题及解决方法出现的问题：在编辑器内编写程序时出现“字符数无效”的警告。解决方法：将输入法切换至英文模式，重新输入。附件A沈阳工业大学实验报告（适用计算机程序设计类）专业班级：机械工程学号：201820065姓名：李奇实验步骤或程序：1、编写牛顿法M文件2、调用newton.m文件，输入求解13、运行结果如下4、编写割线法M文件（1）子程序（2）主程序5、调用secant.m文件，输入方程26、运行结果如下沈阳工业大学实验报告（适用计算机程序设计类）专业班级：机械工程学号：201820065姓名：李奇实验名称：解方程组1.实验目的：掌握高斯消去法，列主元高斯消去法和LU分解法求解方程组能够使用MATLAB编程2.实验内容：(1)学习掌握高斯消元，列主元，LU分解原理。(2)学习MATLAB语言的语法格式。(3)编写MATLAB程序求解下列方程组45.35.633.4100.1652.13.24875.72247.252.498.25382.3357.14.2323.72xyzmnxyzmnxyzmnxyzmnxyn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.实验步骤或程序（经调试后正确的源程序）见附件B5．程序运行结果见附件B6．出现的问题及解决方法附件B沈阳工业大学实验报告（适用计算机程序设计类）专业班级：机械工程学号：201820065姓名：李奇实验步骤或程序：1、编写高斯消元法M文件2、调用Gaussxiaoqu.m文件，输入方程组13、运行结果如下24、编写高斯列主元法M文件5、调用Gauss.m文件，输入方程组36、运行结果如下47、编写LU分解法M文件58、调用LU.m文件，输入方程组9、运行结果如下