太原理工大学数值计算方法题库
- 格式:doc
- 大小:355.00 KB
- 文档页数:16
计算机数值计算方法试题 计算机数值计算方法试题一一、 填空题(每空1分,共17分)1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数,需对分( )次。
2、迭代格式)2(21-+=+k k k x x x α局部收敛的充分条件是α取值在( )。
3、已知⎪⎩⎪⎨⎧≤≤+-+-+-≤≤=31)1()1()1(2110)(233x c x b x a x x x x S 是三次样条函数,则a =( ),b =( ),c =( )。
4、)(,),(),(10x l x l x l n 是以整数点n x x x ,,,10 为节点的Lagrange 插值基函数,则∑==nk kx l0)(( ),∑==nk k jk x lx 0)(( ),当2≥n 时=++∑=)()3(24x l x xk k n k k( )。
,5、设1326)(247+++=x x x x f 和节点,,2,1,0,2/ ==k k x k 则=],,,[10n x x x f和=∆07f。
6、5个节点的牛顿-柯特斯求积公式的代数精度为 ,5个节点的求积公式最高代数精度为 。
7、{}∞=0)(k k x ϕ是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族,其中1)(0=x ϕ,则⎰=14)(dx x x ϕ 。
8、给定方程组⎩⎨⎧=+-=-221121b x ax b ax x ,a 为实数,当a 满足 ,且20<<ω时,SOR 迭代法收敛。
9、解初值问题00(,)()y f x y y x y '=⎧⎨=⎩的改进欧拉法⎪⎩⎪⎨⎧++=+=++++)],(),([2),(]0[111]0[1n n n n n n n n n n y x f y x f h y y y x hf y y 是阶方法。
10、设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=11001a a a a A ,当∈a ( )时,必有分解式T LL A =,其中L 为下三角阵,当其对角线元素)3,2,1(=i l ii 满足( )条件时,这种分解是唯一的。
11级(12/07/03)
一、基础题(40分)
(一)、单项选择(2×5=10分)
1、求解常微分方程的预估—校正法的局部截断误差为( )。
2、过
3
4
5、
(二)
1、是一日插值基函数在节点上的取值是______________。
2、设分段多项式,
,
是以0,2,3为节点的三次样条
函数。
则a =____________,b =____________, c =____________。
3、设,则关于节点,,的二阶向前差分为_________。
4、5个节点的牛顿—科特斯求积公式的代数精度为________,5个节点的求积公式最高代数精度为________。
5、设,则a的取值范围为________A可分解为A = LL T,且当L满足________,分解是唯一的。
6、设是切比雪夫正交多项式系,则的正交区间为________,它的权
7迭8。
1
2
(1)
(2)
3、用二步法求解一阶常微分方程初值问题
,
,问:如何选择、的值,才能使该方法的阶数尽可能高?写出此时的局部截断误差主项。
三、计算题(15×2=30分)
1、(1)设,,是区间[-1,1]上权函数为的最高项系数为1的正交多项式组,其中,,求。
(2)利用,,求函数在[-1,1]上的二次最佳平方逼近多项式。
2、已知求解方程组Ax = b的分量迭代格式:,
,,,;,,,
(1)试求出矩阵格式及迭代矩阵。
(2)证明当A为严格对角占优矩阵,时,该迭代格式收敛。
《数值计算⽅法》试题集及答案要点《数值计算⽅法》复习试题⼀、填空题:1、----=410141014A ,则A 的LU 分解为A ?=。
答案:--??--=15561415014115401411A 2、已知3.1)3(,2.1)2(,0.1)1(===f f f ,则⽤⾟普⽣(⾟⼘⽣)公式计算求得?≈31_________)(dx x f ,⽤三点式求得≈')1(f 。
答案:2.367,0.25 3、1)3(,2)2(,1)1(==-=f f f ,则过这三点的⼆次插值多项式中2x 的系数为,拉格朗⽇插值多项式为。
答案:-1,)2)(1(21)3)(1(2)3)(2(21)(2--------=x x x x x x x L4、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字;5、设)(x f 可微,求⽅程)(x f x =的⽜顿迭代格式是();答案)(1)(1n n n n n x f x f x x x '---=+6、对1)(3++=x x x f ,差商=]3,2,1,0[f (1 ),=]4,3,2,1,0[f ( 0 );7、计算⽅法主要研究( 截断 )误差和( 舍⼊ )误差; 8、⽤⼆分法求⾮线性⽅程f (x )=0在区间(a ,b )内的根时,⼆分n 次后的误差限为(12+-n a b );9、求解⼀阶常微分⽅程初值问题y '= f (x ,y ),y (x 0)=y 0的改进的欧拉公式为()],(),([2111+++++=n n n n n n y x f y x f hy y);10、已知f (1)=2,f (2)=3,f (4)=5.9,则⼆次Newton 插值多项式中x 2系数为( 0.15 );11、两点式⾼斯型求积公式?10d )(x x f ≈(?++-≈1)]3213()3213([21d )(f f x x f),代数精度为( 5 );12、解线性⽅程组A x =b 的⾼斯顺序消元法满⾜的充要条件为(A 的各阶顺序主⼦式均不为零)。
数值计算方法试题一一、 填空题(每空1分,共17分) 1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数,需对分( )次。
2、迭代格式)2(21-+=+k k k x x x α局部收敛的充分条件就是α取值在( )。
3、已知⎪⎩⎪⎨⎧≤≤+-+-+-≤≤=31)1()1()1(2110)(233x c x b x a x x x x S 就是三次样条函数,则a =( ),b =( ),c =( )。
4、)(,),(),(10x l x l x l n Λ就是以整数点n x x x ,,,10Λ为节点的Lagrange 插值基函数,则∑==nk kx l0)(( ),∑==nk k jk x lx 0)(( ),当2≥n 时=++∑=)()3(204x l x xk k n k k( )。
5、设1326)(247+++=x x x x f 与节点,,2,1,0,2/Λ==k k x k 则=],,,[10n x x x f Λ 与=∆07f。
6、5个节点的牛顿-柯特斯求积公式的代数精度为 ,5个节点的求积公式最高代数精度为 。
7、{}∞=0)(k kx ϕ就是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族,其中1)(0=x ϕ,则⎰=14)(dx x x ϕ 。
8、给定方程组⎩⎨⎧=+-=-221121b x ax b ax x ,a 为实数,当a 满足 ,且20<<ω时,SOR 迭代法收敛。
9、解初值问题00(,)()y f x y y x y '=⎧⎨=⎩的改进欧拉法⎪⎩⎪⎨⎧++=+=++++)],(),([2),(]0[111]0[1n n n n n n n n n n y x f y x f h y y y x hf y y 就是阶方法。
10、设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=11001a a a a A ,当∈a ( )时,必有分解式T LL A =,其中L为下三角阵,当其对角线元素)3,2,1(=i l ii 满足( )条件时,这种分解就是唯一的。
数值计算方法试题一一、 填空题(每空1分,共17分)1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数,需对分( )次。
2、迭代格式)2(21-+=+k k k x x x α局部收敛的充分条件是α取值在( )。
3、已知⎪⎩⎪⎨⎧≤≤+-+-+-≤≤=31)1()1()1(2110)(233x c x b x a x x x x S 是三次样条函数,则a =( ),b =( ),c =( )。
4、)(,),(),(10x l x l x l n 是以整数点n x x x ,,,10 为节点的Lagrange 插值基函数,则∑==nk kx l0)(( ),∑==nk k jk x lx 0)((),当2≥n 时=++∑=)()3(204x l x xk k nk k ( )。
5、设1326)(247+++=x x x x f 和节点,,2,1,0,2/ ==k k x k 则=],,,[10n x x x f 和=∆07f。
6、5个节点的牛顿-柯特斯求积公式的代数精度为 ,5个节点的求积公式最高代数精度为 。
7、{}∞=0)(k kx ϕ是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族,其中1)(0=x ϕ,则⎰=14)(dx x x ϕ 。
8、给定方程组⎩⎨⎧=+-=-221121b x ax b ax x ,a 为实数,当a 满足 ,且20<<ω时,SOR 迭代法收敛。
9、解初值问题00(,)()y f x y y x y '=⎧⎨=⎩的改进欧拉法⎪⎩⎪⎨⎧++=+=++++)],(),([2),(]0[111]0[1n n n n n n n n n n y x f y x f h y y y x hf y y 是阶方法。
10、设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=11001a a a a A ,当∈a ( )时,必有分解式T LL A =,其中L 为下三角阵,当其对角线元素)3,2,1(=i l ii 满足( )条件时,这种分解是唯一的。
《数值计算⽅法》试题集及标准答案(-)-《数值计算⽅法》试题集及答案(-)-————————————————————————————————作者:————————————————————————————————⽇期:《计算⽅法》期中复习试题⼀、填空题:1、已知3.1)3(,2.1)2(,0.1)1(===f f f ,则⽤⾟普⽣(⾟⼘⽣)公式计算求得≈31_________)(dx x f ,⽤三点式求得≈')1(f 。
答案:2.367,0.252、1)3(,2)2(,1)1(==-=f f f ,则过这三点的⼆次插值多项式中2x 的系数为,拉格朗⽇插值多项式为。
答案:-1,)2)(1(21)3)(1(2)3)(2(21)(2--------=x x x x x x x L3、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字;4、设)(x f 可微,求⽅程)(x f x =的⽜顿迭代格式是( );答案)(1)(1n n n n n x f x f x x x '---=+5、对1)(3++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 );6、计算⽅法主要研究( 截断 )误差和( 舍⼊ )误差;7、⽤⼆分法求⾮线性⽅程 f (x )=0在区间(a ,b )内的根时,⼆分n 次后的误差限为( 12+-n a b );8、已知f (1)=2,f (2)=3,f (4)=5.9,则⼆次Newton 插值多项式中x 2系数为( 0.15 ); 11、两点式⾼斯型求积公式?1d )(xx f ≈(++-≈1)]3213()3213([21d )(f f x x f ),代数精度为( 5 );12、为了使计算32)1(6)1(41310---+-+=x x x y 的乘除法次数尽量地少,应将该表达式改写为11,))64(3(10-=-++=x t t t t y ,为了减少舍⼊误差,应将表达式19992001-改写为199920012+ 。
数值计算方法试题一一、 填空题1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数,需对分(10)次。
2、迭代格式)2(21-+=+k k k x x x α局部收敛的充分条件是α取值在()0,22(-)22,0()。
3、已知⎪⎩⎪⎨⎧≤≤+-+-+-≤≤=31)1()1()1(2110)(233x c x b x a x x x x S 是三次样条函数,则a =(3),b =(3),c =(1)。
4、)(,),(),(10x l x l x l n Λ是以整数点n x x x ,,,10Λ为节点的Lagrange 插值基函数,则∑==nk kx l)((1),∑==nk k jk x lx 0)((j x ),当2≥n 时=++∑=)()3(204x l x xk k nk k (324++x x )。
5、设1326)(247+++=x x x x f 和节点,,2,1,0,2/Λ==k k x k 则=],,,[10n x x x f Λ 6和=∆07f 25.236494526!77==⨯。
6、5个节点的牛顿-柯特斯求积公式的代数精度为9,5个节点的求积公式最高代数精度为9。
7、{}∞=0)(k kx ϕ是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族,其中1)(0=x ϕ,则⎰=14)(dx x x ϕ0。
8、给定方程组⎩⎨⎧=+-=-221121b x ax b ax x ,a 为实数,当a 满足1<a ,且20<<ω时,SOR 迭代法收敛。
9、解初值问题00(,)()y f x y y x y '=⎧⎨=⎩的改进欧拉法⎪⎩⎪⎨⎧++=+=++++)],(),([2),(]0[111]0[1n n n n n n n n n n y x f y x f h y y y x hf y y 是2阶方法。
10、设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=111a aa a A ,当∈a (22,22-)时,必有分解式T LL A =,其中L 为下三角阵,当其对角线元素)3,2,1(=i l ii 满足(0>ii l )条件时,这种分解是唯一的。
数值计算方法试题一一、填空题(每空1分,共17分)1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数,需对分( )次。
2、迭代格式)2(21-+=+k k k x x x α局部收敛的充分条件是α取值在( )。
3、已知⎪⎩⎪⎨⎧≤≤+-+-+-≤≤=31)1()1()1(2110)(233x c x b x a x x x x S 是三次样条函数,则a =( ),b =( ),c =( )。
4、)(,),(),(10x l x l x l n 是以整数点n x x x ,,,10 为节点的Lagrange 插值基函数,则∑==nk kx l)(( ),∑==nk k jk x lx 0)(( ),当2≥n 时=++∑=)()3(204x l x xk k nk k ( )。
5、设1326)(247+++=x x x x f 和节点,,2,1,0,2/ ==k k x k 则=],,,[10n x x x f 和=∆07f 。
6、5个节点的牛顿-柯特斯求积公式的代数精度为 ,5个节点的求积公式最高代数精度为 。
7、{}∞=0)(k kx ϕ是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族,其中1)(0=x ϕ,则⎰=104)(dx x x ϕ 。
8、给定方程组⎩⎨⎧=+-=-221121b x ax b ax x ,a 为实数,当a 满足 ,且20<<ω时,SOR迭代法收敛。
9、解初值问题00(,)()y f x y y x y '=⎧⎨=⎩的改进欧拉法⎪⎩⎪⎨⎧++=+=++++)],(),([2),(]0[111]0[1n n n n n n n n n n y x f y x f h y y y x hf y y 是阶方法。
10、设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=11001a a a a A ,当∈a ( )时,必有分解式T LL A =,其中L 为下三角阵,当其对角线元素)3,2,1(=i l ii 满足( )条件时,这种分解是唯一的。
枯藤老树昏鸦,小桥流水人家,古道西风瘦马。
夕阳西下,断肠人在天涯。
《数值计算方法》复习试题一、填空题:1、⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=410141014A ,则A 的LU 分解为A ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦。
答案:⎢⎢⎢⎣⎡--=15401411A 3、)3(,2)2(,1)1(=-=f f f 答案:-1,2(21)(2-=x x L 4* 2 )位有效数字; 5( );答案6 ),=]4,3,2,1,0[f ( 0 );7舍入 )误差;8a ,b )内的根时,二分n 次后的误差限为( 12+n );10、已知f (1)=2,f (2)=3,f (4)=5.9,则二次Newton 插值多项式中x 2系数为( 0.15 ); 11、 解线性方程组A x =b 的高斯顺序消元法满足的充要条件为(A 的各阶顺序主子式均不为零)。
12、 为了使计算32)1(6)1(41310---+-+=x x x y 的乘除法次数尽量地少,应将该表达式改写为11,))64(3(10-=-++=x t t t t y ,为了减少舍入误差,应将表达式19992001-改写为199920012+ 。
13、 用二分法求方程01)(3=-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间为 0.5,1 ,进行两步后根的所在区间为 0.5,0.75 。
14、 求解方程组⎩⎨⎧=+=+042.01532121x x x x,该迭代格式的迭代矩阵的谱半径)(M ρ15、 设46)2(,16)1(,0)0(===f f f ,则1l 插值多项式为 (716)(2+=x x x x N 16、 求积公式⎰∑=≈ba k nk k x f A x x f )(d )(0的代数精度以( 高斯型 )求积公式为最高,具21内的根精确到三位小数,需对分( 10 )次。
22≤≤≤≤3110x c x 是三次样条函数,则a =( )。
上机练习题一、求非线性方程的根。
1、 求方程()cos 0f x x x =-=在0 1.5x =附近的是根,要求精度满足3110k k x x -+-<.(牛顿切线法)>> NewtonIterationx0=1.5 del=0.001 N=20k x(k)0 1.500000 1 0.784472结果:0.7395192、 求方程32()0.80f x x x =--=在01x =附近的是根,求出具有四位有效数字的根近似值..(简单迭代法))(1n n x x ϕ=+312)8.0()(+=x x ϕ程序clear clcphi=inline('(x^2+0.8)^(1/3)'); %迭代函数 x0=input('x0='); del=input('del='); N=input('N='); n=1;fprintf('\n %2d %f ',0,x0); while n<N x=phi(x0); if abs(x-x0)<delfprintf('\n \n 近似解=%f \n',x); return endfprintf('\n %2d %f ',n,x); n=n+1; x0=x; endfprintf('\n \n %f d 次迭代后未达到精度要求. \n',N);运行结果 x0=1del=0.00001N=200 1.000000 1 1.216440 2 1.316116 3 1.363004 4 1.385180 5 1.395688 6 1.400671 7 1.403034 8 1.404155 9 1.404687 10 1.404939 11 1.405059 12 1.405116 13 1.405143 14 1.405155近似解=1.405162二、求解线性方程组(直接法或迭代法)1、⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-----22118118344108318311231224321x x x x x使用高斯-赛德尔迭代法求解代码clear clcn=input('n=');%矩阵的阶数 A=input('A=');%系数矩阵 b=input('b='); x=input('x=');%自变量 epsilon=input('\n 精度='); N=input('\n 最大迭代次数N=');fprintf('\n %d:',0);for i=1:nfprintf('%f',x(i));end%以下是迭代过程for k=1:N%这是第k步迭代,迭代前的向量在x0[]中,迭代后的向量在x[]中; normal=0;for i=1:nt=x(i);x(i)=b(i);for j=1:nif j~=ix(i)=x(i)-A(i,j)*x(j);endendx(i)=x(i)/A(i,i);temp=abs(x(i)-t);% 求范数于迭代在同一个循环中;if temp>normalnormal=temp; %这里用的是无穷范数endend%第i不迭代结束;fprintf('\n %d: ',k);for i=1:nfprintf('%f',x(i));%输出迭代过程endif normal<epsilonreturnendendfprintf('\n \n 迭代% d 次后仍未求得满足精度的解\n',N);结果n=4A=[2,2,1,-3;-2,1,-1,-3;8,-1,3,8;10,4,4,3]b=[8,1,-1,8]x=[1,-1,2,-2]精度=0.001最大迭代次数N=100:1.000000-1.0000002.000000-2.0000001: 1.000000-1.0000002.000000-2.000000>>故原方程的解为⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--=2211x。
数值计算方法考试试题一、选择题(每小题4分,共20分)1. 误差根据来源可以分为四类,分别是( A )A. 模型误差、观测误差、方法误差、舍入误差;B. 模型误差、测量误差、方法误差、截断误差;C. 模型误差、实验误差、方法误差、截断误差;D. 模型误差、建模误差、截断误差、舍入误差。
2. 若132)(356++-=x x x x f ,则其六阶差商=]3,,3,3,3[6210 f ( C ) A. 0; B. 1; C. 2; D. 3 。
3. 数值求积公式中的Simpson 公式的代数精度为 ( D )A. 0;B. 1;C. 2;D. 3 。
4. 若线性方程组Ax = b 的系数矩阵A 为严格对角占优矩阵,则解方程组的Jacobi 迭代法和Gauss-Seidel 迭代法 ( B )A. 都发散;B. 都收敛C. Jacobi 迭代法收敛,Gauss-Seidel 迭代法发散;D. Jacobi 迭代法发散,Gauss-Seidel 迭代法收敛。
5. 对于试验方程y y λ=',Euler 方法的绝对稳定区间为( C )A. 02≤≤-h ;B. 0785.2≤≤-h ;C. 02≤≤-h λ;D. 0785.2≤≤-h λ ; 二、填空题(每空3分,共18分)1. 已知⎪⎪⎭⎫⎝⎛--='-=4321,)2,1(A x ,则 =2x 5,=1Ax 16 ,=2A 22115+2. 已知3)9(,2)4(==f f ,则 f (x )的线性插值多项式为)6(2.0)(1+=x x L ,且用线性插值可得f (7)= 2.6 。
3. 要使20的近似值的相对误差界小于0.1%,应至少取 4 位有效数字。
三、利用下面数据表,1. 用复化梯形公式计算积分dxx f I )(6.28.1⎰=的近似值;解:1.用复化梯形公式计算 取2.048.16.2,4=-==h n 1分分分分7058337.55))6.2()2.08.1(2)8.1((22.04))()(2)((231114=+++=++=∑∑=-=f k f f b f x f a f hT k n k k10.466758.030146.042414.425693.12014f (x ) 2.6 2.4 2.2 2.0 1.8 x2. 用复化Simpson 公式计算积分dxx f I )(6.28.1⎰=的近似值。
枯藤老树昏鸦,小桥流水人家,古道西风瘦马。
夕阳西下,断肠人在天涯《数值计算方法》复习试题一、填空题:3得,f(x)dx --------------------- ,用三点式求得f (1)答案:2.367,0.25 3、f (1)1, f (2)2, f(3) J 则过这三点的二次插值多项式中x 2的系数为 ______ ,拉格朗日插值多项式为 __________________________4、 近似值x * 0.231关于真值x 0.229有(2 )位有效数字;5、设f (x )可微,求方程x f (x)的牛顿迭代格式是();X n 1 x n xn f(xn) 答案 1 f(x n )6、对 f(x) x 3 x j 差商 f[0,1,2,3]( 1 ), f[0,1,2,3,4]( 0 );7、 计算方法主要研究(截断)误差和(舍入)误差; 8用二分法求非线性方程f (x)=0在区间(a,b)内的根时,二分n 次后b a的误差限为(莎 );4 10 AA1 411、 0 14则A 的LU 分解为1410 A1 4 1 15 41 0答案:4 15 156 15 2、已知 f (1)1.0,f ⑵ 1.2, f (3) 仁,则用辛普生(辛卜生)公式计算求答案:-1,L 2(X ) ;(x 2)(x 3)2(x 1)(x 3) £(x 1)(x 2)用二分法求方程f(x) x 3 x 1 0在区间[0,1]内的根,进行一步后根的所在区间为 0.5,1,进行两步后根的所在区间为0.75 o1I —15、计算积分0L xdx ,取4位有效数字。
用梯形公式计算求得的近似值为 0.4268 ,用辛卜生公式计算求得的近似值为 0.4309,梯 形公式的代数精度为1,辛卜生公式的代数精度为3 。
3x 1 5x 2116、求解方程组0.2x1 4x2 0的高斯一塞德尔迭代格式为—x 1(k 1} (1 5x 2k))/31(k 1)(k 1)x2x1/20_,该迭代格式的迭代矩阵的谱半径(M)= —12 — o9、求解一阶常微分方程初值问题y= f (x,y), y(x o )=y o 的改进的欧拉公式为(h yn1yn2[f(xn 'yn)f(xn1'yn1)]10、已知f(1) = 2, f(2) = 3, f ⑷=5.9,则二次Newton 插值多项式中x 2系数为(0.15 );11、1113131“ 两点式高斯型求积公式。
数值计算方法试题一一、 填空题1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数,需对分(10)次。
2、迭代格式)2(21-+=+k k k x x x α局部收敛的充分条件是α取值在()0,22(-)22,0()。
3、已知⎪⎩⎪⎨⎧≤≤+-+-+-≤≤=31)1()1()1(2110)(233x c x b x a x x x x S 是三次样条函数,则a =(3),b =(3),c =(1)。
4、)(,),(),(10x l x l x l n 是以整数点n x x x ,,,10 为节点的Lagrange 插值基函数,则∑==nk kx l)((1),∑==nk k jk x lx 0)((j x ),当2≥n 时=++∑=)()3(204x l x xk k nk k (324++x x )。
5、设1326)(247+++=x x x x f 和节点,,2,1,0,2/ ==k k x k 则=],,,[10n x x x f 6和=∆07f 25.236494526!77==⨯。
6、5个节点的牛顿-柯特斯求积公式的代数精度为9,5个节点的求积公式最高代数精度为9。
7、{}∞=0)(k kx ϕ是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族,其中1)(0=x ϕ,则⎰=14)(dx x x ϕ0。
8、给定方程组⎩⎨⎧=+-=-221121b x ax b ax x ,a 为实数,当a 满足1<a ,且20<<ω时,SOR 迭代法收敛。
9、解初值问题00(,)()y f x y y x y '=⎧⎨=⎩的改进欧拉法⎪⎩⎪⎨⎧++=+=++++)],(),([2),(]0[111]0[1n n n n n n n n n n y x f y x f h y y y x hf y y 是2阶方法。
10、设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=111a aa a A ,当∈a (22,22-)时,必有分解式T LL A =,其中L 为下三角阵,当其对角线元素)3,2,1(=i l ii 满足(0>ii l )条件时,这种分解是唯一的。
二、选择题1、解方程组b Ax =的简单迭代格式g Bx xk k +=+)()1(收敛的充要条件是(2)。
(1)1)(<A ρ, (2) 1)(<B ρ, (3) 1)(>A ρ, (4) 1)(>B ρ2、在牛顿-柯特斯求积公式:⎰∑=-≈ba ni i n i x f C a b dx x f 0)()()()(中,当系数)(n iC 是负值时,公式的稳定性不能保证,所以实际应用中,当(1)时的牛顿-柯特斯求积公式不使用。
(1)8≥n , (2)7≥n , (3)10≥n , (4)6≥n ,(1)二次; (2)三次; (3)四次; (4)五次4、若用二阶中点公式)),(4,2(1n n n n n n y x f hy h x hf y y +++=+求解初值问题1)0(,2=-='y y y ,试问为保证该公式绝对稳定,步长h 的取值范围为(3)。
(1)20≤<h , (2)20≤≤h , (3)20<<h , (4)20<≤h三、12bx a y +=解:},1{x span =Φ ⎥⎦⎤⎢⎣⎡=2222383125191111TA []3.730.493.320.19=Ty 解方程组y A AC A T T = 其中⎥⎦⎤⎢⎣⎡=3529603339133914A A T ⎥⎦⎤⎢⎣⎡=7.1799806.173y A T 解得:⎥⎦⎤⎢⎣⎡=0501025.09255577.0C 所以 9255577.0=a , 0501025.0=b 2、用8=n 的复化梯形公式(或复化 Simpson 公式)计算dxe x ⎰-10时,(1)试用余项估计其误差。
(2)用8=n 的复化梯形公式(或复化 Simpson 公式)计算出该积分的近似值。
解:001302.0768181121)(12][022==⨯⨯≤''--=e f h a b f R T η])()(2)([2)8(71∑=++=k k b f x f a f hT ]36787947.0)41686207.047236655.05352614.060653066.07788008.08824969.0(21[161++++++⨯+=6329434.0=四、1、方程013=--x x 在5.1=x 附近有根,把方程写成三种不同的等价形式(1)31+=x x 对应迭代格式311+=+n n x x ;(2)xx 11+=对应迭代格式nn x x 111+=+;(3)13-=x x 对应迭代格式131-=+n n x x 。
判断迭代格式在5.10=x 的收敛性,选一种收敛格式计算5.1=x 附近的根,精确到小数点后第三位。
选一种迭代格式建立Steffensen 迭代法,并进行计算与前一种结果比较,说明是否有加速效果。
解:(1)321(31)(-+=')x x ϕ,118.05.1<=')(ϕ,故收敛; (2)xx x 1121)(2+-='ϕ,117.05.1<=')(ϕ,故收敛; (3)23)(x x ='ϕ,15.135.12>⨯=')(ϕ,故发散。
选择(1):5.10=x ,3572.11=x ,3309.12=x ,3259.13=x ,3249.14=x ,32476.15=x ,32472.16=x Steffensen 迭代:k k k k k k k x x x x x x x +---=+)(2))(())((21ϕϕϕϕ11211)1(33323++-++-+-=k k k k k x x x x x计算结果:5.10=x ,324899.11=x ,324718.12=x 有加速效果。
2、已知方程组f AX =,其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=4114334A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=243024f(1) 列出Jacobi 迭代法和Gauss-Seidel 迭代法的分量形式。
(2)求出Jacobi 迭代矩阵的谱半径,写出SOR 迭代法。
解:Jacobi迭代法:⎪⎪⎪⎩⎪⎪⎪⎨⎧=+-=+-=-=+++ ,3,2,1,0)24(41)330(41)324(41)(2)1(3)(3)(1)1(2)(2)1(1k x x x x x x x k k k k k k kGauss-Seidel迭代法:⎪⎪⎪⎩⎪⎪⎪⎨⎧=+-=+-=-=+++++ ,3,2,1,0)24(41)330(41)324(41)1(2)1(3)(3)1(1)1(2)(2)1(1k x x x x x x x k k k k k k k⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--=+-=-0430430430430)(1U L D B J ,790569.0)410(85)(==或J B ρSOR迭代法:⎪⎪⎪⎩⎪⎪⎪⎨⎧=+-+-=+-+-=-+-=+++++ ,3,2,1,0)24(4)1()330(4)1()324(4)1()1(2)(3)1(3)(3)1(1)(2)1(2)(2)(1)1(1k x x x x x x x x x x k k k k k k k k k k ωωωωωω五、1、取步长1.0=h ,求解初值问题⎪⎩⎪⎨⎧=+-=1)0(1y y dxdy用改进的欧拉法求)1.0(y 的值;用经典的四阶龙格—库塔法求)1.0(y 的值。
解:改进的欧拉法:⎪⎩⎪⎨⎧+=++=+=+=++++095.0905.0)],(),([21.09.0),()0(111)0(1n n n n n n n n n n n n y y x f y x f h y y y y x hf y y所以1)1.0(1==y y ;经典的四阶龙格—库塔法:⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧++=++=++==++++=+),()2,2()2,2(),(]22[6342312143211hk y h x f k k hy h x f k k hy h x f k y x f k k k k k hy y n n n n n n n n n n 04321====k k k k ,所以1)1.0(1==y y 。
2、求一次数不高于4次的多项式)(x p 使它满足)()(00x f x p =,)()(11x f x p =,)()(00x f x p '=',)()(11x f x p '=',)()(22x f x p =解:设)(3x H 为满足条件⎩⎨⎧='='=1,0)()()()(33i x f x H x f x H i i i i 的Hermite 插值多项式,则 21203)()()()(x x x x k x H x p --+= 代入条件)()(22x f x p =得:212202232)()()()(x x x x x H x f k ---=六、(下列2题任选一题,4分) 1、数值积分公式形如⎰'+'++=≈1)1()0()1()0()()(f D f C Bf Af x S dx x xf 1,试确定参数D C B A ,,,使公式代数精度尽量高;2,设]1,0[)(4C x f ∈,推导余项公式⎰-=1)()()(x S dx x xf x R ,并估计误差。
解:将32,,,1)(x x x x f =分布代入公式得:201,301,207,203-====D B B A构造Hermite 插值多项式)(3x H 满足⎩⎨⎧='='=1,0)()()()(33i x f x H x f x H i i i i 其中1,010==x x则有:⎰=103)()(x S dx x xH , 22)4(3)1(!4)()()(-=-x x f x H x f ξdx x x f dx x S x f x x R 213)4(10)1(!4)(])()([)(-=-=⎰⎰ξ1440)(60!4)()1(!4)()4()4(1023)4(ηηηf f dx x x f=⨯=-=⎰2、 用二步法)],()1(),([111101---+-+++=n n n n n n n y x f y x f h y y y θθαα求解常微分方程的初值问题⎩⎨⎧=='00)(),(y x y y x f y 时,如何选择参数θαα,,10使方法阶数尽可能高,并求局部截断误差主项,此时该方法是几阶的。