《数值分析》实验报告
1
f实验一、误差分析
误差问题是数值分析的基础,又是数值分析中一个困难的课题。在实际计算中,如果选用了不同的算法,由于舍入误差的影响,将会得到截然不同的结果。因此,选取算法时注重分析舍入误差的影响,在实际计算中是十分重要的。同时,由于在数值求解过程中用有限的过程代替无限的过程会产生截断误差,因此算法的好坏会影响到数值结果的精度。
一、实验目的
1、通过上机编程,复习巩固以前所学程序设计语言及上机操作指令;2、通过上机计算,了解误差、绝对误差、误差界、相对误差界的有关概念;
3、通过上机计算,了解舍入误差所引起的数值不稳定性。
二、实验任务
对
01220,计算定积分
算法1:利用递推公式
y
10
x
dx
x5
y
1
5y
1
取
y0
10
x
1dx5
l
6l
5
0182322
1220
算法2:利用递推公式
注意到
y
1
15
15
y
20191
111x20dx1x20dx11x20dx1
12660
0x5
50
105
2
f取
y20
1120105
1126
0008730
思考:从计算结果看,哪个算法是不稳定的,哪个算法是稳定的。
算法1:
tlog60log50
0yzeros121y1tfork221yk1k5yk1
1e
dy16y721
运行结果:
a
s0182304116
239141170695873432935049
算法2:
yzeros211
1y11105112620fork2112yk115kyk5
1e
d
运行结果:
y00884005800043100343002850024300212001880016900154001410013000120001120010500099000930008900081000950
由数据对比可知,算法2较为稳定。
3
f实验二、插值法
插值法是函数逼近的一种重要方法,它是数值积分、微分方程数值解等数值计算的基础与工具,其中多项式插值是最常用和最基本的方法。拉格朗日插值多项式的优点是表达式简单明确,形式对称,便于记忆,它的缺点是如果想要增加插值节点,公式必须整个改变,这就增加了计算工作量。而牛顿插值多项式对此做了改进,当增加一个节点时只需在原牛顿插值多项式基础上增加一项,此时原有的项无需改变,从而达到节省计算次数、节约存储单元、应用较少节点达到应有精度的目的。
一、实验目的
1、理解插值的基本概念,掌握各种插值方法,包括拉格朗日插值和牛顿插值等,注意其不同特点;2、通过实验进一步理解并掌握各种插值的基本算法。
二、实验任务
1、已知函数表
xi056160056280056401056521yi082741082659082577082495用二次拉格朗日插值多项式求x05635时的函数近似值。2、已知函数表
xi04
055
065
08
yi041075057815069675088811
用牛顿插值r