天津大学化工数学课件--第四章边值问题.pptVIP

打靶法的关键点 如何将边值问题转化为初值问题 解非线性方程，寻找合适的m值使初值问题与边值问题等价 其次对方程在xi处离散化处理： 用差商近似导数： 整理得： 其中： 第二步 用追赶法求解三对角线方程组： 总 结 1、 1、离散化,确定节点及差分方程 2、用追赶法解差分方程 有限差分法的关键点 例2： 解： 取 h=0.2 N+1=(1-0)/0.2＝5 N=4, 有四个内节点，即 设 ui= yi，求 整理得： 将 代入 得 * 4.3 常微分方程边值问题 1、第一类边界条件 4.3.1 边界条件及边界问题的分类 2、第二类边界条件 3、第三类边界条件 r = a r = b 建立微分方程： 建立边界条件： 例1 ：传热问题 恒温 ●第一类边界条件 －给定边界处函数满足的条件 ●第一类边值问题 例 2：传热问题 r = b 建立微分方程： 建立边界条件： 绝热 ●第二类边界条件 －给定边界处导函数满足的条件 ●第二类边值问题 例 3：传热问题 r = b 建立微分方程： 建立边界条件： 对流传热 ●第三类边界条件 －给定边界处函数和导数共同满足的条件 ●第三类边值问题 基本思想：将边值问题转化为等价的初值问题来求解。 一、 打靶法 4.3.2 求解方法 对于二阶常微分方程第一类边值问题 （Ⅱ） （Ⅰ） 解等价吗？ ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● y(x) ¤ 打靶法的几何说明 x =a y(a)=α ● x =b y(b)=β ● 对于初值问题 合适的 m 值应满足： 即： 用弦截法求解非线性方程： ●打靶法的求解步骤： 边值问题 初值问题 y (b)m 求 解 寻找m值 YES 得到边值 问题的解 NO 例题4：试用打靶法求解二阶边值问题,ε=10 - 4 解：首先将边值问题化为初值问题， （ⅰ） 设 （ⅱ） 取h = 0.1,用四阶龙格－库塔法求解 解（ⅱ） 解（ⅱ） 再解初值问题（ⅱ） X Y1 Y2 0.10000 1.03958 10.78717 9.70000 7.79941 -1.4477E+001 9.80000 6.38256 -1.3854E+001 9.90000 5.02974 -1.3197E+001 10.00000 3.74433 -1.2506E+001 ● ● ● ● 返回 ● ● ● ● ● ● ● ● 四阶龙格-库塔法求解一阶常微分方程组的初值问题 -------------------------------------------- 四阶龙格-库塔法求解一阶常微分方程组的初值问题 --------------------------------------------- X Y1 Y2 0.10000 1.234731 2.73708 9.60000 5.38702 -1.4893E+001 9.70000 3.93162 -1.4209E+001 9.80000 2.54638 -1.3490E+001 9.90000 1.23475 -1.2737E+001 10.00000 0.00002 -1.1952E+001 y(10)= ● ● ● ● ● ● ● ● ● ● ● ● ε=10-4 思考:如何求解第二及第三类边值问题 转化的初值问题为： 合适的 m 值应满足： 即： 用弦截法求解非线性方程： 转化的初值问题为： 合适的 m 值应满足： 即： 用弦截法求解非线性方程：