交大数值分析 题库 课后题目 大全 期末

合集下载

山东交通学院数值分析期末复习题

山东交通学院数值分析期末复习题

数值分析A 卷复习题一、填空题(每题3分,共30分)1.已知 3.201,0.57a b ==是经过四舍五入后得到的近似值,则a b +有( 2 )位有效数字。

2.设()i l x 是以,(0,1,,9)k x k k ==为节点的Lagrange 插值基函数,则9()==∑ik kl k ( i )。

3.n 个求积节点的插值型求积公式的代数精度至少为( 1n - )。

4.设有矩阵2304A -⎡⎤=⎢⎥⎣⎦,则1A =( 7 )。

5.三次样条函数是在各个子区间上的( 3 )次多项式。

6.用牛顿下山法求解方程303x x -=的下山条件是( 1()()k k f x f x +≤ )。

7.已知3n =时的Newton-Cotes 系数(3)018C =,则(3)3C =( 1/8 )。

8.若线性代数方程组Ax b =的系数矩阵A 为严格对角占优阵,则雅可比迭代和高斯-塞德尔迭代都 ( 收敛 )。

9.用Gauss-Seidel 迭代法解方程组1212423x ax ax x +=⎧⎨+=-⎩,其中a 为实数,方法收敛的充要条件是a 满足a <<。

二、单选题(每题3分,共30分)1.用1+x 近似表示xe 所产生的误差是( C )误差。

A.模型B.观测C.截断D.舍入 2.已知数21234721,0.721,0.700,710x x x x -====⨯是由四舍五入得到的,则它们的有效数字的位数应分别为( A )。

A.3,3,3,1B.3,3,3,3C.3,3,1,1D.3,3,3,2 3.设(1)1,(0)3,(2)4-===f f f ,则抛物线插值多项式中2x 的系数为( A )。

A.-0.5B.0.5C.2D.-24. 为求方程3210x x --=在区间[]1.3,1.6内的一个根,把方程改写成下列形式,并建立相应的迭代公式,迭代公式不收敛的是( A )。

A.2111,:k x x x +==-迭代公式1221111,:k kx x x x +=+=+迭代公式C.32213111/,:()k kx x x x +=+=+迭代公式 D.23212111,:k k kk xx x x x x +-==+++迭代公式5.线性方程组AXB =能用高斯消元法求解的充分必要条件为( D )。

数值分析课后习题及答案

数值分析课后习题及答案

