第六章解偏微分方程.PDFVIP

第六章 解偏微分方程 1 差分法解偏微分方程 2 热传导方程 ut = a 2uxx 的差分公式 2.1 显式格式 热传导方程可以写成差分形式(右边取t时刻的值计算) u(x, t + ∆t)− u(x, t) ∆t ≈ a2u(x + ∆x, t)− 2u(x, t) + u(x−∆x, t) (∆x)2 即 u(x, t+∆t) ≈ u(x, t)+ ∆t (∆x)2 a2 [u(x + ∆x, t)− 2u(x, t) + u(x−∆x, t)] 令x = i4x, t = j4t, i, j = 0, 1, 2, · · ·n − 1, r = ∆t (∆x)2 a2, 上式可写 成显式差分公式 u(i, j + 1) = (1− 2r)u(i, j) + r[u(i + 1, j)− 2u(i, j) + u(i− 1, j)] 稳定条件为∆t ≤ (∆x) 2 2a2 ,截断误差为O((∆x)2, ∆t). 例 细杆传热问题 定解问题是 ut = a 2uxx u(0, t) = 0, u(l, t) = 0 u(x, t = 0) = ϕ(x) 其中0 ≤ x ≤ 20, a = 10且 ϕ(x) = { 1, (10 ≤ x ≤ 11) 0, (x 10, x 11) 根据上面的公式,编写如下程序 x=0:20; t=0:0.01:1; a2=10;