《《计算方法例题分析》word版.doc》由会员分享，可在线阅读，更多相关《《计算方法例题分析》word版.doc（30页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、例题分析一例1设准确值x*= =31415926，当分别取近似值x=314和x=31416和x=31415时，求绝对误差、绝对误差限及有效数字位数。解：近似值x=3140.314101，即m=1，它的绝对误差是00015926，有x-x*=0.00159260.5101-3即n=3，故x=314有3位有效数字。x=314准确到小数点后第2位，又近似值x=3.1416，它的绝对误差是0.0000074，有x-x*=0.00000740.5101-5即m=1，n5，x=3.1416有5位有效数字。而近似值x=3.1415，它的绝对误差是0.0000926，有x-x*=0.00009260.5101
2、-4即m=1，n4，x=31415有4位有效数字。这就是说某数有s位数，若末位数字是四舍五入得到的，那么该数有s位有效数字；例2指出下列各数具有几位有效数字，及其绝对误差限和相对误差限：200040002009000900000解：因为x1=20004020004101，它的绝对误差限000005=05101-5，即m=1，n=5，故x=20004有5位有效数字。 a1=2，相对误差限；x2=000200，绝对误差限0000005，因为m=2，n=3，x2=000200有3位有效数字。 a1=2，相对误差限x3=9000，绝对误差限为05100，因为m=4, n=4，x3=9000有4位有效数
3、字，a=9，相对误差限x4=900000，绝对误差限0005，因为m=4，n=6，x4=900000有6位有效数字，相对误差限为由x3与x4可以看到小数点之后的0，不是可有可无的，它是有实际意义的。例3ln2=069314718，精确到103的近似值是多少？解：精确到10-30001，意旨两个近似值x1，x2满足x1-x20001，由于近似值都是四舍五入得到的，要求满足x1-x20001，近似值的绝对误差限应是00005，故至少要保留小数点后三位才可以。故In20693。 例4试利用f(x)的数据表计算积分，并估计计算误差分析在f(x)的表达式不知道的情况下，如何去求f(x)的积分值呢?若利用
4、本章的知识，即可利用已知的f(x)的数据表构造f(x)的二次插值多项式p2(x)，以作为f(x)的近似函数，并进而以p2(x)的积分值作为所求积分值的近似。至于误差的计算，也可由误差f(x)-p2(x)出发进行估计。解：根据拉格朗日插值公式，利用给定的数据表，可构造出f(x)的二次插值多项式插值余项为由此得积分近似值积分值的误差为其中例5给定f(x)在节点a，b上的函数值与导数值f(a)，f(b)，f(a)。试求一个二次多值式H2(x)，使之满足插值条件H2(x)=f(a)，H2(x)=f(b)，(1)分析 构造插值多项式的基本方法是基函数法，即对每一个插值条件建立一个与之相应的插值基函数。基
5、函数的形式要与所求的插值函数相一致。然后用给定的插值数据与基函数作线性组合，就可得到所求的插值函数。解：法一 与(1)中三个插值条件相应，依次建立三个插值基函数，是二次多项式且满足标准的基函数插值条件利用待定系数法容易求得则所求的二次插值多项式为法二可先根据给定条件H2(x)=f(a)，H2(b)=f(b)作出牛顿插值(或拉格朗日插值)多项式，然后再加带有待定系数的一项，所加项自然应保证在a，b处取值为零，故而可取k(x-a)(x-b)，再由条件确定待定系数k。设H2(x)=f(a)+fa，b(x-a)+k(x-a)(x-b)。于是所以注由于二次多项式由H2(a)，f(b)，f(a)三个条件所
6、唯一确定，所以本题由各种方法所求得的解，实质上是相同的。例题分析二例6已知函数y=f(x)的观察数据为试构造f(x)的拉格朗日多项式Pn(x)，并计算f(1)。解：先构造基函数所求三次多项式为：例7已知函数y=f(x)的数据如表中第2，3列。计算它的各阶均差。解：依据均差计算公式，结果列表中。计算公式为：一阶均差二阶均差例8设x0，x1，x2，xn是n+1个互异的插值节点，lk(x)(k=0，1，2，n)是拉格朗日插值基函数，证明：证明当f(x)=1时， 由于，故有.例9已知数据表如下：用最小二乘法求拟合这组数据的曲线。分析首先根据已知数据，在坐标平面上画出相应的点，然后再画出曲线的粗略图形。
7、如图31。由图形确定拟合函数的类型。在具体问题中也可结合考虑问题的物理意义和经验。最后由最4乘法建立法方程组，求出待定的参数，即可得拟合曲线的方程。并可比较拟合值、实验值算出各点的误差图31解：根据图31，取幂函数y=axb作拟合函数，其中a，b 为待定参数。根据曲线拟合的思想，令由 (a，b)求出a，b由极值的必要条件得方程组这是关于a、b的非线性方程组，求解很困难。于是，将问题转化为线性问题求解。为此，将y=axb两边取对数有lgy=lga+blgx.令=lgy，z=lgx，c=lga。上式化为=c+bz由(xi，yi)可得到相应的(zi，i)，于是得如下数据表：这样，待定系数c，b即为内
8、容提要中所述的和的线性组合系数。建立c，b所满足的法方程组其中由方程组(1)解得C=01624，b=20150。从而a=10c=14534，Y=14534，x2.0150。比较拟合值、实验值并算出各点的误差如下表注：通常针对一组数据的图形，可以选择不同的拟合函数类进行求解，最后按误差大小决定取舍。例题分析三例10满足条件p(0)=p(0)=0，p(1)=1，p(2)=2的插值多项式p(x)=_解：设所求的为p(x)=a0+a1x+a2x2+a3x3由插值条件知解之得a2=3/2a3=1/2所求的插值多项式为p(x)=1/2x3+3/2x2例11选择填空题1.通过四个互异节点的插值多项式P(x)
9、，只要满足()，则P(x)是不超过一次的多项式。(A)初始值y0=0(B)一阶均差为0(C)二阶均差为0(D)三阶均差为0解答：因为二阶均差为0，那么牛顿插值多项式为N(x)=f(x0)+f(x0，x1)(x-x0)它是不超过一次的多项式。故选择(C)正确。2. 拉格朗日插值多项式的余项是()，牛顿插值多项式的余项是()(A) (B) f(x，x0，x1，x2，xn)(xx1)(xx2)(xxn1)(xxn)(C)(D)f(x，x0，x1，x2，xn)(xx0)(xx1)(xx2)(xxn1)(xxn)解答：(A)，(D)。例12证明方程1xsinx0在区间0，1内有一个根，使用二分法求误差不
10、超过0510-4的根要迭代多少次？证明令f(x)1xsinx，f(0)=10，f(1)=sin10(x0.1)，故f(x)0在区间0，1内有唯一实根。给定误差限0510-4，有只要取n14。例13证明：方程f(x)=x3+4x2-10=0在1，2内有一个根，并用对分法求此根；若要求误差xn-x*=10-5，估计至少需要对分多少次?分析若连续函数f(x)满足f(a)f(b)0，则取a2=a1，b2=x1，否则取a2=a1，b2=x1。继续下一步，计算，这样得到为要使xn-x*，只需有。解：易见f(x)在1，2上连续，f(1)=-5，f(2)=14，且f(x)=3x2+8x0，x1，2，故f(x)
11、=0在1，2内有唯一根。为使误差xn-x*10-5，只需，即166，所以只需对分17次就能达到给定的精度。具体计算结果列于表1。注x12=1364990234，而x*=136523001。例14用迭代法求方程x54x20的最小正根。计算过程保留4位小数。分析 容易判断1，2是方程的有根区间。若建立迭代格式，即，此时迭代发散。建立迭代格式，此时迭代收敛。解：建立迭代格式，取初始值x0=1取x*1.5185例题分析四例15用弦截法求方程x3x210，在x=1.5附近的根。计算中保留5位小数点。分析 先确定有根区间。再代公式。解：f(x)=x3x21，f(1)=1，f(2)=3，有根区间取1，2。取
12、x1=1，迭代公式为取x*1.46553，f(1.46553)0.000145例16选择填空题1. 设函数f(x)在区间a，b上连续，若满足_，则方程f(x)=0在区间a，b一定有实根。答案：f(a)f(b)0解答：因为f(x)在区间a，b上连续，在两端点函数值异号，由连续函数的介值定理，必存在c，使得f(c)=0，故f(x)=0一定有根.2. 用简单迭代法求方程f(x)=0的实根，把方程(x)=0表成x=(x)，则f(x)=0的根是()(A)y=x与y=(x)的交点(B) y=x与y=(x)交点的横坐标(C)y=x与x轴的交点的横坐标(D) y=(x)与x轴交点的横坐标答案：(B)解答：把f
13、(x)=0表成x=(x)，满足x=(x)的x是方程的解，它正是y=x与y=(x)的交点的横坐标。3.为求方程x3x21=0在区间1.3，1.6内的一个根，把方程改写成下列形式，并建立相应的迭代公式，迭代公式不收敛的是()(A) (B) (C) (D) 答案：(A)解答：在(A)中，故迭代发散。在(B)中，故迭代收敛。在(C)中，故迭代收敛。在(D)中，类似证明，迭代收敛。例17设x*为方程x=g(x)的根，g(x)在x*附近连续，且g(x*)0，使对任意x0x*-，x*+，迭代格式xn+1=g(xn)(n=0，1，)收敛于x*。分析根据迭代法收敛条件，只须证明存在0，使g(x)在x*-，x*+
14、上满足；x*-g(x)x*+及g(x)L1。证因g(x*)1，由g(x)在x*附近的连续性知，存在0L0，使得当xx*-，x*+时有g(x)L1)，讨论用牛顿法求解方程f(x)=0时，迭代格式是否收敛?分析易见x*=0为，f(x)=0的根，且不是单根(因f(x*)=0)。但牛顿法收敛条件中的f(x*)0只是充分条件之一，不是必要条件，因此应从收敛定义具体讨论。解法f(x)=xa，f(x)=axa-1相应的牛顿迭代公式为对任意初值x0，即此格式收敛于f(x)=0的根。又所以此格式为线性收敛。这表明，当f(x*)=0时，牛顿法仍可收敛，但未必为二阶收敛。法二可设，在x*=0附近，由一般迭代法的收敛
15、条件，显然此格式是收敛的。例19证明：迭代公式是计算的三阶方法(ao)，这里初值x0是充分靠近x*=的。分析首先说明f(x)=x2-a=0和方程x=g(x)=x(x2+3a)/(3x2+a)解之间的关系，再验证迭代法收敛条件是否满足，以及利用收敛阶的概念证明是三阶的。证 令f(x)=x2-a=0， (1)(2)显然，方程(1)和方程x=g(x)在不包含x=0点的区间上是同解的。又由g(x)在x*=附近连续，且 (3)知g(x*)=g()，满足g(x*)6(B)=6(C)6解答：当a6时，线性方程组的系数矩阵是严格对角占优矩阵，由教材第3章定理知，迭代解一定收敛。应选择(A)。例23分别用高斯消
16、去法、列主元素消去法和直接三角分解法，解线牲方程组(要求在计算中取四位有效数字，最后结果舍入成三位有效数字)。解(1)高斯消去法解得x=(-0.399，-0.0998，0.400)T(2)高斯列主元素消去法解得x=(0.490，0.0513，0.368)T(3)直接三角分解法由Ly=b，解得y=(1.001，1002，2006)T再由Ux=y,得x=(0.400，0.0998，0.400)T注此题用四位浮点数进行计算，精确解舍入到三位有效数字为x=(0.490，0.0510，0.368)T可见列主元素消去法精度比较高。例24分别用雅可比迭代法、超松弛法(=0.8，1，1.2)求线性方程组的前三
17、次迭代解(已知(x0)=(0，0，0)T，在运算中要求取五位有效数字)。分析为了使所求解的线性方程组能够使用迭代法(要求aii0，i=1，2，n)和所使用的迭代法收敛得更快，通常都在使用迭代法之前对所给的线性方程组作一定的预处理，使得所得到的同解方程组的迭代矩阵的范数或谱半径变得比较小。一般应使所得同解方程组系数阵的对角元素在同行元素中绝对值最大(最好是使其系数阵为对角占优阵)。为此，须把第一、三两个方程对调，即得到与原方程组同解的线性方程组，再对此方程组使用迭代法求解。解：(1)雅可比迭代法，迭代公式为计算结果为(2)SOR法，迭代公式为依次取=0.8，1，1.2；其计算结果如下注(1)可以
18、证明，如果所给方程组不作如上的预处理，则解原方程组的雅可比迭代法、高斯一赛德尔迭代法均不收敛。(2)对大型(高阶)线性方程组，使用某种迭代法是否收敛，在迭代过程中就能发现，一般都不作事先论证。但是适当的预处理是必要的。(3)松弛因子的值对解具体方程组的SOR法是否收敛或收敛快慢影响很大，在实际计算中，通常是选择几个不同的值进行试算，比较后选出适合具体题目的较好的松弛因子。(4)=1时的SOR法就是高斯一赛德尔迭代法本题的精确解为x*=(3，2，1)T。例题分析七例25用迭代改善法解病态方程组分析由于方程组的系数阵的两行几乎成线性关系，因此本方程组可能是病态的。事实上迭代改善法的作法是，先用某种
19、方法(最好是三角分解法)求得Ax=b的近似解x(0)，然后按下面公式进行迭代改善：其中r(k)要用双字长进行计算。解：由及Ly=b，Ux=y，得y=(34.97，0.0181)T，x(0)=(1.949，3.051)T，及r(0)=(0.00051，0)T.由Ly(0)=r(0)及Uz(0)=y(0)，得y(0)=(0.00051，-0.0002914)，z(0)=(0.0493，-0.0494)T，x(1)=x(0)+z(0)=(1.998，3.002)T，及r(1)=(0.0002，0)T.由Ly(1)=r(1)及Uz(1)=y(1)得y(1)=(0.00002，-0.00001143)T
20、，z(1)=(0.00193，-0.00194)T，x(2)=x(1)+z(1)=(1.9999，3.0001)T.本题精确解为x*=(2，3)T.注迭代改善法不仅适用于病态方程组的求解，也可以改善非病态方程组的解。例26试确定求积公式的代数精度。解：当f(x)取1，x，x2，计算求积公式何时精确成立。(1) 取f(x)=1，有：左边，右边2(2) 取f(x)=x，有：左边，右边0(3)类似导出，取f(x)=x2，x3，有左边=右边(4) 取f(x)=x4，有：左边=2/5，右边=2/9当k3求积公式精确成立，而x4公式不成立，可见该求积公式具有3次代数精度。例27确定求积公式(1)中的系数A
21、-1，A0，A1，使其代数精度尽可能高，并指出求积公式的代数精度。分析按照使求积公式具有较高的代数精度的原则来确定求积系数，是构造求积公式的一种方法。本题中有三个系数需要确定，因此可建立三个方程使其满足。根据代数精度的定义，可设求积公式(1)对，f(x)=1，x，x2精确成立，从而建立关于A-1，A0，A1的三个方程，并求解。解：令求积公式(1)对，f(x)=1，x，x2精确成立，得解得此时求积公式(1)为(2)通过直接验证，求积公式(2)对f(x)=x3也精确成立，而对f(x)=x4不能精确成立，所以它的代数精度为3。 例28试用梯形公式、科茨公式和辛卜生公式计算定积分(计算结果取5位有效数
22、字)(1)用梯形公式计算(2)用柯特斯公式系数为(3)如果要求精确到10-5，用复化辛卜生公式，截断误差为只需把0.5，14等分，分点为0.5，0.625，0.75，0.875，1例题分析八例29用三点高斯勒让德求积公式计算积分解：做变量替换，有。查表得节点为0.774 596 669 和0；系数分别为0.555 555 5556和0.888 888 8889例30已知函数值f(1.0)=0.250 000，f(1.1)=0.226757，f(1.2)=0.206 612，用三点公式计算f(x)在x=1.0，1.1，1.2处的导数值。解：三点导数公式为k=1，2，3，n1本例取x0=1.0，x
23、1=1.1，x2=1.2，y0=0.250 000，y1=0.226757，y2=0.206 612，h=0.1。于是有例31选择填空题1. 如果用复化梯形公式计算定积分,要求截断误差不超过0.510-4，试问n()(A)41(B)42(C)43(D)40解答：复化的梯形公式的截断误差为，n=40.8,取n41。故选择(A)。2.已知n=3时，柯特斯系数，那么_解答：由柯特斯系数的归一性质，.例32用尤拉法解初值问题，取步长h=0.2。计算过程保留6位小数。解：h=0.2，f(x)=-y-xy2。首先建立尤拉迭代格式当k=0，x1=0.2时，已知x0=0，y0=1，有y(0.2)y1=0.21
24、(401)0.8当k1，x2=0.4时，已知x1=0.2，y1=0.8，有y(0.4)y2=0.20.8(40.20.8)0.6144当k=2，x3=0.6时，已知x2=0.4，y2=0.6144，有y(0.6)y3=0.20.6144(40.40.4613)=0.8例33用尤拉预报校正公式求解初值问题，取步长h=0.2，计算y(0.2)，y(0.4)的近似值，小数点后至少保留5位。解：步长h=0.2，此时f(x,y)=-y-y2sinx尤拉预报校正公式为：有迭代公式：当k=0，x0=1，y0=1时，x1=1.2，有当k=1，x1=1.2，y1=0.71549时，x2=1.4，有y(1.4)y
25、2=0.71549(0.9-0.10.71549sin 1.2)-0.1(0.47697+0.476972sin 1.4)=0.52608例题分析九例34写出用四阶龙格库塔法求解初值问题的计算公式，取步长h=0.2计算y(0.4)的近似值。至少保留四位小数。解：此处f(x，y)=83y，四阶龙格库塔法公式为其中k1=f(xk,yk)；k2=f(xn+h，yk+hk1)；k3=f(xk+h，yn+hk2)；k4=f(xk+h，yk+hk3)本例计算公式为：其中k1=83yk；k2=5.62.1yk；k3=6.322.37yk；k4=4.2081.578yk当x0=0，y0=2，y(0.2)y1=
26、1.2016+0.5494y0=1.2016+0.54942=2.3004y(0.4)y2=1.2016+0.5494y1=1.2016+0.54942.3004=2.4654.例35证明：改进欧拉方法 (1) 可精确求解常微分方程初值问题y=at+b，y(t0)=y0分析此题除可直接求解验证外，也可以从此方法的局部截断误差来证明。用归纳法，在假设yn=y(tn)的前提下去证明yn+1=y(tn+1)。利用泰勒展开式，将y(tn+1)和(1)式中的yn+1进行展开，从而可得误差y(tn+1)-yn+1=0证y(t0)=y0，并假设y(yn)=yn而以代入，得即y(tn+1)-yn+1=0.例36用欧拉折线法求解初值问题 (1) 问当步长取什么值时，此方法是稳定的?分析通常用典型问题y=y来考察方法的稳定性。在欧拉法中，因f(t，y)=y，故yn+1=yn+hyn=(1+h)yn因此，当要求前一步的误差在传给后一步时不扩大，那么必须且只须1+h1(2)由此给出步长h的取值范围。解：根据分析，本题中欧拉公式(3)当且仅当，即时，欧拉公式稳定，由于在所考虑的范围0.1t3，01，所以的取值范围应为0h0.2.