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

1 第六章 物理学中非线性方程的求根问题 2 §6.1 根的搜索和二分法 高于四次的代数方程无精确求根公式。一般的超 越方程更无法求其精确解，只能求其近似解。 3 由于曲线最低点和最高点相差10m，有 求解，需数值求解。 4 例：水槽由半圆柱体水平放置而成，如下图所示。圆柱体长L，半径r，当给定水槽内盛水的体积V后，要求计算从水槽边沿到水面的距离H。 已知L=25m，r=2m，求V分别为10，50，100m3的H。 5 建立直角坐标系如图所示： 可得：x2 + (r-y)2 = r2 又因为：dV = 2·x·L·dy 所以： 积分可得： 注意：上述方程中含有参数r，L均为已知常数，通过上式，当V=10、50、100时，要求变量H，这显然是一个非线性方程。问如何求解？ 6 例：在我方前沿阵地1000米处有一座高为50米的山丘，山丘上建有敌方一座碉堡，求我方的大炮在什么角度下以最小的速度发射炮弹就能摧毁敌军的这座碉堡？ 解： 由抛体运动的轨道方程可得： 整理可得： 对上述方程求导以后得到的是一超越方程，可用数值方法求解。 7 例：静电除尘器由半径为ra=0.84m的金属圆筒（阳极）和半径为rb的同轴圆细线（阴极）组成。当它们加上一高电压V=50kv时，圆筒内就产生了一强大电场，圆筒内的空气被电离，混浊的空气通过这个圆筒时灰尘粒子与离子碰撞而带电，于是在电场的作用下奔向电极，并落下沉积在圆筒底部而被扫出，达到了清洁空气的目的。为了在中心轴线处产生6.0MV/m的电场强度而击穿空气，试求静电除尘器中心线的粗细。 解： 设静电除尘器中心细线上所带电荷 线密度为λ，则在距中心轴线r位置 处的电场强度为： 8 中心细线与金属圆筒的电势差为： 两式联立消去λ可得： 令r=rb可得在中心线表面处的电场强度Eb与其半径rb的关系： 将△V，Eb，ra的值代入上式，可得关于 rb的一超越方程，可用数值方法求解。 9 一、根的搜索 1．作图法 2．逐步搜索法 有根区间。 当 10 ∴ 方程在[-1, 2]内至少有一个实根。 11 若事先确定实根上、下界，关于上面方程根的绝对值的 上下界有以下结论： 12 13 二、二分法 14 以上过程，得到一系列有根区间： 前一个区间长度的一半， 15 二分法算法实现 16 编程步骤 17 求该方程在区间[0.5, 1.25]内根的近似值。 ∴[0.5, 1.25]为有根区间。 故所求根的近似值为： 二分法优点:算法简单，收敛性有保障。 缺点:收敛速度慢。 18 输入：1.0，5.0 输出结果： 2.528919219970703 19 解： 物理分析： 20 21 例6：水槽由半圆柱体水平放置而成，如下图所示。圆柱体长L，半径r，当给定水槽内盛水的体积V后，要求计算从水槽边沿到水面的距离H。 已知L=25m，r=2m，求V分别为10，50，100m3的H。 二分法求解本例题 由前推导得出的结果： 即： 将不同的V值代入上述方程，调用二分法程序计算： V= 10m3，H=1.71351242065430m V= 50m3，H=1.13508224487305m V=100m3，H=0.57899093627929m H必然在[0，2]区间内 输入x1=0，x2=2，可得： 22 作 业 23 §6.2 函数迭代法 迭代公式可能收敛， 可能发散。 24 从而得迭代公式： 25 相应的迭代公式为： 有多种，迭代公式可以有多个，如： 迭代发散! 26 反复递推得 27 定理 28 解：取小球为研究对象，根据已知条件有 该微分方程的解为 29 采用迭代法求解，将把方程化成 30 例9：在我方前沿阵地1000米处有一座高为50米的山丘，山丘上建有敌方一座碉堡，求我方的大炮在什么角度下以最小的速度发射炮弹就能摧毁敌军的这座碉堡？ 解： 由前面分析可知： 若将方程整理为： 31 此方程在迭代过程中，θ的取值范围(0, 90°)，而 1 30.000000000000000 2 46.240880181691520 3 46.431649042200080 4 46.431203734034090 32 加速算法－埃特金(Aitken)法 33 则有 不断进行下去，归纳有： 34 35 作 业 。 36 §6.3 牛顿迭代法 37 以上迭代方法－牛顿法（切线法） 38 由于牛顿法仍是迭代法，其迭代函数为 牛顿法收敛性分析 则牛顿法收敛。 上节定理可知，若 由 39 40 故牛顿法是收敛的。 凡满足关系式 结 论 41 定理 42 利用一阶泰勒公式及零阶泰勒公式有