文档介绍：作彷2丈大拿
课程名称：
数值分析
实验项目：
解方程/解方
程组
专业班级：
机械工程
姓名：一
李奇学号：
201820065
实验室号：_
实验组号：
实验时间：
批阅时间：
指导教帅：_
曲绍波一成绩：
沈阳工业大学实验报告
(适用计算机程序设计类)
专业班级：机械工程学号： 201820065 姓名：李奇
实验名称：解方程
实验目的：
掌握牛顿法和割线法求解方程
能够使用MATLAB^程求解
实验内容：
学****掌握牛顿法，割线法原理。
学****MATLAB!言的语法格式。
x
编写MATLAB序，求解万程cosx-xe =0的最小正根，
取初值X0=°,当|S-对＜〔"时，迭代终止。
实验方案(程序设计说明)
牛顿法是一种线性化的近似方法，其基本思想是将非线性方程转化
为线性方程来进行迭代求解，当计算结果满足精度时，计算中止。
割线法具有超线性的收敛速度，求解速度快。根据所掌握知识将牛顿法和割线法用计算机语言表达出来，并运行求解，比较两种求解方法的优劣。
实验步骤或程序(经调试后正确的源程序)
见附件A
程序运行结果
见附件A
出现的问题及解决方法
出现的问题：在编辑器内编写程序时出现“字符数无效”的警告。
解决方法：将输入法切换至英文模式，重新输入。
附件A 沈阳工业大学实验报告
（适用计算机程序设计类）
专业班级：机械工程学号： 201820065 姓名：李奇
实验步骤或程序：
1、编写牛顿法M文件
区荣疆器-C:\Users\liqi\Desktop\

syns f m:
f= input。清ffiAf1 xl=，):
df>diff<f);
x0= input C清输入迭代初值*。=);
e=input 请输入根的误差眼I'):
N=input (* iH输入迭代次数限N=‘):
k=l：
8
g
10
11
12
13
14
while (k<N)
k=k0 ;
x1-xO-eval(f)/eval(df);
if abs (xl-xO) <e
fprintf「：<=%. 6f \n迭代次数为:%d\n?, k)
break
else
xO=x 1:
k=k+l;
end
end
- if k>^M
- fprintf (T 败5’)
- end
2、，输入求解
命令打窗口
» newton
fx清输入f（N）=
3、运行结果如下
命令行窗口 ®
>> nevton
请输入 f (x)=cos (x)-x*«cp (z)
请输入迭代初值成=。
请输入粮的误差限定1『-6
请输入迭代次数限股S
x=0i517757
迭代次数为：6
4、编写割线法M文件
子程序
届宾嚎器-C:\U s e rs\l i q i\ D e s kto p\ki n, m
fun,m ：■: +
□ function y = fun(x)
- 笋=cos (x)-x*exp(1：);
主程序
囱器-Ct\U rs\l i q i\D es kto p\^e ca nt. m
K +
function x = secant (x0? xl,tol)
- if nargin < 10
- tol - L 0l6 :
- end
- x = xO - fun(xOJ * (xO - xl) / (fun(sO) - fun(xl));
- n = 1 ；
一三while (abs(xO-xl) > tol)蛙(n <= 1000)
S - si = x0;
- ^0 = x；
- k = xO - fun(xO) * (xO - si) / (funfsO) - funtxl)):
- n - n + k
- end
- n
5、，输入方程
命令行窗口
>> secant (0, 0. 2, la-6)
命令行分口
>> secant (0, 0. 2, 1l6)
n =
e
ans =
0,5178
沈阳工业大学实验报告
(适用计算机程序设计类)
专业班级：机械工程学号： 201820065 姓名：李奇
实验名称：解方程组
实验目的：
掌握高斯消去法，列主元高斯消去法和 LU分解法求解方程组
能够使用MATLAB^程
实验内容：
学****掌握高斯消元，列主元，LU分解原理。
学****MATLAB!言的语法格式。