数值分析实验二分法.pdfVIP

石家庄经济学院

实验报告

(学院)系:信息工程学院

专业:计算机科学与技术

姓名:000

班级:计科二班

学号:5111091111112

任课教师:00

《数值分析》实验报告

姓名学号日期

实验室152机房设备编号

实验题目用二分法求f(x)=x^3+x^2-3x-3=0的正根（精确到１e-3）

牛顿法求x-cos(x)=0的实根，精确到小数后六位

一实验目的

使用matlab编程，熟悉matlab编程环境。

分别使用二分法和牛顿法求方程的根，掌握二分法和牛顿法求方程根的过程及其算法。

深入分析二分法的优缺点和牛顿法与弦割法进行比较。

二实验内容

1二分法实验f(x)=x^3-x^2-3*x-3=0的正根（精度精确到０.００１）

算法实现

步1：构造区间[a,b],取a=1,b=2.先把c=(a+b)/2代入方程f(x),如果fb*fc0说明根的

区间在[a,c]，及令b=c;

步2：如果fb*fc不大于0说明根的区间在[c,b]，及令a=c;继续在新的区间内二分，直

到达到要求的精度0.001；

步3：循环执行步骤1,2，直到步骤2的条件成立；

步4：取区间中点及c为所求方程的根；

参数说明

a,b:区间端点

c:每次所取区间的中间值

x：最后所求的近似值

ya,yb:对应的所求函数值

abs（）:求绝对值函数

源代码

clc

closeall

formatlong;

f=input(f(x)=);

qujian=input(输入区间=);

err=input(误差范围;);

a=qujian(1);

b=qujian(2);

yc=1;

while((b-a)err)(yc~=0);

c=(a+b)/2;%设定中点

x=a;

ya=eval(f);

x=b;

yb=eval(f);

x=c;

yc=eval(f);

ifya*yc0%当根不在左边的子区间时，把右边的子区间设为迭代范围

b=c;

else%当根不在左边的子区间时，把右边的子区间设为迭代范围

a=c;

end

x0=c

End

运行结果

erfenfa

fx=

4.596249527821783e-005

x=

1.73205566406250

2牛顿法实验f(x)=x-cosx的正根（精度精确到0.000001）

算法实现

1：构造区间[-1,1],代入方程f(x),如果f(-1)*f(1)0说明根的区间在[-1,1];

2：利用fzero（’x-cos(x)’），取得临界点0.5;

3：取f(x)的一介导数，构造牛顿方法;

4：在构造区间[-1,1]中计算得到所求的方程根；

5：循环次数k为10次;

6:误差的大小为1E-6;

参数说明

Error:误差限