本文参考的是来自mooc上北京师范大学彭芳麟老师的计算物理基础
基础知识
偏微分方程的三种类型
- 椭圆型
-
- 初始条件:无
-
- 抛物型
-
- 初始条件:初始温度分布
-
- 双曲型
-
- 初始条件:初始位移与初始速度
-
边界条件
- Dirchlet边界条件
- 区域边界的函数值
- Neumann边界条件
- 给出边界上函数的法向导数
- 混合边界条件
- 给出边界上函数及其法向导数的线性组合
差分法解热传导方程
热传导方程:
令
可以获得显式公式
由此递推公式可以得出下列矩阵
例1
求解程序
x=0:20;
a2=10;
r=a2*0.01;
u=zeros(21,25);
u(10:11,1)=1;
for j=1:24
u(2:20,j+1)=(1-2*r)*u(2:20,j)+r*(u(1:19,j)+u(3:21,j));
plot(x,u(:,j));
axis([0 21 0 1]);
pause(0.1)
end
同时上述的传导方程也有隐式公式
利用
由此递推公式可以得到
引入算符
则上述矩阵可以化为
<