index previous next

Step-by-Step Instructions for using iPeak

1. Use the Matlab load or File > Import data command to get your data into the Matlab workspace as single vector, a pair of vectors (x and y), or a matrix with the independent variable (x) in the first column and the dependent variable (y) in the second column (e.g. "mydata").

2. Make sure the latest version of ipeak.m is in the Matlab path. Type ipeak(mydata);

3. Matlab displays the data in the Figure window, with the entire signal in the lower panel and a zoomed-in section in the top panel. In the figure below, the data are the historical sunspot counts by year from
this data set. Note: Make sure you don't click on the “Show Plot Tools” button in the toolbar above the figure; that will disable normal program functioning. If you do; close the Figure window and start again.
 
4. Use the left and right cursor keys to pan back and forth across the signal...

To pan faster, use the < and > keys. To pan by one point, use the [ and ] keys.
 
Use
the up and down cursor keys to zoom in and out.

To zoom faster, use the / and ' keys.

5. Use the A and Z keys to adjust the Amplitude Threshold (AmpT) to eliminate peaks that are too small to be of interest.

The current value of AmpT is displayed at the bottom of the window.

6. Use the S and X keys to adjust the Slope Threshold (SlopeT) to discriminate on the basis of peak width. Larger values of this parameter will neglect broad features of the signal. (Using a different signal here).

The current value of SlopeT is displayed at the bottom of the window.

7. Use the D and C keys to adjust the Smooth Width (SmoothW), the width of the smooth function that is applied to data before the slope is measured. Larger values of SmoothWidth will neglect narrow, sharp features. (In this figure the data are from ipeakdemo.m.)

The current value of SmoothW is displayed at the bottom of the window.

8. Use the F and V keys to adjust the Fit Width (FitW) to adjust the number of points around the "top part" of the (unsmoothed) peak that are taken to estimate the peak heights, positions, and widths.

The current value of FitW is displayed at the bottom of the window.

9. Use the Space and Tab keys to jump to the next and previous detected peaks. The measured peak parameters for the selected peak are displayed on the left side of the upper panel.



10. Use the Shift-G key to cycle between Gaussian, Lorentzian, and flat-top shape modes. The current mode is shown at the top of the window. The measured peak parameters change accordingly.  In this example, the peak is much more nearly Gaussian, so the Gaussian mode gives a better fit and more accurate values.


11. Use the T key to cycle between baseline correction (Autozero) modes: None, Linear, Quadratic, and Flat. The current mode is displayed at the top of the window. The measured peak parameters change accordingly.  In this example, the baseline seems to be linear, tilting down from left to right, so the Linear mode gives a better fit and more accurate values.


12. Use the Y key to switch the Y axis between linear and log mode.  In this example, the dynamic range of the Y values is so large that it's much easier to see on a log Y axis.

The y-axis label of the bottom panel changes accordingly.

13. Press the R key to display the peak report table (for the historical sunspot data in sunspots.txt):

Gaussian shape mode (press Shift-G to change)
Amplitude Threshold (AmpT) = 19.02
Slope Threshold (SlopeT) = 0.0048225
Smooth Width (SmoothW) = 5 points
Fit Width (FitW) = 5 points
No baseline correction

    Peak#   Position    Height    Width      Area
      1     1704.9      43.374    4.1307     190.74
      2     1717.5      61.043    4.6045     299.23
      3     1727.3      101.73    5.6906     616.28
      4     1738.1      106.81    4.9657     564.67
      5     1749.4      73.725    6.3741     500.28
      6     1760.9      71.524    6.0076     457.45
      7     1769.8      103.35    5.1104     562.29
..... etc

14. Press the E key to display the peak statistics table:

Peak Summary Statistics

26 peaks detected
No baseline correction
          Interval   Height    Width     Area
Maximum    15.153   195.8527   8.356   1025.5604
Minimum     8.2455   41.5679   4.1307   190.7358
Mean       11.0096   96.5382   5.6021   569.2662
% STD      13.9588   38.2597  15.4761    36.6246

This shows the well-known 11-year sunspot period as the mean of the interval. Plots of the histograms are also shown in Figure window 2.

15. Press the U key to switch between peak and valley modes.

16.
Press the N key to perform a (Normal) curve fit all numbered peaks in the upper window. Alternatively, press the M key to perform a (Multiple) curve fit all numbered peaks in the entire signal. The curve fit results are displayed in Figure 2 just as they are for the peakfit function and printed out in the command window.

17. Press K to display the complete set of keystroke commands..




This page is part of "A Pragmatic Introduction to Signal Processing", created and maintained by Prof. Tom O'Haver , Department of Chemistry and Biochemistry, The University of Maryland at College Park. Comments, suggestions and questions should be directed to Prof. O'Haver at toh@umd.edu.

Last updated,  November, 2016

Website created with SeaMonkey.


index previous next

Number of unique visits since May 17, 2008:

counter create hit