一维扩散方程的有限差分方法.pdfVIP

2010年第35期 O高校讲坛。 科技信息 一维扩散方程的有限差分方法 杨晓君潘艳冯建军 (云南师范大学云南昆明650092 J 【摘要】本文用有限差分方法求解一雏扩散方程。通过构造两种成分格式讨论其收敛性。和稳定性。利用Matalb程序做差分解和误差的 图像 【关键词】向前差分格式；向后差分格式；收敛性：稳定性 本文将研究以下一维扩散方程初边值问题 令A=；； I警≈嚣“如”“o，”“o．叫· K 1u(x，O)=uo(x)=sin(仃x)．￡Ee0，1】， (1．I) ‘ (1．10) lu(O，I)2妒(I)=l-cost，tE[0。∞]， 差分格式改为 Iu(1，1)却(t)爿，IE[0，∞]。 吩=(1-2aA)吩+珏^(％I+I)， (1．11) 其中Ⅱ是正常数． 本文使用的基本方法是通过构造向前和向后差分格式。并讨论差 令‘=茸一u(丐，tI)，由(1．Io)，(1．11)可得 分格式的性质．借助Matalb软件做出差分解。解的误差l墨I像． 々=(1面)气诎(％l+分1)呵r(丐，‘I)， 若拟≤1．那么 1讨论向前差分格式 IFlI≤(1，d)IfI“认IfI“r(丐，k)。 (1．12) 第一步：建立向前差分格式． 假设“扛，t)是(1．1)充分光滑。 威断误差可知，I7．(丐，l-)l 在对(1．1)建立差分网格时．取h=l∥，J为一个正整数，作出网格线 2 ≤J|If(．r“)． x=jh，(j=O。1，…，．，)，忙盯(0≤n≤[所])．设解在所考虑的区域R：{0≤ 茗《1，O≤l≤7-}，内部格点为(『^。盯)，(1吲≤．『-1，1≤，I≤[T／r]一I)．那 令E=supIfI，代入(1．22)得翌 f1≤(1一拟)E诎(E城)卅)《E撕(枞2)，(1．13) 么r(¨-)-坐血掣呻亟出丝争业亟世， 于是有 (1．2) E一．1《E+朋吁(件_11)， (1．14) 其中，Ir(丐，气)l≤(手+各)m肛I咛I． 由递推式。可得 对于区域R的边界格点．由所给定解条件。有 E≤‰+材耵(m)， (1．15) 吩=uju驰)；q=O，1．…，刀 (1．3) 在初始时间层Io上，有‘=l‘q i=o，我们有岛～pJ-