有限差分声波方程正演

这是一个用matlab编写的有限差分的声波方程正演的程序。希望对大家有用~

% CREWES elastic finite-difference seismic modelling
% The top boundary is a free surface suitable for propagating surface
% （Rayleigh） waves.
%
% script drivers
% fdModDown ..... build a geological model from top to bottom
% fdModUp ....... build a geological model from bottom to top
% fdlayers ...... run a forward model select parameters by modifying code
%                    The control menu code is in this module
%
% Modelling subroutines
% fdModUD ....... Make Vp Vs and density vectors X/Z boundary matrices
% fdLayDoc ...... Carry f-d model through several time steps
% plotPack ...... plot a snapshot （state） with a choice of plot types
% plotcolourk ... Plot vector data with direction coded around the spectrum
% plotcoloura ... Plot vector data with direction coded around the spectrum
%                   with absolute amplitude （consistent scaling）
% plotptwist2 ... Plot pressure and twist functions of the displacement field
%
% Private subroutines
% fdModUD
% fdInitGen
% fdInitMod2
% fdLoadModel
% fdInitPad
% fdInitTrace
% fdInitArray1
% fdInitArray
% fdInitBound
% fdInitForce
% fdLoadPrior
% fdSEGY2
% plotTraces
% padright6
% zbounds
%
% Dot mat files
% The two files wncEdge.mat and wncUnity.mat （wnc: wave number correction）
% are generated to work with surface waves. The generating software which
% accounts for body waves is not yet ready to be released.
%
% The main script fdlayers writes ASCII files and a parameter file for
% command ascii2segy which is also part of the 2008/9 CREWES software release.

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1588  2010-05-27 14:48  fd\contents.m

     文件        415  2010-05-27 14:48  fd\fdCREWES.m

     文件       2515  2010-05-27 14:48  fd\fdInitArray.m

     文件       2027  2010-05-27 14:48  fd\fdInitArray1.m

     文件       2675  2010-05-27 14:48  fd\fdInitBound.m

     文件       2325  2010-05-27 14:48  fd\fdInitForce.m

     文件       2590  2010-05-27 14:48  fd\fdInitGen.m

     文件       3348  2010-05-27 14:48  fd\fdInitMod2.m

     文件       3897  2010-05-27 14:48  fd\fdInitPad.m

     文件       2282  2010-05-27 14:48  fd\fdInitTrace.m

     文件      15426  2010-05-27 14:48  fd\fdLayDoc.m

     文件       8750  2010-05-27 14:48  fd\FDlayers.m

     文件       2337  2010-05-27 14:48  fd\fdLoadModel.m

     文件       2282  2010-05-27 14:48  fd\fdLoadPrior.m

     文件       2782  2010-05-27 14:48  fd\fdModDown.m

     文件       4256  2010-05-27 14:48  fd\fdModUD.m

     文件       2826  2010-05-27 14:48  fd\fdModUp.m

     文件       4884  2010-05-27 14:48  fd\fdSEGY2.m

     文件       4293  2010-05-27 14:48  fd\padright6.m

     文件       3658  2010-05-27 14:48  fd\pedgeQuad2b.m

     文件       6041  2010-05-27 14:48  fd\plotcoloura.m

     文件       5367  2010-05-27 14:48  fd\plotcolourk.m

     文件       4514  2010-05-27 14:48  fd\plotPack.m

     文件       6254  2010-05-27 14:48  fd\plotptwist2.m

     文件       2532  2010-05-27 14:48  fd\plotTraces.m

     文件       2180  2010-05-27 14:48  fd\plotwiglx.m

     文件        318  2010-05-27 14:48  fd\suppress.m

     文件        218  2010-05-27 14:48  fd\wncEdge.mat

     文件        176  2010-05-27 14:48  fd\wncUnity.mat

     文件       2793  2010-05-27 14:48  fd\zbounds.m

............此处省略4个文件信息