Practise: MUSIC teleseismic Back-Projection Multiple Signal - - PowerPoint PPT Presentation

▶

Jun 02, 2023 445 likes •590 views

Practise: MUSIC teleseismic Back-Projection Multiple Signal Classification (MUSIC) q MUSICBP is a MATLAB implementation of the teleseismic MUSIC back-projections. q Open wrapper General_BP.m q Initialization (creating project folder, etc) q Read

SLIDE 1

Practise: MUSIC teleseismic Back-Projection

SLIDE 2

qMUSICBP is a MATLAB implementation of the teleseismic MUSIC back-projections. qOpen wrapper General_BP.m qInitialization (creating project folder, etc) qRead seismograms in SAC format qAlign Seismograms by initial P arrival qRun either MUSIC or beamforming solver qThe code produces movies, summary plot and distance vs time plot.

Multiple Signal Classification (MUSIC)

SLIDE 3

Cd Cp –r /home/seismology_lecture/lmeng/musicbp-master ./ Cd ~/musicbp-master/MUSICBP Matlab & Open General_BP.m

SLIDE 4

General_BP.m

SLIDE 5

Initialization

%Set Initial_Flag=1 %run General_BP.m

SLIDE 6

Cp ./PaluAUdata/*SAC ./Palu2018/Data/ %Set ReadBP_Flag=1 %run General_BP.m

Reading Seismograms

SLIDE 7

%Set AlignBP_Flag=1 % Set bandChoice from 1 to 4 %run General_BP.m 4 times. The code progressively aligns the sesimograms at increasing frequencies.

Aligning Seismograms

SLIDE 8

First Alignment

SLIDE 9

Second Alignment

SLIDE 10

Third Alignment

SLIDE 11

Fourth Alignment

SLIDE 12

Final Station Map

SLIDE 13

Running MUSIC back-projection

Movie/snapshot Summary Plot Distance vs Time Normalized Power

SLIDE 14

Use the beamforming solver and see if there’s any difference. Try your own alignment and station selections If you messed up the parameters in General_BP.m, check General_BP_Backup.m for the set of parameters I use.