第一章 绪论(12) 第二章 插值法(40-42)2、当2,1,1-=x 时,4,3,0)(-=x f ,求)(x f 的二次插值多项式。

[解]372365)1(34)23(21)12)(12()1)(1(4)21)(11()2)(1()3()21)(11()2)(1(0))(())(())(())(())(())(()(2221202102210120120102102-+=-++--=+-+-⨯+------⨯-+-+-+⨯=----+----+----=x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L 。

3、给出x x f ln )(=的数值表用线性插值及二次插值计算54.0ln 的近似值。

X 0.4 0.5 0.6 0.7 0.8 x ln -0.916291 -0.693147 -0.510826 -0.357765 -0.223144[解]若取5.00=x ,6.01=x ,则693147.0)5.0()(00-===f x f y ,510826.0)6.0()(11-===f x f y ,则604752.182321.1)5.0(10826.5)6.0(93147.65.06.05.0510826.06.05.06.0693147.0)(010110101-=---=--⨯---⨯-=--+--=x x x x x x x x x y x x x x y x L ,从而6202186.0604752.19845334.0604752.154.082321.1)54.0(1-=-=-⨯=L 。

若取4.00=x ,5.01=x ,6.02=x ,则916291.0)4.0()(00-===f x f y ,693147.0)5.0()(11-===f x f y ,510826.0)6.0()(22-===f x f y ,则 217097.2068475.404115.2)2.09.0(5413.25)24.0(3147.69)3.01.1(81455.45)5.06.0)(4.06.0()5.0)(4.0()510826.0()6.05.0)(4.05.0()6.0)(4.0()693147.0()6.04.0)(5.04.0()6.0)(5.0(916291.0))(())(())(())(())(())(()(22221202102210120120102102-+-=+--+-⨯++-⨯-=----⨯-+----⨯-+----⨯-=----+----+----=x x x x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L ,从而61531984.0217097.21969765.259519934.0217097.254.0068475.454.004115.2)54.0(22-=-+-=-⨯+⨯-=L补充题:1、令00=x ,11=x ,写出x e x y -=)(的一次插值多项式)(1x L ,并估计插值余项。

数值分析期末考试题

数值分析期末考试题

数值分析期末考试题一、选择题1. 在数值分析中,用于求解线性方程组的雅可比方法属于以下哪种迭代法?A. 直接迭代法B. 间接迭代法C. 外推法D. 松弛法2. 插值法中,拉格朗日插值多项式的主要特点是?A. 适用于多项式插值B. 适用于函数值已知的情况C. 只适用于单点插值D. 适用于分段插值3. 在数值积分中,辛普森法则是一种?A. 单区间求积公式B. 双区间求积公式C. 三区间求积公式D. 多区间求积公式4. 误差分析中,截断误差通常与以下哪个概念相关?A. 舍入误差B. 舍入误差的补偿C. 条件数D. 病态条件5. 非线性方程求解中,牛顿法的收敛速度通常?A. 较慢B. 较快C. 与初始值有关D. 与方程的性质有关二、填空题1. 在求解三对角线性方程组时,托马斯算法是一种________方法。

2. 多项式插值中,牛顿插值多项式可以通过________法来构建。

3. 数值积分中,高斯求积法是一种________方法。

4. 误差传递的估计通常通过________公式来进行。

5. 非线性方程的求解中,二分法是一种________方法。

三、简答题1. 请简述数值分析中的条件数概念及其在解方程中的应用。

2. 描述线性方程组迭代法中的收敛性判断方法,并给出收敛域的计算公式。

3. 解释插值和拟合的区别,并举例说明各自的应用场景。

4. 阐述数值积分中梯形法则的原理及其误差估计方法。

5. 讨论非线性方程求解中不动点理论和收敛性的关系。

四、计算题1. 给定线性方程组如下,请使用高斯消元法求解未知数x、y、z的值: \[\begin{cases}2x + y + z = 6 \\x + 3y + 2z = 11 \\3x + y + 4z = 17\end{cases}\]2. 假设有一个函数f(x) = sin(x),给定插值节点如下,请使用拉格朗日插值法构造一个三次插值多项式,并计算在x=π/4处的插值误差。

《数值分析》A卷期末考试试题及参考答案

《数值分析》A卷期末考试试题及参考答案

一、单项选择题(每小题3分,共15分) 1、用Simpson 公式求积分1401x dx +⎰的近似值为 ( ).A.2924 B.2429C.65D. 562、已知(1)0.401f =,且用梯形公式计算积分2()f x dx ⎰的近似值10.864T =,若将区间[0,2]二等分,则用递推公式计算近似值2T 等于( ). A.0.824 B.0.401 C.0.864 D. 0.8333、设3()32=+f x x ,则差商0123[,,,]f x x x x 等于( ).A.0B.9C.3D. 64的近似值的绝对误差小于0.01%,要取多少位有效数字( ). A.3 B.4 C.5 D. 25、用二分法求方程()0=f x 在区间[1,2]上的一个实根,若要求准确到小数 点后第四位,则至少二分区间多少次( ).A.12B.13C.14D. 15二、填空题(每小题4分,共40分)1、对于迭代函数2()=(3)ϕ+-x x a x ,要使迭代公式1=()ϕ+k k x x则a 的取值范围为 .2、假设按四舍五入的近似值为2.312,则该近似值的绝对误差限为 .3、迭代公式212(3)=,03++>+k k k k x x a x a x a收敛于α= (0)α>. 4、解方程4()530f x x x =+-=的牛顿迭代公式为 . 5、设()f x 在[1,1]-上具有2阶连续导数,[1,1]x ∀∈-,有1()2f x ''≤,则()f x 在[1,1]-上的线性插值函数1()L x 在点0处的误差限1(0)R ≤______.6、求解微分方程初值问题2(0)1'=-⎧⎨=⎩y xy yy ,0x 1≤≤的向前Euler 格式为 .7、设310131013A -⎛⎫⎪=-- ⎪ ⎪-⎝⎭,则A ∞= .8、用梯形公式计算积分112-⎰dx x 的近似值为 . 9、设12A 21+⎡⎤=⎢⎥⎣⎦a 可作Cholesky 分解,则a 的取值范围为 . 10、设(0)1,(0.5) 1.5,(1)2,(1.5) 2.5,(2) 3.4f f f f f =====,若1=h ,则用三点公式计算(1)'≈f .三、解答题(共45分) 1、给定数据用复化Simpson 公式计算1.381.30()f x dx ⎰的近似值,并估计误差,小数点后保留3位. (8分)2、用直接三角分解法求线性代数方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡432631531321321x x x 的解. (8分) 3、求()λx ,使得迭代公式1()()λ+=+k k k k f x x x x 求方程2()31=+-f x x x 的根的相应迭代序列{}k x 具有平方收敛. (5分)4、已知数据试对数据用最小二乘法求出形如=+y x b的拟合曲线. (8分) 5、已知(2)8f -=,(0)4f =,(2)8=f ,试求二次拉格朗日插值多项式. (8分) 6、设矩阵A 如下,根据谱半径判断用Jacobi 迭代法求解方程组Ax b =的敛散性.(8分)1102111221012A ⎡⎤-⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦一、单项选择题(每小题3分,合计15分) 1、A 2、D 3、C 4、C 5、D 二、填空题(每小题3分,合计30分) 1、0<<a ; 2、31102-⨯; 3;4、4135345++-=-+k k k k k x x x x x ; 5、14; 6、1(2)+=+-n n n n n y y h x y y ; 7、5;8、34-; 9、3>a ;10、1.2;三、计算题(合计55分) 1、给定数据用复化Simpson 公式计算 1.381.30()f x dx ⎰的近似值,并估计误差,小数点后保留3位. (8分)解: 401024S [()4()()]6-=++x x f x f x f x ………… 1分 1.38 1.30(3.624 4.20 5.19)6-=+⨯+ 0.341= ………… 2分20422012234S [()4()()][()4()()]66--=+++++x x x xf x f x f x f x f x f x =0.342 ………… 6分2211[]15-≈-I S S S =-⨯40.6710 ………… 8分 2、用直接三角分解法求线性代数方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡432631531321321x x x 的解. (8分) 解:设111213212223313233u u u 123100135l 100u u 136l l 100u ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=*⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦………… 1分 111=u ,212=u ,313=u ,121=l ,131=l 122=u ,223=u ,132=l133=u ,133=l …………6分所以⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=111011001L ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=100210321U …………7分 由b Ly =得Ty )1,1,2(=;由y Ux =得Tx )1,1,1(-=. ………… 8分3、求()λx ,使得迭代公式1()()λ+=+k k k k f x x x x 求方程2()31=+-f x x x 的根的相应迭代序列{}k x 具有平方收敛.(6分)解:要使迭代序列具有平方收敛,则()0ϕ'*=x ………… 2分 而()()()ϕλ=+f x x x x ,即 ………… 3分 2()()()()10()λλλ''**-**+=*f x x x f x x …………4分 而()0*=f x 则有()1()λ'*=-*f x x ………… 5分所以()()23λ'=-=--x f x x ………… 6分4、已知数据试对数据用最小二乘法求出形如=+ay x b的拟合曲线. (8分) 解:因为11=+b x y a a ,令0111,,,====b a a y x x a a y……2分 则有法方程01461061410⎛⎫⎛⎫⎛⎫= ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭a a ……5分解出014,1==-a a ,则1,4=-=-a b ……7分 所以1=4-y x……8分5、已知(2)8f -=,(0)4f =,(2)8=f ,试求二次拉格朗日插值多项式. (7分)解:01()(2)8l x x x =- …………2分 211()(4)4l x x =-- …………4分21()(2)8l x x x =+ …………6分 2012()()(2)()(0)()(2)L x l x f l x f l x f =-++24=+x …………7分6、设矩阵A 如下,根据谱半径判断用Jacobi 迭代法求解方程组Ax b =的敛散性.(8分)1102111221012A ⎡⎤-⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦解:100010001D ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,00010021002L ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦,10021002000U ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦…………3分1100211()0221002J B D L U -⎡⎤⎢⎥⎢⎥⎢⎥=+=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦…………5分 2102111()0222102J E B λλλλλλ⎡⎤-⎢⎥⎢⎥⎢⎥-=--=-=⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦…………6分()2J B ρ=…………7分 所以用Jacobi 迭代法求解方程组Ax b =收敛 …………8分。

数值分析考试卷及详细答案解答汇总

数值分析考试卷及详细答案解答汇总

姓名 __________ 班级 ___________ 学号 _____________一、选择题i.F (2,5,-3,4)表示多少个机器数(C ).A 64B 129C 257D 256 2. 以下误差公式不正确的是(D )A ・ £(迎 *一七 *)« 5(Xj*)+£(£ *) c ,£(“*•£ *)«|^2 *k (-'l*) + |时住2 *)3. 设° =(、任_1)6,从算法设计原则上定性判断如下在数学上等价的表达式,哪一个在数值计算上将给出°较好的近似值? (D )A ———B 99-70V2C (3-2V2)3D —— (V2 +1)6 (3 + 204. 一个30阶线性方程组,若用Crammer 法则来求解,则有多少次乘法?(A ) A31X29X30! B 30X30X30! C31X30X31! D 31X29X29!5. 用一把有亳米的刻度的米尺来测量桌子的长度,读出的长度1235mm,桌子的精确长度 记为(D ) A 1235mm B 1235-0.5mm C 1235+0.5nun D 1235±0.5mm二、填空1. 构造数值算法的基本思想是 近似替代、离散化、递推化 。

2. 十进制123.3转换成二进制为1111011.0而1。

3. 二进制110010.1001转换成十进制为 50.5625 。

4. 二进制o.ioi 转换成十进制为-o75.已知近似数X *有两位有效数字,则其相对误差限 5%。

6.1112=0.69314718...,精确到 10一’的近似值是 0.693。

* *7. x = ;r = 3.1415926・・・,则“ =3.1416 , =3.141的有效数位分别为5 和 3 __________ o8. 设卅=2.001,严=-0.8030是由精确值x 和y 经四舍五入得到的近似值,则兀* +y *的误差限____________________ o9.设x = 2.3149541•…,取5位有效数字,则所得的近似值卅二2.3150 。

数值分析课后习题与解答

数值分析课后习题与解答

课后习题解答第一章绪论习题一1.设x>0,x*的相对误差为δ,求f(x)=ln x的误差限。

解:求lnx的误差极限就是求f(x)=lnx的误差限,由公式(已知x*的相对误差满足,而,故即2.下列各数都是经过四舍五入得到的近似值,试指出它们有几位有效数字,并给出其误差限与相对误差限。

解:直接根据定义和式(有5位有效数字,其误差限,相对误差限有2位有效数字,有5位有效数字,3.下列公式如何才比较准确?(1)(2)解:要使计算较准确,主要是避免两相近数相减,故应变换所给公式。

(1)(2)4.近似数x*=0.0310,是 3 位有数数字。

5.计算取,利用:式计算误差最小。

四个选项:第二、三章插值与函数逼近习题二、三1. 给定的数值表用线性插值与二次插值计算ln0.54的近似值并估计误差限. 解:仍可使用n=1及n=2的Lagrange插值或Newton插值,并应用误差估计(5.8)。

线性插值时,用0.5及0.6两点,用Newton插值误差限,因,故二次插值时,用0.5,0.6,0.7三点,作二次Newton插值误差限,故2. 在-4≤x≤4上给出的等距节点函数表,若用二次插值法求的近似值,要使误差不超过,函数表的步长h 应取多少?解:用误差估计式(5.8),令因得3. 若,求和.解:由均差与导数关系于是4. 若互异,求的值,这里p≤n+1.解:,由均差对称性可知当有而当P=n+1时于是得5. 求证.解:解:只要按差分定义直接展开得6. 已知的函数表求出三次Newton均差插值多项式,计算f(0.23)的近似值并用均差的余项表达式估计误差.解:根据给定函数表构造均差表由式(5.14)当n=3时得Newton均差插值多项式N3(x)=1.0067x+0.08367x(x-0.2)+0.17400x(x-0.2)(x-0.3) 由此可得f(0.23) N3(0.23)=0.23203由余项表达式(5.15)可得由于7. 给定f(x)=cosx的函数表用Newton等距插值公式计算cos 0.048及cos 0.566的近似值并估计误差解:先构造差分表计算,用n=4得Newton前插公式误差估计由公式(5.17)得其中计算时用Newton后插公式(5.18)误差估计由公式(5.19)得这里仍为0.5658.求一个次数不高于四次的多项式p(x),使它满足解:这种题目可以有很多方法去做,但应以简单为宜。

数值分析课后部分习题答案



x * = 2.00021 = 0.200021 × 101 ,即 m = 1
1 1 × 10m − n = × 10−3 , 2 2
由有效数字与绝对误差的关系得 即
m − n = −3 ,所以, n = 2 ; y* = 0.032 = 0.32 × 101 ,即 m = 1
由有效数字与绝对误差的关系得 即
m − n = −3 ,所以, n = 4 ; z * = 0.00052 = 0.52 × 10−3 ,即 m = −3
1 1 × 10m − n = × 10−3 , 2 2
由有效数字与绝对误差的关系得 即
m − n = −3 ,所以, n = 0 .
1 1 × 10m − n = × 10−3 ,Fra bibliotek2 2=
f [x1 , x2 ,⋯ , x n ]-f [ x0 , x1 ,⋯ , x n−1 ] g[ x1 , x2 ,⋯ , x n ] − g[ x0 , x1 ,⋯ , x n−1 ] + x n − x0 x n − x0
( x − 1)( x − 2)( x − 3) 1 =- ( x − 1)( x − 2)( x − 3) , (0 − 1)(0 − 2)(0 − 3) 6
x ( x − 2)( x − 3) 1 = x ( x − 2)( x − 3) , (1 − 0)(1 − 2)(1 − 3) 2 x( x − 1)( x − 3) 1 =- x( x − 1)( x − 3) , (2 − 0)(2 − 1)(2 − 3) 2 x( x − 1)( x − 2) 1 = x ( x − 1)( x − 2) , (3 − 0)(3 − 1)(3 − 2) 6

数值分析期末考卷

数值分析期末考卷一、选择题(每题4分,共40分)A. 插值法B. 拟合法C. 微分法D. 积分法A. 高斯消元法B. 高斯赛德尔迭代法C. 共轭梯度法D.SOR方法3. 下列哪个算法不是求解非线性方程的方法?A. 二分法B. 牛顿法C. 割线法D. 高斯消元法A. 梯形法B. 辛普森法C. 高斯积分法D. 复化求积法A. 欧拉法B. 龙格库塔法C.亚当斯法D. 高斯消元法A. 幂法B. 反幂法C. 逆迭代法D. QR算法A. 梯度下降法B. 牛顿法C. 共轭梯度法D. 高斯消元法A. 拉格朗日插值法B. 牛顿插值法C. 埃尔米特插值法D. 分段插值法A. 前向差分法B. 后向差分法C. 中心差分法D. 拉格朗日插值法A. 牛顿法B. 割线法C. 雅可比迭代法D. 高斯消元法二、填空题(每题4分,共40分)1. 数值分析的主要任务包括数值逼近、数值微积分、数值线性代数和______。

2. 在求解线性方程组时,迭代法的收敛速度与______密切相关。

3. 牛顿法的迭代公式为:x_{k+1} = x_k f(x_k)/______。

4. 在数值积分中,复化梯形公式的误差为______。

5. 求解常微分方程初值问题,龙格库塔法的阶数取决于______。

6. 矩阵特征值的雅可比方法是一种______方法。

7. 梯度下降法在求解无约束优化问题时,每次迭代的方向为______。

8. 拉格朗日插值多项式的基函数为______。

9. 数值微分中的中心差分公式具有______阶精度。

10. 在求解非线性方程组时,牛顿法的迭代公式为:x_{k+1} =x_k J(x_k)^{1}______。

三、计算题(每题10分,共60分)1. 给定数据点(1,2),(2,3),(3,5),(4,7),求经过这四个数据点的拉格朗日插值多项式。

2. 用牛顿迭代法求解方程x^3 2x 5 = 0,初始近似值为x0 = 2,计算前三次迭代结果。

数值分析-课后习题答案

(2)如果A为对称正定矩阵,则Cond2(A)=1/n,1和n 分别为A的最大和最小特征值.
证明 (1)A正交,则ATA=AAT=E,Cond2(A)=A2A-12=1. (2)A对称正定,ATA=A2, A2=1. A-12=1/n.
精选课件
12
三.习题3 (第75页)
3-2.讨论求解方程组Ax=b的J迭代法和G-S迭代法的收
计算结果如下:
x x 1 2 ( (k k 1 1 ) ) 3 2 1 2 .x 5 2 (x k ) 1 (k 1 )
k
J法x1(k)
0
1.01
1
0.98
2
2.03
3
1.94
4
5.09
5
4.82
6
14.27
J法x2(k) 1.01 0.485 0.53 -1.045 -0.91 -5.635精选课件 -5.23
1.01
1.01
1
0.66
0.995
0.66
1.17
2
0.67
1.17
0.553333
1.223333
3
0.553333
1.165
0.517778
1.241111
4
0.556667
1.223333
0.505926
1.247037
5
0.517778
1.221667
0.501975
1.249012
6
0.518889
3 4精1选 课件
1
1
5
2-5.对矩阵A进行LDLT分解和GGT分解,并求解方程组
Ax=b,其中
16 4 8
1

数值分析课后习题全解

* x8 − x8 = 0.3 ≤
x7 = 0.96 × 10 −3
* x7 − x7 = 0
x7 精确
x8 = −8700.3
1 × 10 0 2
x8 具有 4 位有效数字, x8 = −8700 精确
2.以下各数均为有效数字: (1) 0.1062 + 0.947; (3)2.747 × 6.83; (2)23.46―12.753; (4)1.473 / 0.064 。
R R er ( R ) = 2πh ⋅ er ( R ) = er ( R ) S 2πRh 0.05 = 0.0025 20
er ( S ) ≈ er ( R ) ≤
答 计算体积的相对误差限为 0.005,计算侧面积的相对误差限为 0.0025
9.试改变下列表达式,使计算结果比较精确: (1) (2) (3) (4)
(3) x1 = 2.747
e( x1 ) ≤
x2 = 6.83
x1 x2 = 18.76201,
1 1 × 10 − 3 , e( x2 ) ≤ × 10 − 2 2 2
e( x1 x2 ) ≈ x2 e( x1 ) + x1e( x2 ) ≤ x2 e( x1 ) + x1 e( x2 )
1 1 1 ≤ 6.83 × × 10 − 3 + 2.747 × × 10 − 2 = × 10 − 2 × (0.683+2.747)=0.01715 2 2 2
e( A1 ) 10 −2 =1 er ( A1 ) = ≤ 0.01 A1
不能肯定所得结果具有一位有效数字。
2 ) A* = 0.01 ( 2.01 + 2.00 ) ,
A2 = 0.01 (1.42 + 1.41) = 0.01 2.83 = 0.00353356
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

1.填空
1) 设()3()1f x x =-,则()f x 关于[0,1]C
的f ¥= 1 ,1f = 14, 2f = 17。

2) 设A 是正定矩阵,则A 的cholesky 的分解 唯一 (唯一,不唯一).
3) 用梯形公式计算积分23
2x e dx -≈⎰ 9.219524E-003:此值比实际值 大 (大,小)
4) 用Euler 方法解初值问题 0(0)1
y y y ì¢-=ïïíï=ïî的近似解的最终表达式n y =()1n h + (取步长x h n =);当n 时,lim n n y = x e 。

