20150318数值分析学生版作业
数值分析作业及参考答案
数值分析第一次作业及参考答案1. 设212S gt =，假定g 是准确的，而对t 的测量有0.1±秒的误差，证明当t 增加时S 的绝对误差增加，而相对误差却减少。
解：2**22211()0.122()0.10.2()1122,(),().r r e S S S gt gt gt e S gt e S t gt gt t e S e S =-=-====∴↑↑↓2. 设2()[,]f x C a b ∈且()()0f a f b ==，求证2''1max ()()max ().8a x ba xb f x b a f x ≤≤≤≤≤-解：由112,0),(,0)()()0()00.a b L x l x l x =⨯+⨯=（两点线性插值 插值余项为"111()()()()()()[,]2R x f x L x f x a x b a b ξξ=-=--∈ [,].x a b ∴∀∈有12211()()"()()()max "()[()()]221()()1max "()[]()max "().228a x ba xb a x b f x R x f x a x b f x x a b x x a b x f x b a f x ξ≤≤≤≤≤≤==--≤---+-≤=-21max ()()max "()8a xb a x b f x b a f x ≤≤≤≤∴≤-3. 已测得函数()y f x =的三对数据：（0，1），（－1，5），（2，－1），（1）用Lagrange 插值求二次插值多项式。
（2）构造差商表。
（3）用Newton 插值求二次插值多项式。
解：(1)Lagrange 插值基函数为0(1)(2)1()(1)(2)(01)(02)2x x l x x x +-==-+-+-同理 1211()(2),()(1)36l x x x l x x x =-=+ 故2202151()()(1)(2)(2)(1)23631i i i p x y l x x x x x x x x x =-==-+-+-++=-+∑（2）令0120,1,2x x x ==-=，则一阶差商、二阶差商为0112155(1)[,]4,[,]20(1)12f x x f x x ---==-==-----0124(2)[,,]102f x x x ---==-22()1(4)(0)1*(0)(1)31P x x x x x x =+--+-+=-+4. 在44x -≤≤上给出()xf x e =的等距节点函数表，若用二次插值求x e 的近似值，要使截断误差不超过610-，问使用函数表的步长h 应取多少？解：()40000(),(),[4,4],,,, 1.x k x f x e f x e e x x h x x h x x th t ==≤∈--+=+≤考察点及(3)200044343()()[(()]()[()]3!(1)(1)(1)(1)3!3!.(4,4).6f R x x x h x x x x h t t t e t h th t h e h e ξξ=----+-+≤+⋅⋅-=≤∈-则436((1)(1)100.006.t t t h --+±<< 在点 得5. 求2()f x x =在［a,b ］上的分段线性插值函数()h I 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位有效数字。
32031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +，b a ⨯有几位有效数字.〔有效数字的计算〕 解：3*1021-⨯≤-aa ，2*1021-⨯≤-b b ，而1811.2=+b a ，1766.1=⨯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 x x 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 v6设x 的相对误差为%a ,求nx y =的相对误差。
数值分析第三版课本习题及答案
第一章 绪 论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 cs a b c ∆∆∆∆≤++第二章 插值法1. 根据(2.2)定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n n n 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. 给出f (x )=ln x 的数值表用线性插值及二次插值计算ln 0.54 的近似值.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 jj x l x x k 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 ba xb f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()xf x e =的等距节点函数表,若用二次插值求xe 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少?9. 若2n n 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 lf 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()()()bbbba a a a f 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 可用(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. 试证{}*()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次逼近多项式,使误差不超过0.005.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n nF 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()sin 2f x x=在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.26. 用最小二乘法求一个形如2y a bx =+的经验公式,使它与下列数据拟合,并求均方误差.27.28. 在某化学反应里,根据实验所得分解物的浓度与时间关系如下:用最小二乘拟合求.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)()hh f x dx A f h A f A f h --≈-++⎰; (2)21012()()(0)()hh f x dx A f h A f A fh --≈-++⎰;(3)[]1121()(1)2()3()/3f x dx f f x fx -≈-++⎰;(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.用辛普森公式求积分1xedx-⎰并计算误差.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.用复化梯形公式求积分()b af x dx⎰,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍入误差)?8.1xedx-,要求误差不超过510-.9. 卫星轨道是一个椭圆,椭圆周长的计算公式是S a =θ,这里a 是椭圆的半长轴,c是地球中心与轨道中心(椭圆中心)的距离,记h 为近地点距离,H 为远地点距离,6371R =公里为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第一颗人造卫星近地点距离439h =公里,远地点距离2384H =公里,试求卫星轨道的周长.10. 证明等式3524sin3!5!n nnnππππ=-+-试依据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相比较。
数值分析习题含答案
x1 )
f (x0)
(x
x 0 )( x x0 x1
x1 )
f ' ( x0 )
(x ( x1
x0)
2 2
x0 )
f ( x1 )
R ( x)
其中 R(x) 由以下计算得到： 构造辅助函数：
(t ) f (t ) N 2 (t ) (t (x x0 ) (t x0 ) ( x
2 2
x1 ) x1 )
f [ 2 ,2 ] =-2089 ，
0 1 2 7
0 1 7
f (x)
M ,
x
[ a , b ] ，证明：在任意相邻两节点间
R1 ( x )
1 8
Mh
2
。
x xi x xi M
1
f ( ) R1 i ( x ) 2 M 8 h 2,
h ,
2
x
8 ,n
[ xi , xi
1
]
R1 ( x )
max R1 i ( x )
1 2
s
2
[( x
xi
1
))( x
x
i
1 2
)( x
x i )]
e
4
h
3
[ s( s
1)( s
1)] 24
3 9
e h
4
3
10
6
3!
8
h
1 . 317
则用二次插值的步长应：
h
0 .6585
10
2
2-6 对区间 [a,b] 作步长为 h 的剖分，且 做线性插值，其误差限为 证明：区间上的误差限： 误差限： 2-7 设 f ( x ) 解： 自变量 1 2
数值分析习题（含答案）
数值分析习题（含答案）第一章绪论姓名学号班级习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。
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π=的绝对误差限与相对误差限。
数值分析作业题(1)
第一章 误差与算法1. 误差分为有__模型误差___, _观测误差___, __方法误差____，___舍入误差____, Taylor 展开式近似表达函数产生的误差是_方法误差 .2. 插值余项是插值多项式的 方法误差。
3. 0.2499作为1/4的近似值，有几位有效数字？00.24990.249910,0m =⨯=即，031|0.2499|0.00010.5100.510,34m n n ---=<⨯=⨯=即223.1428751...,7=作为圆周率的近似值，误差和误差限分别是多少，有几位有效数字？2133.142875 3.14159260.00126450.5100.510---=<⨯=⨯有3位有效数字.* 有效数字与相对误差的关系 4. 利用递推公式计算积分110,1,2,...,9n x n I x e dx n -==⎰错误！未找到引用源。
, 建立稳定的数值算法。
1111111110011,n 2,...,9n x n x n x n x n n I x e dx x dex en x e dx nI ------===-=-=⎰⎰⎰该算法是不稳定的。
因为：11()()...(1)!()nn n I n I n I εεε-=-==- 111n n I I n n-=-， 10110I =5. 衡量算法优劣的指标有__时间复杂度,__空间复杂度_.6. 时间复杂度是指：.算法需耗费时间的度量， 两个n 阶矩阵相乘的乘法次数是 3n ,则称两个n 阶矩阵相乘这一问题的时间复杂度为3()O n .二 代数插值1.根据下表数据建立不超过二次的Lagrange 和Newton 插值多项式，并写出误差估计式，以及验证插值多项式的唯一性。
x 0 1 4 f(x) 1 9 3 Lagrange:设0120120,1,4;()1()9()3x x x f x f x f x ======则，， 对应i x 的标准基函数)(x l i 为：1200102()()(1)(x 4)1()(1)(x 4)()()(01)(04)4x x x x x l x x x x x x ----===------1()...l x =2()...l x =因此，所求插值多项式为：220()()()....i i i P x f x l x ===∑(3)2()()(0)(1)(x 4)3!f R x x x ξ=---Newton:构造出插商表：xi f(xi ) 一 二 三 0 1 1 9 8 4 3 -2 -5/2 所以, 所求插值多项式为：2001001201()()[,]()[,,]()()518(0)(0)(1)2...P x f x f x x x x f x x x x x x x x x x =+-+--=+----=插值余项：2()[0,1,4,](0)(1)(x 4)R x f x x x =---2. 已知函数f(0)=1,f(1)=3,f(2)=7,则f[0,1]=___2________, f[0,1,2]=____1______)('],[000x f x x f =3. 过0,1两节点构造三次Hermite 插值多项式，使得满足插值条件：f (0)=1, f ’(0)=0 , f (1) =2, f ’(1)=1设0101010,1,()1()2'()0,'()1x x f x f x f x f x ======则，，写出插商表：xi f(xi) 一 二 三 0 1 0 1 0 1 a 1 1 1 a 1 0 a-1因此, 所求插值多项式为：222000000100011012232()()[,]()[,,]()[,,,]()()10(0)1(0)1(0)(1)21P x f x f x x x x f x x x x x f x x x x x x x x x x x x x x =+-+-+--=+-+----=-++插值余项：222()[0,0,1,1,](1)R x f x x x =-4. 求f (x)=sin x 在[a,b]区间上的分段线性插值多项式，并写出误差估计式。
数值分析大作业
数值分析大作业数值分析大作业姓名：黄晨晨学号：S1*******学院：储运与建筑工程学院学院班级：储建研17-2实验3.1 Gauss消去法的数值稳定性实验实验目的：理解高斯消元过程中出现小主元即很小时引起方程组解数值不定性实验内容：求解方程组Ax=b，其中（1）A1=0.3×10?1559.14315.291?6.130?1211.29521211，b1=59.1746.7812;（2）A2=10?7013 2.099999999999625?15?10102，b2=85.90000000000151;实验要求：(1)计算矩阵的条件数，判断系数矩阵是良态的还是病态的(2)用Gauss列主元消去法求得L和U及解向量x1,x2∈R4(3)用不选主元的高斯消去法求得L和U及解向量x1,x2∈R4(4)观察小主元并分析对计算结果的影响（1）计算矩阵的条件数，判断系数矩阵是良态的还是病态的代码：format longeformat compactA1=[0.3*10^-15,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1] b1=[59.17;46.78;1;2]n=4C1=cond(A1,1) %C1为A1矩阵1范数下的条件数C2=cond(A1,2) %C2为A1矩阵2范数下的条件数C3=cond(A1,inf) %C3为1矩阵谱范数下的条件数结果：C1 =1.362944708720448e+02C2 =6.842955771253409e+01C3 =8.431146*********e+01显然A1矩阵为病态矩阵将矩阵A2，b2输入上述代码中求得A2矩阵的条件数为：C1 =1.928316831682894e+01C2 =8.993938090170119e+00C3 =1.835643564356072e+01显然A2矩阵也为病态矩阵(2)用Gauss列主元消去法求得L和U及解向量x1,x2∈R4代码：for k=1:n-1a=max(abs(A1(k:n,k)))if a==0returnendfor i=k:nif abs(A1(i,k))==ay=A1(i,:)A1(i,:)=A1(k,:)A1(k,:)=yx=b1(i,:)b1(i,:)=b1(k,:)b1(k,:)=xbreakendendif A1(k,k)~=0A1(k+1:n,k)=A1(k+1:n,k)/A1(k,k)A1(k+1:n,k+1:n)=A1(k+1:n,k+1:n)-A1(k+1:n,k)*A1(k,k+1:n) elsebreakendendL=tril(A1,0);for i=1:nL(i,i)=1;endLU=triu(A1,0)y1=L\b1x1=U\y1得到如下结果：x1 =3.845714853511634e+001.609517394778522e+00-1.547605454206655e+011.041130489899787e+01将A2=[10,-7,0,1;-3,2.0999********,6,2;5,-1,5,-1;0,1,0,2]b2=[8;5.900000000001;5;1]代入上述代码求得结果如下：X2 =4.440892098500626e-16-9.999999999999993e-019.999999999999997e-011.000000000000000e+00(3)用不选主元的高斯消去法求得L和U及解向量x1,x2∈R4代码：format longeformat compactA1=[0.3*10^-15,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1] b1=[59.17;46.78;1;2][L,U]=lu(A1)y1=L\b1x1=U\y1求得如下结果：x1=3.845714853511634e+001.609517394778522e+00-1.547605454206655e+011.041130489899787e+01将A2=[10,-7,0,1;-3,2.0999********,6,2;5,-1,5,-1;0,1,0,2] b2=[8;5.900000000001;5;1]代入上述代码，求得结果如下：x 2 =4.440892098500626e-16 -9.999999999999993e-01 9.999999999999997e-01 9.999999999999999e-01（2）（3）求得结果相同，可验证结果正确。
数值分析作业
数值分析课后作业：习题一1.在字长为3的十进制计算机上计算f （3.33）和g （3.33），其中f(x)=x 4-x 3+3x 2+x-2,g(x)=(((x-1)x+3)x+1)x-2解： m=3; f=@(x)digit(digit(x^4,m)- digit(x^3,m)+ digit(3*x^2,m)+ digit(x-2,m),m); g=@(x)digit(digit(digit( digit(digit(digit( (x-1)*x,m)+3,m)*x,m)+1,m)*x,m)-2,m); f(3.33) g(3.33) 有ans = 121 ans =121 2.下列各近似值的绝对误差限都是1021⨯-3，试指出它们各有几位有效数字：x=1.00052， y=0.05， z=0.00052.解：当 x=1.00052时， 由丨X*—X 丨 ≤0.5×10-3 得 x=1.00052 有四位有效数字； 同理 y=0052 有两位有效数字 Z=0.00052有零位有效数字 3,计算圆的面积，要使其相对误差限为1%，问测量半径r 允许的相对误差限是多少？ 解：设圆的面积为S ， 由题意有|e(S)|≤1%。
又S=πr 2 dS=2πr dr 所以 dS/S=(2πrdr)/(πr 2)=2(dr/r)∴|e(r)|≈21|e(S)|≤0.5×1%=0.5% 11.数组与矩阵是Matlab 编程的基础，试学习Matlab 的数组与矩阵的表示方法，并举例介绍数组、矩阵的常见运算. 解：>> syms a b c d; >> a=[1 2 3];>> b=[4 5 6];>> a+bans =5 7 9>> b-aans =3 3 3>> a.*bans =4 10 18 >> a.^2 ans = 1 4 9>> c=[1 2 3;1 2 3;1 2 3];>> d=[4 5 6;4 5 6;4 5 6];>> cc = 1 2 3 1 2 3 1 2 3d = 4 5 6 4 5 6 4 5 6 >> c+dans =5 7 9 5 7 9 5 7 9>> d-cans = 3 3 33 3 33 3 3 12.学习使用Matlab 命令help 和doc 学习自己感兴趣的Matlab 的运算、函数或命令的用法，并对于任意给定的实数a,b,c,编写Matlab 程序求方程ax 2+bx+c=0的根. 解：x 1=a ac b b b 24)sgn(2---, x 2=1ax c1 x>0 其中 sgn = 0 x=0 -1 x<0 disp('Please input the coefficients of');disp('quadratic equation ax^2+bx+c=0, respectively') a=input('a='); b=input('b='); c=input('c=');m=3; if abs(a)<eps & abs(b)<eps error End if abs(a)<eps disp('Since a=0, quadrtic equation degen erates into a linear equation.') disp('The only solution of the linear equtio n is')x=digit(-c/b,m) return Enddelta=b^2-4*a*c; temp=sqrt(delta); x 1=(-b+temp)/(2*a) ; x 2=(-b-temp)/(2*a) ;err1=abs(a*x 1^2+b*x 1+c) ; err2=abs(a*x 2^2+b*x 2+c) ; if b>0x 1=(-b-temp)/(2*a) End if b<0x 1=(-b+temp)/(2*a) End if b=0x 1=temp/(2*a) Endx 2=c/(a*x 1)err1=abs(a*x 1^2+b*x 1+c) err2=abs(a*x 2^2+b*x 2+c) if abs(a)<epsdisp('Since a=0, quadrtic equation degen erates into a linear equation.')disp('The only solution of the linear equtio n is')x=digit(-c/b,m) return Enddelta=digit(digit(b^2,m)-digit(4*digit(a*c,m),m),m);temp=digit(sqrt(delta),m);x 1=digit(digit(-b+temp,m)/digit(2*a,m),m); x 2=digit(digit(-b-temp,m)/digit(2*a,m),m); err1=abs(a*x 1^2+b*x 1+c); err2=abs(a*x 2^2+b*x 2+c); if b>0x 1=digit(digit(-b-temp,m)/digit(2*a,m),m) ; End if b<0x 1=digit(digit(-b+temp,m)/digit(2*a,m),m); End if b=0x 1=digit(temp/digit(2*a,m),m); Endx 2=digit(digit(c/a,m)/x1,m) ; err1=abs(a*x 1^2+b*x 1+c) ; err2=abs(a*x 2^2+b*x 2+c) ; 14分别利用ln (1+x)=11,)1(11≤<--+∞=∑x nx nn n 和ln11...),12...53(2111253<<-++++++=-++x n x x x x x x n ，给出计算ln2的近似方法，编写相应的Matlab 程序，并比较算法运行情况. 解：方法一： x=1; s=0;for k=1:100s=s+(-1)^(k+1)*(x^k)/k; end sq=log(2)err=abs(t-q) ans= t =0.6882 q =0.6931 err = 0.0050方法二x=1/3; s=0;for k=1:2:100 s=s+(x^k)/k; end t=2*s q=log(2)err=abs(t-q) Ans= t =0.6931 q =0.6931 err =2.2204e-16所以方法二较方法一好。
数值分析习题与答案.doc
第一章绪论习题一1.设x>0,x*的相对误差为δ，求f(x)=ln x的误差限。
解：求lnx的误差极限就是求f(x)=lnx的误差限，由公式(1.2.4)有已知x*的相对误差满足，而，故即2.下列各数都是经过四舍五入得到的近似值，试指出它们有几位有效数字，并给出其误差限与相对误差限。
解：直接根据定义和式(1.2.2)(1.2.3)则得有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),使它满足解：这种题目可以有很多方法去做，但应以简单为宜。
数值分析作业答案
数值分析作业答案【篇一：《数值分析》作业参考答案 2】>一. 选择题1． a； 2.b； 3.b； 4.d; 5.c; 6.d; 7.c; 8.b; 9.d; 10.c； 11.b；12.a； 13.a; 14.c; 15.a; 16.b; 17.d; 18.a19.d,20.c,21.a,22.d,23.c,24.c. 二. 填空题1. 3，3，3 ；2. 1，2/3；3. 100！2^100 ；4. (-1 ,1)；5.x?cosx1； 6. g(x)?x?，2；1?sinx(x0?x1)(x0?x2)?(x0?xn)7. (-1, 1)； 8. x； 9. 4 ； 10. 5，9 ； 11. xn?1?31c(xn?), 2； 2xnb?a1x2?112. x?x?1 ； 13. 10/9, 4； 14. 10, 55, 550; 15. b?16.17.33318．?三.1. p(x)?x2?2x?12. p(x)?f(x)3. 2.a?c?a?b11?i1?3i219.2x?x20. a? . ,,324429121016, 代数精确度为5 ,b?,a?9953.证明：cond(aa)?||aa||?||(aa)?1|| 设??max{aa的特征值的模}，??max{(a?1)a?1的特征值的模}，则上式=????||a||2?||a?1||2?(cond(a))2?1???24.1. (12分)l???????2?1??1???????1????0?2u?, , x??????4?1?1???????0??1??????1?31?2. (8分)seidel收敛，因为a 实正定对称阵. 迭代格式(k)?x1(k?1)?(2?x2)/2?(k?1)(k)?(x1(k?1)?x3)/2?x2?(k?1)(k?1)(k)x?(?2?x?x24)/2?3(k?1)(k?1)?x4?(1?x3)/2?12??35. p(x)??2x?1，余项cosx?p2(x)?|x(x?|?6.?6254*6426. 证明：当a?0时，结论显然成立；|ytax|当a?0时，因|yax|?||y||2||a||2||x||2，故 sup?||a||2；x?0,y?0||x||2||y||2t又aa是实对称矩阵，故存在正交阵p?(p1,p2,?,pn)t??1???tt?使得paap?d???， pi是特征值?i对应的特征向量。
数值分析第2次作业
数值分析第二次作业工院雷祺学号：S201503003主程序%p95第一题主程序%雷祺2015.11.5close all;clear;clc;syms tfun=sym(1/(1+25*t^2));%最佳平方逼近fn=Squ_appr(fun,3,-1,1);%勒让德多项式最佳平方逼近fnl=MyLegendre(fun,3);i=0:1:10;x=-1+0.2.*i;y=1./(1+25.*x.^2);xx=linspace(-1,1,100);yy=1./(1+25.*xx.^2);%最小二乘A=MyLeast_squ2(x,y,3);%正交多项式最小二乘A2=MyLeast_squ2(x,y,3);% 三次样条[fl,yys]=mysplineN(x,y,0,0,xx);% % 多项式% [f,yyL]=myLanguage(x,y,xx);yynh=subs(fn,'t',xx);yynh1=subs(fnl,'t',xx);AA=subs(A,'t',xx);AA2=subs(A2,'t',xx);% %帕德拟合% yPade=MyPade_2(fun,1,1)subplot(2,1,1);plot(x,y,'*',xx,yy,xx,yys,'--',xx,yynh,'-.',xx,yynh1,'o',xx,AA,xx,AA2,'+','Linewidth',0.5);legend('插值节点','f(x)','三次样条插值','最佳平方逼近','勒让德多项式最佳平方逼近','最小二乘','正交多项式最小二乘');title('插值与拟合的比较');xlabel('x');ylabel('y');box on;grid on;subplot(2,1,2);plot(xx,yy-yys,'--',xx,yy-yynh,'-.',xx,yy-yynh1,'.',xx,yy-AA,xx,yy-AA2,'+','Linewidth',0.5);legend('三次样条插值','最佳平方逼近','勒让德多项式最佳平方逼近','最小二乘','正交多项式最小二乘');title('插值与拟合误差');xlabel('x');ylabel('y');box on;grid on;% figure% plot(xx,AA-AA2)运行结果：%p96第2题主程序%雷祺2015.11.5close all;clear;clc;syms tx=[0.0 0.1 0.2 0.3 0.5 0.8 1.0];y=[1.0 0.41 0.50 0.61 0.91 2.02 2.46];%最小二乘A3=MyLeast_squ(x,y,3);%三次A4=MyLeast_squ(x,y,4);%四次%%%%勒让德多项式最小二乘%%%%%%%%A6=MyLegendre2(x,y,4);%三角插值A5=mytriange(y,3);xx=2*pi*x-pi;%转换区间为(-pi,pi)%计算数值AA3=subs(A3,'t',x);AA4=subs(A4,'t',x);AA6=subs(A6,'t',x);AA5=subs(A5,'t',xx);subplot(2,1,1)plot(x,y,x,AA3,'o',x,AA4,'diamond',x,AA5,'pentagram',x,AA6,'+','Linewidth',1.5);title('第二题画图');legend('原节点','三次最小二乘','四次最小二乘','四次三角插值','四次勒让德多项式','Location','NorthWest');xlabel('x');ylabel('y');box on;grid on;subplot(2,1,2)plot(x,zeros(1,7),x,y-AA3,'o',x,y-AA4,'diamond',x,y-AA5,'pentagram',x,y-AA6,'+','Linewidth',2); title('误差');legend('参考线','三次最小二乘','四次最小二乘','三次三角插值','四次勒让德多项式','Location','NorthWest');xlabel('x');ylabel('E');box on;grid on;运行结果：%p96第3题主程序%雷祺2015.11.5close all;clear;clc;j=linspace(0,31,32);x=-pi+j*pi/16;%取32个点xx=linspace(-pi,pi,500);%500f=x.^2.*cos(x);ff=xx.^2.*cos(xx);S =mytriange(f,16);%三角逼近SS=subs(sym(S),'t',xx);plot(xx,ff,xx,SS,'-.','Linewidth',1.3);%绘图legend('原函数','三角插值');title('快速傅里叶变换的三角插值');xlabel('x');ylabel('y');grid on;box on;axis on;vpa(S,6)% plot(xx,ff-SS,'Linewidth',1.3);%绘图运行结果：ans =0.066436*cos(9.0*t) - 0.315508*cos(4.0*t) - 0.0801571*cos(8.0*t) - 0.0389271*cos(16.0*t) - 2.2352*cos(2.0*t) + 5.43643e-16*sin(2.0*t) - 1.81507e-16*sin(4.0*t) + 8.83477e-17*sin(8.0*t) - 5.72546e-17*sin(16.0*t) + 0.194081*cos(5.0*t) - 0.0570435*cos(10.0*t) - 7.96259e-17*sin(9.0*t) + 0.0504667*cos(11.0*t) + 5.16803e-17*sin(5.0*t) + 4.66853e-16*sin(10.0*t) + 0.638104*cos(3.0*t) - 0.134645*cos(6.0*t) - 0.0458343*cos(12.0*t) - 4.34987e-16*sin(11.0*t) + 0.0426232*cos(13.0*t) - 7.26949e-18*sin(3.0*t) + 6.68525e-17*sin(6.0*t) + 3.1259e-16*sin(12.0*t) + 0.101006*cos(7.0*t) - 0.040514*cos(14.0*t) - 3.22878e-16*sin(13.0*t) + 0.0393158*cos(15.0*t) - 1.24978e-16*sin(7.0*t) + 3.84891e-16*sin(14.0*t) - 2.99081e-16*sin(15.0*t) + 3.80277*cos(t) - 1.15345e-15*sin(t) - 2.00644逼近与拟合函数最佳平方逼近function coff=Squ_appr(func,n,a,b)% 最佳平方逼近%输入依次为：原函数，逼近函数次数，范围H = zeros(n+1,n+1);var = findsym(sym(func));%以x^k作基func = func/var;for k=1:n+1func = func*var;%f(x)*x^kd(k,1)=int(sym(func),var,a,b); %D向量% funcc=funcc*var;for j=1:n+1H(k,j)=(power(b,k+j-1)-power(a,k+j-1))/(k+j-1);endendcoff = H\d; %求解逼近多项式的系数C\dcoff=poly2sym(flipud(coff),var);% coff=vpa(coff,6);勒让德多项式n次最佳平方逼近function f = MyLegendre(y,n)%勒让德多项式n次最佳平方逼近。
数值分析作业答案(第2章)
2.2.给出x x f ln )(=的数值表：用线性插值及二次插值计算54.0ln 的近似值。
解 线性插值。
由于54.0=x ，介于0.5和0.6之间，故取5.00=x ，6.01=x ，这时插值余项中的))(()(10x x x x x w --=的绝对值最小，于是147693.00-=y ，826510.01-=y ，代入拉格朗日线性插值多项式，得219620.0)826510.0(5.06.05.054.0)147693.0(6.05.06.054.0)54.0(11001011-≈-⨯--+-⨯--=⋅--+⋅--=y x x x x y x x x x L所以219620.0)54.0(54.0ln 1-≈≈L 。
当然还可以按其他方式取0x ，1x ，但近似程度可能差些。
二次插值。
由于54.0=x ，与0.5，0.6及0.4距离较近，故取4.00=x ，5.01=x ，6.02=x ，这时插值余项中的))()(()(210x x x x x x x w ---=的绝对值最小，于是291916.00-=y ，147693.01-=y ，826510.02-=y ，代入拉格朗日二次插值多项式，得320615.0)826510.0()5.06.0)(4.06.0()5.054.0)(4.054.0()147693.0()6.05.0)(4.05.0()6.054.0)(4.054.0()291916.0()6.04.0)(5.04.0()6.054.0)(5.054.0())(())(())(())(())(())(()54.0(2120210121012002010212-≈-⨯----+-⨯----+-⨯----=⋅----+⋅----+⋅----=y 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 L所以320615.0)54.0(54.0ln 2-≈≈L 。
(完整word版)数值分析上机作业1-1解析
考虑一个高次的代数多项式
（E1-1）
显然该多项式的全部根为l，2，…，20，共计20个，且每个根都是单重的（也称为简单的）。现考虑该多项式方程的一个扰动
(E1-2)
其中 是一个非常小的数。这相当于是对（E1-1）中 的系数作一个小的扰动。我们希望比较（E1-1）和（E1-2）根的差别，从而分析方程（E1-1）的解对扰动的敏感性。
ve=zeros(1,21);
ve(21-Numb)=ess;
root=roots(poly(1:20)+ve);
x0=real(root); y0=imag(root);
plot(x0',y0','*');
disp(['对扰动项 ',num2str(Numb),'加扰动',num2str(ess),'得到的全部根为:']);
ess分别为1e-6,1e-8.1e-10,1e-12的图像如下：
从实验的图形中可以看出，当ess充分小时，方程E.1.1和方程E.1.2的解相差很小，当ess逐渐增大时，方程的解就出现了病态解，这些解都呈现复共轭性质。
（2）将扰动项加到x18上后，ess=1e-009时方程的解都比较准确，没有出现复共轭现象。ess=1e-008时误差与x19(ess=1e-009)时相当，即扰动加到x18上比加到x19小一个数量级。对x8的扰动ess=1000时没有出现复共轭，误差很小；对x的扰动ess=10e10时没有出现复共轭，误差很小。因此，扰动作用到xn上时，n越小，扰动引起的误差越小。
if((Numb>20)|(Numb<0))errordlg('请输入正确的扰动项:[0 20]之间的整数!');return;end
数值分析作业(1,2)
数值分析作业(1)1：思考题（判断是否正确并阐述理由）（a）一个问题的病态性如何，与求解它的算法有关系。
(b)无论问题是否病态，好的算法都会得到它好的近似解。
(c)计算中使用更高的精度，可以改善问题的病态性。
(d)用一个稳定的算法计算一个良态问题，一定会得到他好的近似解。
(e)浮点数在整个数轴上是均匀分布。
(f)浮点数的加法满足结合律。
(g)浮点数的加法满足交换律。
(h)浮点数构成有效集合。
(i)用一个收敛的算法计算一个良态问题，一定得到它好的近似解。
√2: 解释下面Matlab程序的输出结果t=0.1;n=1:10;e=n/10-n*t3：对二次代数方程的求解问题20++=ax bx c有两种等价的一元二次方程求解公式22b x ac x -±== 对a=1，b=-100000000，c=1，应采用哪种算法？4：函数sin x 的幂级数展开为： 357sin 3!5!7!x x x x x =-+-+ 利用该公式的Matlab 程序为function y=powersin(x)% powersin. Power series for sin(x)% powersin(x) tries to compute sin(x)from a power seriess=0;t=x;n=1;while s+t~=s;s=s+t;t=-x^2/((n+1)*(n+2))*tn=n+2;end（a ） 解释上述程序的终止准则；（b ） 对于x=/2π、x=11/2π、x =21/2π，计算的精度是多少？分别需要计算多少项？5:指数函数的幂级数展开2312!3!x x x e x =++++ 根据该展开式，编写Matlab 程序计算指数函数的值，并分析计算结果（重点分析0x <的计算结果）。
数值分析作业（2）思考题1：判断下面命题是否正确并阐述理由（a）仅当系数矩阵是病态或奇异的时候，不选主元的Gauss消元法才会失败。
数值分析课程设计大作业
课程设计2013年07月20日设计题目 《数值分析》课程设计学生姓名 ****学 号 ####专业班级指导教师1.1水手、猴子和椰子问题算法分析:设椰子起初的数目为0p ,第一至第五次猴子在夜里藏椰子后，椰子的数目分别为0p ，1p ，2p ，3p ，4p ，再设最后每个人分得x 个椰子，由题意得：15541(1),0,1,2,3,4.(1)=5155k k p p k x p p x +=-==-+所以利用逆向递推方法求解：n=input('n='); for x=1:n p=5*x+1; for k=1:5 p=5*p/4+1; endif p==fix(p) break end enddisp([x,p])执行代码后得： n= 1023 15621 （输入n=1000000） 即最后每个人分得1023个椰子，椰子总数为156211.2当0,1,2,,100n = 时，选择稳定的算法计算积分10d 10nxn xe I x e --=+⎰ 由1100(1)1110010101,1010110(1)10x x n x nxnxn n n x e I I dx e e e I I dx e dx e e n---+---+-++==+++===-+⎰⎰⎰ 得0111(1)1011[(1)],100,99,...,1.10n n n I I I e I n n -+⎧=-⎪⎪⎨⎪=--=⎪⎩由上式可知求n I 时，1n I +的误差的影响被缩小了。
n=100时100I 的近似值为0。
matlab 代码为fprintf('稳定算法:\n')y0=0;n=100;plot(n,y0,'r*');hold onfprintf('y[100]=%10.6f',y0);while(1)y1=1/10*[(1-exp(-n))/n-y0];fprintf('y[%10.0f]=%10.6f',n-1,y1);plot(n-1,y1,'r*') if(n<=1) break;endy0=y1;n=n-1;if mod(n,3)==0,fprintf('\n'),end,end(具体值已省略)编程实现得下图。
数值分析上机作业(2)
一、数值求解如下正方形域上的Poisson 方程边值问题 2222(,)1,0,1(0,)(1,)(1),01(,0)(,1)0,01u u f x y x y x y u y u y y y y u x u x x ⎧⎛⎫∂∂-+==<<⎪ ⎪∂∂⎪⎝⎭⎨==-≤≤⎪⎪==≤≤⎩二、用椭圆型第一边值问题的五点差分格式得到线性方程组为2,1,1,,1,10,1,,0,141,?,?,?,?0,1i j i j i j i j i j ijj N j i i N u u u u u h f i j N u u u u i j N -+-+++----=≤≤====≤≤+， 写成矩阵形式Au=f 。
其中1．三 、编写求解线性方程组Au=f 的算法程序, 用下列方法编程计算, 并比较计算速度。
2．用Jacobi 迭代法求解线性方程组Au=f 。
3．用块Jacobi 迭代法求解线性方程组Au=f 。
4． 用SOR 迭代法求解线性方程组Au=f,用试算法确定最佳松弛因子。
1122N N v b v b u f v b ⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪== ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭4114114ii A -⎛⎫ ⎪- ⎪= ⎪- ⎪-⎝⎭11,12,1,121,22,2,21,2,,2211,12,1,121,22,2,221,2,,(,,...,),(,,...,),......,(,,...,)(,,...,)?,(,,...,)?,......,(,,...,)?1,999,0.10.011T T N N TN N N N N T T N N T N N N N N v u u u v u u u v u u u b h f f f b h f f f b h f f f h N h N ====+=+=+===+取或则或,1,,1,2,...,i j f i j N== 1122NN A I I A A I I A -⎛⎫ ⎪- ⎪= ⎪- ⎪-⎝⎭5．用块SOR 迭代法求解线性方程组Au=f,用试算法确定最佳松弛因子。
20150318数值分析学生版作业.
2014-2015(2)计算机与信息工程学院数值分析作业计科专业_______ 级____ 班姓名: _____________ 学号: ____________第一章绪论一、单项选择题1.用3.1415作为二的近似值时具有( )位有效数字。
(A) 3 ( B) 4 ( C 5(D) 62.已知数X1=721 X2=0.721 X3=0.700 X4=7*10-2是由四舍五入得到的，则它们的有效数字的位数应分别为()。
(A)3，3，3，1 ( B) 3，3，3，3(C) 3, 3, 1, 1 ( D) 3, 3, 3, 2二、填空题1.在一些数值计算中,对数据只能取有限位表示,如2 : 1.414 ,这时所产生的误差称为_______ 差.(填误差的类型)2.为尽量避免有效数字的严重损失，当x • • 1时，应将表达式..x , x改写为____________ 保证计算结果比较精确.3.在数值计算中，通常取e = 2.71，此时产生的误差为___________ 差(填误差的类型).4.__________________________________________________________ 设x=0.231是精确值x*=0.229的近似值，贝U x有_____________________________ 有效数字。
三、计算题221、(本题5分)试确定—作为…的近似值具有几位有效数字，并确定其相对误差限。
第二章插值法一、单项选择题1.通过点(X。
, y°),(x 1, yj的拉格朗日插值基函数l°(x),l 1(x)满足().(A ) l°(x°) =0」1(X1)=0 ( B)(C ) 10(X0) =1,1心1)=0 (D)2.冷点.广宀是给定的互异节点，:「是一个( ).(A)n+1次多项式(C) 次数小于n的多项式二、填空题10 (x 0) = 1,11(X 1 ) = 1〔0 (x0)= 0,11(xJ = 1：是以它们为插值节点的插值多项式，则(B) n次多项式(D)次数不超过n的多项式1.设有节点x0,x lr x 2，其对应的函数y = f(x)的值分别为y0,y 1,y 2则二次拉格朗日插值基函数l o(X)= ______________ .2._____________________________ 已知 f (x) =x2+1,则 f [1,2,3] = .2.已知f(1)=1,f(2) =3,那么y二f(x)以1，2为节点的拉格朗日线性插值多项式为_________ .3.当 x=1，-1，2 时，对应的函数值分别为 f(-1)=0，f(0)=2, f(4)=10,则 f(x)的拉格朗日插值多项式是___________________________________ .4.设f(x) =x2，则滋关于节点x o = 0，x1 =仃2 = 3的二阶向前差分为5.当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商公式的______ ，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的______________ ;如果要估计结果的舍入误差，应该选用插值公式中的_____________ .6.设f(0)=0, f(1)=10, f(2)=20,则fb,1】=—,fb,1,2】=___, f(x)的二次牛顿插值多项式为______________________________ .7.设Ln (X)为f( x)的n次拉格朗日插值多项式，则其插值余项为5 38 已知f(x)=2x +4x -5x,则f[-1,1,0]= _________ f[-3,-2,-1,1,2,3]=9.设f (x) =3x2 +5,X k =kh,(k = 0,1,2,…),则差商f [X n,X n 也X n 卷,X n4d = _____ .10.设l j(x)(j =0,1,2|l| n)是n次拉格朗日插值多项式的插值基函数，则nl j(x)二___________ (i, j =0,1,2川n) ; 'Tj(x)二—。
数值分析习题与答案,DOC
第一章绪论习题一1.设x>0,x*的相对误差为δ，求f(x)=lnx的误差限。
解：求lnx的误差极限就是求f(x)=lnx的误差限，由公式即有2位有效数字，有5位有效数字，3.下列公式如何才比较准确？（1）（2）解：要使计算较准确，主要是避免两相近数相减，故应变换所给公式。
，利用：式计算误差最小。
四个选项：给定的数值表解：仍可使用n=1及n=2的Lagrange插值或Newton插值,并应用误差估计（5.8）。
线性插值时，用0.5及0.6两点，用Newton插值误差限，因，故二次插值时，用0.5，0.6，0.7三点，作二次Newton插值上给出的等距节点函数表，若用二次插值法求，函数表的步长令因得3.若，求和.解：由均差与导数关系于是4.若互异，求可知当有于是得求证.6.已知的函数表求出三次Newton均差插值多项式，计算f(0.23)的近似值并用均差的余项表达式估计误差.解：根据给定函数表构造均差表由于值并估计误差解：先构造差分表计算，用其中计算时用Newton后插公式（5.18)误差估计由公式（5.19）得这里仍为可先造，显然，再令令称为第二类求的表达式，并证明］上带权的正交解：因，即，故法方解得均方程为11.填空题(1)满足条件的插值多项式p(x)=().(2),则f［1,2,3,4］=()，f［1,2,3,4,5］=().(3)设为互异节点，为对应的四次插值基函数，则＝()，＝().是区间［其中，则＝＝））））第4章数值积分与数值微分习题41.分别用复合梯形公式及复合Simpson公式计算下列积分.解本题只要根据复合梯形公式（6.11）及复合Simpson公式（6.13）直接计算即可。
对，取n=8,在分点处计算f(x)的值构造函数表。
按式（6.11）求出，按式（6.13）求得，积分）式估计误差，因，故(1)(2)(3)（1）令代入公式两端并使其相等，得解此方程组得，于是有再令，得故求积公式具有3次代数精确度。
数值分析作业
数值分析作业数值分析作业姓名：王森学号：Z14030271学院：化学工程学院班级：化工8班实验3.1 Gauss消去法的数值稳定性实验实验目的：理解高斯消元过程中出现小主元即很小时引起方程组解数值不稳定性实验题：求解线性方程组实验要求1计算矩阵的条件数判断系数矩阵是良态的还是病态的2用高斯列主元消去法求得L和U及解向量3用不选主元的高斯消去法求得L和U及解向量4观察小主元并分析对计算结果的影响解11.1判断矩阵A1 是否病态程序如下:>>A1=[0.3*10^-15,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1]; cond(A1 =,1)结果如下:ans =136.2945>>A1=[0.3*10^-15,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1]; cond(A1 ,inf)结果如下ans =84.3115>>A1=[0.3*10^-15,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1];cond(A 1,2)结果如下ans =68.4296因为cond(A1 ,1) =136.2945》1；cond(A1 ,2) = 68.4296》1；cond(A1 ,inf) =84.3115》1所以该矩阵A1 是病态矩阵。
1.2判断矩阵A2是否病态程序如下:>> A2=[10,-7,0,1;-3,2.0999********,6,2;5,-1,5,-1;0,1,0,2]; cond(A2,1) ans =19.2832>> A2 =[10,-7,0,1;-3,2.0999********,6,2;5,-1,5,-1;0,1,0,2]; cond(A2,2) ans =8.9939>> A2=[10,-7,0,1;-3,2.0999********,6,2;5,-1,5,-1;0,1,0,2]; cond(A2,inf) ans =18.3564因为cond(A2,1) =19.2832 》1；cond(A2,2) = 8.9939》1；cond(A2,inf)= 18.3564 》1所以该矩阵A2是病态矩阵。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2014-2015(2)计算机与信息工程学院数值分析作业计科专业_______级_____班 姓名:___________学号:____________第一章 绪论 一、单项选择题1.用3.1415作为π 的近似值时具有（ ）位有效数字。
（A ）3 （B ）4 （C ） 5 （D ）62.已知数x 1=721 x 2=0.721 x 3=0.700 x 4=7*10-2是由四舍五入得到的，则它们的有效数字的位数应分别为( )。
(A) 3，3，3，1 ( B) 3，3，3，3 (C) 3，3，1，1 ( D) 3，3，3，2 二、填空题1.在一些数值计算中,对数据只能取有限位表示,如2 1.414≈ ,这时所产生的误差称为_______误差.（填误差的类型）2. 为尽量避免有效数字的严重损失，当1>>x 时，应将表达式x x -+1改写为_________以保证计算结果比较精确.3.在数值计算中，通常取e 2.71= ，此时产生的误差为_________误差（填误差的类型）.4.设x =0.231是精确值x *=0.229的近似值，则x 有_________位有效数字。
三、计算题1、（本题5分）试确定722作为π的近似值具有几位有效数字，并确定其相对误差限。
第二章 插值法 一、单项选择题1. 通过点0011(x ,y ),(x ,y ) 的拉格朗日插值基函数01l (x),l (x)满足 ( ). (A ) 0011l (x )0,l (x )0== ( B) 0011l (x )1,l (x )1== (C )0011l (x )1,l (x )0== (D) 0011l (x )0,l (x )1==2.是给定的互异节点，是以它们为插值节点的插值多项式，则是一个( ).(A) n +1次多项式 (B) n 次多项式(C) 次数小于n 的多项式 (D) 次数不超过n 的多项式 二、填空题1. 设有节点012x ,x ,x ，其对应的函数=y f (x) 的值分别为012y ,y ,y ， 则二次拉格朗日插值基函数0l (x)___________= .2.已知2()1,=+f x x 则[1,2,3]____=f .2. 已知f (1)1,f (2)3,== 那么y f (x)=以1，2为节点的拉格朗日线性插值多项式为_________.3. 当x =1，-1，2时，对应的函数值分别为f (-1)=0，f (0)=2，f (4)=10，则f (x )的拉格朗日插值多项式是 .4. 设2f (x)x = ，则f (x)关于节点012x 0,x 1,x 3=== 的二阶向前差分为___________.5. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商公式的 _____，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的 ___；如果要估计结果的舍入误差，应该选用插值公式中的 ___.6. 设20)2(,10)1(,0)0(===f f f ,则[]___,,=10f []___,2,1,0=f )(x f 的二次牛顿插值多项式为___________________________.7. 设)(x L n 为)(x f 的n 次拉格朗日插值多项式，则其插值余项为_________________.8. 已知53()245,f x x x x =+-则[1,1,0]f -= ,[3,2,1,1,2,3]f ---=_____. 9. 设),,2,1,0(,,53)(2==+=k kh x x x f k 则差商123[,,,]___n n n n f x x x x +++=.10. 设()(0,1,2)j l x j n =是n 次拉格朗日插值多项式的插值基函数，则()j i l x =____________(,0,1,2)i j n =；0()nj j l x ==∑ 。
三、计算题 1（1）写出f (x) 的3次Lagrange 插值多项式3L (x) ； （2）写出f (x) 的3次Newton 插值多项式3N (x) .2. 已知-12 4 5-2 4 5 7 (1) 用拉格朗日插值法求的三次插值多项式；(2) 求x ， 使=0。
3. 给定数据,)(,)(,)(,)(143521100====y y y y 求三次拉格朗日插值多项式)(x L 3.4.已知函数()y f x =在如下节点处的函数值x-1 0 1 2 y143（1） 建立以上数据的差分表；（2） 根据后三个节点建立二阶牛顿后插公式2()P x ，并计算(1.1)y 的近似值；5.已知y=x ，0x =4，1x =9，用线性插值求7的近似值。
6.已知x1 2 3 4 F(x) 021512计算三阶差商f ［1,3,4,7］。
7.已知i x1 3 4 7 f(i x ) 021512求满足以上插值条件的牛顿型插值多项式。
8.设)(x f 为k 次多项式，n x x x x ,,,210为1+n 个互异点，)(x L n 为)(x f 的n 次插值多项式。
若n k <，试证)()(x f x L n ≡。
第三章 函数逼近于计算 一、填空题1.用二次多项式2012(x)a a x a x ,ϕ=++ 其中012a ,a ,a 是待定参数，拟合点1122n n (x ,y ),(x ,y ),,(x ,y ) ，那么参数012a ,a ,a 是使误差平方和____________________取最小值的解。
2.已知数据对=k k (x ,y )(k 1,2,,n) ,用直线y a bx =+拟合这n 个点,则参数a,b 满足的法方程组是__________________. 二、计算题1.已知一组实验数据如下i x 1 23 4 5 i f (x )44.5688.5求它的拟合曲线（直线）.2、已知一组试验数据如下i x 20 40 60 80 100 i f4.35 7.55 10.40 13.80 16.80求它的拟合曲线（直线）。
3．求32f (x)x =在[0,1]上求关于{}span 1,x φ=的一次最佳平方逼近多项式.4.已知如下数据表，试用最小二乘法求它的二次最小平方逼近多项式。
x -1 0 1 2 y125()[]{}245.,111,f x x span x x ϕ=-=在，上求关于的最佳平方逼近。
6. 求x x f =)(在区间[1/4,1]上的关于权函数1=)(x ρ的一次最佳平方逼近多项式.7. 求353323-++=x x x x f )(在区间],[11-上的最佳二次逼近多项式. 8. 已知-2 -1 0 1 2 42135求的形如的二次拟合曲线，并求的近似值。
9.已知n+1个数据点(,)(0,1,2,,)i i x y i n =，请用多种方法建立这些数据点之间的函数关系，并说明各种函数的适用条件。
10.用最小二乘法解下列超定线性方程组：⎪⎩⎪⎨⎧=-=+=+2724212121x x x x x x 11.求3x f(x )=在[0，1]上的一次平方逼近多项式。
第四章 数值积分与数值微分 一、单项选择题1.已知求积公式21121f (x)dx f (1)Af ()f (2)636≈++⎰ ，则A =（ ）.16(A) 13(B) 12(C) 23(D)2.已知n 4= 时牛顿-科特斯求积公式，科特斯系数(4)(4)01716C ,C ,9045== (4)22C 15= ，那么(4)3C =（ ）. (A) 790 (B) 1645 (C) 215 (D) (4)3716239C 190451590=---=3. 已知节点k k (x ,y ),(k 0,1,2,,n),= 插值型两点求导公式是（ ）.1011y (x x )h '≈--(A) 1101y (x x )h '≈--(B) 1011y (y y )h '≈--(C) 1011y (y y )h'≈-(D) 4.求积分公式11f (x)dx f (1)f (1)-≈-+⎰ 是( )次代数精度.( A) 1 (B) 2 (C) 3 (D) 4二、填空题1.求积分公式10211123f (x)dx f ()f ()f ()343234≈-+⎰ 具有_____次代数精度.2.设求积公式nbk k ak 0f (x)dx A f (x )=≈∑⎰ ，若对_______________的多项式积分公式精确成立，而至少有一个m 1+ 次多项式不成立，则称该求积公式具有m 次代数精度.3.已知n 3= 时，科特斯系数(3)(3)(3)01213C ,C C 88=== ，那么(3)3C _____= . 4. 求初值问题00y f (x,y)y(x )y '=⎧⎨=⎩近似解的梯形公式是k 1y +≈___________.5. n 个求积节点的插值型求积公式的代数精确度至少为______次，n 个求积节点的高斯求积公式的代数精度为 .6. 5个节点的牛顿-柯特斯公式代数精度是 .7.1+n 个节点的Gauss 型求积公式具有______次的代数精度.8.为使求积公式1123133()()(0)()33f x dx A f A f A f -≈-++⎰的代数精度尽量高，应使1A =，2A =，3A =，此时公式具有 次的代数精度。
9.数值微分公式)(a f '≈hh a f h a f 2)()(--+的代数精度为_______.三、计算题1. 试用n 1,2,4=的牛顿-科特斯求积公式计算定积分101I dx 1x=+⎰ . 2.已知012113,,,424===x x x(1)推导以这三点为求积节点在[0,1]上的插值型求积公式10120113()()()()424f x dx A f A f A f ≈++⎰；(2)指明求积公式所具有的代数精度；(3)用所求公式计算120⎰x dx .3. 试求使求积公式的代数精度尽量高，并求其代数精度。
4.确定下列求积公式中的待定参数，使其代数精确度尽量高，并指明求积公式所具有的代数精确度.5.已知)(x f 的函数值如下：x1.82.0 2.2 2.4 2.6 )(x f 3.14.46.08.01.00用复合梯形公式和复合辛普森公式求dx x f ⎰6.28.1)(的近似值.6.已知)(x f 的函数值如下表25.15.001)(15.005.01---x f x用复合梯形公式和复合Simpson 公式求dx x f ⎰-11)(的近似值.第五章 常微分方程数值解法一、单项选择题1.解常微分方程初值问题的平均形式的改进欧拉法公式是k 1p c 1y (y y )2+=+,那么p c y ,y 分别为( ).p k k k c k k 1p y y hf (x ,y )y y hf (x ,y )+=+⎧⎨=+⎩(A) p k k k c k k 1k y y hf (x ,y )y y hf (x ,y )+=+⎧⎨=+⎩(B)p k k 1k c k k p y y hf (x ,y )y y hf (x ,y )+=+⎧⎨=+⎩(C) p k k k c k k p y y f (x ,y )y y f (x ,y )=+⎧⎨=+⎩(D)2. 求解常微分方程的二阶R-K 方法的局部截断误差为( ). 1O(h )(A) 2O(h )(B) 3O(h )(C) 4O(h )(D).3.解微分方程初值问题的方法，（ ）的局部截断误差为)(3h O . (A) 欧拉法 (B) 改进欧拉法(C) 三阶龙格—库塔法 (D) 四阶龙格—库塔法 二、计算题1. 写出四阶经典龙格-库塔法求解初值问题y 83yy(0)2'=-⎧⎨=⎩的计算公式,并取步长h 0.2= ,计算y(0.4) 的近似值,小数点后至少保留4位.2.用Euler 方法求解初值问题'(0)0⎧=-⎨=⎩y x yy ，取0.1=h 在区间[0,0.3]计算,结果保留到小数点后4位.3.初值问题 ⎩⎨⎧=>+='0)0(0,y x b ax y 有精确bx ax x y +=221)(，试证明： 用Euler法以h 为步长所得近似解n y 的整体截断误差为n n n n ahx y x y 21)(=-=ε 4.写出用四阶经典的龙格—库塔方法求解下列初值问题的计算公式：（无需计算）⎩⎨⎧=+=1)0(')1y y x y ,10<<x ⎪⎩⎪⎨⎧=+=1)0(13')2y x y y ,10<<x5.用改进欧拉法求解⎩⎨⎧=+=1)0('y yx y )10(≤≤x ，2.0=h ，取两位小数。