数值分析习题(含标准答案)
数值分析试题及答案

数值分析试题及答案一、选择题1. 下列哪个方法不适合用于求解非线性方程的根?A. 二分法B. 牛顿法C. 弦截法D. 正割法2. 当使用二分法求解非线性方程的根时,需要满足的条件是:A. 函数f(x)在区间[a, b]上连续B. 函数f(x)在区间[a, b]上单调递增C. 函数f(x)在区间[a, b]上存在根D. 函数f(x)在区间[a, b]上可导3. 数值积分是通过将定积分转化为求和的方法来近似计算积分值的过程。
下列哪个方法是常用的数值积分方法?A. 矩形法则B. 辛普森规则C. 梯形规则D. 高斯-勒让德法则4. 龙格-库塔法是常用于求解常微分方程的数值解法。
以下哪个选项是描述龙格-库塔法的特点?A. 该方法是一种多步法B. 该方法是一种多项式插值法C. 该方法是一种单步法D. 该方法是一种数值积分法5. 用有限差分法求解偏微分方程时,通常需要进行网格剖分。
以下哪个选项是常用的网格剖分方法?A. 多边形剖分法B. 三角剖分法C. 矩形剖分法D. 圆形剖分法二、解答题1. 将函数f(x) = e^x 在区间[0, 1]上用复化梯形规则进行数值积分,分为6个子区间,求得的近似积分值为多少?解:将区间[0, 1]等分为6个子区间,每个子区间的长度为h = (1-0)/6 = 1/6。
根据复化梯形规则的公式,近似积分值为:I ≈ (1/2) * h * [f(0) + 2f(1/6) + 2f(2/6) + 2f(3/6) + 2f(4/6) + 2f(5/6) +f(1)]≈ (1/2) * (1/6) * [e^0 + 2e^(1/6) + 2e^(2/6) + 2e^(3/6) + 2e^(4/6) +2e^(5/6) + e^1]2. 使用二分法求解方程 x^3 - 3x + 1 = 0 在区间[1, 2]上的根。
要求精确到小数点后三位。
解:首先需要判断方程在区间[1, 2]上是否存在根。
数值分析课后习题答案

0 1
0 10 1 1 0 0 0 1
0 0 12 1 1 2 0 0 0
1 2
0 0 0 1 1 0
1 2
1 2
1 2
1
0 0 0 1 0
1 2
1 2
0
1 2
1 2
0
0
0
341 1 1
2-5.对矩阵A进行LDLT分解和GGT分解,并求解方程组
Ax=b,其中
16 4 8
1
A 4 5 4 , b 2
8 4 22
3
解
16 A 4
4 5
84
44 11
2-3(1).对矩阵A进行LU分解,并求解方程组Ax=b,其中
2 1 1 A1 3 2
4 ,b6
1 2 2
5
解
2 A 1
1 3
1 2
2 11
22
1
5 2
1
3 21来自,所以 A12
1
2 1 1
5 3
2-2(1).用列主元Gauss消元法解方程组
3 2 6x1 4 10 7 0x2 7 5 1 5x3 6
解
3 2 6 4 10 7 0 7 10 7 0 7
r1r2
消元
10 7 0 7 3 2 6 4 0 0.1 6 6.1
r=0.5101-n/3.162…<0.5101-n/3<0.01% 因此只需n=5.即取101/2=3.1623
数值分析试题与答案

