当前位置:文档之家› 数值分析试卷及其答案1

数值分析试卷及其答案1

1. 已知325413.0,325413*2*1==X X 都有6位有效数字,求绝对误差限。

(4分) 解:由已知可知65.01021,0,6,10325413.0016*1=⨯==-=⨯=ε绝对误差限n k k X2分620*21021,6,0,10325413.0-⨯=-=-=⨯=ε绝对误差限n k k X2分2. 已知⎢⎢⎢⎣⎡=001A 220- ⎥⎥⎥⎦⎤440求21,,A A A ∞ (6分)解:{},88,4,1max 1==A1分{},66,6,1max ==∞A1分()AA A T max 2λ=1分⎢⎢⎢⎣⎡=01A A T42⎥⎥⎥⎦⎤-420⎢⎢⎢⎣⎡001220-⎥⎥⎥⎦⎤440=⎢⎢⎢⎣⎡00180⎥⎥⎥⎦⎤3200 2分{}3232,8,1max )(max ==A A T λ1分24322==A3. 设32)()(a x x f -= (6分) ① 写出f(x)=0解的迭代格式② 当a 为何值时,)(1k k x x ϕ=+ (0,1……)产生的序列{}k x 收敛于2解:①迭代格式为:xa x x x ax a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(22321+=+=---=-=+ϕ 3分 ②时迭代收敛即当222,11210)2(',665)('2<<-<-=-=a ax a x ϕϕ 3分4. 给定线性方程组,其中:⎢⎣⎡=13A ⎥⎦⎤22,⎥⎦⎤⎢⎣⎡-=13b 用迭代公式)()()()1(k k k Ax b x x -+=+α(0,1……)求解,问取什么实数α,可使迭代收敛 (8分) 解:所给迭代公式的迭代矩阵为⎥⎦⎤--⎢⎣⎡--=-=ααααα21231A I B 2分其特征方程为0)21(2)31(=----=-αλαααλλB I2分 即,解得αλαλ41,121-=-=2分要使其满足题意,须使1)(<B ρ,当且仅当5.00<<α 2分5. 设方程,其中⎢⎢⎢⎣⎡=211A 212 ⎥⎥⎥⎦⎤-112,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=765b 试讨论解此方程的迭代法的收敛性,并建立迭代格式 (9分) 解:U D L A ++=⎢⎢⎢⎣⎡--=+-=-21)(1U L D B J22--⎥⎥⎥⎦⎤-012 3分,03213=====-λλλλλJ B I2分即10)(<=J B ρ,由此可知迭代收敛1分迭代格式:⎪⎩⎪⎨⎧--=--=+-=++++++)1(2)1(1)1(3)(3)1(1)1(2)(3)(2)1(12276225k k k k k k k k k x x x x x x x x x (0,1,2,3……)3分6. 用分解计算下列3个线性代数方程组:i ib Ax =(1,2,3)其中⎢⎢⎢⎣⎡=222A 331 ⎥⎥⎥⎦⎤421,23121,,974x b x b b ==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡= (12分)解:①11b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=9741x ⎢⎢⎢⎣⎡11111⎥⎥⎥⎦⎤100⎢⎢⎢⎣⎡00221⎥⎥⎥⎦⎤211 3分 由1,即⎢⎢⎢⎣⎡111110 ⎥⎥⎥⎦⎤100⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡974 得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡234 1分 由1,即⎢⎢⎢⎣⎡00221⎥⎥⎥⎦⎤211x1=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡234 得x1=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 2分 ②22b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111由21,即⎢⎢⎢⎣⎡11111⎥⎥⎥⎦⎤100⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡001 1分 由2,即⎢⎢⎢⎣⎡00221⎥⎥⎥⎦⎤211x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡001 得x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0 2分 ③33b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0 由32,即⎢⎢⎢⎣⎡111110 ⎥⎥⎥⎦⎤100⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0 得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-05.05.0 1分由3,即⎢⎢⎢⎣⎡00221⎥⎥⎥⎦⎤211x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-05.05.0 得x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-025.0375.0 2分7. 已知函数(x)有关数据如下:要求一次数不超过3的H 插值多项式,使'11'33)(,)(y x H y x H i i ==(6分) 解:作重点的差分表,如下:3分21021101011001003))(](,,,[))(](,,[)](,[][)(x x x x x x x x f x x x x x x x f x x x x f x f x H --+--+-+= 1+(1)(1)+2(1) =232x x +3分8. 有如下函数表:试计算此列表函数的差分表,并利用前插公式给出它的插值多项式 (7分) 解:由已知条件可作差分表,3分 i ih x x i=+=0 (0,1,2,3)为等距插值节点,则向前插值公式为: 033210022100003!3))()((!2))((!1)()(f h x x x x x x f h x x x x f h x x f x N ∆---+∆--+∆-+==4+5(1) =442++x x4分9. 求f(x)在[-1,1]上的二次最佳平方逼近多项式)(2x P ,并求出平方误差 (8分) 解: 令22102)(x a x a a x P ++=2分取1, , 2x ,计算得:()dx ⎰-1110 ()= dx x ⎰-11=1 ()= dx x ⎰-112=0 ()= dx x⎰-113=0.5 ()= dx x ⎰-114=0 ()= dx x ⎰-11=1()= dx x ⎰-112=0 ()= dx x ⎰-113=0.5 得方程组:⎪⎩⎪⎨⎧==+=5.05.005.011201a a a a 3分 解之得c a a c a 2,1,210-=== (c为任意实数,且不为零)即二次最佳平方逼近多项式222)(cx x c x P -+=1分 平方误差:32),(22222222=-=-=∑=i i i y a fp f ϕδ2分10. 已知如下数据:用复合梯形公式,复合公式计算⎰+=10214dx x π的近似值(保留小数点后三位) (8分)解:用复合梯形公式: )}1()]87()43()85()21()83()41()81([2)0({1618f f f f f f f f f T ++++++++==3.139 4分用复合公式: )}1()]43()21()41([2)]87()85()83()81([4)0({2414f f f f f f f f f S ++++++++==3.142 4分11. 计算积分⎰=2sin πxdx I ,若用复合公式要使误差不超过51021-⨯,问区间]2,0[π要分为多少等分?若改用复合梯形公式达到同样精确度,区间]2,0[π应分为多少等分? (10分)解: ①由公式余项及x x f x x f sin )(,sin )()4(==得544)4(2041021)1()4(360)(max )4(1802)(-≤≤⨯≤=≤n x f n f R x n πππππ2分 即08.5,6654≥≥n n ,取62分即区间]2,0[π分为12等分可使误差不超过51021-⨯1分②对梯形公式同样1)(''max20≤≤≤x f x π,由余项公式得51021)2(122)(-⨯≤≤n f R n ππ2分即255,2.254=≥n n 取 2分 即区间]2,0[π分为510等分可使误差不超过51021-⨯ 1分 12.用改进格式求解初值问题:⎩⎨⎧==++1)1(0sin 2'y x y y y 要求取步长h为0.1,计算y (1.1)的近似值 (保留小数点后三位)[提示:1=0.841.1=0.89] (6分) 解:改进格式为:⎪⎩⎪⎨⎧++=+=+-++-+)],(),([2),(1111n n n n n n n n n n y x f y x f hy y y x hf y y2分于是有⎪⎩⎪⎨⎧+++-=+-=+-++-+-+)sin sin (05.0)sin (1.012112121n n n n n n n n n n n n n x y y x y y y y x y y y y (0,1,2……) 2分由y(1)0y 1,计算得⎪⎩⎪⎨⎧=≈=+-=-838.0)1.1(816.0)1sin 11(1.01121y y y2分即y(1.1)的近似值为0.838 13.][],[],,[lim ],[),,(],,[)(0'000000'0x f x x f x x f x x f b a x b a C x f x x ==∈∈→证明:定义:设(4分) 证明: ]['],[],[],[lim ][][lim]['00000000000x f x x f x x f x x f x x x f x f x f x x x x ===--=→→故可证出4分14. 证明:设n n R A ⨯∈,⋅为任意矩阵范数,则A A ≤)(ρ (6分) 证明:设λ为A 的按模最大特征值,x 为相对应的特征向量,则有λ 1分且λρ=)(A ,若λ是实数,则x 也是实数,得Ax x =λ 1分而x x ⋅=λλ x A x ,⋅≤⋅⋅≤λ故x A Ax 2分由于A x 0x ≤≠λ得到,两边除以 1分故A A ≤)(ρ 1分当λ是复数时,一般来说x 也是复数,上述结论依旧成立。

相关主题