微分方程数值解(学生复习题)

合集下载

数值分析复习题(下)

数值分析复习题(下)

线性代数方程组直接法题1.⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫⎝⎛-----710413221232321x x x 程组用直接三角分解法解方2.21,,4321A A A A ∞⎥⎦⎤⎢⎣⎡--=,求设 3.用顺序高斯消去法解⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛565331743532321x x x4.用列主元高斯消去法解⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛-20111.0310********x x x5.用直接三角分解法解⎪⎪⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫ ⎝⎛7173530103421101002014321x x x x6.用追赶法求解方程组⎪⎪⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫ ⎝⎛022112111131124321x x x x7.试用平方根法解下列对称正定方程组⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛--103422484548416321x x x8.用高斯-若当消去法求下列矩阵的逆⎪⎪⎪⎭⎫ ⎝⎛---=112221111A作业1.用追赶法解三对角方程组Ax=b,其中A=⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--------2100012100012100012100012,b=⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛00001 作业2.用改进的平方根法解方程组⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛---654131*********x x x 作业3.设A=⎪⎪⎭⎫⎝⎛3.01.05.06.0,计算21A A A ,,∞ 作业4.设A 为非奇异矩阵,求证:∞∞≠∞-=yAy Ay 01min1作业5.下述矩阵能否分解为LU (其中L 为单位下三角阵,U 为上三角阵)?若能分解,那么分解是否是唯一?⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=461561552621,133122111,764142321C B A解线性方程组的迭代法1.设方程组Ax=b,⎪⎪⎪⎭⎫⎝⎛-=122111221A ,试讨论解此方程组的J 法和GS 法的收敛性。

《微分方程数值解法》复习、练习题

《微分方程数值解法》复习、练习题

《微分方程数值解法》复习、练习题第一章复习题1、建立差分格式的三个主要步骤(三个离散化)。

2、差分格式的相容性、收敛性概念。

3、Poisson 方程的5点菱形差分格式,矩形、非矩形区域情形边界条件的处理(离散化)。

4、对长方形区域作正方形网格剖分,求解Poisson方程边值问题的五点菱形差分格式,按什么顺序对节点编号,可使差分方程带宽更窄?(按短方向排)5、差分方程有哪些共同特性,求解选用哪类方法?(大型稀疏,带状,主对角占优等,一般采用迭代法)多重网格等略。

6、极值原理。

7、5点菱形差分格式求解Poisson 方程第一边值问题的收敛性。

第一章练习题1、设有边值问题取h=0.1的正方形网格。

(1)用5点菱形格式在内点建立差分格式;(2)用截断误差为的方法离散化第三边界条件(有两种方式);(3)写出整理后的差分方程的矩阵形式2、定义方形算子如下:试讨论5点方形差分方程逼近微分方程的截断误差是几阶?3、设有,取h=1/3,列出5点方形差分格式所得的差分方程。

第二章复习题1、差分格式稳定性与收敛性的定义。

2、有关求特征值的几个结论。

3、判断稳定性的矩阵法和Fourier分析法(Von-Neumann条件)的应用。

4、显隐格式在一般情况下的优缺点。

5、熟悉古典显、隐格式,六点对称隐格式(C-N格式)。

6、叙述Lax等价定理。

7、高维抛物型方程的ADI格式的优点。

8、了解非线性方程差分格式的建立,讨论稳定性的冻结系数法。

第二章练习题1、设有求解抛物型方程组的初值问题的差分格式试写出用Fourier分析法讨论稳定性时的增长矩阵。

2、对上题考虑另一个差分格式试讨论该格式的稳定性。

3、对抛物型方程,考虑著名的Du Fort-Frankel(1953)格式(1)推导该格式是否满足稳定性的Von-Neumann条件?(2)该格式与Richardson格式有什么关系?4、讨论求解的古典显格式的稳定性。

5、写出逼近的古典显格式。

【复习题】数值分析复习题及答案65177

【复习题】数值分析复习题及答案65177

