第5周实践：隐式与Crank-Nicoon方法求解热方程及其收敛性分析.pdfVIP

实践问题

第5周：隐式方法和Crank‑Nicolson方法求解热方程N

请创建一个practical2.m文件来完成以下任务！

实现隐式方法和Crank‑Nicolson方法来求解下面函数u(x,τ)的热方程，并通过

观察误差的变化和双对数图来研究时间方向和空间方向上的收敛阶：

uτ−uxx=0

对于0≤x≤1和0≤τ≤0.1，初始条件为

u(x,0)=2sin(2πx)

边界条件为

u(0,τ)=0

u(1,τ)=0

提示：

•创建一个名为Tridiag的函数，以实现讲座中讨论的三对角求解器，并

按照讲座幻灯片中的伪代码进行实现。

•创建两个新函数ImplicitFDM_Heat.m和CrankNicolsonFDM_Heat.m，并通过

修改ExplicitFDM_Heat.m中的代码来编写这些新函数。