《计算机数值方法》测试题二

《计算机数值方法》测试题

.判断题（1分X 10=10分）（对打V,错打X）

1. 数值方法是指解数值问题的计算机上可执行的系列计算公式。（）

2. ……计算R二e-2.71828〜0.00000182是截断误差。（）

3. 不同的矩阵三角

分解对应着不同的解法，但在本质上，都是经过A=LU的分解计算, 再解Ly=b和Ux=y的线性方程组。（）

4. 一般不用n次多项式做插值函数。（）

5. FUnge现象说明并非插值多项式的次数越高其精度就越高。（）

6. FOmberg算法是利用加速技术建立的。（）

7. 从复合求积的余项表达式看，计算值的精度与步长无关。（）

8. 可用待定系数法和函数值或公式的线性组合构造新的数值函数求解微分方程。

（）

9. 局部截断误差e k （h）与y （X k）的计算值y k有关。（）

10. 对大型线性方程组和非线性方程采用逐次逼近更为合适。（）

.填空题（2分X 5=10分）

1. 设x € [a,b] , x丰X0, 则一阶均差f （x）= ____________________ 。

2. 矩阵A的F-范数||A|| F=___________________________ 。

3. Euler公式为____________________________ 。

4. 矩阵A的条件数Co nd （A）® = ________________。

5._______________________________________________________________________ 设x为准确值，x*为x的一个近似值，近似值x*的相对误差E （x*）=

_________________________________________________________________________ .选择题（2分X 5=10分）

1. 设x=Pi ;则x*=3.1415有（）位有效数字。

（A） 4 位（B）5 位（C）6 位

2 .顺序主元a ii工0 (i=1,2 k）的充要条件是A的顺序主子式D （i=1,2n-1 )

(A)不全为0 (B) 全不为0 (C) 全为0

3 .若存在实数P> 1和c >0,则迭代为P阶收敛的条件是(

(A) k lim 熙二c (B) O(h p) (C) O(h p+1)

k le k |p

4 .方程x3-x2-1=0在x°=1.5附近有根，则迭代格式X k+i=在x°=1.5附近( )

(A)不收敛(B) 局部收敛(C) 不确定

5.下面哪个公式的局部截断误差为0( h3)0 ()

(A) Euler公式(B)三阶Runge— Kutta公式(C)梯形公式

四.计算题(7分X 6=42分)

1. 要使18的近似值的相对误差限小于0.1 %要取几位有效数字？

2. 用Gauss列主元素消去法求解方程组

.11 12x1-3x2+3x3=15

-18x 汁3X2-X3=-15

X1+ x 2+ x 3=6

3 .已知结点如下：不用幵方的办法求117的值

4. x3-2x 2-4x-7=0在区间［3 , 4］内有根，自选迭代法求解方程的根，精确到10-3

1 o

5 .用复合公式求解定积分：° 1/ (1+x2) dx (n=8)

6 .在［0,1］上求解初值问题，取步长h=0.2 , y =x+1, y(0)=1

五.算法设计(7分X 2=14分)

1. Lagrange插值公式为:

n

printf （" 迭代%d 次之后，没有解.\n",m ）;

2. 给出用二分法解X 2-X +2=0的算法框图

六.编程填空（2分X 7=14分）

1 .用牛顿迭代法解方程：e X -3-x=0

#in clude #in clude #defi ne X 0 2 #defi ne m 1000 #defi ne eps 0.000001 mai n()

{int i;double X 1=X 0,X 2=X 0; for(i=O;i< ;i++

)

{pri ntf("%d %f\n",i,x2); x2=(x1-(exp(x1)-3-x1)/(exp(x1)-1)); if(fabs(x2-x1) eps)

{pri ntf("the root is x=%f,k=%d\n",x2 return; }

X 1=X 2;

}

1 i

( X )

i 0

(X-X j ) /

i 0

(X i

-X j

)

给出算法框图

2. 用列主元素消去法解方程组：

严

X i+2X2-X3=3

x仁x 2+5X3=0

4x i+X2-2x 3=0

#in clude

#in clude

#defi ne n 3

static double aa [n][n+1]={{1,2,-1,3},{1,-1,5,0},{4,1,-2,2}};

mai n()

{int i,j,det,k,c;

double a [n+1][ n+2],x [n+1],r,t,m;

for(i=1;i<= ;i++)

for(j=1;j<= ;j++)

a[i][j]=aa[i-1][j-1];

for (k=1;k<=n-1;k++)

{r=a[k][k];c=k;

for(i=k;i< 二n ;i++)

if(fabs(a[i][k]) fabs(r))

{r=a[i][k];c=i;}

if(c!=k)

for(j=k;j< 二n+1;j++)

{t=a[k][j]; =acja[c][j]=t;}