PML二维声波方程有限差分matlab程序

带PML条件的二维声波方程有限差分正演模拟程序，采用matlab编写。

function [ vp ] = create_mod（ model ）
% CREATE_MOD 
% Create Vp-model on-the-fly

     % Example of a simple on-the-fly Vp model
     % - be more creative
     for j=1:model.ny
         for i=1:model.nx
             
             vp（ji） = 4000.0; 
             
         end
     end

end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\
     目录           0  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\
     文件        2099  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\PML.m
     文件         514  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\RHS_AC.m
     文件         305  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\create_mod.m
     文件         702  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\def_acq.m
     文件         992  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\extend_model.m
     文件         375  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\extract_model.m
     文件        4551  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\germaine.m
     文件        1768  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\init_A_AC_5p.m
     文件        4790  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\init_A_AC_9p_mixed_pml.m
     文件        1522  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\inp_model.m
     目录           0  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\input_files\
     文件        1652  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\input_files\inp_model_marmousi.m
     文件        1218  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\read_acq.m
     文件         267  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\readmod.m
     目录           0  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\receiver\
     文件        7190  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\receiver\receiver_OBC.dat
     文件         354  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\shift_acq.m
     目录           0  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\source\
     文件        6802  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\source\source_OBC_VSP.dat
     目录           0  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\start\
     文件      348000  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\start\marmousi_II_marine.vp
     目录           0  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\tmp\
     文件        1239  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\tmp\SRM.m
     目录           0  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\util\
     文件        9800  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\GERMAINE-AC-Matlab-Edition\util\seismic.map
     文件        1925  2018-07-12 03:59  10_2D_acoustic_FDFD_PMLs\QUICKSTART.txt