数值分析试卷及其答案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 也是复数，上述结论依旧成立。