北京大学数值分析试题2015
- 格式:doc
- 大小:326.00 KB
- 文档页数:10
数值分析复习题一、选择题1. 3.142和3.141分别作为π的近似数具有〔 〕和〔 〕位有效数字.A .4和3B .3和2C .3和4D .4和42. 求积公式()()211211()(2)636f x dx f Af f ≈++⎰,那么A =〔 〕A . 16B .13C .12D .233. 通过点()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足〔 〕A .()00l x =0,()110l x = B .()00l x =0,()111l x = C .()00l x =1,()111l x = D . ()00l x =1,()111l x =4. 设求方程()0f x =的根的牛顿法收敛,那么它具有〔 〕敛速。
A .超线性B .平方C .线性D .三次5. 用列主元消元法解线性方程组1231231220223332x x x x x x x x ++=⎧⎪++=⎨⎪--=⎩ 作第一次消元后得到的第3个方程〔 〕.A .232x x -+= B .232 1.5 3.5x x -+= C .2323x x -+= D .230.5 1.5x x -=-二、填空1. 设2.3149541...x *=,取5位有效数字,那么所得的近似值x=.2.设一阶差商()()()21122114,321f x f x f x x x x --===---,()()()322332615,422f x f x f x x x x --===-- 那么二阶差商()123,,______f x x x =3. 设(2,3,1)TX =--, 那么2||||X = ,=∞||||X 。
4.求方程 21.250x x --= 的近似根,用迭代公式x =01x =, 那么 1______x =。
5.解初始值问题 00'(,)()y f x y y x y =⎧⎨=⎩近似解的梯形公式是 1______k y +≈。
求f (x)的Lagrange二次插值多项式和Newton二次插值多项式。
解:(1)由题可知插值基函数分别为故所求二次拉格朗日插值多项式为(2)—阶均差、二阶均差分别为均差表为故所求Newton二次插值多项式为例2、设f(x) x2 3x 2 , x [0,1],试求f(x)在[0,1]上关于(x)1, spa n 1 ,x 的最佳平方逼近多项式。
解:若span 1,x,贝卩o(x) 1, i(x) x,且(x) 1,这样,有所以,法方程为123 4 123■——i2 a o6,经过消兀得 2 a o61 1 ai9 1 6 U |1234123再回代解该方程,得到印4, ao-6 故,所求最佳平方逼近多项式为S; (x) 口4x6例3、设f(x) e x, x [0,1],试求f(x)在[0, 1 ]上关于(x)1 ,方逼近多项式。
解:若span 1,x,贝卩o(x) 1, i(x) x,这样,有所以,法方程为解法方程,得到ao 0.8732, a ( 1.6902,故,所求最佳平方逼近多项式为9例4、用n 4的复合梯形和复合辛普森公式计算积分.xdx。
1解:(1) 用n4的复合梯形公式由于h 2, f x xk 1 2k k(2) 用n4的复合辛普森公式由于h 2, f x x, xk 1 2k k1,2,3, Xi 2 2k k 0,1,2,3,所以,有k 一例5、用列主元消去法求解下列线性方程组的解。
span 1,x的最佳平1,2,3,所以,有解:先消元再回代,得到X3 3,X2 2, 1所以,线性方程组的解为XJ, X22, X33例6、用直接三角分解法求下列线性方程组的解。
解:则由A LU的对应元素相等,有1u ii u12 , u135 6因此,yi 解Ly b,即一 13 0 y29% 得yi 9, y2 4, y3 1544 解Uxy,即00 60xi177.69, x2 476.92,为227.08451315154'21 '31u11 131121u12 u22 u2260 121u13 u23 u2345I31U12 I32U22 132 36 '31u13 '32u23 u33 u33 13 15一成立O2、当n 8时, Newton — cotes 型求积公式会产生数值不稳定性。
一. 填空题(本大题共4小题,每小题4分,共16分)1.设有节点012,,x x x ,其对应的函数()y f x =的值分别为012,,y y y ,则二次拉格朗日插值基函数0()l x 为 。
2.设()2f x x =,则()f x 关于节点0120,1,3x x x ===的二阶向前差分为 。
3.设110111011A -⎡⎤⎢⎥=--⎢⎥⎢⎥-⎣⎦,233x ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,则1A = ,1x = 。
4. 1n +个节点的高斯求积公式的代数精确度为 。
二.简答题(本大题共3小题,每小题8分,共24分)1. 哪种线性方程组可用平方根法求解?为什么说平方根法计算稳定?2. 什么是不动点迭代法?()x ϕ满足什么条件才能保证不动点存在和不动点迭代序列收敛于()x ϕ的不动点?3. 设n 阶矩阵A 具有n 个特征值且满足123n λλλλ>≥≥≥,请简单说明求解矩阵A 的主特征值和特征向量的算法及流程。
三.求一个次数不高于3的多项式()3P x ,满足下列插值条件:i x 1 2 3 i y 2 4 12 i y '3并估计误差。
(10分)四.试用1,2,4n =的牛顿-科特斯求积公式计算定积分1011I dx x=+⎰。
(10分) 五.用Newton 法求()cos 0f x x x =-=的近似解。
(10分) 六.试用Doolittle 分解法求解方程组:12325610413191963630x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥----⎣⎦⎣⎦⎣⎦ (10分)七.请写出雅可比迭代法求解线性方程组123123123202324812231530x x x x x x x x x ++=⎧⎪++=⎨⎪-+=⎩ 的迭代格式,并判断其是否收敛?(10分)八.就初值问题0(0)y yy y λ'=⎧⎨=⎩考察欧拉显式格式的收敛性。
(10分)《数值分析》(A )卷标准答案(2009-2010-1)一. 填空题(每小题3分,共12分) 1. ()1200102()()()()x x x x l x x x x x --=--; 2.7;3. 3,8;4. 2n+1。
第一章 绪论姓名 学号 班级习题主要考察点:有效数字的计算、计算方法的比较选择、误差和误差限的计算。
1 若误差限为5105.0-⨯,那么近似数0.003400有几位有效数字?(有效数字的计算) 解:2*103400.0-⨯=x ,325*10211021---⨯=⨯≤-x x 故具有3位有效数字。
2 14159.3=π具有4位有效数字的近似值是多少?(有效数字的计算) 解:10314159.0⨯= π,欲使其近似值*π具有4位有效数字,必需41*1021-⨯≤-ππ,3*310211021--⨯+≤≤⨯-πππ,即14209.314109.3*≤≤π即取(3.14109 , 3.14209)之间的任意数,都具有4位有效数字。
3 已知2031.1=a ,978.0=b 是经过四舍五入后得到的近似值,问b a +,b a ⨯有几位有效数字?(有效数字的计算)解:3*1021-⨯≤-aa ,2*1021-⨯≤-b b ,而1811.2=+b a ,1766.1=⨯b a 2123****102110211021)()(---⨯≤⨯+⨯≤-+-≤+-+b b a a b a b a故b a +至少具有2位有效数字。
2123*****10210065.01022031.1102978.0)()(---⨯≤=⨯+⨯≤-+-≤-b b a a a b b a ab 故b a ⨯至少具有2位有效数字。
4 设0>x ,x 的相对误差为δ,求x ln 的误差和相对误差?(误差的计算) 解:已知δ=-**xx x ,则误差为 δ=-=-***ln ln xx x x x则相对误差为******ln ln 1ln ln ln xxx x xxx x δ=-=-5测得某圆柱体高度h 的值为cm h 20*=,底面半径r 的值为cm r 5*=,已知cm h h 2.0||*≤-,cm r r 1.0||*≤-,求圆柱体体积h r v2π=的绝对误差限与相对误差限。
北航2010-2015年研究生数值分析报告期末模拟试卷与真题数值分析模拟卷A一、填空(共30分,每空3分)1 设-=1511A ,则A 的谱半径=)(a ρ______,A 的条件数)(1A cond =________. 2 设 ,2,1,0,,53)(2==+=k kh x x x f k ,则],,[21++n n n x x x f =________, ],,[321+++n n n n x x x x f ,=________.3 设≤≤-++≤≤+=21,1210,)(2323x cx bx x x x x x S ,是以0,1,2为节点的三次样条函数,则b=________,c=________.4 设∞=0)]([k k x q 是区间[0,1]上权函数为x x =)(ρ的最高项系数为1的正交多项式族,其中1)(0=x q ,则?=10)(dx x xq k ________,=)(2x q ________.5 设=11001a a a a A ,当∈a ________时,必有分解式,其中L 为下三角阵,当其对角线元素)3,2,1(=i L ii 满足条件________时,这种分解是唯一的.二、(14分)设49,1,41,)(21023====x x x x x f , (1)试求)(x f 在]49,41[上的三次Hermite 插值多项式)(x H 使满足2,1,0),()(==i x f x H i i ,)()(11x f x H '='.(2)写出余项)()()(x H x f x R -=的表达式.三、(14分)设有解方程0cos 2312=+-x x 的迭代公式为n n x x cos 3241+=+,(1)证明R x ∈?0均有?∞→=x x n x lim (?x 为方程的根);(2)取40=x ,用此迭代法求方程根的近似值,误差不超过,列出各次迭代值;(3)此迭代的收敛阶是多少?证明你的结论.四、(16分) 试确定常数A ,B ,C 和,使得数值积分公式有尽可能高的代数精度. 试问所得的数值积分公式代数精度是多少?它是否为Gauss 型的?五、(15分)设有常微分方程的初值问题=='00)(),(y x y y x f y ,试用Taylor 展开原理构造形如)()(11011--++++=n n n n n f f h y y y ββα的方法,使其具有二阶精度,并推导其局部截断误差主项.六、(15分)已知方程组b Ax =,其中= ??=21,13.021b A ,(1)试讨论用Jacobi 迭代法和Gauss-Seidel 迭代法求解此方程组的收敛性.(2)若有迭代公式)()()()1(b Ax a x x k k k ++=+,试确定一个的取值围,在这个围任取一个值均能使该迭代公式收敛.七、(8分)方程组,其中,A 是对称的且非奇异.设A 有误差,则原方程组变化为,其中为解的误差向量,试证明 .其中1λ和2λ分别为A 的按模最大和最小的特征值.数值分析模拟卷B填空题(每空2分,共30分)1. 近似数231.0=*x 关于真值229.0=x 有____________位有效数字;2. 设)(x f 可微,求方程)(x f x =根的牛顿迭代格式是_______________________________________________;3. 对1)(3++=x x x f ,差商=]3,2,1,0[f _________________;=]4,3,2,1,0[f ________;4. 已知???? ??-='-=1223,)3,2(A x ,则=∞||||Ax ________________,=)(1A Cond ______________________ ;5. 用二分法求方程01)(3=-+=x x x f 在区间[0,1]的根,进行一步后根所在区间为_________,进行二步后根所在区间为_________________;6. 求解线性方程组=+=+04511532121x x x x 的高斯—赛德尔迭代格式为_______________________________________;该迭代格式迭代矩阵的谱半径=)(G ρ_______________;7. 为使两点数值求积公式:?-+≈111100)()()(x f x f dx x f ωω具有最高的代数精确度,其求积节点应为=0x _____ , =1x _____,==10ωω__________.8. 求积公式)]2()1([23)(30f f dx x f +≈?是否是插值型的__________,其代数精度为___________。
注:1、教师命题时题目之间不留空白; 2、考生不得在试题纸上答题,教师只批阅答题册正面部分,
8,2⎤⎦为( A )(C )、1;,)n 是n 个互异节点,)n 的拉格朗日
插值基函数,则下列选项中正确的是( C ); )、(0
n
i i i x l =∑
注:1、教师命题时题目之间不留空白;2、考生不得在试题纸上答题,教师只批阅答题册正面部分,
注:1、教师命题时题目之间不留空白;2、考生不得在试题纸上答题,教师只批阅答题册正面部分,
注:1、教师命题时题目之间不留空白;2、考生不得在试题纸上答题,教师只批阅答题册正面部分,
注:1、教师命题时题目之间不留空白;2、考生不得在试题纸上答题,教师只批阅答题册正面部分,。
数值分析复习试题第一章绪论一.填空题1.为精确值的近似值;为一元函数的近似值;*xx ()**x f y =()x f y =1为二元函数的近似值,请写出下面的公式::()**,*y x f y =()y x f y ,2=**e x x =-***r x xe x -=()()()*'1**y f x x εε≈⋅()()()()'***1**r r x f x y x f x εε≈⋅()()()()()**,**,*2**f x y f x y y x y x yεεε∂∂≈⋅+⋅∂∂()()()()()****,***,**222r f x y e x f x y e y y x y y y ε∂∂≈⋅+⋅∂∂2、计算方法实际计算时,对数据只能取有限位表示,这时所产生的误差叫舍入误差。
3、分别用2.718281,2.718282作数e 的近似值,则其有效数字分别有 6 位和7(三位有效数字)。
1.73≈-211.73 10 2-≤⨯4、设均具有3位有效数字,则的相对误差限为 0.0055 。
121.216, 3.654x x ==12x x 5、设均具有3位有效数字,则的误差限为 0.01 。
121.216, 3.654x x ==12x x +6、已知近似值是由真值经四舍五入得到,则相对误差限为0.0000204 .2.4560A x =T x 7、递推公式如果取作计算,则计算到时,误差为,⎧⎪⎨⎪⎩0n n-1y =y =10y -1,n =1,2,0 1.41y =≈10y ;这个计算公式数值稳定不稳定 不稳定 .8110 2⨯8、精确值,则近似值和分别有 3 位和14159265.3*=π141.3*1=π1415.3*2=π4 位有效数字。
9、若,则x 有 6 位有效数字,其绝对误差限为1/2*10-5 。
*2.71828x e x =≈=10、 设x*的相对误差为2%,求(x*)n 的相对误差0.02n11、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字;12、计算方法主要研究( 截断 )误差和( 舍入 )误差;13、为了使计算 的乘除法次数尽量地少,应将该表达式()()2334610111y x x x =++----改写为11,))64(3(10-=-++=x t t t t y ,为了减少舍入误差,应将表达式19992001-改写为199920012+。
北大考研详解与指导北京大学601高等数学考研分析和复习方法指导【介绍】北大考研科目里的“高数601”是针对理科部分专业设置的考试科目(环境科学、地理学、生态学等),主要考察内容为高等数学(微积分)(一般不包括三角级数、换流量、通量、方向导数,对于格林公式、斯托克斯、高斯公式考察也不多。
通过分析真题,育明教育考研专业课咨询师发现,北京大学601高数考查知识面并不很宽,但题目很有难度和深度)。
【题型】目前,从育明教育收集到的近10多年的真题来看,能搜集到的最早的试卷为1994年的版本,而后题型不断变化,从2004年开始趋于稳定,总分150分,包括8道填空题(64分)、4道单选题(24分)、4道解答题(62=16*3+14)。
2010年的分值分布稍有变化:填空8*8+选择7*4+大题14*3(不等式和函数性质、平面几何+函数极值、三角函数不等式)+16(微分方程+级数和函数)=150。
【难度】试题的难度和数一高数部分难度相当,但风格更加灵活多变,讲究方法和技巧,充分的复习可以保证及格,要想冲击高分则要广泛浏览,发散思路。
相对于数一、数三来说,虽然难度偏低一些,但是由于题型的灵活性高,因此,还是需要下一番功夫进行准备的。
【参考资料】官方指定的参考书目为《高等数学》(上、下册),樊映川主编,高等教育出版社,貌似已经买不到了。
育明教育为大家推荐以下参考资料,包括:1、同济六版《高等数学》(上、下册).................掌握基本知识点2、二李《考研数学复习全书(理工类、数一)》高等数学部分.......例题+练习掌握基本解题方法3、陈文灯《考研数学复习全书(数一)》高等数学部分............例题+练习熟悉常用技巧、扩充思路4、《北京大学高等数学复习指导》................巩固知识点、辨析基本概念5、高数601考研真题.....................熟悉题型和试题风格6、《育明教育考研专业课新攻略——北大成规》,含历年真题,重点笔记总结等。
北京大学数值分析试题2015北京大学2014--2015学年第一学期 研究生期末考试试题A (闭卷考试)课程名称:数值分析 注:计算题取小数点后四位 一、填空题(每空3分,共24分) (1) 设1222A ⎡⎤-=⎢⎥-⎢⎥⎣⎦,则A 的奇异值为 。
(2) 设0.00013753x =为真值0.00013759Tx=的近似值,则x有 位有效数字。
(3) 设数据123,,x x x 的绝对误差为0.002,那么123x x x -+的绝对误差约为 ____ _。
(4) )x (l ,),x (l ),x (ln 10是以01,,,,(2)n x x x n ≥为节点的拉格朗日插值基函数,则 20(2)()nkk k xl x =+=∑ 。
(5) 插值型求积公式22=≈∑⎰()()nk k k x f x dx A f x 的求积系数之和0nkk A ==∑ 。
其中2x 为权函数,1≥n 。
(6)已知(3,4),(0,1)TTx y ==,求Householder 阵H使Hx ky =,其中k R ∈。
H= 。
(7) 数值求积公式112()((0)3f x dx f f f -⎡⎤≈++⎢⎥⎣⎦⎰的代数精度为___。
(8) 下面Matlab 程序所求解的数学问题是 。
(输入向量x , 输出S )x =input('输入x :x ='); n=length(x ); S=x (1); for i=2:nif x (i)<S ,S=x (i);else,continue;end end S二、(12分) (1)证明对任何初值 0x R∈,由迭代公式124cos ,0,1,2,...3k k xx k +=+=所产生的序列{}0kk x ∞=都收敛于方程1232cos 0x x -+=的根。
(2)证明它具有线性收敛性。
三、(12分)(1)用辛浦生公式计算积分40x e dx⎰的近似值;(2)若用复化辛浦生公式计算积分40x e dx⎰,问至少应将区间[0,4]多少等分才能保证计算结果有五位有效数字?四、(12分) 已知数据表2102230.510.5i iix y w --(1)构造关于点集和权的正交函数组01{(),()}x x ϕϕ;(2)利用01{(),()}x x ϕϕ拟合已知数据点,并求最小二乘拟合误差2δ。
五、(12分) 利用Gauss 变换阵,求矩阵2113113112A ⎡⎤⎢⎥⎢⎥=⎢⎥-⎢⎥-⎣⎦的LU 分解。
(要求写出分解过程)六、(10分) 已知求解线性方程组Ax=b 的分量迭代格式1(1)()(1)()1,1,2,,i nk k k k iii ij jij j j j iii x x b a xa x i na ω-++===+--=∑∑()(1)试导出其矩阵迭代格式及迭代矩阵; (2)证明当A 是严格对角占优阵,1ω=时此迭代格式收敛。
七、(10分) 用插值极小化方法求 xe f -=)x ( 在[1,2]上的二次插值多项式)x (2P ,并在[1,2]上估计误差。
(已知Chebyshev 多项式)(t T 3的三个零点86600t 0t 86600t 210.,,.==-=)八、(8分)已知求解常微分方程初值问题00'()()()y x f x y y x y=+⎧⎨=⎩的数值格式为2100()'()[1()]2()n n n n n n n n h y y hf x y f x y f x y y x y +⎧⎪=++++++⎨⎪=⎩问此数值格式是几阶格式?北京大学 2014--2015 学年第 一 学期研究生期末考试试题标准答案A (闭卷考试)课程名称: 数值分析一、 填空题(每空3分,共24分)(1) 3 (2)3 (3)0.006 (4)22x+(5) 83 (6)4343--555534345555H H ⎡⎤⎡⎤⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦或(7)3(8)求向量x 的最小值二、(12分) 记2()4cos 3x x ϕ=+,则2'()sin 3x x ϕ=-。
(1)先考虑区间[3,5],当[3,5]x ∈时,2()4cos [3,5]3x x ϕ=+∈ ,22'()sin 133x x ϕ=-<< 。
故对任意初值0[3,5]x∈,由迭代公式124cos ,0,1,2,...3k k xx k +=+=产生的序列{}0kk x ∞= 都收敛于方程1232cos 0x x -+=的根。
(6分)(2)对任意初值0xR∈,有1024cos [3,5]3xx =+∈,将此1x 看成新的迭代初值,则由(1)可知,由迭代公式124cos ,0,1,2,...3k k xx k +=+=产生的序列{}0kk x ∞= 都收敛于方程 1232cos 0x x -+=的根。
(2分)(3)****1**11**22(cos cos )sin ()33222sin ,lim lim sin 1333k k k k k k x k kx x x x x x x x x x x x x x ξξξξ+++→∞→-=-=----=-=-≤<-- (4分)此格式线性收敛性 三、(12分)(1)402404(4) 56.10296x e dx e e e =++=⎰(5分)(2) (4)(),(),xx f x ef x e ==由54(4)(4)45434()4|()||()||()|28802880411028802n b a R S h f f n e n ηη--=-=-≤≤⨯ (5分)14.0371n >至少将区间[0,4] 15等分才能保证计算结果有五位有效数字. (2分)四、(12分)(1)首先构造关于点集和权的首一正交多项式(),0,1,ix i ϕ=显然0()1x ϕ=,设1()()x x a x ϕϕ=+,由1()()x x ϕϕ与正交得000((),)21((),())2x x a x x ϕϕϕ-=-=-=故有 1()1x x ϕ=+。
(4分)(2)设211()()()p x a x a x ϕϕ=+,则01010011((),)((),)9/291/21,((),())24((),())12x y x y a a x x x x ϕϕϕϕϕϕ======191()(1)42p x x ∴=++ (4分)2222000111||||((),())((),())Y a x x a x x δϕϕϕϕ=-- 2299116()2()10.1252428=+-⨯-⨯=≈ (4分)五、(12分)(2)11100021015100010,22001001310001012L L A A ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥-===⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥-⎢⎥⎢⎥⎣⎦⎣⎦(3分)(2)(3)221000210001005010,220100013/51500120001L L A A⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥===⎢⎥⎢⎥-⎢⎥⎢⎥-⎢⎥⎢⎥-⎢⎥⎢⎥⎣⎦⎣⎦(3分)(2)32100210001005010,200100013/51500100021/1313L L A U ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥===⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥--⎢⎥⎢⎥⎣⎦⎣⎦(3分)111122100011002,20105500113L L L L A LU ---⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥===⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦(3分)六、(10分)(1) 1(1)()(1)()1,1,2,,i n k k k k ii iii ii ij jij j j j ia x a xb a xa x i nω-++===+--=∑∑()(1)()(1)()(1)()(1)1()1())()((1))()((1))()k k k k i k k k k Dx Dx b Lx U D x D L x D U x bx D L D U x D L bωωωωωωωωωω++++--=+++--=-++=--++-((1)()11()((1))()k k x B x g B D L U D g D L bωωωωωωω+--=+=-+-=-迭代迭矩阵右端向式量代法的矩阵形 (6分)21A ω=()时,迭代格式为Gauss-seidel 迭代格式,当严格对角占优时,Gauss-seidel 迭代格式收敛。
(4分)七、(10分) 已知Chebyshev 多项式)(t T 3的三个零点86600t 0t 86600t210.,,.==-=,作变量代换)(x 3t 21+=,得三个插值节点210k 3t 21k k,,),(x=+=1.9330x 1.5x 1.0670x 210===,,0.1447x f 0.2231x f 0.3440f(x 210===)(,)(,)构造差商表()ii x f x 一阶差商 二阶差商1.06700.34401.50000.22310.27921.93300.14470.18110.1133--牛顿插值多项式22P (x)0.34400.2792(x 1.0670) 3.5863(x 1.0670)(x 1.5)0.11330.57010.8234x x =--+--=-+ ( 6分)001902216e t t t t t t max 216e x x x x x x 6f x R 2312101t -13121032.)())()(()())()(()()()(=⨯≤---≤---ξ=--≤≤-( 4分)八、(8分)()2121112233()'()[1()]2()'()''()(4)2()'()''()()()'()''()22()(4)n n n n n n n n n n n n n n n n n n n n h y y hf x y f x y f x y hy x hy x y x E y x y h h y x hy x y x O h y x hy x y x O h ++++=++++++=++=-⎛⎫⎛⎫=+++-++ ⎪ ⎪⎝⎭⎝⎭=分此格式二阶精度。
分。