视频1 视频21 视频41 视频61 视频文章1 视频文章21 视频文章41 视频文章61 推荐1 推荐3 推荐5 推荐7 推荐9 推荐11 推荐13 推荐15 推荐17 推荐19 推荐21 推荐23 推荐25 推荐27 推荐29 推荐31 推荐33 推荐35 推荐37 推荐39 推荐41 推荐43 推荐45 推荐47 推荐49 关键词1 关键词101 关键词201 关键词301 关键词401 关键词501 关键词601 关键词701 关键词801 关键词901 关键词1001 关键词1101 关键词1201 关键词1301 关键词1401 关键词1501 关键词1601 关键词1701 关键词1801 关键词1901 视频扩展1 视频扩展6 视频扩展11 视频扩展16 文章1 文章201 文章401 文章601 文章801 文章1001 资讯1 资讯501 资讯1001 资讯1501 标签1 标签501 标签1001 关键词1 关键词501 关键词1001 关键词1501 专题2001
数值分析答案
2025-09-25 23:10:39 责编:小OO
文档


2-1 已知y=f(x)的数值如下:

(1) 

x0123
y2312147
(2)

x-2-101
y154524
求Lagrange插值多项式并写出截断误差。

解:

(1)

(2)

2-2 已知函数lnx的如下数据

x8101214
y2.079442.302592.484912.63906
试分别用Lagrange线性插值和二次插值计算ln(11.85)的近似值,并估计它的截断误差。

解:线性插值公式:

当x=11.85时,

二次插值:

 

误差估计:。

2-3 设为任意给定的n+1个互不相同的节点,证明:

(1)若f(x)为不高于n次的多项式,则f(x)关于这组节点的n次插值多项式就是它自己。

(2)若是关于这组节点的Lagrange基函数,则有恒等式

证明:

(1)

  因为f(x)是n次多项式,所以它的n+1阶导数为零。故f(x)关于这组节点的n次插值多项式就是它自己。

(2) 取,在处进行n次拉格朗日插值,则有

由于,故有  。

(3) 将按二项式展开,得 ,

则 

由上题的结论得:

   。

2-4 已知函数表

x0.10.20.40.60.9
y0.99500.98010.92110.82530.6216
试构造四次Newton插值多项式,计算cos0.47的近似值并估计截断误差。

解:

自变量函数值一阶差商二阶差商三阶差商四阶差商
0.1

0.9950
0.2

0.9801-0.149
0.40.9211-0.295-0.4867
0.60.8253-0.479-0.460.0534
0.90.6216-0.679-0.40.08570.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.16

2-5 在区间[-4,4]上给出f(x)=ex在等距节点下的函数表,若用二次插值求ex的近似值,要使截断误差不超过10-6,问所用函数表的步长应怎样选取?

解:在区间[xi-1,xi]上,记

误差

则用二次插值的步长应:

2-6 对区间[a,b]作步长为h的剖分,且,证明:在任意相邻两节点间做线性插值,其误差限为。

证明:区间上的误差限:

误差限:

2-7 设,计算差商,及.

解:

自变量函数值一阶差商
1-886
2-2975-20
=-20,

2-8 设在有三阶导数,,证明:当

证明:根据已知条件可得到如下表所示的插值条件:

xx0

x1

yf(x0)

f(x1)

y’f’(x0)

建立差商表:

自变量函数值一阶差商二阶差商
x0

f(x0)

x0

f(x0)

f’(x0)

x1

f(x1)

则由newton 插值公式可得:

整理得:

其中R(x)由以下计算得到:

构造辅助函数:   

有,,三个零点,有,,三个零点,则至少有一个零点,记作。

2-9 用下列函数值表构造不超过3次的插值多项式,并建立误差估计式。

x012
f(x)129
f’(x)

3
解:

建立差商表:

自变量函数值一阶差商二阶差商三阶差商
01
121
1232
29741
则由newton 插值公式可得:

误差估计式:。

2-10 求满足下列条件的Hermite插值多项式

xi

12
yi

23
y’i

1-1
解:

       

2-11 求一个不高于4次的插值多项式P4(x),使得

解:根据已知条件可得到如下表所示的插值条件:

x012
P011
P’

01
建立差商表:

自变量函数值一阶差商二阶差商三阶差商四阶差商
00
000
1111
1110-1
210-1

-0.50.25
则由newton 插值公式可得:

2-12 根据下表建立三次样条插值函数

x123
f(x)242
f’1(x)

1-1
解:, , 

列方程:

则三次样条插值函数为:

          =8-16x+13x2-3x3,        。

          =-40+56x-23x2+3x3,      。

2-13 已知y=f(x)的如下数值

x01234
y-8-701956
求三次样条插值函数S(x),满足边界条件

(1)S’(0)=0,S’(4)=48

(2)S”(0)=0,S”(4)=24

解:用三转角算法计算:

(1), , 

, , 

, , 

列方程组:

则三次样条插值函数为:

=x3-8,  。

=x3-8,  。

=x3-8,  。

=x3-8,  。

(2) 

   

列方程组:

则三次样条插值函数为:

=x3-8,  。

=x3-8,  。

=x3-8,  。

=x3-8,  。

用三弯矩算法计算:

(1) , , 

, , 

, , 

列方程组:

(2) 列方程组:

第三章 最佳逼近

3-1 求下列函数在指定区间上得一次最佳平方逼近多项式并估计平方误差

(1) ,

解: 设

法方程组为:

基函数为:,,

得到:,,,

,。

于是法方程组为:

解之得:,。

所以,最佳平方逼近一次多项式为:。

误差估计:

由误差估计式:

(2) ,

解: 设

法方程组为:

基函数为:,,

得到:,,,

,。

于是法方程组为:

