数值分析第一章检测题
- 格式:doc
- 大小:29.00 KB
- 文档页数:2
数值分析习题集及答案数值分析习题集及答案篇一:数值分析习题与答案第一章绪论习题一1.设x>0,x*的相对误差为δ,求f(x)=ln x的误差限。
解:求lnx 的误差极限就是求f(x)=lnx的误差限,由公式()有已知x*的相对误差,故即2.下列各数都是经过四舍五入得到的近似值,试指出它们有几位有效数字,并给出其误差限与相对误差限。
解:直接根据定义和式()()则得有5位有效数字,其误差限有2位有效数字,有5位有效数字,3.下列公式如何才比较准确?(1)(2),相对误差限满足,而解:要使计算较准确,主要是避免两相近数相减,故应变换所给公式。
(1)(2)4.近似数x*=,是 3位有数数字。
5.计算四个选项:取,利用:式计算误差最小。
第二、三章插值与函数逼近习题二、三 1.给定的数值表用线性插值与二次插值计算的近似值并估计误差限. 解:仍可使用n=1及n=2的Lagrange插值或Newton插值,并应用误差估计()。
线性插值时,用及两点,用Newton插值误差限,因,故二次插值时,用,,三点,作二次Newton 插值误差限,故2. 在-4≤x≤4上给出的等距节点函数表,若用二次,函数表的步长h插值法求的近似值,要使误差不超过应取多少? 解:用误差估计式(),令因得3. 若,求和.解:由均差与导数关系于是4. 若的值,这里p≤n+1.解:可知当而当P=n+1时于是得有互异,求,由均差对称性5.求证.解:解:只要按差分定义直接展开得6.已知的函数表求出三次Newton均差插值多项式,计算f()的近似值并用均差的余项表达式估计误差. 解:根据给定函数表构造均差表由式()当n=3时得Newton均差插值多项式N3(x)=+()+()() 由此可得f() N3()= 由余项表达式()可得由于7. 给定f(x)=cosx的函数表篇二:数值分析试题1参考答案参考答案 1 一、1.22.xn?1?xn?3.1, 0 4.7,f(xn)(n?0,1,?) ?f(xn)25 7(k?1)15(k)x2x11336. ? ,1(k?1)x2??x1(k?1)1220?2003??10?2?4二、(1) L?0?13?00?1??(2)1?0?120??,U??0100?5??4000?23100??0?? 3??4?1??l65?a65?(l61u15?l62u25?l63u35?l64u45);u55u56?a55?(l51u16?l52u26?l53u356?l54u46)三、先造差分表如下:(1)选x1?,x2?,x3?,x4?为节点,构造三次向前Newton插值多项式2y1?3y1N(x?th)?y1??y1?t(t?1)?t(t?1)(t?2) 31 2!3!将x1和h代入上式,则有N3(?)?25?2t?1/2*t(t?1)?5/6*t(t?1)(?2)由??解得t?,所以f()?N()?(2) 选x3?,x4?,x5?为节点,构造二次向前Newton插值式N2(x3?th)?y3??y3t?t(t?1)2!将x3和h代入上式,则有N2(?)?20?t?t(t?1) 由+=解得t=,所以 f()?N2()?(3)由f(?)3ht(t?1)(t?2)3!(,0?t?2)R2(x0?th)?f(?)3600有R(2(xi?)?(t?1)(t?2)?**maxt(t?1)(t?2)0?t?23!3!可知f(x)有两位整数,故能保证有两位有效数字。
《数值分析》第⼀章答案习题11.以下各表⽰的近似数,问具有⼏位有效数字?并将它舍⼊成有效数。
(1)*1x =451.023, 1x =451.01;(2)*2x =-0.045 113, 2x =-0.045 18;(3)*3x =23.421 3, 3x =23.460 4;(4)*4x =31, 4x =0.333 3;(5)*5x =23.496, 5x =23.494;(6)*6x =96×510, 6x =96.1×510;(7)*7x =0.000 96, 7x =0.96×310-;(8)*8x =-8 700, 8x =-8 700.3。
解:(1) =*1x 451.023 =1x 451.01=-1*1x x 0.01311021-?≤,1x 具有4位有效数字。
→1x 451.0(2) -=*2x 0.045 113 -=2x 0.045 18=-241021x x 0.045 18045113.0-=0.000 06731021-?<2x 具有2位有效数字,045.02-→x(3)=*3x x =-4604.234213.23=-4213.234604.231 10210391.0-?≤3x 具有3位有效数字,4.233→x (不能写为23.5) (4) =*4x 31 ,=4x 0.3333=-4*4x x 41021000033.0-?<,4x 具有4位有效数字,=4x 0.3333(5) =*5x 23.496,=5x 23.494=-5*5x x =-494.23496.2321021002.0-?<5x具有4位有效数字,→5x 23.50 (不能写为23.49)(6) =*6x 51096?710961.0?==-6*6x x 710001.0-?72101021--??≤6x 具有2位有效数字,57610961096.0?=?=x(7) =*7x 0.00096 371096.0-?=x3*71096.0-?=x =-7*7x x 0 7x 精确(8) 8700*8-=x 8x 3.8700-=8*8x x -010213.0?≤=8x 具有4位有效数字,8x 8700-=精确2.以下各数均为有效数字: (1) 0.1062 + 0.947; (3)2.747?6.83; (2)23.46―12.753; (4)1.473 / 0.064 。
习题11. 以下各表示的近似数,问具有几位有效数字?并将它舍入成有效数。
(1)*1x =451.023, 1x =451.01; (2)*2x =-0.045 113, 2x =-0.045 18; (3)*3x =23.421 3, 3x =23.460 4;(4)*4x =31,4x =0.333 3;(5)*5x =23.496, 5x =23.494; (6)*6x =96×510, 6x =96.1×510; (7)*7x =0.000 96, 7x =0.96×310-; (8)*8x =-8 700, 8x =-8 700.3。
解:(1) =*1x 451.023 =1x 451.01=-1*1x x 0.01311021-⨯≤,1x 具有4位有效数字。
→1x 451.0(2) -=*2x 0.045 113 -=2x 0.045 18=-<⨯-2*241021x x 0.045 18045113.0-=0.000 06731021-⨯<2x 具有2位有效数字,045.02-→x(3)=*3x 23.4213 =3x 23.4604=-3*3x x =-4604.234213.23=-4213.234604.23110210391.0-⨯≤3x 具有3位有效数字,4.233→x (不能写为23.5)(4) =*4x 31,=4x 0.3333=-4*4x x 41021000033.0-⨯< ,4x 具有4位有效数字,=4x 0.3333 (5) =*5x 23.496,=5x 23.494=-5*5x x =-494.23496.2321021002.0-⨯<5x 具有4位有效数字, →5x 23.50 (不能写为23.49)(6) =*6x 51096⨯71096.0⨯= =6x 5101.96⨯710961.0⨯==-6*6x x 710001.0-⨯72101021--⨯⨯≤6x 具有2位有效数字,57610961096.0⨯=⨯=x(7) =*7x 0.00096 371096.0-⨯=x 3*71096.0-⨯=x =-7*7x x 0 7x 精确 (8) 8700*8-=x 8x 3.8700-=8*8x x -010213.0⨯≤= 8x 具有4位有效数字,8x 8700-=精确 2.以下各数均为有效数字: (1) 0.1062 + 0.947; (3)2.747⨯6.83;(2)23.46―12.753; (4)1.473 / 0.064 。
数值分析复习试题第一章 绪论 一. 填空题 1.*x为精确值x 的近似值;()**x f y =为一元函数()x f y =1的近似值;()**,*y x f y =为二元函数()y x f y ,2=的近似值,请写出下面的公式:**e x x =-:***r x xe x -=()()()*'1**y f x x εε≈⋅ ()()()()'***1**r r x f x y x f x εε≈⋅()()()()()**,**,*2**f x y f x y y x y x yεεε∂∂≈⋅+⋅∂∂()()()()()****,***,**222r f x y e x f x y e y y x y y y ε∂∂≈⋅+⋅∂∂ 2、 计算方法实际计算时,对数据只能取有限位表示,这时所产生的误差叫 舍入误差 。
3、 分别用2.718281,2.718282作数e 的近似值,则其有效数字分别有6 位和7 位;又取 1.73≈-211.73 10 2≤⨯。
4、 设121.216, 3.654x x ==均具有3位有效数字,则12x x 的相对误差限为0.0055 。
5、 设121.216, 3.654x x ==均具有3位有效数字,则12x x +的误差限为0.01 。
6、 已知近似值 2.4560A x =是由真值T x 经四舍五入得到,则相对误差限为0.0000204 .7、递推公式,⎧⎪⎨⎪⎩0n n-1y =y =10y -1,n =1,2,如果取0 1.41y ≈作计算,则计算到10y 时,误差为8110 2⨯;这个计算公式数值稳定不稳定 不稳定 . 8、精确值 14159265.3*=π,则近似值141.3*1=π和1415.3*2=π分别有 3位和 4 位有效数字。
9、若*2.71828x e x =≈=,则x 有 6 位有效数字,其绝对误差限为1/2*10-5。
10、 设x*的相对误差为2%,求(x*)n的相对误差0.02n11、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字;12、计算方法主要研究( 截断 )误差和( 舍入 )误差; 13、为了使计算 ()()2334610111y x x x =++---- 的乘除法次数尽量地少,应将该表达式改写为11,))64(3(10-=-++=x t t t t y ,为了减少舍入误差,应将表达式19992001-改写为199920012+。
习题(一)1. 指出四舍五入得到的下列各数有几位有效数字:x 1∗=7.8673,x 2∗=8.0916,x 3∗=0.06213,x 4∗=0.07800,x 5∗=90×103,x 6∗=2.0×10−4解:由有效数字定义得:x 1∗,x 2∗具有5位有效数字x 3∗,x 4∗具有4位有效数字x 5∗,x 6∗具有2位有效数字.2. 设准确值为x=3.78695,y=10,它们的近似值分别为x 1∗=3.7869,x 2∗=3.7870及y 1∗=9.9999,y 2∗=10.1,y 3∗=10.0001,试分析x 1∗,x 2∗,y 1∗,y 2∗,y 3∗分别具有几位有效数字. 解:x 1∗=3.7869=x 1∗=0.37869×101,k 1=1|x 1∗−x|=|3.7869−3.78695|=0.00005≤0.5×10−4=0.5×101−5, 即x 1∗具有5位有效数字;同理,x 2∗=3.7870=0.37870×101,k 2=1|x 2∗−x|=|3.7870−3.78695|=0.00005≤0.5×101−5,所以x 2∗具有5位有效数字; 将y 1∗,y 2∗,y 3∗分别写成y=±10k ×0.α1α2...αn 的表示形式,有:y 1∗=9.9999=0.99999×101,k 3=1;y 2∗=10.1=y 2∗=0.101×102,k 4=2;y 3∗=10.0001=0.100001×102,k 5=2;|y 1∗−y |=|9.9999−10|=0.0001=0.1×10−3≤0.5×101−4,n=4;|y 2∗−y |=|10.1−10|=0.1≤0.5×102−2,n=2;|y 3∗−y |=|10.0001−10|=0.0001=0.1×10−3≤0.5×102−5,n=5;所以y 1∗,y 2∗,y 3∗分别具有4,2,5位有效数字.8.为了使√11的近似值的相对误差不超过0.1%,问至少应取几位有效数字. 解:√11=0.3316624…=0.α1α2...αn ×10k ,α1=3,设x ∗有n 位有效数字,又因为|E x ∗|比值比较小, 故可用E r ∗(x ∗)= |E(x ∗)x ∗|代替相对误差E r ∗(x ∗),用εr ∗=εx ∗代替相对误差限εr 所以εr ∗≤12α1×10−n+1=16×10−n+1 令16×10−n+1≤0.1%,解得n ≥3.22即至少应取4位有效数字.12.如何计算下列函数值才比较精确.(1)11+2x −11+x ,对|x|≪1; (2)√x +1x −√x −1x ,对x ≫1;(3)∫dx 1+x 2N+1N,其中N 充分大; (4)1−cos xsin x ,对|x|≪1;(5)ln(30−√302−1)(开平方用6位函数表);解:(1)原式=1+x−(1+2x)(1+2x)(1+x)=−x (1+2x)(1+x); (2)原式=x+1x −(x−1x )√x+1x +√x−1x =2x √x+1x +√x−1x ;(3)原式=arc tan x|NN+1=arc tan N +1−arc tan N =arc tan N+1−N 1+N(N+1)=arc tan 11+N(N+1); (4)原式=2sinx 222sin x 2cos x 2=tan x2; (5)原式=30+√302−1=−ln(30+√302−1)令f(x)=ln(x −√x 2−1),则f(30)=ln(30−√302−1)=ln(30−√899),记a=30−√899 若用6位开方函数表,则有a ∗=30−29.9833=0.0167,故有ε(a ∗)=0.5×10−4, 而f(30)≈ln a ∗,于是ε(f (30))=ε(ln a ∗)≈|1a ∗|ε(a ∗)=0.50.0167×10−4≈0.003; 又因为f(x)等价于f(x)=-ln(x +√x 2−1),则f (30)=-ln(30+√899),记b=30+√899 同理b ∗=59.9833,进而ε(b ∗)=(2×10−4)−1,对f (30)≈ln b ∗ε(f (30))=ε(ln b ∗)≈|1b ∗|ε(b ∗)=0.559.9833×10−4≈0.834×10−6。
可编辑修改精选全文完整版数值分析-第1章1.填空题(1)为便于算法在计算机上实现,必须将一个数学问题分解为有限次的四则运算;(2)在数值计算中为避免损失有效数字,尽量避免两个相近数作减法运算;为避免误差的扩大,也尽量避免分母的绝对值远小于分子的绝对值;(3)误差有四大来源,数值分析主要处理其中的截断误差和舍入误差;(4)有效数字越多,相对误差越小;2. 用例1.4的算法计算10,迭代3次,计算结果保留4位有效数字。
//见P4解题思路:假定x0是√a的一个近似值,x0>0,则ax0也是√a的一个近似值,且x0和ax0两个近似值必有一个大于√a,另一个小于√a,设想它们的平均值应为√a的更好的近似值,于是x k+1=1 2(x k+ax k),k=0,1,2,……解:取x0=3,按算法x k+1=12(x k+ax k),k=0,1,2,……迭代3次有:x1=12(x0+10x0)=(3+103)≈3.167x2=12(x1+10x1)=(3.167+103.167)≈3.162x3=12(x2+10x2)=(3.162+103.162)≈3.1623. 推导开平方运算的误差限公式,并说明什么情况下结果误差不大于自变量误差。
//见P8解:已知f(x)=√x,设x∗是准确值,令x是x∗的一个近似值,则相对误差e(f(x))=f(x)−f(x∗),由Taylor公式f(x∗)=f(x)0! +f′(x)1!(x∗−x)+f"(x)2!(x∗−x)2+⋯+f n(x)n!(x∗−x)n+R n(x)其中,R n(x)=f n+1(ξ)(n+1)!(x∗−x)n+1将f(x∗)展开分析有:f(x∗)=√x2√x x∗−x)+⋯+f n(ξ)n!(x∗−x)n+R n(x)∴e(f(x))=f(x)−f(x∗)=− (2√x x∗−x)+⋯+f n(ξ)n!(x∗−x)n+R n(x))∴|e(f(x))|≤ ε(f(x))≤|2√x |ε(x)+⋯+|f n(ξ)n!εn(x)|+|R n(x)|忽略二阶以上无穷小,可得f(x)的误差限公式为ε(f(x))≈2√x(x)。
数值分析习题(含答案)第一章绪论姓名学号班级习题主要考察点:有效数字的计算、计算方法的比较选择、误差和误差限的计算。
1 若误差限为5105.0-?,那么近似数0.003400有几位有效数字?(有效数字的计算)解:2*103400.0-?=x ,325*10211021---?=?≤-x x 故具有3位有效数字。
2 14159.3=π具有4位有效数字的近似值是多少?(有效数字的计算)解:10314159.0?= π,欲使其近似值*π具有4位有效数字,必需41*1021-?≤-ππ,3*310211021--?+≤≤?-πππ,即14209.314109.3*≤≤π即取(3.14109 , 3.14209)之间的任意数,都具有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 ba 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π=的绝对误差限与相对误差限。
数值分析习题集(适合课程《数值方法A 》和《数值方法B 》)长沙理工大学第一章 绪 论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2%,求nx 的相对误差.3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====⨯4. 利用公式(3.3)求下列各近似值的误差限:********12412324(),(),()/,i x x x ii x x x iii x x ++其中****1234,,,x x x x 均为第3题所给的数.5. 计算球体积要使相对误差限为1%,问度量半径R 时允许的相对误差限是多少?6. 设028,Y =按递推公式1n n Y Y -= ( n=1,2,…)计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差?7. 求方程25610x x -+=的两个根,使它至少具有四位有效数字≈27.982).8. 当N 充分大时,怎样求211Ndx x +∞+⎰?9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2?10. 设212S gt =假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加,而相对误差却减小.11. 序列{}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字),计算到10y 时误差有多大?这个计算过程稳定吗?12. 计算61)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好?3--13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式ln(ln(x x =-计算,求对数时误差有多大?14. 试用消元法解方程组{101012121010;2.x x x x +=+=假定只用三位数计算,问结果是否可靠?15. 已知三角形面积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量a ,b ,c 的误差分别为,,.a b c ∆∆∆证明面积的误差s ∆满足.s a b c s a b c ∆∆∆∆≤++第二章 插值法1. 根据(2.2)定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n nn n n n x x x V x V x x x x x x x xx x ----==证明()n V x 是n 次多项式,它的根是01,,n x x -,且 101101()(,,,)()()n n n n V x V x x x x x x x ---=--.2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的二次插值多项式.3.4. 给出cos x ,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界.5. 设0k x x kh =+,k =0,1,2,3,求032max ()x x x l x ≤≤.6. 设jx 为互异节点(j =0,1,…,n ),求证:i) 0()(0,1,,);nk kj jj x l x x k n =≡=∑ii) 0()()1,2,,).nk jj j xx l x k n =-≡0(=∑7. 设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8max max a x ba xb f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()xf x e =的等距节点函数表,若用二次插值求xe 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少? 9. 若2nn y =,求4n y ∆及4n y δ. 10. 如果()f x 是m 次多项式,记()()()f x f x h f x ∆=+-,证明()f x 的k 阶差分()(0)k f x k m ∆≤≤是m k -次多项式,并且()0(m l f x l +∆=为正整数).11. 证明1()k k k k k k f g f g g f +∆=∆+∆. 12. 证明110010.n n kkn n k k k k f gf g f g g f --+==∆=--∆∑∑13. 证明1200.n j n j y y y -=∆=∆-∆∑14. 若1011()n n n n f x a a x a x a x --=++++有n 个不同实根12,,,n x x x ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑15. 证明n 阶均差有下列性质: i)若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x =;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+.16. 74()31f x x x x =+++,求0172,2,,2f ⎡⎤⎣⎦及0182,2,,2f ⎡⎤⎣⎦.17. 证明两点三次埃尔米特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔米特插值的误差限.18. 求一个次数不高于4次的多项式()P x ,使它满足(0)(1)P P k =-+并由此求出分段三次埃尔米特插值的误差限. 19. 试求出一个最高次数不高于4次的函数多项式()P x ,以便使它能够满足以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.20. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造一个台阶形的零次分段插值函数()n x ϕ并证明当n →∞时,()n x ϕ在[],a b 上一致收敛到()f x .21. 设2()1/(1)f x x =+,在55x -≤≤上取10n =,按等距节点求分段线性插值函数()h I x ,计算各节点间中点处的()h I x 与()f x 的值,并估计误差.22. 求2()f x x =在[],a b 上的分段线性插值函数()h I x ,并估计误差.23. 求4()f x x =在[],a b 上的分段埃尔米特插值,并估计误差.试求三次样条插值并满足条件i) (0.25) 1.0000,(0.53)0.6868;S S '='= ii)(0.25)(0.53)0.S S "="=25. 若[]2(),f x C a b ∈,()S x 是三次样条函数,证明i)[][][][]222()()()()2()()()bbbbaaaaf x dx S x dx f x S x dx S x f x S x dx"-"="-"+""-"⎰⎰⎰⎰;ii) 若()()(0,1,,)i i f x S x i n ==,式中i x 为插值节点,且01n a x x x b =<<<=,则[][][]()()()()()()()()()b aS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'⎰.26. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可用(8.7)式的表达式).第三章 函数逼近与计算1. (a)利用区间变换推出区间为[],a b 的伯恩斯坦多项式.(b)对()sin f x x =在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做比较. 2. 求证:(a)当()m f x M ≤≤时,(,)n m B f x M ≤≤. (b)当()f x x =时,(,)n B f x x =. 3. 在次数不超过6的多项式中,求()sin 4f x x =在[]0,2π的最佳一致逼近多项式. 4. 假设()f x 在[],a b 上连续,求()f x 的零次最佳一致逼近多项式. 5. 选取常数a ,使301max x x ax≤≤-达到极小,又问这个解是否唯一?6. 求()sin f x x =在[]0,/2π上的最佳一次逼近多项式,并估计误差.7. 求()xf x e =在[]0,1上的最佳一次逼近多项式.8. 如何选取r ,使2()p x x r =+在[]1,1-上与零偏差最小?r 是否唯一?9. 设43()31f x x x =+-,在[]0,1上求三次最佳逼近多项式.10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x .11. 试证{}*()nTx 是在[]0,1上带权ρ=的正交多项式.12. 在[]1,1-上利用插值极小化求11()f x tg x -=的三次近似最佳逼近多项式.13. 设()xf x e =在[]1,1-上的插值极小化近似最佳逼近多项式为()n L x ,若nf L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ϕ=-----,试将()x ϕ降低到3次多项式并估计误差.15. 在[]1,1-上利用幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过0.005.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数.17. 求a 、b 使[]220sin ax b x dx π+-⎰为最小.并与1题及6题的一次逼近多项式误差作比较.18. ()f x 、[]1(),g x C a b ∈,定义 ()(,)()();()(,)()()()();b baaa f g f x g x dxb f g f x g x dx f a g a =''=''+⎰⎰问它们是否构成内积?19. 用许瓦兹不等式(4.5)估计6101x dx x +⎰的上界,并用积分中值定理估计同一积分的上下界,并比较其结果.20. 选择a ,使下列积分取得最小值:1122211(),x ax dx x ax dx----⎰⎰.21. 设空间{}{}10010121,,,span x span x x 1ϕ=ϕ=,分别在1ϕ、2ϕ上求出一个元素,使得其为[]20,1x C ∈的最佳平方逼近,并比较其结果.22. ()f x x =在[]1,1-上,求在{}2411,,span x x ϕ=上的最佳平方逼近.23.sin (1)arccos ()n n x u x +=是第二类切比雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin2f x x=在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.26.2y a bx =+.27.用最小二乘拟合求.29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出一张记录{}{}4,3,2,1,0,1,2,3k x =,试用改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =第四章 数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度:(1)101()()(0)()hhf x dx A f h A f A f h --≈-++⎰; (2)21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰;(3)[]1121()(1)2()3()/3f x dx f f x f x -≈-++⎰;(4)[][]20()(0)()/1(0)()hf x dx h f f h ah f f h ≈++'-'⎰.2. 分别用梯形公式和辛普森公式计算下列积分:(1)120,84xdx n x =+⎰; (2)1210(1),10x e dx n x --=⎰;(3)1,4n =⎰;(4),6n =.3. 直接验证柯特斯公式(2.4)具有5次代数精度.4. 用辛普森公式求积分10x e dx-⎰并计算误差. 5. 推导下列三种矩形求积公式:(1)2()()()()()2ba f f x dxb a f a b a 'η=-+-⎰; (2)2()()()()()2baf f x dx b a f b b a 'η=---⎰; (3)3()()()()()224baa b f f x dx b a f b a +"η=-+-⎰. 6. 证明梯形公式(2.9)和辛普森公式(2.11)当n →∞时收敛到积分()baf x dx⎰.7. 用复化梯形公式求积分()baf x dx⎰,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍入误差)?8.1x e dx-,要求误差不超过510-.9. 卫星轨道是一个椭圆,椭圆周长的计算公式是S a =θ,这里a 是椭圆的半长轴,c 是地球中心与轨道中心(椭圆中心)的距离,记h 为近地点距离,H 为远地点距离,6371R =公里为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第一颗人造卫星近地点距离439h =公里,远地点距离2384H =公里,试求卫星轨道的周长.10. 证明等式3524sin3!5!n n nnππππ=-+-试依据sin(/)(3,6,12)n n n π=的值,用外推算法求π的近似值.11. 用下列方法计算积分31dyy ⎰并比较结果.(1) 龙贝格方法;(2) 三点及五点高斯公式;(3) 将积分区间分为四等分,用复化两点高斯公式.12. 用三点公式和五点公式分别求21()(1)f x x =+在x =1.0,1.1和1.2处的导数值,并估计误()f x第五章 常微分方程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉方法和改进的尤拉方法的近似解的表达式,并与准确解bx ax y +=221相比较。
数值分析第一章检测题
60分钟
一、填空题(每空4')
1.已知x=6
2.1341是由准确数a经四舍五入得到的a的近似值,试
给出x的绝对误差限_______________.
2.设x和y的相对误差均为0.001,则x*y的相对误差约为
____________.
3.π取四位有效数字是,e取五位有效数字是
,此时π/e的误差限是。
4.根据秦九韶算法,多项式p(x)=5x^4+3x^2+12x+4只需计算
次乘法和次加法。
5.要使30的近似值的相对误差小于0.1%,至少要取位有
效数字。
二、简答题(10'×4)
1.科学计算中的误差来源有哪些?
2.截断误差和舍入误差的区别。
3.用公式表示误差、误差限;相对误差、相对误差限。
4.避免误差危害,防止有效数字损失的方法。
三、解答题(8'×3)
1. 用秦九韶算法求多项式p (x )=2x^6+3x^4-9x^2+x-6在x=3处的
值。
2. 计算球面积要使相对误差限为0.1%,问度量半径R 所允许的相
对误差是多少?
3. 取99的6位有效数字9.94987,则以下两种算法各有几位有效数
字?
05013.094987.9109910=-≈-
0501256399.094987
.19194987.910199101==+≈+。