偏微分方程数值解习题课试卷及答案.docxVIP

偏微分方程数值解习题课试卷及答案

考试时间：______分钟总分：______分姓名：______

一、

简述有限差分法求解偏微分方程的基本思想。

二、

对于一维热传导方程初边值问题：

$$

\frac{\partialu}{\partialt}=\alpha\frac{\partial^2u}{\partialx^2},\quad0\lex\leL,\quadt\ge0,

$$

$$

u(x,0)=f(x),\quad0\lex\leL,

$$

$$

u(0,t)=u(L,t)=0,\quadt\ge0,

$$

尝试推导其显式差分格式（FTCS）。

三、

利用分离变量法求解齐次拉普拉斯方程在矩形域$(0\lex\leL,0\ley\leH)$上的定解问题：

$$

\frac{\partial^2u}{\partialx^2}+\frac{\partial^2u}{\partialy^2}=0,\quad0xL,\quad0yH,

$$

$$

u(0,y)=0,\quadu(L,y)=0,\quad0\ley\leH,

$$

$$

u(x,0)=0,\quadu(x,H)=g(x),\quad0\lex\leL,

$$

导出其数值解的代数方程组（以五点差分格式为例）。

四、

考察用FTCS格式求解一维热传导方程初边值问题时，步长比$r=\alpha\Deltat/\Deltax^2$对其稳定性的影响。请说明该格式稳定存在的条件，并解释为何在实际应用中对此条件需要加以限制。

五、

对于二维波动方程初边值问题：

$$

\frac{\partial^2u}{\partialt^2}=c^2\left(\frac{\partial^2u}{\partialx^2}+\frac{\partial^2u}{\partialy^2}\right),\quad0\lex\leL,\quad0\ley\leH,\quadt\ge0,

$$

$$

u(x,y,0)=f(x,y),\quad\frac{\partialu}{\partialt}(x,y,0)=g(x,y),\quad0\lex\leL,\quad0\ley\leH,

$$

$$

u(0,y,t)=u(L,y,t)=0,\quad0\ley\leH,\quadt\ge0,

$$

$$

u(x,0,t)=u(x,H,t)=0,\quad0\lex\leL,\quadt\ge0,

$$

尝试推导其显式差分格式，并分析其稳定性条件。

六、

已知用雅可比迭代法求解线性方程组$Ax=b$时，若矩阵$A$对角占优，则该迭代法收敛。请简述对角占优的定义，并说明对角占优与雅可比迭代法收敛性之间的关系。

七、

将二维拉普拉斯方程$\nabla^2u=0$离散化，采用五点差分格式近似其内点处的偏导数。若将计算区域划分为$N_x\timesN_y$的网格，请推导出利用迭代法（如雅可比法或高斯-赛德尔法）求解该方程数值解的代数方程组的矩阵形式（不必具体写出所有元素，只需说明系数矩阵的结构特点即可）。

试卷答案

一、

有限差分法求解偏微分方程的基本思想是：首先将偏微分方程中的偏导数用有限差分格式近似表示，然后将偏微分方程的初边值问题转化为相应的差分方程（代数方程或方程组），最后在离散的网格点上进行求解，从而得到原偏微分方程数值解的近似值。

二、

推导显式差分格式（FTCS）：

1.将$\frac{\partialu}{\partialt}$在$(x_i,t_n)$处近似为$\frac{u(x_i,t_{n+1})-u(x_i,t_n)}{\Deltat}$。

2.将$\frac{\partial^2u}{\partialx^2}$在$(x_i,t_n)$处近似为$\frac{u(x_{i+1},t_n)-2u(x_i,t_n)+u(x_{i-1},t_n)}{\Deltax^2}$。

3.将这些近似代入原方程，得到：

$$\frac{u(x_i,t_{n+1})-u(x_i,t_n)}{\Deltat}=\alpha\frac{u(x_{i+1},t_n)-2u(x_i,t_n)