一. 填空题(本大题共4小题,每小题4分,共16分)1.设有节点012,,x x x ,其对应的函数()y f x =的值分别为012,,y y y ,则二次拉格朗日插值基函数0()l x 为 。
2.设()2f x x =,则()f x 关于节点0120,1,3x x x ===的二阶向前差分为 。
3.设110111011A -⎡⎤⎢⎥=--⎢⎥⎢⎥-⎣⎦,233x ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,则1A = ,1x = 。
4. 1n +个节点的高斯求积公式的代数精确度为 。
二.简答题(本大题共3小题,每小题8分,共24分)1. 哪种线性方程组可用平方根法求解?为什么说平方根法计算稳定?2. 什么是不动点迭代法?()x ϕ满足什么条件才能保证不动点存在和不动点迭代序列收敛于()x ϕ的不动点?3. 设n 阶矩阵A 具有n 个特征值且满足123n λλλλ>≥≥≥,请简单说明求解矩阵A 的主特征值和特征向量的算法及流程。
三.求一个次数不高于3的多项式()3P x ,满足下列插值条件:i x 1 2 3 i y 2 4 12 i y '3并估计误差。
(10分)四.试用1,2,4n =的牛顿-科特斯求积公式计算定积分1011I dx x=+⎰。
(10分) 五.用Newton 法求()cos 0f x x x =-=的近似解。
(10分) 六.试用Doolittle 分解法求解方程组:12325610413191963630x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥----⎣⎦⎣⎦⎣⎦ (10分)七.请写出雅可比迭代法求解线性方程组123123123202324812231530x x x x x x x x x ++=⎧⎪++=⎨⎪-+=⎩ 的迭代格式,并判断其是否收敛?(10分)八.就初值问题0(0)y yy y λ'=⎧⎨=⎩考察欧拉显式格式的收敛性。
(10分)《数值分析》(A )卷标准答案(2009-2010-1)一. 填空题(每小题3分,共12分) 1. ()1200102()()()()x x x x l x x x x x --=--; 2.7;3. 3,8;4. 2n+1。
数值分析题库答案(含详细解题步骤)

