matlab画微分方程的矢量场图_利用matlab求解偏微分方程

本文参考的是来自mooc上北京师范大学彭芳麟老师的计算物理基础

基础知识

偏微分方程的三种类型

• 椭圆型
  • 初始条件：无
• 抛物型
  • 初始条件：初始温度分布
• 双曲型
  • 初始条件：初始位移与初始速度

边界条件

• Dirchlet边界条件
  • 区域边界的函数值
• Neumann边界条件
  • 给出边界上函数的法向导数
• 混合边界条件
  • 给出边界上函数及其法向导数的线性组合

差分法解热传导方程

热传导方程：

令

可以获得显式公式

由此递推公式可以得出下列矩阵

例1

求解程序

x

同时上述的传导方程也有隐式公式

利用

得到

由此递推公式可以得到

引入算符

则上述矩阵可以化为