Frequency-Multiplexed Multisource Least-Squares Prestack Phase Shift Migration

Procedure:

1. Simulation starts by simply running
   >> run_LSM_muS_phSh_phEn_CGorMomen_scr
   where >> represents MATLAB prompt. But before running this program of run_LSM_muS_phSh_phEn_CGorMomen_scr.m , you need to edit its first 15 lines of code (not counting the comment), in order to customize the path of your directory and your platform.

2. There are some variables in run_LSM_muS_phSh_phEn_CGorMomen_scr.m controlling the behavior of the simulation.
   const_vel (logical) whether to use constant velocity model or varying velocity model.
   Niter (integer) number of iterations for least-squares migration.
   is_sta_Sinu = [a, b] (logical array). a: whether to use static encoding as opposed to dynamic encoding; b: whether to use frequency-multiplexing as opposed to the standard approach of using Ricker wavelet.
   Lms (integer) number of shots per supergather. The number of supergathers is then calculated as ns/Lms, where ns is the number of sources.
   using_CG (integer) which optimization method to use. The recommended value is 0, representing the steepest descent algorithm.
   fnm (string) the stem of the filenames for the output data files. For example, if fnm = 'LSMSP_', then you can expect ouputfiles with the name such as LSMSP_dynSinu_40x8_ref.bin . Here, dyn represents dynamic encoding; Sinu represents frequency-multiplexing; 8 supergathers consisting of 40 shots each, if such parameter settings have been chosen in the program; ref says that this file is the resulting reflectivity model.