二阶椭圆型第一边值问题的数值解法（五点差分格式和有限体积法）附matlab代码

本文链接：https://blog.csdn.net/Sophiayinqianbao/article/details/128407019

本文介绍了使用五点差分格式和有限体积法解决二阶椭圆型第一边值问题的方法，详细阐述了网格剖分、微分算子离散的过程，并附带MATLAB代码实现。通过计算不同网格的最大误差和收敛阶，验证了两种方法的收敛性，结果显示它们都具有二阶收敛阶。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

二阶椭圆型第一边值问题的数值解法（五点差分格式和有限体积法）附matlab代码

这里我们介绍五点差分格式和有限体积法求椭圆型第一边值问题，
题目：
分别采用矩形网格的五点差分格式和有限体积法求椭圆型第一边值问题，
$-\frac{\partial^2u}{\partial x^2}-\frac{\partial^2u}{\partial y^2}+a\frac{\partial u}{\partial x}=0,a>0,$
$u\left(0,y\right)=sin\left(\pi y\right),u\left(1,y\right)=0,0<y<1,$
$u\left(x,0\right)=u\left(x,1\right)=0,0<x<1.$

取参数 a = 40，得到不同网格最大误差和收敛阶。

问题分析

(i) 矩形网格的五点差分法:
Step1：网格剖分。取定沿 x,y 轴方向的步长分别为： $h_x,h_y$ ,便于计算我们取 $h_x=h_y$ ,令 $=\sqrt{h_x^2+h_y^2}$ ,作两族平行于 x,y 轴的直线 $x=x_i,y=y_i$ ，其中 $x_i=x_0+ih_x,i=0,1,2,\cdots,N,y_j=y_0+jh_y,j=0,1,2,\cdots,N$ ,两族直线的交点 $ih_x,jh_y)$ 称为网点或者节点,记为 $x_i , y_j)$ .
Step2：微分算子离散.
在正则节点 $x_i , y_j)$ 上沿 x,y 方向分别使用中心差商公式对 $u_{xx}$ 和 $u_{yy}$ 作近似，得到方程的五点差分格式为：
$-\left(\frac{u_{i+1,j}-2u_{i,j}+u_{i-1,j}}{h_x^2}+\frac{u_{i,j+1}-2u_{i,j}+u_{i,j-1}}{h_y^2}\right)+a\frac{u_{i+1,j}-u_{i-1,j}}{2h_x}=0$
$\Rightarrow$
$\left(-\frac{1}{h_x^2}-\frac{a}{2h_x}\right)u_{i-1,j}+\left(\frac{2}{h_x^2}+\frac{2}{h_y^2}\right)u_{i,j}+\left(-\frac{1}{h_x^2}+\frac{a}{2h_x}\right)u_{i+1,j}+\left(-\frac{1}{h_y^2}\right)u_{i,j-1}+\left(-\frac{1}{h_y^2}\right)u_{i,j+1}=0,i,j=1,2,\cdots,N-1$
Let
$A=\left[\begin{matrix}\frac{2}{h_x^2}+\frac{2}{h_y^2}&-\frac{1}{h_x^2}+\frac{a}{2h_x}&0&\cdots&0&0\\-\frac{1}{h_x^2}-\frac{a}{2h_x}&\frac{2}{h_x^2}+\frac{2}{h_y^2}&-\frac{1}{h_x^2}+\frac{a}{2h_x}&\cdots&0&0\\\cdots&\cdots&\cdots&\cdots&\cdots&\cdots\\0&0&0&\cdots&-\frac{1}{h_x^2}-\frac{a}{2h_x}&\frac{2}{h_x^2}+\frac{2}{h_y^2}\\\end{matrix}\right]$