数值分析实验报告一

/ 歹傍乂才衣寥

实验报告

课程名称：数值分析

实验项目：解方程/解方程组

专业班级：机械工程

姓名:奇学号：201820065实验室号：实验组号：

实验时间：批阅时间：

指导教师：曲绍波成绩：

.C

工业大学实验报告

(适用计算机程序设计类)

专业班级：机械工程学号： 201820065 : ____ 奇实验名称：解方程

工业大学实验报告附件A

(适用计算机程序设计类)

专业班级：机械工程学号： 201820065 : ___ 奇实验步骤或程序：

1、编写牛顿法M文件

器-C:1U er s\l I q i\ D e & kto p ^newto n. m n

ewton.m ;；+

1 - 5*115 f r:

2 -f-dnput「请输入T (s) = ?);

3 -df-dif f Cf >:

4 -x0-input 请输入J8代初值W ):

5 _e=input『请输入根的邊差限于”);

€ -N= input「请输入送代炖册):

f k=l;

S -'while (k^N)

S -x«xO:

1€ -adpQ-eval (f^/evnl (df};

L1 -辻abs (rl-jcO)

12 -fprintf ( x-%. 6f VdS代次数为：・ sl# k)

13 -treafc

14 -else

15 - xO=xl;

16 -k=h+l:

L：-end

1S -end

1© -if k>=N

2C -fprintf ('矢飯S〉

21 -end

2、调用newton.m文件，输入求解

誇令行扇口

口亡祇an

fx请输入f

3、运行结果如下

®>> nekton

请输入芳(x)=cos 'XJ

i青输Xiiffi初値期丸

请输入很的误差®e-10"-6

诸输入迭代戻数限皓8

K-0.517757

迭代次数为:6

4、编写割线法M文件

(1)子程序

厲- &U齟「班li qi\DeskiopXfun.m

® >

1 \珂+ 1

1 fixncticn y = fun(x)

2 —y = cos (x)-x*e?;p :

(2)主程序

iN 斎器-C;\U s-e r5\ 11 q i\D e^kt o peseta nt, m A secant.m | +

1function 咒 = seeant (x0, xl, t ol)

2 -if narcin < 10

3 -toL m 1. Oe-S：

end

生-

5 —x = M O - fun (M O) * (xO - x 1) / (fuiiGiO) - fun (xl)):

5 -n = 1:

n

vliile (aJbs (xO—xl)、tel)聪(n <—10005

i =

3 -xl = xC；

0 -x0 - X ;

10 -x = xO - funtxO) * (xO - xl> / (-fun

11 -n = n H- 1;

12 -Mid

13 --n

5、调用secant.m文件，输入方程

笳夸行窗口

>> 昌gant | Q, Q. Z 1音一6 | 6运行结果如下

命零行窗口

>)secant(0

T 0.2

h

le-6)

n ■

s axis =

0,5178

工业大学实验报告

(适用计算机程序设计类)

专业班级：机械工程________ 学号： 201820065 : ___ 奇实验名称：解方程组

工业大学实验报告

附件B

（适用计算机程序设计类）

专业班级：机械工程________ 学号： 201820065 : ___ 奇

实验步骤或程序：

1、编写高斯消元法M文件

/ 涓辑器-C; ；U sfir5\ I i q i\D e skto p\G a u 5sxia oq u. m

Gaussxi doqu.m :v +

1□ -function ERA. RB, n, Xl^auss^iaoqu (A f b) ■U为亲数拒蹲为右端顶

2 -> A b]:惭向增广拒阵

3 -nr length^);%右湍】页b的维敎

4 -触寸沁（A）:

9 _J?E=rwk(B);%熠广矩阵日的秩

6 =if RE-RAM

7 ―disp（h提示：@^]RA'=KB,所以此方程期无18'）

3 =return

9 -end

10 —if RA—RB

11 -if RA==n

12 -disp C握示:因対RA>RE=m所以址方琨组有唯一的筒J

13 一X=seros(n, 1);

14 - f or p»1: n- 1

15 -□fox h=p十1: n

16 —ffi-E(k,p)/B(p,pk Bgp: n+1) -B (k n p: n4-l) -m* †B tp, p:

n+1):

17 一卜end

18 —L end

19 -WCl-n,rr+L) ;A=Btl ：n1I n) ;X(n)=b (n)/ACn>n);

20 - □for q-n-1.-1:1

21 一X (q) = (b(q)-snni(A(q, f n) *X fq+1; n)))/A(q. q):

22 -end

23 -else

24 -disp（提尔：凶为RA・RBC%肋以此万程组戸尢好第貼J

25 -end

26 - end

2、调用Gaussxiaoqu.m文件，输入方程组

yy Ijaussz i AD

耒定义磴埶事翌宝A a

†A=：l 5. : -=. 5.-7 -IL 1;-? 1 H. 2 ].-：：;! -l.-T. 2 -5.~2 1 : -：： -E：. 2. ■ 3 1_ 2 -I ::-：：

R BllglE -7%73;9<»3^7.1;3.72]:

J& » Gousmioqy IA b)