第五讲——显式差分和隐式差分5.pptxVIP

回顾有限差分法基础差分格式差分方程边界条件的处理相容性、稳定性和收敛性回顾1. 有限差分法的相容性、稳定性和收敛性相容性：针对差分格式而言，在时间步长和空间步长趋近于零的情况下，如果差分格式的截断误差（差分格式与原有偏微分方程之差）的模趋近于零，则该差分格式与原偏微分方程是相容的，或称该差分方程与原偏微分方程具有相容性。稳定性（stability）：如果偏微分方程的严格解析解有界，差分格式给出的解也有界，称该差分格式是稳定的；如果差分格式给出的解是无界的，则称该差分格式是不稳定的。稳定性反映了差分格式在计算中控制误差传递的能力收敛性（convergence）：如果当时间和空间步长趋于零时，FDE解趋于PDE解，称该差分格式是收敛的。如果则称该差分格式是收敛的。收敛性描述的是当差分网格无限细化时，差分方程的解是否具有无限逼近偏微分方程的解的能力Lax等价定理（Lax equivalence theorem）：如果逼近一个给定问题的差分格式是相容的，那么该差分格式的收敛性与稳定性互为充分必要条件。相容性是比较容易满足的。在此基础上，如果满足了稳定性条件，差分格式的收敛性就自动满足。2.5 有限差分法实例U=100(i-1,j+1)(i,j+1)2(i+1,j+1)j+1h2013jh3h1U=0U=0(i-1,j)(i,j)(i+1,j)h44j-1(i-1,j-1)(i,j)(i+1,j-1)123589134671011141512i-1ii+1U=0局部节点编号(i-1,j+1)(i,j+1)2(i+1,j+1)j+1总体节点编号h2013jh3h1(i-1,j)(i,j)(i+1,j)h44j-1(i-1,j-1)(i,j)(i+1,j-1)for j=2:n-1for i=2:m-1; a((j-1)*m+i,(j-1)*m+i+1)=1; a((j-1)*m+i,(j-1)*m+i-1)=1; a((j-1)*m+i,j*m+i)=1; a((j-1)*m+i,(j-2)*m+i)=1; a((j-1)*m+i,(j-1)*m+i)=-4; endendA矩阵非零系数减少，同时引入第一类边界，方程右端项B向量出现非零元素。i-1ii+1内部节点：边界节点：组建A和B矩阵，求解线性方程组得到Xa(121,106)=-0.25;a(135,134)=-0.25;a(135,120)=-0.25;a(15,14)=-0.25;a(15,30)=-0.25;for i=2:14 a(i,i-1)=-0.25; a(i,i+1)=-0.25; a(i,i+15)=-0.25; endfor i=122:134 a(i,i-1)=-0.25; a(i,i+1)=-0.25; a(i,i-15)=-0.25; endfor i=1:7for j=2:14;a(15*i+j,15*i+j-1)=-0.25;a(15*i+j,15*i+j+1)=-0.25;a(15*i+j,15*i+j+15)=-0.25;a(15*i+j,15*i+j-15)=-0.25;endend%Matlab 2Dclear;clc;figure(color,w); a=zeros(135,135);for i=1:135 a(i,i)=1;end;for i=1:7 a(15*i+1,15*i+2)=-0.25; a(15*i+1,15*i+16)=-0.25; a(15*i+1,15*i-14)=-0.25;endfor i=1:7 a(15*i+15,15*i+14)=-0.25; a(15*i+15,15*i+30)=-0.25; a(15*i+15,15*i)=-0.25; Enda(1,2)=-0.25;a(1,16)=-0.25;a(121,122)=-0.25;b=a^(-1);c=zeros(135,1);for i=121:135 c(i,1)=25;endd=b*c;s=zeros(11,17);for i=2:16 s(11,i)=100; endfor i=1:9for j=1:15;s(i+1,j+1)=d(15*(i-1)+j,1);endend subplot(1,2,1),mesh(s)axis([0,17,0,11,0,100])subplot(1,2,2),contour(s,32)2.5 应用实例南加州一次未来大地震的强地面运动的数值模拟盆地效应Cui, 2013Cui, 2013Cui, 2013Cui, 2013总结：1、有限差分方法给出的数值解的精度取决于所用的差分形式（向前、向后、中心）。2、偏微分方程的显式有限差分格式通常是有条件稳定的，为了保证得到精确的数值解，最关键的是需要根据稳定性条件选