数值计算(二分法、简单迭代法、Newton迭代法、弦截法(割线法、双点弦法))

6
-1.50412
12
-1.50489
18
-1.50505
表2-1
初值-1
k
x
1
-1
2
-1
表2-2
初值1.6结果x=1.69028
k
xk
k
xk
k
xk
1
1.6
8
1.68862
15
1.69023
2
1.65669
9
1.68927
16
1.69025
3
1.66987
10
1.68967
17
1.69027
4
1.6779
要求：讨论求解的全过程，对所用算法的局部收敛性，优缺点等作分析及比较，
第
2.1二分法
思想：在函数的单调有根区间内，将有根区间不断的二分，寻找方程的解。
步骤：1.取中点mid=(x0+x1)/2
2.若f(mid)=0,则mid为方程的根，否则比较与两端的符号，若与f(x0)
异号，则根在[x0,mid]之间，否则在[mid,x1]之间。
11
1.68991
18
1.69027
5
1.68278
12
1.69006
19
1.69028
6
1.68573
13
1.69015
20
1.69028
7
1.68753
14
1.6902
表2-3
牛顿迭代法（表
初值-1.5结果x=-1.50507
k
xk
k
xk
1
-1.5
4
-1.50504
2
-1.50471
5
-1.50506
2
1.68071
-0.151106
3
1.69126
0.0157988
4
1.69027
-0.000313515
5
1.69028
-6.3006e-007
表4-3
从测试结果可以看出二分法和简单迭代法的收敛速度远大于牛顿迭代和弦截法的收敛速度。二分法和简单迭代法的公式易于构造和计算，牛顿迭代法虽然收敛高，但要求导数，计算的复杂度高！双点弦法随稍慢于牛顿跌代法，可以用差商代替牛顿迭代法中的导数，降低了计算的复杂度！
-1.00005
1
-0.975
6
-1.00078
11
-0.999976
2
-1.0125
7
-0.999609
12
-1.00001
3
-0.99375
8
-1.0002
13
-0.999994
4
-1.00312
9
-0.999902
14
-1
表1-2
区间[1.5,1.8]
k
xk
k
xk
k
xk
0
1.65
7
1.69102
4．输出x1,即为方程的近似解。
2.3 Newton迭代法
思想：设r是 的根，选取 作为r的初始近似值，过点 做曲线 的切线L，L的方程为 ，求出L与x轴交点的横坐标 ，称x1为r的一次近似值。过点 做曲线
的切线，并求该切线与x轴交点的横坐标 ，称 为r的二次近似值。重复以上过程，得r的近似值序列，其中， 称为r的
2.计算x2=x1-f(x1)(x1-x0)/f(x1)-f(x0);
3.判断f(x2)是否满足精度要求，若没有则按照上述步骤继续迭代，否则输出x2.x2即为方程的近似解。
第
测试结果
函数图像
函数Y=x5-3x3+x-1
二分法
[-1.6,-1.3]
k
xk
k
xk
k
xk
0
-1.45
5
-1.50156
10
-1.50493
14
1.69029
1
1.725
8
1.69043
15
1.69029
2
1.6875
9
1.69014
16
1.69029
3
1.70625
10
1.69028
17
1.69028
4
1.69687
11
1.69036
18
1.69028
5
1.69219
12
1.69032
6
1.68984
13
1.6903
表1-3
简单迭代法（表
初值-1.5
k
xk
k
xk
k
xk
1
-1.5
7
-1.50435
13
-1.50493
2
-1.50217
8
-1.50453
14
-1.50497
3
-1.50287
9
-1.50466
15
1.50499
4
-1.50341
10
-1.50476
16
-1.50501
5
-1.50381
11
-1.50483
17
-1.50504
次近似值
步骤：1.计算原函数的导数f’(x);构造牛顿迭代公式
2.计算 ,若f’(x0)=0,退出计算，否则继续向下迭代。
3.若|x1-x0|满足精度要求，x1即为方程的近似解。
2.4弦截法
思想:为加速收敛，改用两个端点都在变动的弦，用差商替代牛顿迭代公式的导数f’(x)。
步骤：1.构造双点弦法的公式
1
-1.525
6
-1.50391
11
-1.505
2
-1.4875
7
-1.50508
12
-1.50504
3
-1.50625
8
-1.50449
13
-1.50506
4
-1.49688
9
-1.50479
14
-1.50507
表1-1
区间[-1.2,-0.9]
k
xk
k
xk
k
xk
0
-1.05
5
-0.998437
10
3
-1.50497
6
-1.50507
表3-1
初值-1结果x=-1.50507
k
x
1
-1
2
-1
表3-2
初值1.6结果x=1.69028
k
xk
k
xk
1
1.6
5
1.69024
2
1.68602
6
1.69027
3
1.68893
7
1.69028
4
1.68985
8wk.baidu.com
1.69028
表3-3
双点弦法（表
区间[-1.6,-1.3]结果x=-1.50507
k
xk
f(xk)
k
xk
f(xk)
1
-1.5
0.03125
5
-1.50667
0.0784566
2
-1.66149
0.376502
6
-1.505
-0.010079
3
-1.47175
-1.56322
7
-1.50507
0.000440988
4
-1.492
0.186801
8
-1.50507
2.30387e-006
本科生实验报告
实验课程数值计算方法
学院名称信息科学与技术学院
专业名称计算机科学与技术
学生姓名
学生学号
指导教师
实验地点
实验成绩
二〇一六年五月二〇一六年五月
实验一
1.1
实验目的：掌握非线性方程求根的基本步骤及方法，。
实验内容：试分别用二分法、简单迭代法、Newton迭代法、弦截法(割线法、双点弦法)，求x5-3x3+x-1= 0在区间[-8,8]上的全部实根，误差限为10-6。
3并重复上述步骤，直达达到精度要求，则mid为方程的近似解。
2.2简单迭代法
思想：迭代法是一种逐次逼近的方法，它是固定公式反复校正跟的近似值，使之逐步精确，最后得到精度要求的结果。
步骤：1.构造迭代公式f(x)，迭代公式必须是收敛的。
2.计算x1,x1=f(x0).
3.判断|x1-x0|是否满足精度要求，如不满足则重复上述步骤。
表4-1
区间[-1.2,-0.9]结果x=-1
k
xk
f(xk)
1
-1.01393
0.0415678
2
-1.0002
0.000607777
3
-0.999999
-3.11969e-006
4
-1
2.11001e-010
表4-2
区间[1.5,1.8]结果x=1.69028
k
xk
f(xk)
1
1.64403
-0.676455