摘要：
上一篇实现了一维波动方程数值解，这一篇实现二维波动方程数值解。 二维波动方程如下： 写成差分形式： 整理一下就能得到u(i+1,j,k)。 matlab代码如下： clear all;close all;clc; t = 3; %时间范围，计算到3秒 x = 1;y = 1; %空间范围，0-1米 阅读全文
摘要：
上一篇实现了二维热传导方程数值解，这里我们计算波动方程数值解。 波动方程是一种双曲型偏微分方程。 这里依然用差分法计算。 一维波动方程如下： 写成差分形式： 整理一下就能得到u(i+1,j)。 matlab代码如下： clear all;close all;clc; t = 2; %时间范围，计算到 阅读全文
摘要：
上一篇实现了一维热传导方程数值解，这一篇实现二维热传导方程数值解。 套路是一样的，先列微分方程，再改为差分方程，然后递推求解，不同的是一维热传导需要三维显示，而二维热传导需要四维，因此最后做了个三维动态图。 二维热传导方程如下： 另外四条边界都是0。 写成差分方程为： 整理一下就能得到u(i+1,j 阅读全文