解之得:,。

所以,最佳平方逼近一次多项式为:。

误差估计:

由误差估计式:

(3) ,

解: 设

法方程组为:

基函数为:,,

得到:,,,

,。

于是法方程组为:

解之得:,。

所以,最佳平方逼近一次多项式为:。

误差估计:

由误差估计式:

(4) ,

解: 设

法方程组为:

基函数为:,,

得到:,,,

,。

于是法方程组为:

解之得:,。

所以,最佳平方逼近一次多项式为:。

误差估计:

由误差估计式:

3-2 求,在上的最佳平方逼近多项式,并给出平方误差。

解:设

法方程组为:

基函数为:,,,

得到:,,

,,

,,

于是法方程组为:

解之得:,,。

所以,最佳平方逼近一次多项式为:。

误差估计:

由误差估计式:

=2.8814×10-12

3-3求参数,使达到极小。

解:本题也就是求f(x)=sinx的最佳平方逼近一次多项式。

法方程组为:

基函数为:,,

得到:,,,

,。

于是法方程组为:

解之得:,。

3-4已知一组数据如下:

xi2468
yi2112840
用最小二乘法求拟合这组数据的一条直线,并估计平方误差。

解:线性拟合:

根据基函数给出法方程组

求得  

法方程组为:

解得:c0=-12.5,c1=6.55

求得拟合线性多项式函数

p1(x)=-12.5+6.55x

误差为:

先计算出拟合函数值:

xi1111
P1

0.600        13.7026.8039.90
得到:10.7

3-5 已知函数值表

xi-2-112
yi

01210
试用二次多项式拟合这组数据并给出平方误差。

解:二次拟合:

根据基函数给出法方程组

求得  

法方程组为:

解得:c0=58/35=1.6571,c1=0,c2=-3/7=-0.4286

求得拟合线性多项式函数

p2(x)=1.6571-0.4286x2

误差为:

先计算出拟合函数值:

xi-2-1012
P2

-0.05731.22851.65711.2285-0.0573
得到:0.2286

3-6给出下列数据

xi-3-2-124
yi14.38.34.78.322.7
用最小二乘法求形如y=a+bx2的经验公式。

解:根据基函数给出法方程组

求得  

法方程组为:

解得:a=3.5,b=1.2。

3-7确定经验公式中的参数,使之与下列数据拟合:

xi0.10.20.30.40.50.6
yi0.1720.3230.4840.6901.0001.579
解: 将经验公式转化为:

令 ,,,,则上式转化为:。

上表的的数据变为:

xi0.10.20.30.40.50.6
zi5.814

3.096

2.0661.4491.0000.633

这时 ,,

zT =[5.814 3.096 2.066 1.449 1.000 0.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)的关系式为,现测得一组数据如下:

xi12346810121416
yi4.006.418.018.799.539.8610.3310.4210.5310.61
试确定出参数α、β。

解:将经验公式转化为:

上表的的数据变为:

xi12346810121416
zi0.250.1560.125

0.114

0.1050.1010.09680.096

0.095

0.094

这时  ,

zT =[0.25 0.156 0.125 0.114 0.105 0.101 0.0968 0.096 0.095 0.094],

解得:α= 0.1650,β= 0.07。

3-9 用最小二乘法求下列方程组的解

 (1)   (2) 

解:

(1) 

简化为:

两边同乘以系数矩阵的转置矩阵,就得到所需要的法方程组:

具体计算如下:           

解得最小二乘解:x1=26/11,x2=15/11

(2)

简化为:

两边同乘以系数矩阵的转置矩阵,就得到所需要的法方程组:

具体计算如下:

解得最小二乘解:x=1450/487=2.9774,y=597/487=1.2259

第四章  数值积分与数值微分

4-1 用四节点复化梯形公式计算积分

(1) ,    (2) 

解:(1) 

    (2) 

4-2用四节点复化Simpson公式计算积分

   (1) ,  (2) 

解:(1) 

    (2) 

4-3 分别用复化梯形和复化Simpson公式计算积分

并使绝对误差限不超过,问需要将区间[0, 1]多少等分?

解:复化梯形:

所以区间应该409等分。

复化Simpson公式:

所以区间应该6等分。

4-4 利用积分计算ln2时,若采用复化Simpson公式,问应取多少个节点才能使其误差的绝对值不超过。

解:

        所以应26等分,节点数为:2×26+1=53个。

4-5 直接验证Simpson求积公式

具有3次代数精确度。

证明:

当f(x)=1时,

,等式成立。

当 f(x)=x时,

,等式成立。

当 f(x)=x2时,

,等式成立。

当 f(x)=x3时,

,等式成立。

当 f(x)=x4时,

等式不成立,所以Simpson求积公式具有3次代数精度。

4-6 设函数由下表给出,分别用复化梯形和复化Simpson公式计算积分

xi

0.60.81.01.21.41.61.8
f(xi)

5.74.63.53.74.95.25.5
解:

复化梯形公式:

复化Simpson公式:

4-7 用两点Guass型求积公式计算积分

(1) 

(2) 

(3) 

解:

(1)

(2)

(3) 

4-8 用两点Guass-chebgshev公式计算积分

解:

4-9 如何用两点Guass型求积公式计算下列积分:

(1) ,(2) ,(3) ,(4) 。

解:

(1)

(2)(3) (4) 4-10 确定x1,x2,A1,A2使下式成为Guass型求积公式

解:

因为 ,

上面的求积公式显然是两点Guass型求积公式,其中,,

4-11 已知y=f(x)的如下数据

xi

0.60.80.911.11.21.4
yi

0.73600.83650.909511.11051.24461.6017
分别利用三点数值微分公式计算f’(1)和f”(1)。

解:

    下载本文

显示全文
专题