FORTRAN数值方法及其在物理学中应用.pptxVIP

会计学1FORTRAN数值方法及其在物理学中应用 2§6.1 根的搜索和二分法 高于四次的代数方程无精确求根公式。一般的超的解 称为方程的的零点） 根（或称函数 设 是代数多项式或超越函数，则 称为代数方程或超越方程。 越方程更无法求其精确解，只能求其近似解。第1页/共65页 3由于曲线最低点和最高点相差10m，有这是关于未知数的非线性方程，用解析的方法难以求解，需数值求解。即的参数 。例：在相距100m的两个塔（高度相等的点）上悬挂一根 电缆，允许电缆中间下垂10m。要计算两个塔之间 所用电缆的长度，需确定悬链线方程第2页/共65页 4例：水槽由半圆柱体水平放置而成，如下图所示。圆柱体长L，半径r，当给定水槽内盛水的体积V后，要求计算从水槽边沿到水面的距离H。已知L=25m，r=2m，求V分别为10，50，100m3的H。第3页/共65页 5建立直角坐标系如图所示：可得：x2 + (r-y)2 = r2又因为：dV = 2·x·L·dy所以：积分可得： 注意：上述方程中含有参数r，L均为已知常数，通过上式，当V=10、50、100时，要求变量H，这显然是一个非线性方程。问如何求解？Oyxx rHdyy第4页/共65页 6例：在我方前沿阵地1000米处有一座高为50米的山丘，山丘上建有敌方一座碉堡，求我方的大炮在什么角度下以最小的速度发射炮弹就能摧毁敌军的这座碉堡？解：由抛体运动的轨道方程可得：整理可得：求v0的极值，只需要计算并将x=1000，y=50代入即可。对上述方程求导以后得到的是一超越方程，可用数值方法求解。第5页/共65页 7例：静电除尘器由半径为ra的金属圆筒（阳极）和半径为rb的同轴圆细线（阴极）组成。当它们加上一高电压V=50kv时，圆筒内就产生了一强大电场，圆筒内的空气被电离，混浊的空气通过这个圆筒时灰尘粒子与离子碰撞而带电，于是在电场的作用下奔向电极，并落下沉积在圆筒底部而被扫出，达到了清洁空气的目的。为了在中心轴线处产生的电场强度而击穿空气，试求静电除尘器中心线的粗细。解：设静电除尘器中心细线上所带电荷线密度为λ，则在距中心轴线r位置处的电场强度为：第6页/共65页 8中心细线与金属圆筒的电势差为：两式联立消去λ可得：令r=rb可得在中心线表面处的电场强度Eb与其半径rb的关系：将△V，Eb，ra的值代入上式，可得关于rb的一超越方程，可用数值方法求解。第7页/共65页 9一、根的搜索1．作图法在 内至少有一实根，区间 称为有根区间。 2．逐步搜索法适当选择某一区间 ，从 出发，按事先选好的步长，逐点计算 的函数值 ，与 的值异号时，则 就是方程 的一个有根区间。对 若在 上连续， ，则 ],[ba?当第8页/共65页 10例1：找出方程： 的有根区间。 ∵ 解：设 取 ∴ 方程在[-1, 2]内至少有一个实根。 取 ，步长 ，从 出发，有： 在[-1，-0.25]，[0.5，1.25]，[1.25，2]内各有一个根，逐步搜索，便于计算机实现。第9页/共65页 11定 理对 次代数方程 若事先确定实根上、下界，关于上面方程根的绝对值的(1) 若 ，则方程根的绝对值小于 (2) 若 ，则方程根的绝对值大于 上下界有以下结论：即： 第10页/共65页 12例2：求方程 的根的模的上下界。 解：∵ ∴ 根据定理，上界为： ∴下界为： ∴即 第11页/共65页 13二、二分法 (2) 若 ，则有根区间为之间。 取 计算 。 中点设 为连续函数， 的有根区间为 (1) 若 ，则 就是方程 的根。 0)(=xf(3) 若 ，则有根区间为之间。 第12页/共65页 14所以 的长度为： 以上过程，得到一系列有根区间：将新的有根区间记为 ，再将二分，重复当 时，区间 的长度为必要趋于零， 若事先给定精度要求为 ： ，即可停止计算。 误差： 足够大，就有 实际计算时，只要二分的次数 后一个区间都在前一个区间内，且后一个区间的长度均是?前一个区间长度的一半，第13页/共65页 15 二分法算法实现第14页/共65页 16 编程步骤(1) 计算 (2) 计算 (3) 若 ，则 ，否则 (i) ，则根在 内，以 替代 (4) 若 ，计算终止，此时 ，否则转向(2)。 (ii) ，则根在 内，以 替代 第15页/共65页 17例3：设 ，试用二分法求该方程在区间[0.5, 1.25]内根的近似值。解：∵ ∴[0.5, 1.25]为有根区间。 取 起进行计算: 故所求根的近似值为：所产生的误差为： 二分法优点:算法简单，收敛性有保障。缺点:收敛速度慢。第16页/共65页 18例4：用二分法编程求方程在 附近的一个实根。(要求