数值计算方法第四章参考答案
数值方法课后习题答案第4章
![数值方法课后习题答案第4章](https://img.taocdn.com/s3/m/1a8d722558fb770bf78a55c0.png)
解线性方程组迭代法
第四章 解线性方程组迭代法
习题4-1
第四章
解线性方程组迭代法
Байду номын сангаас
第四章
解线性方程组迭代法
第四章
解线性方程组迭代法
第四章
解线性方程组迭代法
第四章
解线性方程组迭代法
第四章
解线性方程组迭代法
第四章
解线性方程组迭代法
第四章
解线性方程组迭代法
第四章
习题4
习题4
习题4
习题4
8.
设A为严格对角优势阵,证明:
习题4
9. A是n阶非奇异阵,B是n阶奇异阵,试求证:
习题4
习题4
P91
P91.
x0
p0 r0
Ap0
x1
r1
p1
Ap1
x2
r2
0
3
7
30/29=
17/29=
1360/841=
1530/841=
14/9=
0.3
P91
1.034482758 0 1 8 10/29= 0.344827586
0.570796875 0.493315839 0.500166165 0.499999398 1.001438281 0.998173633 1.000074653 1.000013383 -0.49943416 -
0.500558834 0.499923587 0.500003961
w=1.03
10 29 a0 =10/29=0.344827586
2890/841=3.436385254 260100/24389=10.66464388 a1 =8381/26010=0.322222222 -289/29= -9.965517218 b0 =289/841=0.343638524
数值计算课后答案4
![数值计算课后答案4](https://img.taocdn.com/s3/m/c3825a01ee06eff9aef8077f.png)
习 题 四 解 答1、设010,1x x ==,写出()x f x e -=的一次插值多项式1()L x ,并估计插值误差。
解:根据已知条件,有设插值函数为1()L x ax b =+,由插值条件,建立线性方程组为解之得111a eb -⎧=-⎨=⎩则11()(1)1L x e x -=-+ 因为(),()x x y x e y x e --'''=-= 所以,插值余项为 所以010101()max max (1)2111248x r x e x x e ξξ-≤≤≤≤-≤-=⨯⨯=。
2、给定函数表选用合适的三次插值多项式来近似计算f(0.2)和f(0.8)。
解:设三次插值多项式为230123()f x a a x a x a x =+++,由插值条件,建立方程组为 即解之得则所求的三次多项式为23()0.41 6.29 3.489.98f x x x x =--+。
所以3、设(0,1,2,,)i x i n =是 n+1个互异节点,证明:(1)0()(0,1,2,,)n k k i i i x l x x k n ===∑;(2)0()()0(0,1,2,,)nk i i i x x l x k n =-==∑。
证明: (1)由拉格朗日插值定理,以x 0,x 1,x 2,…x n 为插值节点,对y=f(x)=x k 作n 次插值,插值多项式为0()()nn i i i p x l x y ==∑,而y i =x i k ,所以0()()()nnk n i i i i i i p x l x y l x x ====∑∑同时,插值余项 所以()nk k i i l x x x =∑结论得证。
(2)取函数()(),0,1,2,,k f x x t k n =-=对此函数取节点(0,1,2,,)i x i n =,则对应的插值多项式为0()()()nk n i i i p x x t l x ==-∑,由余项公式,得(1)(1)011()()()()()()()()0(1)!(1)!nn kk n ki i i r x x t x t l x f x x t x n n ξξππ++==---==-=++∑所以令t=x ,4、给定数据(()f x =(1)试用线性插值计算f(2.3)的近似值,并估计误差;(2)试用二次Newton 插值多项式计算f(2.15)的近似值,并估计误差。
《数值分析》第四章答案
![《数值分析》第四章答案](https://img.taocdn.com/s3/m/834bea40be1e650e52ea994f.png)
习题41. 给定x x f =)(在144,121,100=x 3点处的值,试以这3点建立)(x f 的2次(抛物)插值公式,利用插值公式115求的近似值并估计误差。
再给13169=建立3次插值公式,给出相应的结果。
解:x x f =)( 2121)(-='x x f ,2341)(--=''x x f ,2583)(-='''x x f ,27)4(1615)(--=x x f,72380529.10)115(=f1000=x , 1211=x , 1442=x , 1693=x 100=y , 111=y , 122=y , 133=y))(())(())(())(())(())(()(1202102210120*********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 ----+----+----= )121144)(100144()121115)(100115(12)144121)(100121()144115)(100115(11)144100)(121100()144115)(121115(10)115(2----⨯+----⨯+----⨯=L=2344)6(1512)23(21)29(1511)44)(21()29)(6(10⨯-⨯⨯+-⨯-⨯⨯+----⨯72276.1006719.190683.988312.1=-+=))()((!3)()()(2102x x x x x x f x L x f ---'''=-ξ ,144100<<ξ )44115()121115()100115()(max 61)115()115(1441002-⨯-⨯-⋅'''≤-≤≤x f L f x 296151083615⨯⨯⨯⨯⨯≤-001631.0101631.02=⨯=- 实际误差 22101045.0)115()115(-⨯=-L f))()(())()(())()(())()(()(312101320130201032103x x x x x x x x x x x x y x x x x x x x x x x x x y x L ------+------= ))()(())()(())()(())()((23130321033212023102x x x x x x x x x x x x y x x x x x x x x x x x x y ------+------+ )169100()144100()121100()169115()144115()121115(10)115(3-⨯-⨯--⨯-⨯-⨯=L )169121()144121()100121()169115()144115()100115(11-⨯-⨯--⨯-⨯-⨯+)169144()121144()100144()169115()121115()100115(12-⨯-⨯--⨯-⨯-⨯+)144169()121169()100169()144115()121115()100115(13-⨯-⨯--⨯-⨯-⨯+)48()23(21)54()29(1511)69()44()21()54()29()6(10-⨯-⨯-⨯-⨯⨯+-⨯-⨯--⨯-⨯-⨯= 254869)29()6(1513)25(2344)54()6(1512⨯⨯-⨯-⨯⨯+-⨯⨯-⨯-⨯⨯+ 723571.10409783.0305138.2145186.11473744.1=+-+= ))()()((!4)()()(3210)4(3x x x x x x x x f x L x f ----=-ξ,169100<<ξ)169115)(144115)(121115)(10115(101615241)115()115(73----⨯⨯⨯≤--L f )54()29()6(151016152417-⨯-⨯-⨯⨯⨯⨯=- 0005505.0105505.03=⨯=-实际误差 321023429.0)115()115(-⨯=-L f 2. 设j x 为互异节点),,1,0(n j =求证: (1)k nj j k j x x l x =∑=)(0),,1,0(n k =;(2)0)()(0=-∑=x l x x j knj j ),,1(n k =。
数值计算方法第2版 第4章 插值法
![数值计算方法第2版 第4章 插值法](https://img.taocdn.com/s3/m/54b903e9da38376baf1faefa.png)
则
l ( x ) 1 ( k i ) , k i l ( x ) 0 ( k i ) , i 、 k 0 , 1 , , n k i
lk (x)称为关于节点xi( i=0,1,…,n)的n次插值基函数。
基函数的特点
1. 基函数的个数等于节点数。 2. n+1个节点的基函数是n次代数多项式。 3. 基函数和每一个节点都有关。节点确定,基函数就唯 一的确定。 4. 基函数和被插值函数无关。 5. 基函数之和为1。
公式的结构:它是两个一次函数的线性组合 线性插值基函数
x x 1 l ( x ) , 0 x x 0 1 x x 0 l ( x ) 1 x x 1 0
3 线性插值的几何意义 用直线 P ( x ) 近似代替被插值函数 f ( x ) 。
例
造数学用表。平方根表
给定函数在100、121两点的平方根如下表,试用线性 插值求115的平方根。 x 100 121
其系数行列式
a0 a1 x0 a2 x02 an x0n y0 2 n a0 a1 x1 a2 x1 an x1 y1 2 n a a x a x a x n n yn 0 1 n 2 n
1 x 0 x 02
x 0n
2 n 1 x x x 1 1 ( x x 0 V ( x , x , , x ) 1 i j) 0 1 n 0 j i n
1 xn
x n2 x nn
,a , ,a 0 1 n ,因此P(x)存在且唯一。 方程组有唯一解 a
唯一性说明不论用那种方法构造的插值多 项式只要满足相同的插值条件,其结果都是互 相恒等的。 推论 当f(x)是次数不超过n的多项式时, 其n次插值多项式就是f(x)本身。
数值计算方法第四章习题部分参考答案
![数值计算方法第四章习题部分参考答案](https://img.taocdn.com/s3/m/ebce181a5a8102d276a22fd7.png)
3) 快速弦截法 相应的弦截迭代公式为:
f x k x x x x k 1 k k k 1 f x f x k k 1
x k e 4 c o s x k x x x k k k 1 x x k k 1 e 4 c o s x e 4 c o s x k k 1
2) 弦截法:取 x0
4
, x1
2
相应的弦截迭代公式为:
f xk xk1 xk xk x0 f xk f x0 xk e 4co sxk xk xk x0 exk 4co sxk (e4 4 2 ) 2
因此有 x 1 . 3 , 1 . 6
2 x3
x 又:
,易知 x 为单调递减函数,所以有
2 x 0 . 9 1 0 3 1 3 ( 1 . 3 )
由压缩影像定理知该迭代式收敛。
2) 对于该迭代式,相应的迭代函数为:
x 3 1x2
利用公式作迭代得:
x 1 e 4 c o s x 1 x x x x 0 . 8 7 7 0 0 3 2 1 1 0 x 1 4 e 4 c o s x ( e 4 2) 1 2
x2 e 4cosx2 x3 x2 x2 x0 0.906360 ex2 4cosx2 (e4 4 2 ) 2 x3 e 4cosx3 x4 x3 x3 x0 0.904701 ex3 4cosx3 (e4 4 2 ) 2
2 2 3 x 1 . 6 1 1 . 6 0 . 5 6 4 1 1 3
黄云清版数值计算方法习题解答
![黄云清版数值计算方法习题解答](https://img.taocdn.com/s3/m/2e309b3c2cc58bd63086bd56.png)
第一章 引论〔习题〕2.证明 : 记 x x f =)( ,则)()(***x x x x x xx x f E r +-=-=)(21**x E x x x x x xr ≈-⋅+=.3.证明: 令: )()()(b a fl b a fl b a **-*=δ可估计: 1|)(|-≥*c b a fl β 〔c 为b a *阶码〕,故: 121||--≤c t c ββδt-=121β 于是:)1)()(δ+*=*b a b a fl .4.解 (1) )21()122x x x++.(2) )11(2x x x x x-++.(3) xx x x x x x cos 1sin )cos 1(sin cos 12+≈+=-.6.解 a 的相对误差:由于31021|)(|-⋅≤-=a x x E .x ax x E r -=)(, 221018110921)(--⋅=⨯≤x E r . 〔1Th 〕)(a f 对于)(x f 的误差和相对误差.|11||)(|a x f E ---==()25.021011321⨯⋅≤-+---ax x a =310-33104110|)(|--⨯=-≤a f E r .9.解 递推关系: 1101.100-+-=n n n y y y (1) 取初值 10=y , 01.01=y 计算 可得: 11001.10022-⨯=-y 10001.1-=410-=6310-=y , 8410-=y , 10510-=y , …(2)取初值 50101-+=y , 2110-=y ,记: n n n y y -=ε,序列 {}n ε ,满足递推关系,且 5010--=ε , 01=ε1101.100-+-=n n n εεε, 于是: 5210-=ε, 531001.100-⨯=ε,55241010)01.100(---⨯=ε,55351002.20010)01.100(--⨯-⨯=ε,可见随着 n ε 的主项5210)01.100(--⨯n 的增长,说明该递推关系式是不稳定的.第二章 多项式插值 (习 题)1.方法一. 由 Lagrange 插值公式)()()()()(332211003x l f x l f x l f x l f x L ⋅+⋅+⋅+⋅=)1)((31)2)()(1()1)(()(2123210---=-----=x x x x x x x l ,))(1(2)1)()(1()(21221211--=--+=x x x x x x l ,x x x x x x l )1()()1()1!()(2382121232--=-⋅⋅-+=,)()1(12)()1()(2121213-+=⋅⋅-+=x x x x x x x l . 可得: )21()(23-=x x x L方法二. 令:)()21()(3B Ax x x x L +-=由 23)1(3-=-L , 21)1(3=L , 定A ,B 〔称之为待定系数法〕2.证明(1) 由于 j i j i x l ,)(δ= 故: =)(x L n ∑=ni i k i x l x)( ,当 j x x = 时有: k j j n x x L =)( , n j ,,1,0 =)(x L n 也即为 k x 的插值多项式,由唯一性,有:∑==ni k i ki x x l x)( , n k ,,1,0 =证明(2):利用Newton 插值多项式)(],[)()(0100x x x x f x f x N n -+=)()(],,[100---++n n x x x x x x f)()()()()()(00101x l x x x x x x x x x f n n =----=差商表:f(x) 一阶 二阶 … n 阶差商0x 11x 0101x x -0 )()(11020x x x x --n x 00)()(1010n x x x x --代入)(*式有:)()()()()(1)(020*******n n n x x x x x x x x x x x x x x x N -----++--+=- .)(0x l 为n 次代数多项式,由插值多项式的唯一性:有 )()(0x N x l n ≡.4.解作)(x f 以b a a ,,ε+为节点的Lagrange 插值多项式,有:)()()(22x R x L x f +=, 其中:)()()()()()()()()()(2εεεεε+-+--+-----=a fb a b x a x a f b a b x a x x L)()()()()(b f a b a b a x a x εε------+, )()()(!3)()(2b x a x a x f x R ----'''=εζ , b a <<ζ令: 0→ε 有)()(6)()()(22b x a x f x R x R --'''=→ζ, 又:)()()()([)()(2a f a b ax a f a b a x x b x L εεεεε----+----= )]()()()()(a f a b a x a f a b a x -------+εεεε )()()()()(b f a b a b a x a x εε------+)()()2()(2a f ab a b x x b --+-→)()()()(a f a b a x x b '---+ )()()()(22x P b f a b a x =--+ 故当0→ε 时,成立公式: )()()(x R x P x f +=.5.解:因为34)(3'-=x x f ,2''12)(x x f =)(x f 为凹函数.又从数值表可见:当]5.0,1.0[∈x 时,)(x f 单调下降.有反函数)(1y fx -=)(y f的Newton 插值多项式:)17440.0)(10810.0)(40160.0)(70010.0(01225.0)10810.0)(40160.0)(70010.0(01531.0)40160.0)(70010.0(0096436.0)70010.0(33500.01.0)(4+---+------+--=y y y y y y y y y y y N.337.0)0(4*≈=N x7.解 1)(37++=x x x f .有:=]2,,2,2[71f !7)()7(ξf =1,!8)(]2,,2,2[)8(810ηf f = 0=.9.证明:(1) =⋅-⋅=⋅∆++i i i i i i g f g f g f 11)(i i i i i i i i g f g f g f g f ⋅-⋅+⋅-⋅++++1111i i i i f g g f ∆+∆=+1.(3) n xnn)1()1(-=∆!)()(nh x h x x hn++此题可利用数学归纳法:设 k n = 成立,证明 1+=k n 成立.又 1=n 时是成立的.10.证明: 记: 2]2/)1([)(+=n n n f ,33321)(n n g +++=有: 3)1()()1()(+=-+=∆n n f n f n f 故: ∑-=∆=10)()(n k k f n g ∑-=-+=1)]()1([n k k f k f2]2/)1([)0()(+=-=n n f n f .13.解 作重节点差商的Newton 插值公式)1(]1,1[)1()(+--+-=x f f x P 22)1(]1,0,1,1[)1(]0,1,1[+--++--+x x f x f )1()1(]1,1,0,1,1[2-+--+x x x f重节点差商表:i x i f 一阶 二阶 三阶 四阶10-=x 110-=x 1 201=x 1 0 -212=x 1 0 0 112=x 1 2 2 1 0得 22)1()1(2)1(21)(+++-++=x x x x x P 13+-=x x .17.证: 取 ,00=x 211=x , 12=x , 21=h00=f , 11=f , 12=f记: )(i i x s M ''= , 2,1,0=i有 h x x M h x x M x S 01101)(-+-=''x M x M 102)21(2+-= )21(2)1(2)(212-+-=''x M x M x S 又三弯矩方程为:(2],,[210-=x x x f )244210-=++M M M , )24(41201M M M ++-=.分段积分:⎰⎰+''=''∆121221)]([)]([dx x s dx x s ⎰''12221)]([dx x s ⎰+-+=210201)]21([4dx x M x M ⎰-+-121221)]21()1([4dx x M x M⎰⎰-+-+-+-=121121221201)]21()1([4)]1()21([4dx x M x M dx x M x M由于⎰=-1212241)21(dx x ,⎰=-1212241)1(dx x ,⎰=--121481)1()21(dx x x ,于是:⎰++++=''∆1022212110202]2[61))((M M M M M M M dx x S 又: )24(41201M M M ++-=记 =),(20M M I ⎰∆''12))((dx x S=)()24(41[6120202220M M M M M M +++-+ ])24(81220M M +++ 由00=∂∂M I , 02=∂∂M I. 得: ⎩⎨⎧=+-=-07072020M M M M 即当: 020==M M 时,),(20M M I 达最小故:⎰=⋅⋅≥''∆12212)24(8161))((dx x S ,由最小模原理: ⎰≥''1212)]([dx x f .20.解 利用三弯矩方法 )(i i x s M ''= , 2,1,0=i10=x , 22=x , 32=x⎪⎩⎪⎨⎧-=+=++=+542364622121010M M M M M M M解得: 70-=M , 201=M , 372-=M]2,1[∈x 72431729)(231-+-=x x x x s ]3,2[∈x 105229367219)(232+-+-=x x x x s .第三章 最佳逼近与其实现 (习 题)2.解(1) ⎰'⋅'=badx x g x f g f )()(),( 不是 ),(b a c '中的内积,事实上容易验证:),(),(f g g f = , ),(),(g f g f λλ=),(),(),(w g w f w g f +=+但是0),(=f f 当且仅当 0)(≡x f . 条件不满足,因为:⎰='⋅'=badx x f x f f f 0)()(),(推出0)(≡'x f ,0)(≠=const x f . 因而 ),(g f 不是 ),(b a C '中的内积.(2) ),(g f 是 =],[10b a C {}],[)(,0)(:)(b a C x f a f x f '∈'=空间的内积,这是因为:0),(=f f 推出 0)(='x f ,C x f =)(,又],[10b a C f ∈ ,故 0)(=x f .4.解:由于 0)(],,[2≠''∈x f b a c f ,则)(x f ''于],[b a 上保号,由定理5的推论2可知:)()(1x P x f -的交错点组恰有三个交错点,且 a x =1,b x =3,即: ⎪⎪⎩⎪⎪⎨⎧=-'='-=+-==+-==+-=0)()(,)()()(,)()()(,)()()(122210223103311011αρααρααρααx f x e x x f x e x x f x e x x f x e 故:ab a f b f x f --='=)()()(21α,2)()(2)()(220x a a b a f b f x f a f +⋅---+=α记 c x =2 ,即证得(1).(2) 若 x x f cos )(= ,]2,0[],[π=b a此时由 ab a f b fc f --=')()()( 得:π2sin =c , )2sin(πarc c =,πα21-=πππα2)4(2120-+=2)/2sin(2ππarc ⋅+)4(212-+=πππππ)2sin(arc +. 误差估计:)()(10b f b f E -+=-=ααρ)4(212-+=πππ1)2sin(-+ππarc5.解:选取 α ,使得:=)(αI ||max 211x x x α-≤≤ ,达到极小,即要求 x x *)(*αϕ= ,于]1,0[上一致逼近于2x ,如图 应选*α ,使得: x x x *)(2αϕ-=,于 ]1,0[ 上有两个轮流为正负偏差点,其中之一为1,另一个假设为 ζ 于是: )()1(ζϕα-=,0)(='ζϕ , 〔 ζ为)(x ϕ的极值点〕 得: αζζα+-=-2102=-αζ解得:ζα2= ,0122=-+ζζ,212,1±-=ζ取 12-=ζ , 222-=α. 又: α 是唯一的.6.证明:由最佳一致逼近的特征定理,)(*x P n 为)(x f 的最佳一致逼近多项式,则存在2+n 个点b x x x a n ≤<<<≤+110使得: )()()(*k n k k x P x f x e -==*)1(n k P f --σ.又由于 ],[)(b a C x f ∈ ,于 ),(1+i i x x 中有一个点 i η ,1+<<i i i x x η ,使得: 0)()()(*=-=i n i i P f e ηηη,n i ,,1,0 =即: )(*x P n 为)(x f 满足插值条件: )()(*i i n f P ηη= , n i ,,1,0 =的插值多项式.7.解:求C*,使得:C x f C I bx a R C -=≤≤∈)(max min *)(记 C x f x e -=)()(, 依最佳一致逼近的特征定理:应取)](min )(max [21*],[],[x f x f C b a b a +=*)()(C x f x e -=于 ],[b a 才有两个轮流正负的偏差点,〔即 )(x f 于],[b a 上的最大值点和最小值点〕1x ,2x)(max )(],[1x f x f b a = , )(min )(],[2x f x f b a =此时: *)(max )1()(],[C x f x e b a ii --=σ即 *C 为)(x f 的零次最佳逼近多项式.8.解: 436)(23+++=x x x x f2)(34)3(62031T T T T +++=014T T ++01232112112323T T T T +++= 因为 )(413x T 与零偏差最小,故:012221121123)(T T T x P ++=421132++=x x .为)(x f 的最佳一致逼近多项式.9. 证明:我们仅证明)(x f 是偶函数时,)(x P n 亦是偶函数.由于)(x P n 为)(x f的最佳一致逼近多项式,有:)()()(max ],[f E x P x f n n a a =--和: [,max ()()()]n n a af x P x E f ----=即: )()()(max ],[f E x P x f n n a a =---)(x P n -亦是)(x f 的最佳一致逼近多项式,由最佳一致逼近多项式的惟一性,有: )()(x P x P n n =-即: )(x P n 为偶函数.11.解: 设 x a a x P 10*1)(+= , 2210*2)(x b x b b x P ++=分别为)(x f 的一次、二次最佳平方逼近多项式。
数值分析答案第四章
![数值分析答案第四章](https://img.taocdn.com/s3/m/277e229951e79b896802260d.png)
令
f (x) = x ,则
0 = −1 + 2 x1 + 3 x2
令 f ( x ) = x 2 ,则
2 2 = 1 + 2 x12 + 3 x2
从而解得
⎧ x1 = −0.2899 ⎧ x1 = 0.6899 或⎨ ⎨ ⎩ x2 = 0.5266 ⎩ x2 = 0.1266
令 f ( x ) = x 3 ,则
∫
1
−1
f ( x)dx = ∫ x3 dx = 0
−1
1
[ f ( −1) + 2 f ( x1 ) + 3 f ( x2 )] / 3 ≠ 0
故
∫
1
−1
f ( x)dx = [ f (− 1) + 2 f ( x1 ) + 3 f ( x2 )] / 3不成立。
h
因此,原求积公式具有 2 次代数精度。 (4)若
7 h T8 = [ f ( a) + 2∑ f ( xk ) + f ( b)] = 0.11140 2 k =1
复化辛普森公式为
7 7 h S8 = [ f ( a) + 4∑ f ( x 1 ) + 2∑ f ( xk ) + f ( b)] = 0.11157 k+ 6 k=0 k =1 2 1
令 f ( x ) = x 2 ,则
b 1 3 3 2 f ( x ) dx = ∫a ∫a x dx = 3 (b − a ) b −a 1 3 3 [7 f ( x0 ) + 32 f ( x1 ) + 12 f ( x2 )+ 32 f ( x (b − a ) 3 )+ 7 f ( x 4 )]= 90 3 b
计算方法各章习题及答案
![计算方法各章习题及答案](https://img.taocdn.com/s3/m/1e9e77f90066f5335b812119.png)
第二章数值分析4^92.1 已知多项式通过下列点:1 3答案:q(x) = p(x) -r(x) X5X4X3-3X 1 .2 22.2观测得到二次多项式2的值:表中p2(x)的某一个函数值有错误,试找出并校正它.答案:函数值表中p2(-1)错误,应有p2(-1) = 0 .2.3利用差分的性质证明12■ 22■川,n2=n(n ■ 1)(2n ■ 1)/6.2.4当用等距节点的分段二次插值多项式在区间[-1,1]近似函数e x时,使用多少个节点能够保证误差不超过丄10-6.2答案:需要143个插值节点.2.5 设被插值函数f (x) • C4[a,b] , H3h)(x)是f (x)关于等距节点b — aa ^Xo :::捲:::川:::x n=b的分段三次艾尔米特插值多项式,步长h .试估计n ||f(x)-H3h)(x)||::.答案:||住)-出5)仪川:乞令人4.384第三章函数逼近3.1求f(x)二sinx,x,[0,0.1]在空间门=span{1,x, x2}上最佳平方逼近多项式,并给出平方误差.答案:f (x) =sin X的二次最佳平方逼近多项式为sin x p2(x) = -0.832 440 7 105 1.000 999 1x - 0.024 985 1x2,二次最佳平方逼近的平方误差为20.12 12■ = 0 (sinx) - P 2(x))2dx =0.989 310 7 10•3.2确定参数a,b 和c ,使得积分1 ---------------------------2 1 I (a,b,c)[ax 2 bx c -1 -x 2]dx 取最小值.J 1 — x 2810答案:a, b = 0, c =3 二3 二3.3 求多项式f (x) =2x 4 x 3 5x 2 1在[-1,1]上的3次最佳一致逼近多项式p(x)-答案:f (x)的最佳一致逼近多项式为p(x) = X ’ 7x2 3.43.4用幕级数缩合方法,求 f(x)=e x (―1兰XW1)上的3次近似多项式 p 6,3(x),并估计 || f(X )-P 6,3(X )II ::.答案:p5,3(x) =0.994 574 65 + 0.997 395 83x+0.542 968 75x 2 十 0.177 083 33x 3, || f (x) - p 6,3 (x) |^<0.006 572 327 71 一3.5 求f (x) -e x ( -1乞x 乞1)上的关于权函数「(X )-的三次最佳平方逼近小-x 2多项式 Q(x),并估计误差 || f(x)-$(x)||2 和 || f(x)-S 3(x) ||::.答案:§3(x) =0.994 571 0.997 308x 0.542 991x 20.177 347x 3,|| f (x) -S 3(x) ||2 = 0.006 894 83, || f (x) - §3(x)||严 0.006 442 575.第四章数值积分与数值微分14.1用梯形公式、辛浦生公式和柯特斯公式分别计算积分 X n dx (n -1,2,3,4),并与精确值比较.答案:计算结果如下表所示I 2 0. 5 0. 333 333 0. 250 000 0. 208 333 I 30. 5 0. 333 333 0. 250 000 0. 200 000 精确值0. 50. 333 3330. 250 0000. 200 0004.2 确定下列求积公式中的待定参数,使得求积公式的代数精度尽量高,并指明所确定的求积公式具有的代数精度.h(])仁 f (x)dx 止 A_i f (-h) + A f (0) + A f (h)11 (2)J(x)dx: 3【f(-1) 2f(X i ) 3f(X 2)]hh2⑴ of(x)dxVf(O) f(h)「h[f g f(h)]答案:(1)具有三次代数精确度 (2)具有二次代数精确度 (3)具有三次代数精确度. 4.3 设h = % - X 0,确定求积公式r (x - x o ) f (x)dx = h 1 2[ Af (x o ) + Bf (x i )] + h 3[C 「(x o ) + Df^)] + R[ f ]xo中的待定参数 A, B,C, D ,使得该求积公式的代数精确度尽量高,并给出余项表达式.37 1 if 4)(叮)6答案:A = —, B— ,C —, D — , R[f]=— _) h ,其中 (x o ,xi).202030 20 14404.4设P 2(x)是以0,h,2h 为插值点的f(x)的二次插值多项式,用F 2(x)导出计算积分3h3 4 5If (x)dx 的数值积分公式I h ,并用台劳展开法证明:I - l h h f (0) O(h ).力83h3答案:I h P 2(x)dx h[ f(0) 3f (2h)].0 4(3)取7个节点处的函数值.1sin x4.6用变步长的复化梯形公式和变步长的复化辛浦生公式计算积分Idx .要x1o 1«求用事后误差估计法时,截断误不超过10和 10 .1(1) 运用复化梯形公式计算上述积分值,使其截断误差不超过丄10」. 2(2) 取同样的求积节点,改用复化辛浦生公式计算时,截断误差是多少?(3) 要求的截断误差不超过10“ ,若用复化辛浦生公式,应取多少个节点处的函数值? 答案:(1)只需n — 7.5,取9个节点,I : 0.9464.5 给定积分I 二1sin xdx|R n [f]耳一孟宀皿盂日中0.271估2 2答案:使用复化梯形公式时,I T^ 0.946满足精度要求;使用复化辛浦生公式时,I s4 =0.946 083满足精度要求.4.7 ( 1 )利用埃尔米特插值公式推导带有导数值的求积公式1 323 1 3>5.2用矩阵的直接三角分解法解方程组 广1 0 2 0、「5、0 10 1 X 2312 4 3X 3仃10 1 0 3丿 g<7;答案: &=2 , x 3 = 2 , x 2 = 1, X| = 1 .ba f(x)dx 二 其中余项为b —a(b 「a)2[f(a)f(b)] — ' 丿[f (b)-f (a)] R[f], 2 12R[f]=U 54!30 f ( 4()),(a,b).其中(2)利用上述公式推导带修正项的复化梯形求积公式h 2 f(x)dx :T^—[ f (X N ) - f (x 。
数值计算方法第四章插值1
![数值计算方法第四章插值1](https://img.taocdn.com/s3/m/95f6200f964bcf84b9d57bd6.png)
代数插值
代数插值
当f(x)是次数不超过n的多项式时,给定n+1个节点,其n次插值多项式就是f(x)本身.
代数插值几何意义
拉格朗日插值 逐次线性插值 牛顿插值 等距节点插值 反插值 埃尔米特插值 分段插值法 三次样条插值
拉格朗日插值 线性插值
格朗日插值 抛物线插值
基函数之和为1.
拉格朗日插值 n次插值
当插值点x∈(a,b)时称为内插,否则称为外插。
内插的精度高于外插的精度。
拉格朗日插值余项
余项 设函数f(x)在包含节点x0 , x1 ,…, xn的区间[a,b]上有n+1阶导数,则
拉格朗日插值
活动14
写出3次拉格朗日插值多项式及余项
拉格朗日插值
拉格朗日插值
作业5
已知函数表
应用拉格朗日插值公式计算f(1.300)的近似值.
数值计算方法
苏 强
江苏师范大学连云港校区
数学与信息工程学院 E-mail: 412707233@
数值计算方法 第四章 插值与曲线拟合
没有明显的解析表达式
使用不便的解析表达式
简单函数代替
插值问题
插值问题
代数插值 插值函数
被插值函数 插值节点
插值区间
三角多项式插值 有理函数插值
代数插值
抛物线插值
三点插值
拉格朗日插值 抛物线插值
抛物线插值
三点插值
拉格朗日插值 抛物线插值
拉格朗日插值 n次插值
称为关于节点
的n次插值基函数.
拉格朗日插值n次插值
基函数的个数等于节点数.
n+1个节点的基函数是n次代数多项式 基函数和每一个节点都有关。节点确定,基函数就唯一的确定。 基函数和被插值函数无关
数值计算方法(宋岱才版)课后答案
![数值计算方法(宋岱才版)课后答案](https://img.taocdn.com/s3/m/350bcd17a76e58fafab003f9.png)
第一章 绪论一 本章的学习要求(1)会求有效数字。
(2)会求函数的误差及误差限。
(3)能根据要求进行误差分析。
二 本章应掌握的重点公式(1)绝对误差:设x 为精确值,x *为x 的一个近似值,称e x x **=-为x *的绝对误差。
(2)相对误差:r e e x***=。
(3)绝对误差限:e x x ε***==-。
(4)相对误差限:r x x xxεε*****-==。
(5)一元函数的绝对误差限:设一元函数()()()0,df f x f x dx εε***⎛⎫==⋅ ⎪⎝⎭则。
(6)一元函数的相对误差限:()()1r df f x dx f εε****⎛⎫=⋅ ⎪⎝⎭。
(7)二元函数的绝对误差限:设一元函数()()(),0,f f x y f y y εε***⎛⎫∂==⋅ ⎪∂⎝⎭则。
(8)二元函数的相对误差限:()()()1r f f f x y x y f εεε******⎡⎤⎛⎫∂∂⎛⎫⎢⎥=⋅+⋅ ⎪ ⎪∂∂⎝⎭⎢⎥⎝⎭⎣⎦。
三 本章习题解析1. 下列各数都是经过四舍五入得到的近似值,(1)试指出它们有几位有效数字,(2)分别估计1123A X X X ***=及224X A X **=的相对误差限。
12341.1021,0.031,385.6,56.430x x x x ****====解:(1)1x *有5位有效数字,2x *有2位有效数字,3x *有4位有效数字,4x *有5位有效数字。
(2)1111123231312123,,,,A A AA x x x x x x x x x x x x ∂∂∂====∂∂∂由题可知:1A *为1A 的近似值,123,,x x x ***分别为123,,x x x 近似值。
所以()()111rA A Aεε***=()()()12311111123A A A x x x A X X X εεε*******⎡⎤⎢⎥=++⎢⎥⎢⎥⎣⎦⎛⎫⎛⎫⎛⎫∂∂∂ ⎪ ⎪ ⎪∂∂∂⎝⎭⎝⎭⎝⎭43123131212311111010100.215222x x x x x x x x x **-**-**-***⎡⎤=⨯⨯+⨯⨯+⨯⨯=⎢⎥⎣⎦()222222424441,,,X A Ax A X x x x x ∂∂===-∂∂则有同理有2A *为2A 的近似值,2x *,4x *为2x ,4x 的近似值,代入相对误差限公式:()()222rA A Aεε***=()()24212224A A X X A X X εε*****⎡⎤⎢⎥=+⎢⎥⎢⎥⎣⎦⎛⎫⎛⎫∂∂ ⎪ ⎪∂∂⎝⎭⎝⎭()33542224411*********X X X X X **--***⎡⎤⎢⎥=⨯⨯+⨯⨯=⎢⎥⎣⎦2. 正方形的边长大约为100cm ,怎样测量才能使其面积误差不超过21cm ? 解:设正方形的边长为x ,则面积为2S x =,2dsx dx=,在这里设x *为边长的近似值,S *为面积的近似值:由题可知:()()1ds s x dx εε***=≤⎛⎫ ⎪⎝⎭即:()21x x ε**⋅≤ 推出:()10.005200xcm ε*≤=。
数值计算方法第4章4-06反插值
![数值计算方法第4章4-06反插值](https://img.taocdn.com/s3/m/465757b88762caaedc33d41b.png)
(0 1)(0 4)
(4 1)(4 0)
x 3 ,代入
p(3) (3 0)(3 4) 0 (3 1)(3 4) 2 (3 1)(3 0) 10 8
(1 0)(1 4)
(0 1)(0 4)
(4 1)(4 0)
(2)由于 f (x) 是单调连续函数,用反插值,将函数表转换成反
函数表
y f (x)
-1
0
2
10
1
x f 1(y) - 3
-1
0
4
?
已知连续函数 f ( x) 在x 1,0,2,3 的值分别是-4,-1,0,3,用牛
顿插值求(1) f (1.5) 的近似值。(2) f ( x) 0.5 时,x 的近似值。
解 (1)根据已知条件列表
x
-3 -1
0
4
3
f (x) - 1
0
2
10
?
取靠近 3 的三个节点- 1,0,4,作拉格朗日二次插值
p(x) (x 0)( x 4) 0 (x 1)( x 4) 2 (x 1)( x 0) 10 将
(1 0)(1 4)
y f (x) 则 x f 1 ( y) ,有函数表。
y
-4
-1
0
3
0.5
x
-1
0
2
3
?
根据已知 x f 1 ( y) 的函数值,构造差商表。
y
x
-4
-1
-1
0
0
2
3
3
牛顿插值多项式
一阶
1/ 3 2
1/ 3
数值计算方法课后习题答案(李庆扬等)
![数值计算方法课后习题答案(李庆扬等)](https://img.taocdn.com/s3/m/ee6dcb310b4c2e3f5727637c.png)
第一章 绪论(12)1、设0>x ,x 的相对误差为δ,求x ln 的误差。
[解]设0*>x 为x 的近似值,则有相对误差为δε=)(*x r ,绝对误差为**)(x x δε=,从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x , 相对误差为****ln ln )(ln )(ln x x x x rδεε==。
2、设x 的相对误差为2%,求n x 的相对误差。
[解]设*x 为x 的近似值,则有相对误差为%2)(*=x r ε,绝对误差为**%2)(x x =ε,从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε,相对误差为%2)()(ln )(ln ***n x x x nr==εε。
3、下列各数都是经过四舍五入得到的近似数,即误差不超过最后一位的半个单位,试指出它们是几位有效数字:1021.1*1=x ,031.0*2=x ,6.385*3=x ,430.56*4=x ,0.17*5⨯=x 。
[解]1021.1*1=x 有5位有效数字;0031.0*2=x 有2位有效数字;6.385*3=x 有4位有效数字;430.56*4=x 有5位有效数字;0.17*5⨯=x 有2位有效数字。
4、利用公式(3.3)求下列各近似值的误差限,其中*4*3*2*1,,,x x x x 均为第3题所给的数。
(1)*4*2*1x x x ++; [解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε;(2)*3*2*1x x x ;[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε;(3)*4*2/x x 。
计算机数值方法 第四章作业
![计算机数值方法 第四章作业](https://img.taocdn.com/s3/m/2974f7c14028915f804dc257.png)
n −1 b−a T8 = f ( a ) + 2∑ f ( xi ) + f (b) 2n i =1
1 = 0 + 2∑ f ( xi ) + 0.679570 ≈ 0.358726 16 i =1
7
习题4(P162):用复合梯形公式和复合 用复合梯形公式和复合simpson公式计算定积分 习题 公式计算定积分
若要该计算式具有2价代数精度,则必须有: 若要该计算式具有 价代数精度,则必须有: 价代数精度
1 2 2 2 (1+ 2α + 3β ) = 3 3
α β
= 1 − 5 6 6
(2) )
= 1 + 5 6 6
至此, )、(2)两式,可求出α 至此,由(1)、( )两式,可求出α、β: )、(
1
α β
解:
xe x = 。 I =∫ dx ,n=8。 0 (1 + x ) 2
1
i+1/2 1/2 3/2 5/2 7/2 复合 1/16 3/16 5/16 7/16 Simpson Xi+1/2 f (xi+1/2 ) 0.058934 0.160385 0.247952 0.327919 公式 i 9/2 11/2 13/2 15/2 xi 9/16 11/16 13/16 15/16 f (xi ) 0.404365 0.480134 0.557355 0.637733
I − Tn = − 8−2 2 1 1 1 h ( × 2η −3 ) ≤ 4 h 2 ≤ 10−5 12 2 2 2
8
1 计算ln2时,若采用复合 时 d x = 2 ln 2 计算 x
b−a 2 I − Tn = − h f "(η ) 12
数值计算课后全部答案(整合)
![数值计算课后全部答案(整合)](https://img.taocdn.com/s3/m/4f7cf8bb0029bd64783e2cc2.png)
目录第一章-----------------------------------------1 第二章-----------------------------------------4 第三章-----------------------------------------9 第四章-----------------------------------------15 第五章-----------------------------------------20 第六章-----------------------------------------27 第七章-----------------------------------------30第一章数值计算中的误差习题一1.1 下列各近似数的绝对误差限是最末位的半个单位,试指出它们各有几位有效数字。
1x =-3.105 , 2x =0.001, 3x =0.100, 4x =253.40, 5x =5000, 6x =5⨯310.答案:4,1,3,6,4,1.1.2 设100>*x >10,x 是*x 的有五位有效数字的的近似数,求x 的绝对误差限。
答案:当10<x<100时,因为有5位有效数字,所以绝对误差限为0.005. 1.3 求下列各近似数的相对误差限和有效数字位数: 1) 123x x x ++,2) 124x x x 3) 24x x 答案:()10.0005e x ≤()20.0005e x ≤()30.0005e x ≤ ()40.005e x ≤ ()50.5e x ≤ ()60.5e x ≤1)()()()()123123e x x x e x e x e x ++=++≤()()()123e x e x e x ++3221.5100.15100.510---≤⨯=⨯≤⨯2123()0.1510x x x ε-++=⨯123123123()()0.0004993...0.0004994r x x x e x x x x x x ε++++==≤++123x x x ++=-3.004 精确到小数点后两位,所以有三位有效数字。
《数值计算方法》复习资料
![《数值计算方法》复习资料](https://img.taocdn.com/s3/m/b62f2b54b0717fd5370cdcc2.png)
实用文档《数值计算方法》复习资料第一章数值计算方法与误差分析第二章非线性方程的数值解法第三章线性方程组的数值解法第四章插值与曲线拟合第五章数值积分与数值微分第六章常微分方程的数值解法自测题课程的性质与任务数值计算方法是一门应用性很强的基础课,在学习高等数学,线性代数和算法语言的基础上,通过本课程的学习及上机实习、使学生正确理解有关的基本概念和理论,掌握常用的基本数值方法,培养应用计算机从事科学与工程计算的能力,为以后的学习及应用打下良好基础。
第一章数值计算方法与误差分析一考核知识点误差的来源类型;绝对误差和绝对误差限,相对误差和相对误差限,有效数字;绝对误差的传播。
二复习要求1.知道产生误差的主要来源。
2.了解绝对误差和绝对误差限、相对误差和相对误差限和有效数字等概念以及它们之间的关系。
3.知道四则运算中的误差传播公式。
实用文档三例题例 1 设x*= =3.1415926⋯近似值 x=3.14 = 0.314× 101,即 m=1,它的绝对误差是- 0.001 592 6 ,⋯有即 n=3,故 x=3.14 有 3 位有效数字 .x=3.14准确到小数点后第 2 位 .又近似值 x=3.1416,它的绝对误差是0.0000074 ⋯,有即 m=1,n= 5, x=3.1416 有 5 位有效数字 .而近似值x=3.1415,它的绝对误差是0.0000926 ⋯,有即 m=1,n= 4, x=3.1415 有 4 位有效数字 .这就是说某数有s 位数,若末位数字是四舍五入得到的,那么该数有s 位有效数字;例 2指出下列各数具有几位有效数字,及其绝对误差限和相对误差限:2.000 4-0.002 009 0009 000.00解因为 x1=2.000 4= 0.200 04× 101, 它的绝对误差限 0.000 05=0.5 × 10 1―5,即m=1,n=5, 故 x=2.000 4 有 5 位有效数字 . a1=2,相对误差限x2=- 0.002 00,绝对误差限0.000 005,因为 m=-2,n=3 ,x2=- 0.002 00 有 3 位有效数字 . a1=2 ,相对误差限r ==0.002 5实用文档x3=9 000 ,绝对误差限为0.5× 100,因为 m=4, n=4, x3=9 000 有 4 位有效数字, a=9 ,相对误差限r== 0.000 056x4=9 000.00 ,绝对误差限0.005,因为 m=4, n=6, x4=9 000.00 有 6 位有效数字,相对误差限为r== 0.000 000 56由 x3与 x4可以看到小数点之后的0,不是可有可无的,它是有实际意义的.例 3 ln2=0.69314718⋯,精确到10-3的近似值是多少?解精确到 10-3= 0.001,意旨两个近似值x1,x2满足,由于近似值都是四舍五入得到的,要求满足,近似值的绝对误差限应是=0.0005,故至少要保留小数点后三位才可以。
数值计算方法第04章矩阵特征值与特征向量的计算
![数值计算方法第04章矩阵特征值与特征向量的计算](https://img.taocdn.com/s3/m/476e812743323968011c92ab.png)
• 计算出k=2时的x和y。 • (保留四位有效数字)
22
二、幂法的加速
因为幂法的收敛速度是线性的,而且依赖 于比值 2 /1 ,当比值接近于1时,幂法收敛 很慢。幂法加速有多种,介绍两种。
23
幂法的加速—原点移位法 应用幂法计算矩阵A的主特征值的收敛速度主要
26
4 14 0 , 2.9, 用原点移位法求矩 例:A 5 13 0 0 1 0 2.8 -4 阵A的按模最大的特征值,要求误差不超过10 。 解:取x (0) (1,1,1)T , 按x ( k 1) ( A pI )x (k )进行计算 0 6.9 14 A 0 I 5 10.1 0 0 0.1 1 (3.1000568, 2.214326, 0.9687661) 4 3.1000568
在一定条件下, 当k充分大时: 相应的特征向量为:
x 1 x
x
( k 1)
( k 1 ) i (k ) i
10
幂法的理论依据 对任意向量x(0), 有 x ( 0 ) i ui , 设1不为零.
i 1 n
x
( k 1 )
Ax
n i 1
(k )
A
k 1
x
(0) n
1 Ak 1 i ui i k i ui i 1
k 1 1
2 k 1 n k 1 1u1 ( ) 2 u2 ( ) n un 1 1
k 1 1 1u1
故 1 xi( k 1) xi( k ) x(k+1)为1的特征向量的近似向量(除一个因子外).
数值计算方法课后习题答案
![数值计算方法课后习题答案](https://img.taocdn.com/s3/m/51ba3849f5335a8102d220c1.png)
第一章 绪论(12)1、设0>x ,x 的相对误差为δ,求x ln 的误差。
[解]设0*>x 为x 的近似值,则有相对误差为δε=)(*x r ,绝对误差为**)(x x δε=,从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x , 相对误差为****ln ln )(ln )(ln x x x x rδεε==。
2、设x 的相对误差为2%,求n x 的相对误差。
[解]设*x 为x 的近似值,则有相对误差为%2)(*=x r ε,绝对误差为**%2)(x x =ε,从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε,相对误差为%2)()(ln )(ln ***n x x x nr==εε。
3、下列各数都是经过四舍五入得到的近似数,即误差不超过最后一位的半个单位,试指出它们是几位有效数字:1021.1*1=x ,031.0*2=x ,6.385*3=x ,430.56*4=x ,0.17*5⨯=x 。
[解]1021.1*1=x 有5位有效数字;0031.0*2=x 有2位有效数字;6.385*3=x 有4位有效数字;430.56*4=x 有5位有效数字;0.17*5⨯=x 有2位有效数字。
4、利用公式(3.3)求下列各近似值的误差限,其中*4*3*2*1,,,x x x x 均为第3题所给的数。
(1)*4*2*1x x x ++; [解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε;(2)*3*2*1x x x ;[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε;(3)*4*2/x x 。
第四章 自测题及答案
![第四章 自测题及答案](https://img.taocdn.com/s3/m/3b497bdd76a20029bd642dd2.png)
22. 已知 xi f ( xi ) -2 -1 0 1 2 4 2 1 3 5 求 f ( x ) 的二次拟合曲线 p 2 ( x) ,并求 f (0) 的近似值。 解: xi yi x i2 x i3 x i4 i 0 -2 4 4 -8 16 1 -1 2 1 -1 1 2 0 1 0 0 0 3 1 3 1 1 1 4 2 5 4 8 16 0 15 10 0 34
解:取 x0 0
l30
x1 2
x2 3
x3 5
l31 ( x x0 )( x x2 )( x x3 ) ( x1 x0 )( x1 x2 )( x1 x3 )
( x x1 )( x x2 )( x x3 ) ( x0 x1 )( x0 x2 )( x0 x3 ) ( x x0 )( x x1 )( x x3 ) ( x2 x0 )( x2 x1 )( x2 x3 )
当 n 4 时, ( xk4 xk2 3)l k ( x) _______________________________________.
k 0 n
答案: 1; x j ;
x4 x2 3
3. 设 li ( x ) 是以 xk k(k 0,1,, 9) 为节点的Lagrange插值基函数,则
1
湖北民族学院理学院《数值计算方法》教学辅导材料
陈以平编写
x3 0 x 1 9. 已知 S ( x) 是三次样条函数,则 1 ( x 1) 3 a ( x 1) 2 b( x 1) c 1 x 3 2 =______ ; b =________; c =__________. a 答案: 3;3; 1
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
f ( xk ) e xk − 4 cos ( xk ) = xk − xk xk +1 = xk − f ′ ( xk ) e + 4sin ( xk )
利用公式作迭代得:
e − 4 cos ( x0 ) x1 = x0 − x0 = 0.911878, e + 4sin ( x0 )
x0
e x1 − 4 cos ( x1 ) x2 = x1 − x1 = 0.904810, e + 4sin ( x1 ) e − 4 cos ( x2 ) x3 = x2 − x2 = 0.904788 e + 4sin ( x2 )
因此可以得到牛顿迭代式为:
2 2 axk − 1 axk + 1 xk +1 = xk − = 2axk 2axk 1 −2 解法2:f ( x ) = 2 − a, f ′ ( x ) = 3 , x x 1 −2 相应的:f ( xk ) = 2 − a, f ′ ( x ) = 3 , xk xk
2 x13 + 1 x2 = 2 = 1.8795 3x1 − 3
3 2 x2 + 1 x3 = 2 = 1.8794, x∗ = 1.879 3x2 − 3
f ( x ) = ax 2 − 1, f ′ ( x ) = 2ax 21、解法1:由题意知 f ( xk ) = axk2 − 1, f ′ ( xk ) = 2axk 相应的:
数值计算方法第四章参考答案
7.、 证明:1) 对于该迭代式,相应的迭代函数为:
1 ϕ ( x ) = 1 + 2 ,易知该函数为单调递减函数, x 而 ϕ (1.3) = 1.59, ϕ (1.3) = 1.39
因此有 ϕ
( x ) ∈ [1.3,1.6] 2 ′( x) = − ϕ 又: ,易知 ϕ ′ ( x ) 为单调递减函数,所以有 x
因此可以得到牛顿迭代式为:
xk +1 = xk −
1 − axk2 −2
2 xk
3 xk
a 3 3 = xk − xk 2 2
x x 23. 解: 1)由题意得 f ( xk ) = e k − 4 cos ( xk ) , f ′ ( xk ) = e k + 4sin ( xk )
相应的牛顿迭代公式为:
x2
2) 弦截法:取 x0
=
π
4
, x1 =
π
2
相应的弦截迭代公式为:
xk +1 = xk −
f ( xk ) − f ( x0 )
f ( xk )
( xk − x0 )
xk e − 4 cos ( xk ) = xk − ( xk − x0 ) π e xk − 4 cos ( xk ) − (e 4 − 4 ⋅ 2 ) 2
(
)
利用公式作迭代得:
x1 e − 4 cos ( x1 ) x2 = x1 − ( x1 − x0 ) = 0.877003 π e x1 − 4 cos ( x1 ) − (e 4 − 4 ⋅ 2 ) 2
x2 e − 4 cos ( x2 ) x3 = x2 − ( x2 − x0 ) = 0.906360 π e x2 − 4 cos ( x2 ) − (e 4 − 4 ⋅ 2 ) 2 x3 e − 4 cos ( x3 ) x4 = x3 − ( x3 − x0 ) = 0.904701 π e x3 − 4 cos ( x3 ) − (e 4 − 4 ⋅ 2 ) 2
相应的牛顿迭代公式为:
3 3 f ( xk ) xk − 3 xk − 1 2 xk + 1 = xk − = 2 xk +1 = xk − 2 3 xk − 3 3xk − 3 f ′ ( xk )
利用公式作迭代得:
2 x + 1 17 x1 = = = 1.8889 3x − 3 9
3 0 2 0
由压缩影像定理知该迭代式收敛。
2 − 2 又: ′ ( x ) = x(1 + x) 3 ϕ 3
2 2 − ϕ ′ ( x ) ≤ ⋅1.6 ⋅ (1 + 1.6 ) 3 = 0.5641 < 1 3
2 3 17、 解:1) 由题意得 f ′ ( xk ) = 3xk − 3, f ( xk ) = xk − 3 xk − 1
(
利用公式作迭代得:
x1 e − 4 cos ( x1 ) x2 = x1 − ( x1 − x0 ) = 0.877003 π e x1 − 4 cos ( x1 ) − (e 4 − 4 ⋅ 2 ) 2 e x2 − 4 cos ( x2 ) ( x2 − x1 ) = 0.898545 x3 = x2 − x x1 e 2 − 4 cos ( x2 ) − ( e − 4 cos ( x1 ) ) e x3 − 4 cos ( x3 ) ( x3 − x2 ) = 0.904866 x4 = x3 − x x2 e 3 − 4 cos ( x3 ) − ( e − 4 cos ( x2 ) )
3
ϕ′( x) ≤
2 = 0.9103 < 1 3 (1.3)
由压缩影像定理知该迭代式收敛。
2) 对于该迭代式,相应的迭代函数为:
ϕ ( x) = 1+ x
3
2
,易知该函数为单调递增函数,
而 ϕ (1.3) = 1.39, ϕ (1.3) = 1.53
因此有 ϕ
( x ) ∈ [1.3,1.6] 易知 ϕ ′ ( x ) 为单调递增函数,所以有
3) 快速弦截法 相应的弦截迭代公式为:
f ( xk ) xk +1 = xk − ( xk − xk −1 ) f ( xk ) − f ( xk −1 ) e xk − 4 cos ( xk ) = xk − x e k − 4 cos ( xk ) − e xk −1 − 4 cos ( xk −1 )