数值分析偏微方程数值解法.pptVIP

阜师院数科院第十章偏微分方程数值解法10-*Poisson方程差分格式的建立（续1） Poisson方程（10-6）在点(k,j)处可表示为:(10-8)(10-7)其中：第29页，共43页，星期日，2025年，2月5日阜师院数科院第十章偏微分方程数值解法10-*在式（10-8）中略去R(k,j)即得与方程（10-6）相近似的差分方程:Poisson方程差分格式的建立（续2）式（10-9）为其截断误差表示式.（10-10）式（10-10）中方程的个数等于正则内点的个数，而未知数uk,j则除了包含正则内点处解u的近似值外，还包含一些非正则内点处u的近似值，因而方程个数少于未知数个数。在非正则内点处Poisson方程的差分近似不能按式（10-10）给出，需要利用边界条件得到。第30页，共43页，星期日，2025年，2月5日阜师院数科院第十章偏微分方程数值解法10-*Poisson方程边界条件的处理边界条件的处理可以有各种方案，下面介绍较简单的两种。（1）直接转移用最接近非正则内点的边界点上的u值作为该点上u值的近似，这就是边界条件的直接转移。如图10-2，点R(k,j)为非正则内点，其最接近的边界点为Q点，则有 （10-11）将式（10-11）代入式（10-10），方程个数即与未知数个数相等。式（10-11）可以看作是用零次插值得到非正则内点处u的近似值，容易求出，其截断误差为O(h+τ)。第31页，共43页，星期日，2025年，2月5日WY阜师院数科院第十章偏微分方程数值解法阜师院数科院第十章偏微分方程数值解法WY阜师院数科院第十章偏微分方程数值解法数值分析偏微方程数值解法阜师院数科院第十章偏微分方程数值解法10-*第1页，共43页，星期日，2025年，2月5日阜师院数科院第十章偏微分方程数值解法10-*第十章目录§1差分方法的基本概念1.1偏微分方程的定解问题1.2差分方法的基本概念§2椭圆型方程第一边值的差分方法2.1差分格式的建立2.2差分格式解的存在唯一性§3抛物型方程的差分解法及其稳定性3.1差分格式的建立3.2差分格式的稳定性§4双曲型方程的差分解法4.1几种简单的差分格式4.2差分格式的收敛性与稳定性第2页，共43页，星期日，2025年，2月5日阜师院数科院第十章偏微分方程数值解法10-*补充知识“高数”中接触了一些简单偏微分，也接触了简单偏微分方程，如： 其中： 1.2.满足：第3页，共43页，星期日，2025年，2月5日阜师院数科院第十章偏微分方程数值解法10-*补充知识（续1）3.2sin(x+2y-3z)=x+2y-3z满足：4.满足：5.满足：6.满足：上面是已知函数，，验证满足等式，反过来，将等式视为方程，则是求解方程，得到解函数。第4页，共43页，星期日，2025年，2月5日阜师院数科院第十章偏微分方程数值解法10-*因此偏微分方程：1.含偏微分的等式，2.求解偏微分方程、求含多个自变量的函数3.带有初值、边界条件。常微分方程的求解已很困难，通过分门别类研究，能求得一些特殊类型方程的解（只含一个变量），即便是一阶方程，也很难求出解析解表达式，也因此，在上一章我们研究了一阶微分方程的数值解法。补充知识（续2）第5页，共43页，星期日，2025年，2月5日阜师院数科院第十章偏微分方程数值解法10-*§1差分方法的基本概念要求解偏微分方程比求解常微分方程更难，因此寻求偏微分方程的数值解更显重要，实际上，绝大部分偏微分方程不可能求到解析函数解，基本上都是数值解法。一般来说，偏微分方程从实际问题抽出后，多是下列几种类型:（1）泊阿松方程（Poisson），又称为椭圆型方程：?：自变量的变化区域，有