椭圆型方程的有限差分法幻灯片.ppt

椭圆型方程的有限差分法第1页，共77页，编辑于2022年，星期六11差分逼近的基本概念差分逼近的基本概念第2页，共77页，编辑于2022年，星期六区间的剖分剖分1 1 区间的剖分区间的剖分第3页，共77页，编辑于2022年，星期六1 微分方程离散微分方程离散(差分方程）差分方程）1 1 微分方程离散微分方程离散(差分方程）差分方程）第4页，共77页，编辑于2022年，星期六 第5页，共77页，编辑于2022年，星期六 第6页，共77页，编辑于2022年，星期六 第7页，共77页，编辑于2022年，星期六定义定义1.11.1第8页，共77页，编辑于2022年，星期六定义定义1.21.2第9页，共77页，编辑于2022年，星期六第10页，共77页，编辑于2022年，星期六定义定义1.31.3第11页，共77页，编辑于2022年，星期六定理定理1.11.1（相容（相容+稳定稳定=收敛）收敛）第12页，共77页，编辑于2022年，星期六2 2 一维差分格式一维差分格式第13页，共77页，编辑于2022年，星期六2.12.1直接差分化直接差分化第14页，共77页，编辑于2022年，星期六ab图1第15页，共77页，编辑于2022年，星期六第16页，共77页，编辑于2022年，星期六第17页，共77页，编辑于2022年，星期六第18页，共77页，编辑于2022年，星期六第19页，共77页，编辑于2022年，星期六2.2 2.2 积分插值法积分插值法第20页，共77页，编辑于2022年，星期六 第21页，共77页，编辑于2022年，星期六 第22页，共77页，编辑于2022年，星期六 第23页，共77页，编辑于2022年，星期六2.3 2.3 变分变分-差分法差分法数值计算中，数值计算中，我们学习过我们学习过Lagrange插值多项式插值多项式公式：公式：第24页，共77页，编辑于2022年，星期六Lagrange插值多项式插值多项式 先从最简单的线性插值先从最简单的线性插值(n=1)开始。这时插值问题就是求一次多项式开始。这时插值问题就是求一次多项式P1(x)=a0+a1x 使它满足条件使它满足条件P1(x0)=y0,P1(x1)=y1,令令P1(x)=l0(x)y0+l1(x)y1,由于由于l0(x0)=1,l0(x1)=0,l0(x0)=0,l1(x1)=1.这样这样l0(x)含有因子含有因子x-x1,令令 l0(x)=(x-x1),再利用再利用 l0(x0)=1确定其中的系数，确定其中的系数，结果得到结果得到x-x1 l0(x)=-,x0-x1类似的可得到类似的可得到 x-x0 l1(x)=-,x1-x0这样这样x-x1 x-x0P1(x)=-y0+-y1,x0-x1 x1-x0 l0(x),l1(x)称为以称为以x0,x1 为节点的为节点的插值基函数插值基函数。第25页，共77页，编辑于2022年，星期六 第26页，共77页，编辑于2022年，星期六 第27页，共77页，编辑于2022年，星期六2.4 2.4 边值条件的处理边值条件的处理第28页，共77页，编辑于2022年，星期六 第29页，共77页，编辑于2022年，星期六第30页，共77页，编辑于2022年，星期六3 3 矩形网的差分格式矩形网的差分格式第31页，共77页，编辑于2022年，星期六3.1 3.1 五点差分格式五点差分格式第32页，共77页，编辑于2022年，星期六第33页，共77页，编辑于2022年，星期六第34页，共77页，编辑于2022年，星期六(i,j)(I,j-1)(i,j+1)(i+1,j)(i-1,j)第35页，共77页，编辑于2022年，星期六第36页，共77页，编辑于2022年，星期六第37页，共77页，编辑于2022年，星期六第38页，共77页，编辑于2022年，星期六第39页，共77页，编辑于2022年，星期六ABCD第40页，共77页，编辑于2022年，星期六第41页，共77页，编辑于2022年，星期六第42页，共77页，编辑于2022年，星期六3.2 3.2 边值条件的处理边值条件的处理第43页，共77页，编辑于2022年，星期六第44页，共77页，编辑于2022年，星期六第45页，共77页，编辑于2022年，星期六第46页，共77页，编辑于2022年，星期六3.3 3.3 极坐标形式的差分格式极坐标形式的差分格式第47页，共77页，编辑于2022年，星期六第48页，共77页，编辑于2022年，星期六第49页，共77页，编辑于2022年，星期六第50页，共77页，编辑于2022年，星期六3 3 三角网的差分格式三角网的差分格式第51页，共77页，编辑于2022年，星期六第52页，共77页，编辑于2022年，星期六第53页，共77页，编辑于2022年，星期六第54页，共77页，编辑于2022年，星期六第55页，共77页，编辑于2022年，星期六第56页，共77页，编辑于2022年，星期六第57页，共77页，编辑于2022年，星期六第58页，共77页，编辑于2022年，星期六例子1,2第59页，共77页，编辑于2022年，星期六第60页，共77页，编辑于2022年，星期六例子3第61页，共77页，编辑于2022年，星期六3 3 极值定理极值定理第62页，共77页，编辑于2022年，星期六第63页，共77页，编辑于2022年，星期六第64页，共77页，编辑于2022年，星期六第65页，共77页，编辑于2022年，星期六第66页，共77页，编辑于2022年，星期六第67页，共77页，编辑于2022年，星期六第68页，共77页，编辑于2022年，星期六5.2 5.2 极值定理极值定理第69页，共77页，编辑于2022年，星期六第70页，共77页，编辑于2022年，星期六第71页，共77页，编辑于2022年，星期六第72页，共77页，编辑于2022年，星期六第73页，共77页，编辑于2022年，星期六第74页，共77页，编辑于2022年，星期六第75页，共77页，编辑于2022年，星期六.差分方程、相容条件、稳定性、差分方程、相容条件、稳定性、LAX等等价定理、先验估计、极值定理等概念；价定理、先验估计、极值定理等概念；.构造差分方程方法构造差分方程方法(直接差分化、积分插直接差分化、积分插值法和变分值法和变分-差分法差分法)，矩形网和三角网的差分格，矩形网和三角网的差分格式，边界条件的处理。式，边界条件的处理。（重点重点）.如何将偏微分方程构造成相应的差分方程、如何将偏微分方程构造成相应的差分方程、对该格式的敛速估计对该格式的敛速估计（难点难点）主要内容主要内容第76页，共77页，编辑于2022年，星期六重点：重点：LAX等价定理，构造矩形网和三角网的等价定理，构造矩形网和三角网的各种差分格式。各种差分格式。难点：难点：如何将偏微分方程构造成相应的差分方程、如何将偏微分方程构造成相应的差分方程、对该格式的敛速估计。对该格式的敛速估计。重点难点第77页，共77页，编辑于2022年，星期六