2008-2009第一学期数值计算方法A答案
- 格式:doc
- 大小:293.00 KB
- 文档页数:8
2008—2009学年第一学期《数值计算方法》试卷答案及评分标准专业班级姓名学号开课系室信息与计算科学系考试日期2009年1月12日说明:1.封面及题目的背面为草稿纸.2.答案必须写在该题下方空白处,不得写在草稿纸中,否则答案无效.3.正文共7页一、选择题(4×5=20分)1.已知:⎪⎪⎭⎫⎝⎛=4321A ,则∞A = D 。
A ) 2;B ) 3;C ) 6;D ) 7。
2.数值*x 的近似值310113.0-⨯=x ,若满足*x x -≤ B ,则称x 有3位有效数字。
A )12; B ) 61102-⨯; C )11102-⨯; D)31102-⨯。
3.若用复化梯形公式,计算10x e dx ⎰,区间至少应分 B 等份,才能保证结果有5位有效数字。
A ) 67;B ) 68;C ) 69;D ) 70。
4.设,1221⎪⎪⎭⎫⎝⎛-=A 则=)(1A cond C 。
A )51; B )54; C )59; D )57。
5.形如∑⎰=≈nk k k bax f A dx x f 0)()(的插值型数值求积公式,它的代数精度最高不会超过 A 阶。
A ) 12+n ;B ) 22+n ;C )n 2;D ) 12-n 。
二、填空题(4×5=20分)1.使用SOR 迭代方法求解线性方程组b Ax =时,若使得方法收敛,则必有迭代因子∈ω (0,2) 。
2.已知()(2)(3)(4)(5)f x x x x x =----,则差商]5,4,3[f = 0 。
3.设序列{}k x 是收敛于方程0)(=x f 的根*x 的迭代序列,记,*k k x x e -=,1,0=k ,表示各步的迭代误差,使得0lim1≠=+C e e pkk 成立,称序列{}k x 是p阶收敛。
4.)(x f 在],[b a 上有连续的n 阶导数,)()1(x f n +在),(b a 上存在,设)(x f 的n 次Lagrange插值多项式为)(x P n ,余项=-=)()()(x P x f x R n n)()!1()(0)1(∏=+-+nj j n x x n f ξ。
5.设m n n m R b R x R A ∈∈∈⨯,,,n R ∈η是方程组b Ax =的最小二乘解的充要条件为η是法方程组b A Ax A T T =的解。
三、计算题(共20分)1、(6分)利用牛顿迭代法求130的近似值,取130=x 为初始值,精确到六位有效数字。
解:设x 为所求,即求01302=-x 的正根,130)(2-=x x f .取130=x ,有迭代公式1k x +=k x -)()(k k x f x f '=),2,1,0(2130221302=+=--k x x x x x k k k k k ……3分1x ==⨯+13213021311.50002x =5000.11213025000.11⨯+=11.40223x =4022.11213024022.11⨯+=11.4018 4018.11213024018.114⨯+=x =4018.11213024018.11⨯+=11.4018x *≈11.4018 ……6分2、(7分)已知:8)3(,5)2(,1)1(===f f f ,构造()y f x =以3,2,1=x 为节点的拉格朗日插值多项式。
解:22110020)()()()()(y x l y x l y x l y x l x P i i i n ++==∑=, ……2分)3)(2(21)31)(21()3)(2())(())(()(2010210--=----=----=x x x x x x x x x x x x x l ,……3分)3)(1()32)(12()3)(1())(())(()(2101201---=----=----=x x x x x x x x x x x x x l ,……4分)2)(1(21)23)(13()2)(1())(())(()(1202102--=----=----=x x x x x x x x x x x x x l ,……5分4211218)2)(1(215)3)(1(1)3)(2(21)(2-+-=⨯--+⨯---⨯--=x x x x x x x x x p n ……7分。
3、(7分)求出改进的Euler 方法的绝对稳定区间。
解:写出改进的Euler 的公式:)),(),((2111+++++=n n n n n n y x f y x f hy y ……2分 将其应用于试验方程得:)(211++++=n n n n y y hy y λλ……4分 整理得:nn n y h E y h hy )(21211λλλ=-+=+……5分要使格式绝对稳定必满足:1)( h E λ ……6分解得:0 h λ∞-,即绝对稳定区间为)0,(-∞。
……7分四、(10分)构造下列方程组收敛的雅可比和高斯-赛德迭代格式⎪⎩⎪⎨⎧=++=++-=+-143101410353102321321321x x x x x x x x x 。
解:由于直接分离321,,x x x 不能满足收敛条件,因此将方程组变形为:⎪⎩⎪⎨⎧=+-=+-=++141035310214310321321321x x x x x x x x x + ……2分 对应系数矩阵⎪⎪⎪⎭⎫ ⎝⎛-=103131021310A 为严格对角占忧,分离得⎪⎪⎪⎩⎪⎪⎪⎨⎧+--=++=+--=5710310*********7101103213312321x x x x x x x x x 4分则收敛的Jacobi 迭代格式为:⎪⎪⎪⎩⎪⎪⎪⎨⎧+--=++=+--=+++57103101211035157101103)(2)(1)1(3)(3)(1)1(2)(3)(211k k k k k k k k k x x x x x x x x x )( 7分 则收敛的G-S 迭代格式为:⎪⎪⎪⎩⎪⎪⎪⎨⎧+--=++=+--=++++++57103101211035157101103)1(2)1(1)1(3)(3)1(1)1(2)(3)(211k k k k k k k k k x x x x x x x x x )( 10分五、(10分)用列主元Gauss 消去法求解线性方程组⎪⎩⎪⎨⎧-=++-=++=++754217743322321321321x x x x x x x x x 。
解:线性方程组的系数增广矩阵为:→⎪⎪⎪⎭⎫ ⎝⎛--=754217743322)|(b A ⎪⎪⎪⎭⎫⎝⎛--754233221774(第一次选主元) 2分⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛---→213217215025212301774(第一次Gauss 消元) 4分 ⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛----→252123021321721501774(第二次选主元) 6分 ⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛--→56560021321721501774 (第二次Gauss 消元) 8分 则得方程组的精确解()Tx 2,2,1*-=。
10分六、(10分)确定待定参数,使其代数精度尽是高,并指明所得公式具有的代数精度。
)()0()()(10122h f A f A h f A dx x f hh++-≈-⎰-解:令2()1,,f x x x =使公式两边准确成立,则有⎪⎪⎩⎪⎪⎨⎧++=++-=++=------⎰⎰⎰120122221012210122001A h A A h dx x hA A hA xdx AA A dx hhhh h h ……3分 即⎪⎪⎩⎪⎪⎨⎧+=-=++=---)(316)(0411*******A A h h A A h A A A h 可求得h A h A A 34,38011-===- ……5分则由该系数所确定的求积公式为)(38)0(34)(38)(22h hf hf h hf dx x f hh+--≈⎰-下面验证该公式的代数精度。
令3()f x x =,则0)(22322==⎰⎰-dx x dx x f hhhh-038)(38)(38)0(34)(3833=+-=+--hh h h h hf hf h hf 即公式左右两边相等。
……7分 令4()f x x =,则522422564)(h dx x dx x f hhhh==⎰⎰-- 54431638)(38)(38)0(34)(38h hh h h h hf hf h hf =+-=+--即公式两边不相等。
……9分 所以此求积公式具有三次代数精度。
……10分七、(10分)设x x x A R A n n δ+≠∈⨯,,0d e t , 分别满足方程组b Ax =,b b x x A A δδδ+=++))((。
其中0≠b ,且A δ适当小,使)(1A cond AAδ。
求证:⎪⎪⎭⎫ ⎝⎛+-≤-b b A A A A A cond xxδδδδ11)(,这里用到的是任何一种向量范数及从属于它的矩阵范数。
证明: 设βδδα=+-==-)(,1A A A A A,由条件)(1A cond AAδ可知11 A A δαβ-= ……1分 由摄动定理,即n n R C A ⨯∈,,若1-A 存在,且1,,1 αββα≤--C A A ,则C 可逆,且αβα-≤-11C 。
可知A A δ+可逆,且()AAA A A δδ1111----≤+, ……3分则扰动方程组b b x x A A δδδ+=++))((存在唯一解,())(1b b A A x x δδδ++=+-。
……4分由此可知()()x A b A A x A A b b A A x )()()()(11δδδδδδδ-+=+-++=-- ……6分⎪⎪⎭⎫ ⎝⎛+-≤++≤---x A b b x A A A A A A x A b A A x δδδδδδδ1111)()( …9分 得出:⎪⎪⎭⎫ ⎝⎛+-≤-b b A A A A A cond x xδδδδ11)( ……10分 定理得证。