13级研究生数值分析习题
数值分析课后习题答案

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
武汉大学研究生数值分析考试试题

武 汉 大 学2013~2014学年第一学期硕士研究生期末考试试题科目名称:数值分析 学生所在院: 学号: 姓名:1、(10分)已知方程 32cos 120x x --=.(1)估计出含根的区间;(2)讨论迭代格式 12cos 43n n x x +=+ 的收敛性; (3)写出求解此方程的牛顿迭代格式,并讨论初值0x 取何值时迭代收敛。
2、(10分)用杜利特尔(Doolittle )分解算法求解方程 b Ax =,其中1514227135A 轾犏犏=犏犏臌 146116b 轾犏犏=犏犏臌3、(14分)设线性方程组Ax b =的系数矩阵()ij n n A a ´=,0,1,,ii a i n ?L .()12,,,Tn b b b b =L (1)分别写出Jacobi 迭代格式及 Gauss-Seidel 迭代格式;(2)证明:Gauss-Seidel 迭代格式收敛的充分必要条件是方程1112121222120nnn n nn a a a a a a a a a l l l l l l =L L M M ML 的根的模1l <.4、(12分)已知 )(x f y = 的数据如下:求)(x f 的Hermite 插值多项式)(3x H 及其余项。
5、(10 求形如 2y a bx cx =++ 的拟合曲线。
6、(10分)确定常数 a ,b 的值,使积分[]220(,)sin I a b a bx x dx p =+-ò 取得最小值。
7、(10分)已知三次Legendre(勒让德)多项式331()(53)2L x x x =-,[1,1]x ? 试确定常数,(1,2,3)i i A x i =,使求积公式 31122333()()()()f x dx A f x A f x A f x -?+ò有尽可能高的代数精度,并判断它是否为高斯型求积公式。
8、(12分)对于下面求解常微分方程初值问题 ⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy 的单步法:112121(3)4(,)22(,)33n n n n n n h y y k k k f x y k f x h y hk +ìïï=++ïïïï=íïïïï=++ïïïî(1)验证它与微分方程相容;(2)确定此单步法的绝对稳定域9、(12分)设初值0x充分接近*x =0a >为常数),证明:迭代格式212(3),0,1,3n n n n x x a x n x a++==+L 三阶收敛于*x,并求lim n x -参考答案(2014-1-10)1、 含根区间:[π,4];因为2()sin 13g x x ¢=<,所以迭代收敛; 在含根区间[π,4],f 的一阶导数恒正,二阶导数恒负,所以取初值为π时牛顿必收敛。
数值分析习题集及答案

数值分析习题集及答案数值分析习题集(适合课程《数值⽅法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. 利⽤公式求下列各近似值的误差限:********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 -=…)计算到100Y .(五位有效数字),试问计算100Y 将有多⼤误差7. 求⽅程25610x x -+=的两个根,使它⾄少具有四位有效数字.8. 当N 充分⼤时,怎样求211Ndx x +∞+?9. 正⽅形的边长⼤约为100㎝,应怎样测量才能使其⾯积误差不超过1㎝210. 设212S gt =假定g 是准确的,⽽对t 的测量有±秒的误差,证明当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 cs a b c ≤++第⼆章插值法1. 根据定义的范德蒙⾏列式,令2000011211121()(,,,,)11n n n n n n n n n x x x V x V x x x x x x x xxx ----==证明()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. 给出f(x )=ln x 的数值表⽤线性插值及⼆次插值计算ln 的近似值.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,,);nkkj j j x l x xk n =≡=∑ii)()()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 b a x b 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 上的分段埃尔⽶特插值,并估计误差. 24. 给定数据表如下:试求三次样条插值并满⾜条件 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 =<<<=,则[][][]()()()()()()()()()baS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'?.26. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可⽤式的表达式).第三章函数逼近与计算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. 试证{}*()nT x 是在[]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次逼近多项式,使误差不超过.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数.17. 求a 、b 使[]22sin 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. ⽤许⽡兹不等式估计601x 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()sin 2f x x=在[]1,1-上按勒让德多项式及切⽐雪夫多项式展开,求三次最佳平⽅逼近多项式并画出误差图形,再计算均⽅误差.25. 把()arccos f x x =在[]1,1-上展成切⽐雪夫级数.26. ⽤最⼩⼆乘法求⼀个形如2y a bx =+的经验公式,使它与下列数据拟合,并求均⽅误差.28. 在某化学反应⾥,根据实验所得分解物的浓度与时间关系如下:⽤最⼩⼆乘拟合求.29. 编出⽤正交多项式做最⼩⼆乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出⼀张记录{}{}4,3,2,1,0,1,2,3k x =,试⽤改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =第四章数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量⾼,并指明所构造出的求积公式所具有的代数精度:101()()(0)()hh f x dx A f h A f A f h --≈-++?;(2)21012()()(0)()hhf 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. 直接验证柯特斯公式具有5次代数精度.4. ⽤⾟普森公式求积分1xedx-?并计算误差.5. 推导下列三种矩形求积公式:(1)2()()()()()2baf f x dx b a f a b a 'η=-+()()()()2baf f x dx b a f b b a 'η=---?;(3)3()()()()()224baa b f f x dx b a f b a +"η=-+-?.6. 证明梯形公式和⾟普森公式当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 nn n π-试依据sin(/)(3,6,12)n n n π=的值,⽤外推算法求π的近似值.11. ⽤下列⽅法计算积分31dyy ?并⽐较结果.(1) 龙贝格⽅法; (2) 三点及五点⾼斯公式;(3) 将积分区间分为四等分,⽤复化两点⾼斯公式.12. ⽤三点公式和五点公式分别求21()(1)f x x =+在x =,和处的导数值,并估计误差.()f x 的值由下表给出:第五章常微分⽅程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉⽅法和改进的尤拉⽅法的近似解的表达式,并与准确解bx ax y +=221相⽐较。
研究生数值分析期末考试试题A答案