【关键字】复习题数值分析复习题一、选择题1. 3.142和3.141分别作为的近似数具有()和()位有效数字.A.4和3 B.3和2 C.3和4 D.4和42. 已知求积公式,则=()A. B.C.D.3. 通过点的拉格朗日插值基函数满足()A.=0,B.=0,C.=1,D.=1,4. 设求方程的根的牛顿法收敛,则它具有()敛速。

A.超线性B.平方C.线性D.三次5. 用列主元消元法解线性方程组作第一次消元后得到的第3个方程().A.B.C. D.二、填空1. 设,取5位有效数字,则所得的近似值x= .2.设一阶差商,则二阶差商3. 设, 则,。

4.求方程的近似根,用迭代公式,取初始值,那么5.解初始值问题近似解的梯形公式是6、,则A的谱半径=。

7、设,则和。

8、若线性代数方程组AX=b 的系数矩阵A为严格对角占优阵,则雅可比迭代和高斯-塞德尔迭代都。

9、解常微分方程初值问题的欧拉(Euler)方法的局部截断误差为。

10、为了使计算的乘除法运算次数尽量的少,应将表达式改写成。

11. 设, 则, .12. 一阶均差13. 已知时,科茨系数,那么14. 因为方程在区间上满足,所以在区间内有根。

15. 取步长,用欧拉法解初值问题的计算公式.16.设是真值的近似值,则有位有效数字。

17. 对, 差商( )。

18. 设, 则。

19.牛顿—柯特斯求积公式的系数和。

20. 若a=2.42315是2.42247的近似值,则a有( )位有效数字.21. 是以为插值节点的Lagrange插值基函数,则( ).22. 设f (x)可微,则求方程的牛顿迭代格式是( ).23. 迭代公式收敛的充要条件是。

24. 解线性方程组Ax=b (其中A非奇异,b不为0) 的迭代格式中的B称为( ). 给定方程组,解此方程组的雅可比迭代格式为( )。

25、数值计算中主要研究的误差有和。

26、设是n次拉格朗日插值多项式的插值基函数,则;。

微分方程数值解习题(李立康)

微分方程数值解习题(李立康)

常微分方程习题 《李立康》习题1.用Euler 方法求初值问题⎩⎨⎧=-='0)0(21u tuu 在1=t 时的近似解(取41=h )。

2.初值问题1300u u u()⎧⎪'=⎨⎪=⎩ 有解3223/u(t )t ⎛⎫= ⎪⎝⎭。

但若用Euler 方法求解,对一切N T ,和HTh =,都只能得到N t u t , (2)1,0==,试解释此现象产生的原因。

3.用Euler 方法计算⎩⎨⎧=='1)0(u uu 在1=t 处的值,取161和41=h ,将计算结果与精确值e =)1(u 相比较。

