五点差分格式课件.pptVIP

微分方程数值解计算科学系杨韧椭圆型方程的五点差分格式

第三章椭圆型方程的差分格式椭圆型方程的五点差分格式

§3.1正方形区域中的Laplace方程Dirichlet的差分模设Ω是xy平面中的具有正方形边界的一个有界区域，考虑Laplace方程的第一边值Dirichlet）问题椭圆型方程的五点差分格式

l,m+1网格节点（l，m）处的二阶中心差商代替二阶微商l–1,ml,ml+1,ml,m–1椭圆型方程的五点差分格式

Laplace方程的五点差分格式(3.6)为-Ul,m+1截断误差为O(h2)。-Ul–1,m4Ul,m-Ul+1,m-Ul,m–1椭圆型方程的五点差分格式

令则Laplac方程的五点差分格式为（3.8）即椭圆型方程的五点差分格式

例1用五点差分格式求解Laplace方程在区域内的近似解，边界值为：取。椭圆型方程的五点差分格式

解网格点如图所示u(1,4)=180u(2,4)=180u(3,4)=180U7U4U8U5U9U6u(4,3)=0u(4,2)=0u(0,3)=80u(0,2)=80u(0,1)=80UUU123u(4,1)=0u(1,0)=20u(2,0)=20u(3,0)=20椭圆型方程的五点差分格式

椭圆型方程的五点差分格式

椭圆型方程的五点差分格式

矩阵方程AU=K，K由边界条件所确定，解得U=[UUUUUUUUU]’=A-1K123456789=[55.714343.214327.142979.642970.000045.3571112.8571111.785784.2857]T椭圆型方程的五点差分格式

椭圆型方程的五点差分格式

加密网格，取h=0.5椭圆型方程的五点差分格式

椭圆型方程的五点差分格式

定义向量为从左到右，自下而上的自然次序排列的未知函数值，则正方形区域Ω中的内部节点上的(M-1)个线性方程2写为矩阵方程AU=K,其中K由边界条件确定.椭圆型方程的五点差分格式

椭圆型方程的五点差分格式

§3.2Neumann的差分模表示函数u沿着边界的外法线方向导数，在正方形的四个顶点上法向量没有定义，取平均值代替。椭圆型方程的五点差分格式

讨论左边界x=0上的导数边值条件的差分模拟0,m+1又由点（0，m）的五点差分格式-1,m0,m1,m0,m-1消去U，得-1,m椭圆型方程的五点差分格式

边界x=0上(3.14)边界x=1上(3.15)边界y=0上(3.16)边界y=1上(3.17)椭圆型方程的五点差分格式

边界x=0边界x=1-U-UM,m+10,m+14U,4U,-2UMm0m-2U1,mM-1,m-U,-U,0m-1Mm-1边界y=1边界y=0-Ul-1,M4U-Ul+1,M-2Ul,1l,M-2U-Ul,M-1l-1,04U-Ul,0l+1,0椭圆型方程的五点差分格式

在顶点（0，0）,取偏导数的平均值作为外法线方向导数0,1用一阶中心差商代替微商-1,00,01,00,-1在顶点（0，0）,五点差分格式为故椭圆型方程的五点差分格式

在四个顶点（0，0）（0，M）（M，0）（M，M）椭圆型方程的五点差分格式

例1在单位正方形区域Ω上解Laplace方程的Nenmann问题U7U4U1U8U顶点9解网格节点如图所示U内点U边界点56U2U3椭圆型方程的五点差分格式

矩阵方程为椭圆型方程的五点差分格式

令则矩阵方程为椭圆型方程的五点差分格式

§3.3混合（Robins）边值条件椭圆型方程的五点差分格式

例1用五点差分格式求解Laplace方程在区域内的近似解，边界值为：取。椭圆型方程的五点差分格式

解网格节点如图所示u(1,4)=180u(2,4)=180u(3,4)=180U10U7U11U8U12U9u(0,3)=80u(4,3)=0u(4,2)=0u(0,2)=80u(0,1)=80u(0,0)=80U4U5U6u(4,1)=0u(4,0)=0UUU123椭圆型方程的五点差分格式

椭圆型方程的五点差分格式

椭圆型方程的五点差分格式

椭圆型方程的五点差分格式

解矩阵方程AU=KU=[71.821856.854332.2342…75.216561.680636.0412…87.363678.610350.2502…115.6276115.146886.3492]T椭圆型方程的五点差分格式

椭圆型方程的五点差分格式

椭圆型方程的五点差分格式

椭圆型方程的五点差分格式

作业:1、P.1