2010年秋研究生数值分析期末考试试题答案一、单选题(4*5=20分)1、D; 2、B ; 3、D ; 4、B ; 5、D 。
二、填空题(4*5=20)1、4; 2、⎪⎪⎪⎭⎫ ⎝⎛323203*⎪⎪⎪⎪⎭⎫⎝⎛320323; 3、)]23()0()23([3f f f ++-∏;4、kk k k x x x x 2221--=+;5、9.605。
三、(10分)由两点三次Hermite 插值多项式公式秋得:)2()(23x x x H -=,设所求多项式223)1()()(-+=x Ax x H x P ,。
(4分) 由P(2)=1,得A=1/4,。
(4分) 故22)3(41)(-=x x x P 。
.。
(2分) 四、(10分)设⎪⎪⎪⎭⎫⎝⎛⎪⎪⎪⎭⎫ ⎝⎛--=1001001*10010021321u u l l l A ,由追赶法公式求得, 15/56,15/4,4/15,4/1,432211=-==-==l u l u l ,。
(4分) 由Ly=d,求得T y )77.0,87.0,25.0(=,(3分) 由Ux=y,求得,T x )5179.0,0714.1,7679.0(=(3分)五、(10分)Jacobi 迭代计算格式:⎪⎩⎪⎨⎧++-=--=--=+++3/)221(5/)327(24)(2)(1)1(3)(3)(1)1(2)(3)(2)1(1k k k k k k k k k x x x x x x x x x 。
(2分) G-S 迭代计算格式: ⎪⎩⎪⎨⎧++-=--=--=++++++3/)221(5/)327(24)1(2)1(1)1(3)(3)1(1)1(2)(3)(2)1(1k k k k k k k k k x x x x x x x x x 。
(2分) 由于016415)(3=-+=-λλλJ B I del ,,11516)(>=J B ρ即Jacobi 迭代发散;。
数值分析课后习题答案

x2 6.6667x2 8.205
再解
1
15 56
x31.785,7得 x35.769
1 25069x4 0.47847x4 1.4872
1 x5 5.3718 x5 5.3718
2-10.证明下列不等式:
(1)x-yx-z+z-y; (2)|x-y|x-y;
证明 (1)x-y=(x-z)+(z-y)x-z+z-y
b.用Gauss消元法
102 x y 1 x y 2
回代得解: y=1, x=0.
102 x Байду номын сангаасy 1
100y 100
再用列主元Gauss消元法
102 x y 1 x y 2
回代得解: y=1, x=1.
x y
y 1
2
2-8.用追赶法求解方程组:
4 1
x1 100
1 4 1
x2 0
3-8.判定求解下列方程组的SOR方法的收敛性.
2 1 0 0 x1 1
1
0 0
2 1 0
1 2 1
0 12
x2 x3 x4
0 00
解 直接可验证系数矩阵A是负定矩阵,所以-A是对称
1-3.为了使101/2的相对误差小于0.01%,试问应取几位 有效数字?
解 因为101/2=3.162…=0.3162…10,若具有n位有效 数字,则其绝对误差限为0.5 101-n ,于是有
r=0.5101-n/3.162…<0.5101-n/3<0.01% 因此只需n=5.即取101/2=3.1623
1 2
0
12 1,
1 2
1 2
0
12
《数值分析》上机习题

《数值分析》上机习题1.用Newton法求方程X7 - 28X4 + 14=0在(0.1,1.9)中的近似根(初始近似值取为区间端点,迭代6次或误差小于0.00001)。
#include<stdio.h>#include<math.h>int main(){ float x1,x,f1,f2;static int count=0;x1=0.1 ;//定义初始值do{x=x1;f1=x*(x*x*x*x*x*x-28*x*x*x)+14;f2=7*x*x*x*x*x*x-4*28*x*x*x;//对函数f1求导x1=x-f1/f2; count++; }while(fabs(x1-x)>=1e-5);printf("%8.7f\n",x1); printf("%d\n",count);return 0;}2.已知函数值如下表:试用三次样条插值求f(4.563)及f’(4.563)的近似值。
include <stdio.h>#include <math.h>#define N 11main(){double B[N+1][N+1],m,x,u[N+1],y[N+1],c[N+1],d[N+1];double e[N+1]={2,0,4.15888308,6.5916738,8.3177664,9.6566268,10.750557,11.675460 6,12.47667,13.1833476,13.8155106,14.0155106};int i;x=4.563;B[0][0]=-2;B[0][1]=-4;B[N][N-1]=4;B[N][N]=2;for(i=1;i<N;i++){B[i][i-1]=1;B[i][i]=4;B[i][i+1]=1;}u[0]=B[0][0];y[0]=e[0];for(i=1;i<N;i++){m=B[i][i-1]/u[i-1];u[i]=B[i][i]-m*B[i-1][i];y[i]=e[i]-m*y[i-1];}c[N]=y[N]/u[N];for(i=N-1;i>=0;i--)c[i]=(y[i]-B[i][i+1]*c[i+1])/u[i];for(i=0;i<12;i++){m=fabs(x-i);if(m>=2)d[i]=0;else if(m<=1)d[i]=0.5*fabs(pow(m,3))-m*m+2.0/3;elsed[i]=(-1.0/6.0)*fabs(pow(m,3))+m*m-2*fabs(m)+4/3.0;} m=0;for(i=0;i<12;i++) m=m+c[i]*d[i];printf("f(%4.3f)=%f\n",x,m);printf("f'(4.563)=%lf\n",(c[4]-c[2])/2); }3.用Romberg 算法求)00001.0(sin )75(32314.1=+⎰ε允许误差dx x x x x . #include "stdafx.h" #include<stdio.h> #include<math.h> float f(float x) {float f=0.0;f=pow(3.0,x)*pow(x,1.4)*(5*x+7)*sin(x*x); return (f); } main() {int i=1,j,k,n=12;float T[12],a=1.0,b=3.0,s=0.0; T[0]=0.5*(b-a)*(f(a)+f(b));for(j=1;j<n-1;j++){ for(k=1;k<=pow(2,j-1);k++) s+=f(a+(2*k-1)*(b-a)/pow(2,j)); T[j]=0.5*(T[j-1]+(b-a)*s/pow(2,j-1)); s=0.0; }T[11]=(4*T[1]-T[0])/(float)3;for(;fabs(T[11]-T[0])>0.00001;i++) {T[0]=T[11];for(j=1;j<n-1-i;j++) T[j]=(pow(4,i)*T[j+1]-T[j])/(pow(4,i)-1); T[11]=(pow(4,i+1)*T[1]-T[0])/(pow(4,i+1)-1); }printf("%f\n",T[11]); }4. 用定步长四阶Runge-Kutta 求解一、程序要求用定步长四阶法求解 y1’=1y2’=y3 y3’=1000-1000y2-100y3(y1(0)=y2(0)=y3(0)=0) h=0.0005,打印yi(0.025),yi(0.045),yi(0.085),yi(0.1),(i=1,2,3)h =0.0005,打印y i (0.025) ,⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧===--===0)0(0)0(0)0(10010001000//1/321323321y y y y y dt dy ydt dy dt dyy i(0.045) ,y i(0.085) ,y i(0.1) ,(i=1,2,3)#include "stdafx.h"#include <stdio.h>#include <math.h>double F(double x,double y[4],double f[4]){f[1]=0*x+0*y[1]+0*y[2]+0*y[3]+1;f[2]=0*x+0*y[1]+0*y[2]+1*y[3]+0;f[3]=0*x+0*y[1]-1000*y[2]-100*y[3]+1000;return(1);}void main(){double F(double x,double y[4],double f[4]);doubleh=0.0005,x=0,Y[4],k[5][4],s[4],f[4],det,m[4]={0.025,0.045,0.085,0.1};int i,j,t; for(t=0;t<=3;t++)/*龙格-库塔算法*/{for(j=0;j<=3;j++)Y[j]=0; //每求一组值后将初值条件还原为0 for(i=1;i<=int(m[t]/h);i++){for(j=1;j<=3;j++)s[j]=Y[j];det=F(x,s,f);for(j=1;j<=3;j++)k[1][j]=h*f[j]; /*四阶古典公式中的k值和求和的计算*/ for(j=1;j<=3;j++)s[j]=Y[j]+0.5*k[1][j];det=F(x+0.5*h,s,f);for(j=1;j<=3;j++)k[2][j]=h*f[j];for(j=1;j<=3;j++)s[j]=Y[j]+0.5*k[2][j];det=F(x+0.5*h,s,f);for(j=1;j<=3;j++)k[3][j]=h*f[j];for(j=1;j<=3;j++)s[j]=Y[j]+k[3][j];det=F(x+h,s,f);for(j=1;j<=3;j++)k[4][j]=h*f[j];for(j=1;j<=3;j++)Y[j]=Y[j]+(k[1][j]+2*k[2][j]+2*k[3][j]+k[4][j])/6;x+=h;} for(j=1;j<=3;j++)printf("y[%d](%f)=%f ",j,m[t],Y[j]); printf("\n"); } } 5.⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=40.00001 4.446782 2.213474- 0.238417 1.784317 0.037585- 1.010103- 3.123124 2.031743- 4.446782 30.719334 3.123789 1.103456- 2.121314 0.71828- 0.336993 1.112348 3.067813 2.213474- 3.123789 14.7138465 0.103458- 3.111223- 2.101023 1.056781- 0.784165- 1.7423820.238417 1.103456- 0.103458- 9.789365 0.431637 3.741856- 1.836742 1.563849 0.718719 1.7843172.1213143.111223- 0.431637 19.8979184.101011 2.031454 2.189736 0.113584-0.037585- 0.71828- 2.101023 3.741856- 4.101011 27.108437 3.123848 1.012345- 1.112336 1.010103- 0.336993 1.056781- 1.836742 2.031454 3.123848 15.567914 3.125432- 1.061074- 3.123124 1.112348 0.784165- 1.563849 2.189736 1.012345- 3.125432- 19.141823 2.115237 2.031743- 3.067813 1.742382 0.718719 0.113584- 1.112336 1.061074- 2.115237 12.38412A Tb )5.6784392- 4.719345 1.1101230 86.612343- 1.784317 0.84671695 25.173417- 33.992318 2.1874369(= 用列主元消去法求解Ax=b 。
数值分析典型习题2

= ( I − AX 0 )
2k
11/16
X k = A [ I − ( I − AX 0 ) ]
2k
−1
ρ ( I − AX 0 ) < 1
−1
lim ( I − AX 0 )
k →∞
2k
2k
=0
lim X k = lim A [ I − ( I − AX 0 ) ] = A −1
k →∞ k →∞
2 1 2 5 − 2 − 2 3 1 1 − 2 x1 4 x 7 3 − 2 2 = 3 5 x3 − 1 2 3 x4 0
13/16
练习5: 练习 求矩阵的 2-范数,
|| x || − || y || ≤|| x − y ||
证: || x || = || (x – y )+ y || ≤|| x – y || + || y || || x || – || y || ≤|| x – y || 同理, 同理 || y || – || x || ≤|| y – x || =|| x – y || || x || – || y ||≥ – || x – y || – || x – y || ≤ || x || – || y || ≤|| x – y ||
|λ(D– L ) – U | = 0
a 12 ⋯
行列式对应的矩阵为
λa 11 λa 21 C (λ ) = ⋮ λ a n1
λa 22
⋮
λa n 2
a1n ⋯ a 2n ⋱ ⋮ ⋯ λa nn
矩阵的主对角占优性质, 当|λ | > 1时,利用 矩阵的主对角占优性质,得 时 利用A矩阵的主对角占优性质
数值分析原理习题答案

数值分析原理习题答案数值分析原理习题答案数值分析原理是计算数学的一个重要分支,它研究如何利用数值方法解决数学问题。
在学习数值分析原理的过程中,我们常常会遇到一些习题,下面我将为大家提供一些数值分析原理习题的答案,希望对大家的学习有所帮助。
1. 习题:使用二分法求方程f(x) = x^3 - 2x - 5 = 0在区间[1, 2]内的根,并给出根的近似值,要求误差不超过0.001。
答案:首先,我们需要确定方程f(x)在区间[1, 2]内是否有根。
根据介值定理,如果f(1)和f(2)异号,则方程在区间[1, 2]内有根。
计算f(1) = -6和f(2) = 1,可知f(1)和f(2)异号,因此方程在区间[1, 2]内有根。
接下来,我们可以使用二分法求解方程的根。
首先,将区间[1, 2]等分为两个子区间[1, 1.5]和[1.5, 2],计算f(1.5) = -1.375。
由于f(1)和f(1.5)异号,根据介值定理,方程在区间[1, 1.5]内有根。
接下来,我们将区间[1, 1.5]等分为两个子区间[1, 1.25]和[1.25, 1.5],计算f(1.25) = 0.796875。
由于f(1)和f(1.25)异号,根据介值定理,方程在区间[1, 1.25]内有根。
重复以上步骤,我们可以得到方程的近似根为1.365234375,误差小于0.001。
2. 习题:使用拉格朗日插值多项式求解以下问题:已知函数f(x)在区间[-1, 1]上的节点值为f(-1) = -1,f(0) = 0,f(1) = 1,求f(0.5)的近似值。
答案:根据拉格朗日插值多项式的定义,我们可以得到以下公式:f(x) ≈ L(x) = f(-1) * L0(x) + f(0) * L1(x) + f(1) * L2(x)其中,L0(x),L1(x),L2(x)分别为拉格朗日基函数,其定义如下:L0(x) = (x - 0) * (x - 1) / ((-1 - 0) * (-1 - 1)) = x^2 - xL1(x) = (x + 1) * (x - 1) / ((0 + 1) * (0 - 1)) = -x^2 + 1L2(x) = (x + 1) * (x - 0) / ((1 + 1) * (1 - 0)) = x^2 + x将已知节点值代入公式,我们可以得到:f(x) ≈ L(x) = -1 * (x^2 - x) + 0 * (-x^2 + 1) + 1 * (x^2 + x)= -x^2 + x + x^2 + x= 2x将x = 0.5代入公式,我们可以得到f(0.5)的近似值为1。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
13级研究生数值分析习题
第一章 误差及相关问题
内容及纲目:
1) 舍入误差和截断误差
2) 绝对误差和相对误差
3) 误差的传播和计算函数值
4) 算法的数值稳定性
5) 计算中需要注意的问题
1. 用x 近似,sin x 即,sin x x ≈δδ],,0[∈x 最大为多少时,该近似计算的截断误差不超过10-7
. 2. 设,0>x x 的相对误差为δ,求x ln 的绝对误差。
3.的相对误差不超过0.1%,应取几位有效数字?
解:知识点:有效数字和相对误差间的关系。
4,设近视数*x 有n 位有效数字,所以有: *11|()|1024n r e x -≤
⨯⨯,令:11100.1%24n -⨯≤⨯,解得: 3.097,n ≥所以有4位有效数字。
4. 227
作为=3.1415926π有几位有效数字? 5. 误差的来源?计算中需要注意的几个问题.
第二章 函数插值
内容及纲目:
1) 插值多项式的存在性与唯一性
2) 插值多项式的构造方法(lagrange 插值,Newton 插值,等距
节点的插值)
3) 带导数的插值函数构造,Hermite 插值,误差估计和构造方法
4) 差分和差商的定义、性质和联系
5) 三次样条插值公式及误差估计
1. ]2,,2,2,2,[]2,,2,2,2,[,13)(72162147 x f x f x x x x f 和求+++=。
2. 已知12144,11121,10100===,分别用线性插值和抛物插值法,求115的近似值。
3. (分三次Hermite 插值),仅给定10,x x 和相应的函数值10,y y 及其微商10,m m ,构造插值函数)(x H ,)(x H 满足条件:1).)(x H 是不超过三次的多项式;
2). ,)(,)(1100y x H y x H ==1100)(,)(m x H m x H ='='。
4. 构造 不超过3次的插值多项式,使其满足:.3)1(;0)2(,2)1(,1)0(='===f f f f
5. 设f(x) ∈C 2[a,b],且)(a f = )(b f =0,求证:b x a x f ≤≤|)(|max )(81a b -≤ 2 b
x a x f ≤≤|)(''|max 。
6. 已知插值节点数据(-2,17),(0,1),(1,2),(2,17),构造三次Lagrange 插值多项式)(3x L .
7. 依据下列函数表建立不超过三次的Lagrange 插值多项式及Newton 插值多项式,并验证插值多项式的唯一性。
8.求一个次数不超过3的多项式3()p x ,满足条件:(1)2,(2)4,(3)12,(2)3f f f f '====,并求插值误差。
9. 求一个次数不高于4的多项式4()P x ,使它满足1)2(,1)1()1(,0)0()0(44444=='=='=P P P P P ,并求插值误差。
10. 求满足条件1)2(,1)1(,3)2(,2)1(-='='==f f f f 的Hermite 插值多项式,并求插值误差。
11.给定数据(0)0,(1)1,(0)0,(1)1,(1)0,f f f f f ''''=====构造次数不高于4的多项式4()P x 及插值余项。
12给出数据(1)2,(0)3,(1)4,(f f f f -====试求三次样条函数()S x ,并满足(1)5,(3)29S S ''-==。
并计算(2)?S =
13. 设x 0, x 1,…… x n 为n +1个相异的插值点,)(x l i (i=0,1,2,……,n)为Lagrange 插值基函数, 证明:(1)∑≡n i i x l ;1)((2) ),,2,1,0(,)(n k x x x l k k
i n
i i ==∑ 第三章 函数与数据的逼近
内容及纲目:
1) 范数、内积 定义
2) 正交多项式的性质及构造
3) 常用的正交多项式
4) 离散数据的最小二乘曲线拟合
1.求,a b 值使得 ⎰--2
0]sin [π
x b ax 2达到最小.
2.给定一组数据,,,2,1),,(m i y x i i =用最小二乘法做拟合直线bx a x p +=)(,求出.,b a
3.求[0,1]区间上关于权函数是x x ln )(-=ρ的正交多项式系的前三项。
4.给出关于点集{-2,-1,0,1,2}的首项系数为1的正交多项式系的前四项。
5、判定函数3
1,,12-
x x 在]1,1[-上两两正交,并求一个三次多项式,使其在]1,1[-上与上述函数两两正交 6. 1) 用最小二乘法求解矛盾方程组:⎪⎩
⎪⎨⎧===+-1y -2x -94y -x 13y 2x
2)求解下列矛盾方程组的最小二乘解
⎪⎪⎩⎪⎪⎨⎧=+-=-=+-=-433
22212
12
12121x x x x x x x x 7. 在]1,25.0[的所有连续函数的集合]1,25.0[C 中,给定,)(x x f =子集},,1{x span =φ对于]1,25.0[)(),(C x g x f ∈,定义内积⎰
=125.0)()(),(dx x g x f g f ,试在φ中寻找一个线性函数bx a +,使它为x 的最佳平方逼近函数 8.求()sin ,[0,0.1]f x x x =∈在空间2(){1,,}x span x x Φ=上的最佳平方逼近多项式,并给出平方误差。
9.利用正交化方法求[0,1]上带权1()ln
x x ρ=的前三个正交多项式012(),(),()p x p x p x 。
10.求函数3f(x)x =在区间[0,1]上的二次最佳平方逼近多项式*2().S x
11. 写出Legendre 多项式的表达式,证明Legendre 多项式的正交性.并用Legendre 多项式构造4f(x)x =在[1,1]-上的二次最佳平方逼近多项式*
2().S x 并估计平方逼近误差2
2||||.δ
第四章 数值积分和数值微分
内容及纲目:
1) 插值型求积公式
2) Guass 求积公式
3) 复合数值求积公式
4) 外推方法
5) 数值微分
1.用复化Simpson 公式计算积分 dx x ⎰21ln 的近似值时,为使结果具有4位有效数字,需要取多少个节点处的函数值。
2.计算积分
⎰10dx e x ,如用复化梯形公式,问区间应分多少等份才能保证计算结果有五位有效数字。
3.确定下列求积公式的参数101,,A A A -,使其代数精度尽量高,并指明求积公式所具有的代数精度。
⎰++-=--h h h f A f A h f A dx x f )()0()()().1(101.
)].()0([)]()0([2)().2(2
h f f h h f f h dx x f h '-'++⎰≈α (3) 1201()[(1)2()3()].3
h
f x dx f f x f x ≈-++⎰求12??x x == (4) 确定下列求积公式的参数D C B A ,,,,使其代数精度尽量高,并指明求积公式所具有的
代数精度。
)]()([)]()([)()(103102010x f D x f C h x Bf x Af h dx x f x x x x
'+'+⎰+≈- 其中01x x h -=
4.建立Guass 求积公式
1
0011()()().f x dx A f x A f x ≈+⎰ 5.求10,A A 及10,x x 使
⎰-+≈1
11100)()()(x f A x f A dx x f 达到最大的代数精度。
6. 利用复合梯形求积公式计算积分:
⎰10sin dx x x ,使截断误差不超过31021-⨯。
取同样步长,改用复合
Simpson 求积公式计算,问截断误差界是多少?
7. 已知 .43,21,41310===
x x x 1)
推导一这3个点作为求积节点在]1,0[上的插值型求积公式。
2)
指明求积公式的代数精度; 3) 用所求公式计算dx x ⎰1
02
8.用三点公式求2)1(1
)(x x f +=在
2.1,1.1,0.1=x 的导数值,并估计误差。
9.推导4点数值微分公式001231()[111892]6f x f f f f h '=-+-+3
(4)().4
h f ξ-. 其中10 3.(0,1,2),i i h x x i x x ξ+=-=<<
10. 用4n =的复化梯形公式和2n =的复化Simpson
公式计算积分9
1的近似值.。