4.设),(u t f 满足定理2.1的条件,对改进Euler 法(2.10)式证明: (1)其局部截断误差为)()(1243h O t u h -'''-;(2)当1<hL 时,其整体截断误差满足:)1(22--≤Lt n lT m e hLRe εε (3)方法具有二阶收敛速度且稳定。

5.导出用改进Euler 法求解⎩⎨⎧=='1)0(u uu 计算公式mmh h u ⎪⎪⎭⎫ ⎝⎛-+=22 取41=h 计算)1(u 的近似值,并与习题3的结果比较。

6.就初值问题⎩⎨⎧=+='0)0(u bat u 分别导出用Euler 方法和改进Euler 法求近似解的表达式,并与真解bt t au +=22相比较。

7.证明改进Euler 法的绝对稳定区域是整个左半平面0)Re(<h 。

8.对初值问题⎩⎨⎧=-='1)0(2u u u 用41=h 的Euler 方法求解,求出实际计算值t u 与真解tu +=11在)1(u 处的误差,并将它与定理2.3的估计式(2.22)式相比较。

9.证明:Runge-Kutta 方法中);,(h u t ϕ关于u 或t 满足Lipschitz 条件的充分条件是),(u t f 关于t 或u 满足Lipschitz 条件。

微分方程数值解法(戴嘉尊_第二版)习题讲解

微分方程数值解法(戴嘉尊_第二版)习题讲解
第四章双曲型方程的差分方法............................................................25
电子文档制作:成都信息工程学院数学学院杨韧吴世良,2010年4月
第二章抛物型方程的差分方法..............................................................8
第三章椭圆型方程的差分方法............................................................16
成都信息工程学院>>精品课程>>微分方程数值解
微分方程数值解
习题解答Байду номын сангаас
杨韧吴世良(编)
成都信息工程学院
数学学院
二O一O年四月编写
电子文档制作:成都信息工程学院数学学院杨韧吴世良,2010年4月
成都信息工程学院>>精品课程>>微分方程数值解
第一章常微分方程数值解......................................................................3

偏微分方程数值解期末试题及参考答案

偏微分方程数值解期末试题及参考答案

《偏微分方程数值解》试卷参考答案与评分标准专业班级信息与计算科学开课系室考试日期 2006.4.14命题教师王子亭偏微分方程数值解试题(06A)参考答案与评分标准信息与计算科学专业一(10分)、设矩阵A 对称正定,定义)(),(),(21)(n R x x b x Ax x J ∈-=,证明下列两个问题等价:(1)求n R x ∈0使 )(min )(0x J x J nRx ∈=;(2)求下列方程组的解:b Ax =解: 设n R x ∈0是)(x J 的最小值点,对于任意的n R x ∈,令),(2),()()()(2000x Ax x b Ax x J x x J λλλλϕ+-+=+=, (3分)因此0=λ是)(λϕ的极小值点,0)0('=ϕ,即对于任意的n R x ∈,0),(0=-x b Ax ,特别取b Ax x -=0,则有0||||),(2000=-=--b Ax b Ax b Ax ,得到b Ax =0. (3分) 反之,若nR x ∈0满足bAx =0,则对于任意的x ,)(),(21)0()1()(00x J x Ax x x J >+==+ϕϕ,因此0x 是)(x J 的最小值点. (4分)评分标准:)(λϕ的表示式3分, 每问3分,推理逻辑性1分二(10分)、 对于两点边值问题:⎪⎩⎪⎨⎧==∈=+-=0)(,0)(),()(b u a u b a x f qu dxdu p dx d Lu 其中]),([,0]),,([,0)(min )(]),,([0min ],[1b a H f q b a C q p x p x p b a C p b a x ∈≥∈>=≥∈∈建立与上述两点边值问题等价的变分问题的两种形式:求泛函极小的Ritz 形式和Galerkin 形式的变分方程。

解: 设}0)()(),,(|{110==∈=b u a u b a H u u H 为求解函数空间,检验函数空间.取),(10b a H v ∈,乘方程两端,积分应用分部积分得到 (3分))().(),(v f fvdx dx quv dxdv dx du pv u a b a ba ==+=⎰⎰,),(10b a H v ∈∀ 即变分问题的Galerkin 形式. (3分)令⎰-+=-=b a dx fu qu dxdup u f u u a u J ])([21),(),(21)(22,则变分问题的Ritz 形式为求),(10*b a H u ∈,使)(min )(1*u J u J H u ∈= (4分)评分标准:空间描述与积分步骤3分,变分方程3分,极小函数及其变分问题4分,三(20分)、对于边值问题⎪⎩⎪⎨⎧=⨯=∈-=∂∂+∂∂∂0|)1,0()1,0(),(,12222G u G y x yux u (1)建立该边值问题的五点差分格式(五点棱形格式又称正五点格式),推导截断误差的阶。

微分方程数值解法答案

微分方程数值解法答案

微分⽅程数值解法答案包括基本概念,差分格式的构造、截断误差和稳定性,这些内容是贯穿整个教材的主线。

解答问题关键在过程,能够显⽰出你已经掌握了书上的内容,知道了解题⽅法。

这次考试题⽬的类型:20分的选择题,主要是基本概念的理解,后⾯有五个⼤题,包括差分格式的构造、截断误差和稳定性。

习题⼀1.略2. y y x f -=),(,梯形公式:n n n n n n y hh y y y h y y )121(),(2111+-+=+-=+++,所以0122)1(01])121[()121()121(y hh y h h y h h y hhn h h n n n +--+--+-+=+-+==+-+= ,当0→h 时,x n e y -→。

同理可以证明预报-校正法收敛到微分⽅程的解.3.局部截断误差的推导同欧拉公式;整体截断误差:++++++-++≤1),())(,(11111n nx x n n n n n n n dx y x f x y x f R εε11)(++-++≤n n n y x y Lh R ε,这⾥R R n ≤ ⽽111)(+++-=n n n y x y ε,所以 R Lh n n +=-+εε1)1(,不妨设1()]11111[1111101---++-+-+-≤≤-+-=n n n n Lh Lh Lh R Lh Lh R Lh εεε ]1[2)(02)(00-+≤--x X L x X L eLh R eε4.中点公式的局部截断误差: dx x y x f hx y h x f x y x f yx y n n x x n n n n n n))](,(2)(,2())(,([)(11*1?+++-=-++dx x y x f hx y h x f h x y h x f h x y x y dxx y x f hx y h x f hx y h x f h x y h x f x y x f n n n n x x n n n n n n n x x n n n n n n n n))](,(2)(,2())2(,2([)]2()([))](,(2)(,2())2(,2())2(,2())(,([11++-++++'-'=++-+++++-=??++所以上式为+--+''=?++dx hx x x y e n nx x n n n )2()(11θdx x y x f h x y h x f h x y h x f n n n n x x n n n n))](,(2)(,2())2(,2([1++-++?+ 3218)(LMh h x y Lh e n n ≤+''≤+?中点公式的整体截断误差:dx y x f hy h x f x y x f y x y y x y n n x x n n n n n n n n)],(2,2())(,([)()(111?+++-+-=-++dxy x f hy h x f x y x f h x y h x f x y x f hx y h x f x y x f y x y n n n n n n n n x x n n n n n n n n))],(2,2()))(,(2)(,2()))(,(2)(,2())(,([)(1++-+++++-+-=?+因⽽n n n L h Lh R εεε)21(1+++≤+,R L h Lh n n +++≤-122)21(εε≤≤])21()21(1[2)21(1222222022-+++++++--+++n nL h Lh L h Lh Lh Lh RL h Lh ε )1(00-+≤--x X L x X L e LhR eε 5.略 6.略 7.略8.(1)欧拉法:2.0≤h ;四阶Runge-Kutta ⽅法:278.0≤h (2)欧拉法:3 54≤h ;四阶Runge-Kutta ⽅法:3556.5≤h(3)欧拉法:1≤h ;四阶Runge-Kutta ⽅法:278.0≤h 9.略 10.略习题21.略 2.略 3.略4.差分格式写成矩阵形式为:n n M n M n n n M n M n n e u u u u r t r r r t r r r t r r r t u u u u +?--------= --+-+-++12211221121212121 αβαααβαααβαααβ矩阵的特征值为:)cos(221Mj r r t j πααβλ+-?-=,要使格式稳定,则特征值须满⾜t c j ?+≤1λ,即21≤r α5.利⽤泰勒展式可以得到古典隐式差分格式的截断误差为)(2h t O +?。

微分方程数值解(学生复习题)

微分方程数值解(学生复习题)

一.填空1.Euler 法的一般递推公式为,整体误差为 ,局部截断误差为:.,改进Euler 的一般递推公式 整体误差为,局部截断误差为:。

2.线性多步法绝对稳定的充要条件是。

3.当,则单步法1(,,)0,1,2,,n n n n Tu u h t u h n hϕ+=+=,稳定。

4. 一个相容,稳定的多步法若绝对稳定,则绝对稳定域在。

5. 若,则多步法是相容的。

6.所有内点,界点的差分方程组成一个封闭的线性代数方程组,其系数矩阵是。

7.刚性方程是:8.Runge-Kutta 法的特征值为 ,相容的充要条件为:8.二阶常微分方程边值问题:22,(), ()d uLu qu f a x bdxu a u b αβ⎧=-+=<<⎪⎨⎪==⎩ 的中心差分格式为:P i 的四个相邻点均属于h G ,则称P i 为。

10.逼近泊松方程的五点差分格式的截断误差的阶为。

逼近泊松方程的九点差分格式的截断误差的阶为。

12.SOR 收敛当且仅当松弛因子0,2ω∈(),且Jacobi 迭代收敛。

最佳松弛因子是。

二.判断τ和空间步长h 无限缩小时,差分格式的解是否逼近到微分方程问题的解,这就是差分格式的收敛性问题。

2.单参数的PR 迭代格式的收敛速度与SOR 最佳超松弛法的收敛速度同阶。

3、对称矩阵的普条件数与条件数相同。

4、一级Runge-Kutta 法的绝对稳定域(-2,0)5、若差分方程满足相容条件,且按右端稳定,则差分解收敛至波动方程的解。

6、Euler 法非A 稳定。

7.对任意网比0r >,六点对称格式的解有收敛阶22()O h τ+ 8.对任意网比12r ≤,向前差分格式的解有收敛阶2()O h τ+。

9、相容,稳定的多步法一定绝对稳定。

三.选择1.抛物型方程的加权隐式差分格式的稳定性为()A 绝对稳定B 无条件稳定C 条件稳定D 非条件稳定 2.von Neumann 条件是差分格式稳定的()A 充分条件B 必要条件C 充要条件D 既非充分也非必要条件 3.实系数二次方程20b c λλ--=的根按模小于或者等于1的充要条件是() A 12b c ≤-≤ B 1+2b c ≤≤ C 12c b ≤-≤ D 12c b ≤+≤ 4.若线性多步法A 稳定,则有( ),其中1,2,,i i k λ=()为()()0h ρλσλ-=的根。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

一.填空
1. Euler 法的一般递推公式为 ,整体误差为 ,局部截断误差为: .,改进Euler 的一般递推公式 整体误差为 ,局部截断误差为: 。

2. 线性多步法绝对稳定的充要条件是 。

3.当 ,则单步法1(,,)0,1,2,,n n n n T u u h t u h n h
ϕ+=+= ,稳定。

4. 一个相容,稳定的多步法若绝对稳定,则绝对稳定域在 。

5. 若 ,则多步法是相容的。

6.所有内点,界点的差分方程组成一个封闭的线性代数方程组,其系数矩阵是 。

7.刚性方程是:
8.Runge-Kutta 法的特征值为 ,
相容的充要条件为:
8.二阶常微分方程边值问题:22,(), ()d u Lu qu f a x b dx u a u b αβ⎧=-+=<<⎪⎨⎪==⎩
的中心差分格式为:
9.若内点P i 的四个相邻点均属于h G ,则称P i 为 。

10.逼近泊松方程的五点差分格式的截断误差的阶为 。

逼近泊松方程的九点差分格式的截断误差的阶为 。

11.线性多步法A 稳定的充要条件是 。

12. SOR 收敛当且仅当松弛因子0,2ω∈()
,且Jacobi 迭代收敛。

最佳松弛因子是 。

二.判断
1.当时间步长τ和空间步长h 无限缩小时,差分格式的解是否逼近到微分方程问题的解,这就是差分格式的收敛性问题。

2.单参数的PR 迭代格式的收敛速度与SOR 最佳超松弛法的收敛速度同阶。

3、对称矩阵的普条件数与条件数相同。

4、一级Runge-Kutta 法的绝对稳定域(-2,0)
5、若差分方程满足相容条件,且按右端稳定,则差分解收敛至波动方程的解。

6、Euler 法非A 稳定。

7.对任意网比0r >,六点对称格式的解有收敛阶22()O h τ+
8. 对任意网比12
r ≤,向前差分格式的解有收敛阶2()O h τ+。

9、相容,稳定的多步法一定绝对稳定。

三.选择
1.抛物型方程的加权隐式差分格式的稳定性为()
A 绝对稳定
B 无条件稳定
C 条件稳定
D 非条件稳定
2.von Neumann 条件是差分格式稳定的()
A 充分条件
B 必要条件
C 充要条件
D 既非充分也非必要条件
3.实系数二次方程20b c λλ--=的根按模小于或者等于1的充要条件是() A 12b c ≤-≤ B 1+2b c ≤≤ C 12c b ≤-≤ D 12c b ≤+≤
4.若线性多步法A 稳定,则有( ),其中1,2,,i i k
λ= ()为()()0h ρλσλ-=的根。

A Re 01,1,2,,i h i k λ<⇒>= B 1Re 0i h λ≥⇒≥ C Re 01,1,2,,i h i k λ≤⇒≤= D 1Re 0i h λ<⇒<
5.一个相容,稳定的多步法若绝对稳定,则绝对稳定域在()
A 下半平面
B 上半平面
C 左半平面
D 右半平面
6.线性多步法稳定的充要条件是()
A 第一特征式()ρλ满足根条件
B 第一特征式()ρλ严格满足根条件
C ()()0h ρλσλ-=满足根条件
D ()()0h ρλσλ-=严格满足根条件
7. P 阶K 步法的局部截断误差的阶为( )
A p O h ()
B 1p O h +()
C 1k O h +()
D 1k O h +()
8. 线性多步法绝对稳定的充要条件是( )
A 第一特征式()ρλ满足根条件
B 第一特征式()ρλ严格满足根条件
C ()()0h ρλσλ-=满足根条件
D ()()0h ρλσλ-=严格满足根条件
9.Euler 法的整体误差为( )
A O h ()
B 2O h ()
C 1O h -()
D 1O ()
四.计算
1.试求差分方程初值问题:
21012320n n n u u u u u ++--=⎧⎨==⎩
的解。

2.已知显式方法
[]2110110n n n n n u u u h f f ααββ+++++=+
(1) 取1α为参数,确定001αββ,,,使方法至少是二阶的;
(2) 当1α取何值时,方法满足根条件;
3. k 步线性法:[]2
n k n n k n hk u u f f ++=++,证明其A 稳定。

4.证明11n n n u u hf ++=+对所有的(),0h ∈-∞都绝对稳定。

5.由待定系数法构造边值问题:
,()()0u f a x b u a u b ''=<<⎧⎨==⎩ 的中心差分格式。

6.求正三角网上的差分格式。

7.用有限体积法推导五点格式。

8.写出扩散方程22u u a t x
∂∂=∂∂的向前,向后差分方程(中心差分格式,用第n 层计算第n+1层),并把有限差分方程改写成便于计算的迭代格式(矩阵形式),2
a r h τ
=为网比。

9.计算差分格式()11-1n n n n j j j j u u r u u ++=--,(其中,0a r a h
τ=>)的增长因子,并根据von Neumann 条件给出差分格式稳定性条件。

10. 已知线性多步法:
212412+333
n n n n h u u u f +++-= 试求它的阶及误差常数。

11.计算向前,向后等差分格式的增长因子,并给出稳定性条件。

12. Adams 二步外插法:2113122n n n n u u h f f +++⎡⎤-=-⎢⎥⎣⎦
,试求其绝对稳定域。

五.证明题
1.将三层差分格式改写为改写成等价的二层差分格式,写出其增长矩阵,并由 von Neumann 条件证明该格式是否稳定。

其他例子关于证明差分格式稳定或者不稳定(参考书上的课后习题及例题)。

2. 求N 阶三角阵:
01101101C=10110⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦ 或者111-111-11C=1-1111-⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦
的特征值和特征向量,并证明矩阵是病态的。

3. 证明Euler 向后公式A 稳定:11n n n u u hf ++=+。

4. 证明:梯形公式:[]112n n n n h u u f f ++=+
+,证明其A 稳定。

相关文档
最新文档