5) 令f(x)=3x 7+ x 4+3x+1, 则f[20, 21,…,27]= 3 ;f[20, 21,…,28]= 0
2.用牛顿法求解Leonardo 方程
02010223=-++x x x 要求 6110-+<-k k x x 。

解答 由上题知,0)(=x f 在(1,2)内有一个根,且0)2(,0)(>>''f x f ,故应取20=x ,利用牛顿迭代公式 ,2,1,0),(/)(1='-=+k x f x f x x k k k k 计算结果如下: k k x k k x 0
1 3 1.368869419 1
1.6 4 1.368808109 2
1.383388704 5 1.368808108 684510101.0--<⨯=-x x ,故取368808108.1*5=≈x x 。

3.给定方程组
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--111 211111112321x x x 证明Jacobi 方法发散。

证明 对Jacobi 方法,迭代矩阵为 ⎥⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎢⎣⎡---=0212110121210B 设其特征值为λ,则 i B I 25,0,045
3,213±
===+=-λλλλλ 12/5)(>=B ρ,故Jacobi 方法发散。

4.若{ϕ0(x ), ϕ1(x ),…, ϕn (x )}是[a,b]上的正交族。

∑==n k k k
x a x 0)()(ϕϕ为f(x)的最佳平方逼近。

证明: 10 n k f a k k k k ,,,)
,(),( ==ϕϕϕ 证:法方程系数矩阵为Q T Q=()()()()()()()()()⎪⎪⎪⎪⎪⎭⎫ ⎝⎛n n n n n n ϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕ,,,,,,,,,10
1110101000 =()()()⎪⎪⎪⎪⎪⎭⎫ ⎝
⎛n n ϕϕϕϕϕϕ,000,00
0,1100 此时法方程为 ⎪⎪⎪⎪⎪⎭
⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫ ⎝⎛),(),(),(),(),(),(n n n n f f f a a a ϕϕϕϕϕϕϕϕϕ 10101100故 ,,1,0 ),(),(n k y a k k k k ==ϕϕϕ 5.设l k (x )是关于互异节点x 0, x 1,…, x n , 的Lagrange 插值基函数,证明
(1) m n
k k m k
x x l x =∑=0)( m=0,1,…,n
(2) ∑=-n k k m k x l x x 0)()(≡0 m=1,2,…,n
证明:由插值唯一性定理知(1)。

展开知(2)
6、 试确定常数A ,B ,C 和a ,使得数值积分公式
22()()(0)()f x dx Af a Bf Cf a -≈-++⎰
有尽可能高的代数精度。

试问所得的数值积分公式代数精度是多少?它是否为Gauss 型? 解 由()1f x = 得 4A B C ++=
由()f x x = 得 0aA aC -+=
由2()f x x = 得 22163a A a C +=
由3()f x x = 得 330a A a C -+=
由4()f x x = 得 44645
a A a C += 可 得101612,,995A C B a ====? 代数精度是5, 是Gauss 型积分公式。

相关文档
最新文档