数值分析答案Word文档格式.doc
《数值分析答案Word文档格式.doc》由会员分享,可在线阅读,更多相关《数值分析答案Word文档格式.doc(18页珍藏版)》请在冰豆网上搜索。
2.63906
试分别用Lagrange线性插值和二次插值计算ln(11.85)的近似值,并估计它的截断误差。
线性插值公式:
当x=11.85时,
二次插值:
误差估计:
。
2-3设为任意给定的n+1个互不相同的节点,证明:
(1)若f(x)为不高于n次的多项式,则f(x)关于这组节点的n次插值多项式就是它自己。
(2)若是关于这组节点的Lagrange基函数,则有恒等式
证明:
因为f(x)是n次多项式,所以它的n+1阶导数为零。
故f(x)关于这组节点的n次插值多项式就是它自己。
(2)取,在处进行n次拉格朗日插值,则有
由于,故有。
(3)将按二项式展开,得,
则
由上题的结论得:
。
2-4已知函数表
0.1
0.2
0.4
0.6
0.9
0.9950
0.9801
0.9211
0.8253
0.6216
试构造四次Newton插值多项式,计算cos0.47的近似值并估计截断误差。
自变量
函数值
一阶差商
二阶差商
三阶差商
四阶差商
-0.149
-0.295
-0.4867
-0.479
-0.46
0.0534
-0.679
-0.4
0.0857
0.040375
P4(x)=0.9950-0.149(x-0.1)-0.4867(x-0.1)(x-0.2)+0.0534(x-0.1)(x-0.2)(x-0.4)
+0.040375(x-0.1)(x-0.2)(x-0.4)(x-0.6)
当x=0.47时,P4(x)=0.8916
2-5在区间[-4,4]上给出f(x)=ex在等距节点下的函数表,若用二次插值求ex的近似值,要使截断误差不超过10-6,问所用函数表的步长应怎样选取?
在区间[xi-1,xi]上,记
误差
则用二次插值的步长应:
2-6对区间[a,b]作步长为h的剖分,且,证明:
在任意相邻两节点间做线性插值,其误差限为。
区间上的误差限:
误差限:
2-7设,计算差商,及.
-886
-2975
-2089
=-2089,
,
2-8设在有三阶导数,,证明:
当
根据已知条件可得到如下表所示的插值条件:
x0
x1
f(x0)
f(x1)
y’
f’(x0)
建立差商表:
则由newton插值公式可得:
整理得:
其中R(x)由以下计算得到:
构造辅助函数:
有,,三个零点,有,,三个零点,则至少有一个零点,记作。
则
2-9用下列函数值表构造不超过3次的插值多项式,并建立误差估计式。
f(x)
9
f’(x)
7
误差估计式:
2-10求满足下列条件的Hermite插值多项式
xi
yi
y’i
2-11求一个不高于4次的插值多项式P4(x),使得
P
P’
-0.5
0.25
2-12根据下表建立三次样条插值函数
f’1(x)
,
列方程:
则三次样条插值函数为:
=8-16x+13x2-3x3,。
=-40+56x-23x2+3x3,。
2-13已知y=f(x)的如下数值
-8
-7
19
56
求三次样条插值函数S(x),满足边界条件
(1)S’(0)=0,S’(4)=48
(2)S”(0)=0,S”(4)=24
用三转角算法计算:
(1),,
列方程组:
=x3-8,。
(2)
用三弯矩算法计算:
(1),,
(2)列方程组:
第三章最佳逼近
3-1求下列函数在指定区间上得一次最佳平方逼近多项式并估计平方误差
(1),
设
法方程组为:
基函数为:
,,
得到:
,,,
,。
于是法方程组为:
解之得:
所以,最佳平方逼近一次多项式为:
由误差估计式:
(2),
(3),
(4),
3-2求,在上的最佳平方逼近多项式,并给出平方误差。
设
,,。
=2.8814×
10-12
3-3求参数,使达到极小。
本题也就是求f(x)=sinx的最佳平方逼近一次多项式。
3-4已知一组数据如下:
6
11
28
40
用最小二乘法求拟合这组数据的一条直线,并估计平方误差。
线性拟合:
根据基函数给出法方程组
求得,
解得:
c0=-12.5,c1=6.55
求得拟合线性多项式函数
p1(x)=-12.5+6.55x
误差为:
先计算出拟合函数值:
P1
0.600
13.70
26.80
39.90
10.7
3-5已知函数值表
0
试用二次多项式拟合这组数据并给出平方误差。
二次拟合:
c0=58/35=1.6571,c1=0,c2=-3/7=-0.4286
p2(x)=1.6571-0.4286x2
P2
-0.0573
1.2285
1.6571
0.2286
3-6给出下列数据
-3
14.3
8.3
4.7
22.7
用最小二乘法求形如y=a+bx2的经验公式。
a=3.5,b=1.2。
3-7确定经验公式中的参数,使之与下列数据拟合:
0.3
0.5
0.172
0.323
0.484
0.690
1.000
1.579
将经验公式转化为:
令,,,,则上式转化为:
上表的的数据变为:
zi
5.814
3.096
2.066
1.449
0.633
取
这时,,
zT=[5.8143.0962.0661.4491.0000.633],
a0=-1.9674,a1=0.9761,a2=0.5034。
则a=1.939,b=-3.908,c=1.987。
3-8在某化学反应里,生成物的质量浓度y(10-3g/cm3)与时间t(min)的关系式为,现测得一组数据如下:
16
4.00
6.41
8.01
8.79
9.53
9.86
10.33
10.42
10.53
10.61
试确定出参数α、β。
将经验公式转化为:
0.156
0.125
0.114
0.105
0.101
0.0968
0.096
0.095
0.094
这时,
zT=[0.250.1560.1250.1140.1050.1010.09680.0960.0950.094],
α=0.1650,β=0.0789。
3-9用最小二乘法求下列方程组的解
(1)
(2)
简化为:
两边同乘以系数矩阵的转置矩阵,就得到所需要的法方程组:
具体计算如下:
解得最小二乘解:
x1=26/11,x2=15/11
x=1450/487=2.9774,y=597/487=1.2259
第四章