第 1 页/共 22 页1. 正方形的边长大约为100cm ,应怎样测量才干使面积误差不超过1cm 22. 已测得某场地长l 的值为110=*l m ,宽d 的值为80=*d m ,已知 2.0≤-*l l m,1.0≤-*d d m, 试求面积ld s =的绝对误差限与相对误差限.3.为使π的相对误差小于0.001%,至少应取几位有效数字?4.设x的相对误差界为δ,求n x的相对误差界.5.设有3个近似数a=2.31,b=1.93,c=2.24,它们都有3位有效数字,试计算p=a+bc的误差界和相对误差界,并问p的计算结果能有几位有效数字?第 3 页/共 22 页6. 已知333487.034.0sin ,314567.032.0sin ==,请用线性插值计算3367.0sin 的值,并预计截断误差.7. 已知sin0.32=0.314567, sin0.34=0.333487, sin0.36= 0.352274,用抛物插值计算sin0.3367的值, 并预计误差.8. 已知16243sin ,sin πππ===请用抛物插值求sin50的值,并预计误差9. . .6,8,7,4,1)(,5,4,3,2,1求四次牛顿插值多项式时设当==i i x f x第 5 页/共 22 页10. 已知4)2(,3)1(,0)1(=-=-=f f f , 求函数)(x f 过这3点的2次牛顿插 值多项式.11. 设x x f =)(,并已知483240.1)2.2(,449138.1)1.2(,414214.1)0.2(===f f f ,试用二次牛顿插值多项式计算(2.15)f 的近似值,并研究其误差12. 设],[)(b a x f 在上有四阶延续导数,试求满意条件)2,1,0()()(==i x f x P i i 及)()(11x f x P '='的插值多项式及其余项表达式.13. 给定3201219(),,1,,44f x x x x x ====试求()f x 在1944⎡⎤⎢⎥⎣⎦,上的三次埃尔米特插值多项式()P x ,使它满意11()()(0,1,2),()(),i i P x f x i P x f x ''===并写出余项第 7 页/共 22 页表达式.14. 设],1,0[,23)(2∈++=x x x x f 试求)(x f 在]1,0[上关于,,1{,1)(x span x =Φ=ρ}2x 的最佳平方逼近多项式15.已知实验数据如下:用最小二乘法求形如y=a+bx2的拟合曲线,并计算均方误差.16.已知数据表如下第 9 页/共 22 页x i 1 2 3 4 5 y iωi4 4.56 8 8.5 2 1 3 1 1试用最小二乘法求多项式曲线与此数据组拟合17. .1)(},1{span ,1]41[)(的最佳平方逼近多项式中的关于上的在在求==Φ=x x x x f ρ18. 决定求积公式⎰++≈10110)1()(32)0()(f A x f f A dx x f 中的待定参数110,,A x A , 使其代数精度尽量高,并指出所决定的求积公式的代数精度.19. 用复化辛普森公式计算积分⎰=10dx e I x , 问区间[0,1]应分多少等分才干使截断误差不超过?10215-⨯第 11 页/共 22 页20. 利用下表中给出的数据,分离用复化梯形公式和复化辛甫生公式计算定积分dx x I ln 21⎰=的近似值(要求结果保留到小数点后六位)21. 用复化梯形公式和复化辛甫生公式计算积分⎰=6.28.1)(dx x f I ,函数)(x f 在某些节点上的值如下图:(本题共14分)22. 决定公式⎰+≈101100)()()(x f A x f A dx x f x 的系数1010,,,x x A A ,使其具有最高代数精度23. 决定求积公式⎰++≈1110)1()(32)0()(f A x f f A dx x f 中的待定参数110,,A x A ,使其代数精度尽量高,并指出所决定的求积公式的代数精度第 13 页/共 22 页24.用LU 分解法求解以下方程组 (10分)123123142521831520x x x ⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪= ⎪⎪ ⎪ ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭25.用LU 分解法求解以下方程组⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛8892121514131615141321x x x26. 用LU 分解法求解以下方程组⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫⎝⎛542631531321321x x x27. 设方程组b Ax =,其中⎪⎪⎪⎭⎫⎝⎛-=220122101A ,Tb ⎪⎭⎫ ⎝⎛-=32,31,21, 已知它有解Tx ⎪⎭⎫⎝⎛-=0,31,21,若右端有小扰动61021-∞⨯=bδ,试预计由此引起的解的相对误差.第 15 页/共 22 页28. 设方程组b Ax =,其中212 1.0001A -⎛⎫= ⎪-⎝⎭,11.0001b -⎛⎫= ⎪⎝⎭,当右端向量b 有误差00.0001δ⎛⎫= ⎪⎝⎭b 时,试预计由此引起的解的相对误差(用∞范数计算)29. 给定b Ax =,其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=111a a a a a a A 证实:(1) 当121<<-a 时,A 对称正定,从而GS 法收敛. (2) 惟独当2121<<-a 时,J 法收敛.30. 对于线性方程组⎪⎩⎪⎨⎧-=+-=-+=+1242043 16343232121x x x x x x x ,列出求解此方程组的Jacobi 迭代格式,并判断是否收敛。
数值分析期末试题及答案

数值分析期末试题及答案一、选择题(每题5分,共20分)1. 在数值分析中,下列哪个算法不是用于求解线性方程组的?A. 高斯消元法B. 牛顿法C. 雅可比法D. 追赶法答案:B2. 插值法中,拉格朗日插值法属于:A. 多项式插值B. 样条插值C. 线性插值D. 非线性插值答案:A3. 以下哪个选项不是数值分析中的误差来源?A. 截断误差B. 舍入误差C. 计算误差D. 测量误差答案:C4. 在数值积分中,梯形法则的误差项是:A. O(h^2)B. O(h^3)C. O(h)D. O(1)答案:A二、填空题(每题5分,共20分)1. 牛顿插值法中,插值多项式的一般形式为:______。
答案:f(x) = a_0 + a_1(x-x_0) + a_2(x-x_0)(x-x_1) + ...2. 牛顿迭代法求解方程的根时,迭代公式为:x_{n+1} = x_n -f(x_n) / __________。
答案:f'(x_n)3. 在数值分析中,______ 用于衡量函数在区间上的近似积分值与真实积分值之间的差异。
答案:误差4. 线性方程组的解法中,______ 法是利用矩阵的LU分解来求解。
答案:克兰特三、解答题(每题10分,共60分)1. 给定函数f(x) = e^(-x),使用拉格朗日插值法,求x = 0.5时的插值值。
解答:首先选取插值节点x_0 = 0, x_1 = 0.5, x_2 = 1,对应的函数值分别为f(0) = 1, f(0.5) = e^(-0.5), f(1) = e^(-1)。
拉格朗日插值多项式为:L(x) = f(0) * (x-0.5)(x-1) / (0-0.5)(0-1) + f(0.5) * (x-0)(x-1) / (0.5-0)(0.5-1) + f(1) * (x-0)(x-0.5) / (1-0)(1-0.5)将x = 0.5代入得:L(0.5) = 1 * (0.5-0.5)(0.5-1) / (0-0.5)(0-1) + e^(-0.5) * (0.5-0)(0.5-1) / (0.5-0)(0.5-1) + e^(-1) * (0.5-0)(0.5-0.5) / (1-0)(1-0.5)计算得L(0.5) = e^(-0.5)。
数值分析期末考试题及答案

数值分析期末考试题及答案一、选择题(每题2分,共20分)1. 在数值分析中,下列哪个算法用于求解线性方程组?A. 牛顿法B. 高斯消元法C. 插值法D. 傅里叶变换答案:B2. 以下哪个选项不是数值分析中的误差类型?A. 舍入误差B. 截断误差C. 测量误差D. 累积误差答案:C3. 多项式插值中,拉格朗日插值法的特点是:A. 插值点必须等距分布B. 插值多项式的次数与插值点的个数相同C. 插值多项式是唯一的D. 插值多项式在插值点处的值都为1答案:B4. 在数值分析中,下列哪个方法用于求解非线性方程?A. 辛普森法则B. 牛顿迭代法C. 欧拉法D. 龙格-库塔法答案:B5. 以下哪个是数值稳定性的指标?A. 收敛性B. 收敛速度C. 条件数D. 误差传播答案:C二、简答题(每题10分,共20分)1. 简述高斯消元法求解线性方程组的基本原理。
答案:高斯消元法是一种直接解法,通过行变换将增广矩阵转换为上三角形式,然后通过回代求解线性方程组。
它包括三个基本操作:行交换、行乘以非零常数、行相加。
2. 解释什么是数值稳定性,并举例说明。
答案:数值稳定性是指数值解对输入数据小的扰动不敏感的性质。
例如,某些数值方法在计算过程中可能会放大舍入误差,导致结果不可靠,这样的方法就被认为是数值不稳定的。
三、计算题(每题15分,共30分)1. 给定线性方程组:\[\begin{align*}x + 2y - z &= 4 \\3x - y + 2z &= 1 \\-x + y + z &= 2\end{align*}\]使用高斯消元法求解该方程组,并给出解。
答案:首先将增广矩阵转换为上三角形式,然后回代求解,得到\( x = 1, y = 2, z = 1 \)。
2. 给定函数 \( f(x) = x^2 - 3x + 2 \),使用拉格朗日插值法在\( x = 0, 1, 2 \) 处插值,并求出插值多项式。
《数值分析》练习题及答案解析

《数值分析》练习题及答案解析一、单选题1. 以下误差公式不正确的是( D )A .()1212x x x x ∆-≈∆-∆B .()1212x x x x ∆+≈∆+∆C .()122112x x x x x x ∆≈∆+∆D .1122()x x x x ∆≈∆-∆ 2. 已知等距节点的插值型求积公式()()352kkk f x dx A f x =≈∑⎰,那么3kk A==∑( C )A .1 B. 2 C.3 D. 4 3.辛卜生公式的余项为( c )A .()()32880b a f η-''-B .()()312b a f η-''-C .()()()542880b a f η--D .()()()452880b a f η--4. 用紧凑格式对矩阵4222222312A -⎡⎤⎢⎥=-⎢⎥⎢⎥--⎣⎦进行的三角分解,则22r =( A ) A .1 B .12C .–1D .–25. 通过点()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足( D ) A .()00l x =0,()110l x = B . ()00l x =0,()111l x = C .()00l x =1,()111l x = D . ()00l x =1,()111l x =6. 用二分法求方程()0f x =在区间[],a b 上的根,若给定误差限ε,则计算二分次数的公式是n ≥( D )A .ln()ln 1ln 2b a ε-++ B. ln()ln 1ln 2b a ε-+-C.ln()ln 1ln 2b a ε--+ D. ln()ln 1ln 2b a ε--- 7.若用列主元消去法求解下列线性方程组,其主元必定在系数矩阵主对角线上的方程组是( B )A .123123123104025261x x x x x x x x x -+=⎧⎪-+=⎨⎪-+=-⎩ B 。
数值分析考试题和答案

数值分析考试题和答案一、单项选择题(每题2分,共20分)1. 在数值分析中,插值法的主要目的是()。
A. 求解线性方程组B. 求解非线性方程C. 构造一个多项式来近似一个函数D. 求解微分方程答案:C2. 线性方程组的高斯消元法中,主元为零时,应采取的措施是()。
A. 停止计算B. 回代求解C. 转置矩阵D. 行交换答案:D3. 以下哪种方法不是数值积分方法()。
A. 梯形规则B. 辛普森规则C. 牛顿法D. 复合梯形规则答案:C4. 以下哪种方法用于求解非线性方程的根()。
A. 欧几里得算法B. 牛顿迭代法C. 高斯消元法D. 线性插值法答案:B5. 在数值分析中,最小二乘法主要用于()。
A. 求解线性方程组B. 求解非线性方程C. 曲线拟合D. 微分方程数值解答案:C6. 以下哪种方法不是数值微分方法()。
A. 前向差分B. 后向差分C. 中心差分D. 欧拉方法答案:D7. 以下哪种方法用于求解常微分方程的初值问题()。
A. 欧拉方法B. 龙格-库塔方法C. 牛顿迭代法D. 高斯消元法答案:B8. 在数值分析中,矩阵的特征值问题可以通过()方法求解。
A. 高斯消元法B. 幂迭代法C. 牛顿迭代法D. 梯形规则答案:B9. 以下哪种方法不是数值稳定性分析中的方法()。
A. 绝对稳定性B. 相对稳定性C. 条件数D. 牛顿法答案:D10. 在数值分析中,条件数用于衡量()。
A. 算法的效率B. 算法的稳定性C. 算法的准确性D. 算法的复杂度答案:B二、填空题(每题2分,共20分)1. 在数值分析中,插值多项式的次数最高为______,其中n是插值点的个数。
答案:n-12. 线性方程组的高斯消元法中,如果某行的主元为零,则需要进行______。
答案:行交换3. 梯形规则的误差与被积函数的______阶导数有关。
答案:二4. 牛顿迭代法中,每次迭代需要计算______。
答案:函数值和导数值5. 最小二乘法中,残差平方和最小化时,对应的系数向量是______。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
]第一章 绪论姓名 学号 班级习题主要考察点:有效数字的计算、计算方法的比较选择、误差和误差限的计算。
1若误差限为5105.0-⨯,那么近似数有几位有效数字(有效数字的计算) 解:2*103400.0-⨯=x ,325*10211021---⨯=⨯≤-x x 故具有3位有效数字。
2 14159.3=π具有4位有效数字的近似值是多少(有效数字的计算) 解:10314159.0⨯= π,欲使其近似值*π具有4位有效数字,必需!41*1021-⨯≤-ππ,3*310211021--⨯+≤≤⨯-πππ,即14209.314109.3*≤≤π即取( , )之间的任意数,都具有4位有效数字。
3已知2031.1=a ,978.0=b 是经过四舍五入后得到的近似值,问b a +,b a ⨯有几位有效数字(有效数字的计算)解:3*1021-⨯≤-aa ,2*1021-⨯≤-b b ,而1811.2=+b a ,1766.1=⨯b a 2123****102110211021)()(---⨯≤⨯+⨯≤-+-≤+-+b b a a b a b a故b a +至少具有2位有效数字。
2123*****10210065.01022031.1102978.0)()(---⨯≤=⨯+⨯≤-+-≤-b b a a a b b a ab 故b a ⨯至少具有2位有效数字。
4设0>x ,x 的相对误差为δ,求x ln 的误差和相对误差(误差的计算)~解:已知δ=-**xx x ,则误差为 δ=-=-***ln ln xx x x x则相对误差为******ln ln 1ln ln ln xxx x xxx x δ=-=-5测得某圆柱体高度h 的值为cm h 20*=,底面半径r 的值为cm r 5*=,已知cm h h 2.0||*≤-,cm r r 1.0||*≤-,求圆柱体体积h r v2π=的绝对误差限与相对误差限。
(误差限的计算)解:*2******2),(),(h h r r r h r r h v r h v -+-≤-ππ绝对误差限为πππ252.051.02052)5,20(),(2=⨯⋅+⨯⋅⋅⋅≤-v r h v相对误差限为%420120525)5,20()5,20(),(2==⋅⋅≤-ππv v r h v 6设x 的相对误差为%a ,求nx y =的相对误差。
(函数误差的计算)解:%**a x x x =-,)%(******na xx x nxx x yy y nnn =-≤-=-【7计算球的体积,为了使体积的相对误差限为%1,问度量半径r 时允许的相对误差限为多大(函数误差的计算) 解:球体积为334)(r r v ⋅⋅=π,3**34)(r r v ⋅⋅=π 欲使%13344)()()(**3**2***=-=⋅⋅-⋅⋅=-r r r r r r r r v r v r v ππ,必须%31**=-r r r 。
8设⎰-=11dx e x eI x n n ,求证: (1))2,1,0(11 =-=-n nI I n n(2)利用(1)中的公式正向递推计算时误差逐步增大;反向递推计算时误差逐步减小。
(计算方法的比较选择)解:11111011011111][-------=-=-==⎰⎰⎰n x n xn x nx n n nI dx e x nedx e xn ex e de x eI 111101)1(----=-==⎰e e e dx eeI x[如果初始误差为*000I I -=ε,若是向前递推,有0221*11*!)1()1()1()1()1(εεεεn n n n nI nI I I nn n n n n n n -==--=-=---=-=----可见,初始误差0ε的绝对值被逐步地扩大了。
如果是向后递推n n I nn I 111-=-,其误差为n n n I I εεεε!)1(211)1(11)1111()1111(221*110-==⋅-=-=---= 可见,初始误差n ε的绝对值被逐步减少了。
第二章 插值法姓名 学号 班级^习题主要考察点:拉格朗日插值法的构造,均差的计算,牛顿插值和埃尔M 特插值构造,插值余项的计算和应用。
1已知1)2(,1)1(,2)1(===-f f f ,求)(x f 的拉氏插值多项式。
(拉格朗日插值) 解法一(待定系数法):设c bx ax x L ++=2)(,由插值条件,有⎪⎩⎪⎨⎧=++=++=+-12412c b a c b a c b a 解得:3/4,2/1,6/1=-==c b a 。
故 342161)(2+-=x x x L 。
解法二(基函数法):由插值条件,有1)12)(12()1)(1(1)21)((11()2)(1(2)21)(11()2)(1()(⋅-+-++⋅-+-++⋅------=x x x x x x x L·)1)(1(31)2)(1(21)2)(1(31-++-+---=x x x x x x3421612+-=x x 2已知9,4,10===x x x y ,用线性插值求7的近似值。
(拉格朗日线性插值)解:由插值节点与被插函数,可知,240==y ,391==y ,其线性插值函数为565134942949)(+=⋅--+⋅--=x x x x L 7的近似值为6.25135657)7(≈=+=L 。
3若),...1,0(n j x j =为互异节点,且有)())(())(()())(())(()(11101110n j j j j j j j n j j j x x x x x x x x x x x x x x x x x x x x x l ----------=+-+-,试证明),...1,0()(0n k x x l xnj k jk j =≡∑=。
(拉格朗日插值基函数的性质) 解:考虑辅助函数∑=-=nj k jkj x x l xx F 0)()(,其中,n k ≤≤0,),(∞-∞∈x 。
)(x F 是次数不超过n 的多项式,在节点i x x =(n i ≤≤0)处,有)()()(0=-=-=-=∑=ki k i k i i i k i nj k i i j k j i x x x x l x x x l x x F 这表明,)(x F 有n+1个互异实根。
故0)(≡x F ,从而∑=≡nj k jkj x x l x)(对于任意的n k ≤≤0均成立。
4已知352274.036.0sin ,333487.034.0sin ,314567.032.0sin ===,用抛物线插值计算3367.0sin 的值并估计截断误差。
(拉格朗日二次插值) 解:由插值条件,其抛物线插值函数为 {314567.0)36.032.0)(34.032.0()36.0)(34.0()(⋅----=x x x L333487.0)36.034.0)(32.034.0()36.0)(32.0(⋅----+x x352274.0)34.036.0)(32.036.0()34.0)(32.0(⋅----+x x将3367.0=x 代入,计算可得:3304.0)3367.0(≈L 。
其余项为:)36.0)(34.0)(32.0(!3sin )(----=x x x x r ξ其中,36.032.0<<ξ )36.0)(34.0)(32.0(61)(---≤x x x x r 故误差的上界为:71014.2)36.03367.0)(34.03367.0)(32.03367.0(61)3367.0(-⨯≤---≤r 。
|5用余弦函数x cos 在00=x ,41π=x ,22π=x 三个节点处的值,写出二次拉格朗日插值多项式, 并近似计算6cos π及其绝对误差与相对误差,且与误差余项估计值比较。
(拉格朗日二次插值)解:由插值条件,二次拉格朗日插值多项式为0)4/2/)(02/()4/)(0(21)2/4/)(04/()2/)(0(1)2/0)(4/0()2/)(4/()(⋅----+⋅----+⋅----=ππππππππππππx x x x x x x L22)2/(28)2/)(4/(8πππππ----=x x x x8508.09242)2/6/(6/28)2/6/)(4/6/(8)6(22≈+=----=ππππππππππL 绝对误差为:0153.01828439924223)6(6cos≈--=+-=-ππL 相对误差为:0179.028428439)6()6(6cos≈+--=-πππL L余项为: 、)2/)(4/(!3sin )(ππξ--=x x x x r ,其中,2/0πξ<<其余项的上界为:)2/)(4/(61)(ππ--≤x x x x r 0239.06)26)(46(661)6(43≈=--≤πππππππr 比较可知,实际计算所得的绝对误差较余项公式所估计出的值要小一些。
6已知函数值212)6(,82)4(,46)3(,10)1(,6)0(=====f f f f f ,求函数的四阶均差]6,4,3,1,0[f 和二阶均差]3,1,4[f 。
(均差的计算)解:采用列表法来计算各阶均差,有从表中可查得:]6,4,3,1,0[=f 。
故6]3,1,4[=f 。
其实,根据均差的对称性,6]4,3,1[]3,1,4[==f f ,该值在第一个表中就可以查到。
7设)())(()(10n x x x x x x x f ---= 求][1,0p x x x f 之值,其中1+≤n p ,而节点)1,1,0(+=n i x i 互异。
(均差的计算)解:由均差可以表示成为函数值的线性组合,有∑=-+-------=pi p i p i i i i i i i i p x x x x x x x x x x x xx f x x x f 0111101,0))(())(())(()(][而0)(=i x f p i ≤≤0,故0][1,0=p x x x f 。
;8如下函数值表建立不超过三次的牛顿插值多项式。
(牛顿插值多项式的构造) 解:故 )2)(1(4)1(381)(----++=x x x x x x x N 。
、9求一个次数小于等于三次多项式)(x p ,满足如下插值条件:2)1(=p ,4)2(=p ,3)2(='p ,12)3(=p 。
(插值多项式的构造)解法一(待定系数法):设d cx bx ax x p +++=23)(,则c bx ax x p ++='23)(2,由插值条件,有⎪⎪⎩⎪⎪⎨⎧=+++=++=+++=+++123927341242482d c b a c b a d c b a d c b a 解得:6,15,9,2-==-==d c b a 。