Wavelet Analysis of a Seismic Event

The seismographic data for an earthquake event are contained in the variable seismicData.
We will analyze this data with a Daubechies least asymmetric filter of order 8.

wtdata  = 
      WaveletPacketTransform[seismicData,  LeastAsymmetricFilter[8] ];

Phase Space Results

This shows the graph of the seismic event above the Phase Space Plot of the result from WaveletPacketTransform. The higher frequency components appear at the onset of the seismic shock, as we would expect.

\(seismicplot\  = ListPlot[seismicData, \ \ plotoptions1]; 

    psplot\  = PhaseSpacePlot[GrayCodeOrder[wtdata], \ plotoptions2]; 


    Show[GraphicsArray[{{seismicplot}, {psplot}}], 
      DisplayFunction -> $DisplayFunction\ ]; \)