数值模拟1(FDM)教材课程.ppt

地球物理数值模拟方法;一、 课程性质、目的和任务;二、 教学要求、形式和方法;三、 教学内容;三、 教学内容;四、教材及参考书目;第一讲：概论;第一讲：概论;第一讲：概论;物理模拟法：最原始和最直接的正演方法，理论上说能够解决所有正演问题，因为只要能构建出实验模型，就能得到结果。;数学模拟法：大部分的地球物理过程都能以数学的形式加以描述。 ;解析法;数值模拟方法;有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。 该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 有限差分法的缺点是必需进行整个区域的剖分，并且要求网格比较规则，空间网格最好为直角网格。;有限元法的特点;积分方程法（边界元法）的特点;数学补充：关于偏微分方程的一些基本概念;Laplace方程;几个典型方程;定解条件;二阶方程;第二讲：有限差分法;第1章 差分;第1章差分;第2章 一维波动方程的有限差分法;这个有限差分方程的截断误差是;如果已知在时间 和 ，空间 和 处的 位移值，则可以根据差分波动方程中时间 和空间 处的位移 值的显式表达式，通过简单数学运算实现对整个弹性波波场 的时间－空间逐层推进计算。这种有限差分格式称为显式格式。; ;举例;网格剖分;规则网格和不规则网格;clear N=4010; dx=0.0024; dt=0.0005; c=dt*dt/dx/dx; x=linspace(0,1,420); u(1:420,1)=0; u(181:240,1)=0.05*sin(pi*x(181:240)*7); u(2:419,2)=u(2:419,1)+c/2*(u(3:420,1)-2*u(2:419,1)+u(1:418,1)); h=plot(x,u(:,1),linewidth,3); axis([0,1,-0.05,0.05]); set(h,EraseMode,xor,makersize,18) for k=2:N set(h,Xdata,x,Ydata,u(:,2)); drawnow; u(2:419,3)=2*U(2:419,2)-U(2:419,1)+C*(U(3:420,2)... -2*U(2:419,2)+U(1:418,2)); u(2:419,1)=u(2:419,2); u(2:419,2)=u(2:419,3); end;其???，;作业;有限差分格式相容性、收敛性和稳定性;;;;第2章 一维有限差分;第2章 一维有限差分;有限差分离散（网格剖分）;有限差分近似（代替导数）;3、有限差分方程（代数方程）;3、有限差分方程（代数方程）;第2章 一维有限差分;例题：;例题：;有限差分收敛性？;第2章 一维有限差分;有限差分收敛性？——舍位误差;第2章 一维有限差分;第2章 一维有限差分;第2章 一维有限差分;第2章 一维有限差分;第2章 一维有限差分;第3章 二维有限差分;第3章 二维有限差分;第3章 二维有限差分;第3章 二维有限差分;例：;第3章 二维有限差分;第3章 二维有限差分;第3章 二维有限差分;第3章 二维有限差分;第3章 二维有限差分;第3章 二维有限差分;第3章 二维有限差分;第3章 二维有限差分;第3章 二维有限差分;第3章 二维有限差分;第3章 二维有限差分