椭圆型方程的五点差分格式微分方程数值解计算科学系杨韧椭圆型方程的五点差分格式第三章椭圆型方程的差分格式椭圆型方程的五点差分格式§3.1正方形区域中的Laplace方程Dirichlet边值问
题
快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题
的差分模拟设Ω是xy平面中的具有正方形边界的一个有界区域,考虑Laplace方程的第一边值Dirichlet)问题椭圆型方程的五点差分格式网格节点(l,m)处的二阶中心差商代替二阶微商椭圆型方程的五点差分格式Laplace方程的五点差分格式(3.6)为截断误差为O(h2)。椭圆型方程的五点差分格式令 则Laplac方程的五点差分格式为(3.8) 即椭圆型方程的五点差分格式例1用五点差分格式求解Laplace方程 在区域 内的近似解,边界值为: 取。椭圆型方程的五点差分格式解网格点如图所示u(1,0)=20u(2,0)=20u(3,0)=20u(1,4)=180u(2,4)=180u(3,4)=180u(0,3)=80u(0,2)=80u(0,1)=80u(4,3)=0u(4,2)=0u(4,1)=0 U7 U8 U9 U4 U5 U6 U1 U2 U3椭圆型方程的五点差分格式椭圆型方程的五点差分格式椭圆型方程的五点差分格式矩阵方程AU=K,K由边界条件所确定,解得U=[U1U2U3U4U5U6U7U8U9]’=A-1K=[55.714343.214327.142979.642970.000045.3571112.8571111.785784.2857]T椭圆型方程的五点差分格式椭圆型方程的五点差分格式加密网格,取h=0.5椭圆型方程的五点差分格式椭圆型方程的五点差分格式定义向量 为从左到右,自下而上的自然次序排列的未知函 数值,则正方形区域Ω中的内部节点上的(M-1)2 个线性方程 写为矩阵方程AU=K,其中K由边界条件确定.椭圆型方程的五点差分格式椭圆型方程的五点差分格式§3.2Neumann边值问题的差分模拟
表
关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf
示函数u沿着边界的外法线方向导数,在正方形的四个顶点上法向量没有定义,取平均值代替。椭圆型方程的五点差分格式讨论左边界x=0上的导数边值条件的差分模拟又由点(0,m)的五点差分格式消去U-1,m,得0,m+1-1,m0,m1,m0,m-1椭圆型方程的五点差分格式边界x=0上(3.14) 边界x=1上(3.15) 边界y=0上(3.16) 边界y=1上(3.17)椭圆型方程的五点差分格式边界x=0-2U1,m-U0,m+14U0,m-U0,m-1-2UM-1,m-UM,m+14UM,m-UM,m-1-Ul-1,04Ul,0-Ul+1,0-2Ul,1-Ul-1,M4Ul,M-Ul+1,M-2Ul,M-1边界x=1边界y=0边界y=1椭圆型方程的五点差分格式在顶点(0,0),取偏导数的平均值作为外法线方向导数用一阶中心差商代替微商 在顶点(0,0),五点差分格式为 故-1,00,01,00,10,-1椭圆型方程的五点差分格式在四个顶点(0,0)(0,M)(M,0)(M,M)椭圆型方程的五点差分格式例1在单位正方形区域Ω上解Laplace方程的Nenmann问题解网格节点如图所示椭圆型方程的五点差分格式矩阵方程为椭圆型方程的五点差分格式令则矩阵方程为椭圆型方程的五点差分格式§3.3混合(Robins)边值条件椭圆型方程的五点差分格式例1用五点差分格式求解Laplace方程 在区域 内的近似解,边界值为: 取。椭圆型方程的五点差分格式解网格节点如图所示U1U2U3u(1,4)=180u(2,4)=180u(3,4)=180u(0,3)=80u(0,2)=80u(0,1)=80u(0,0)=80u(4,3)=0u(4,2)=0u(4,1)=0u(4,0)=0 U10 U11 U12 U7 U8 U9 U4 U5 U6椭圆型方程的五点差分格式椭圆型方程的五点差分格式椭圆型方程的五点差分格式椭圆型方程的五点差分格式解矩阵方程AU=KU=[71.821856.854332.2342…75.216561.680636.0412…87.363678.610350.2502…115.6276115.146886.3492]T椭圆型方程的五点差分格式椭圆型方程的五点差分格式椭圆型方程的五点差分格式椭圆型方程的五点差分格式作业:1、P.115例3.1取h=1/3,利用五点差分格式写出求解节点上的Ul,m值的线性方程组及矩阵方程。2、P.159习题三2