is an interactive Matlab function that performs smoothing for
time-series signals using the fastsmooth algorithm, with keystrokes
that allow you to adjust the smoothing parameters continuously while observing the
effect on your signal dynamically. Other functions include
differentiation, peak sharpening, and least-squares peak
measurement. View the code here or dowload
the ZIP file with sample data for
you can right-click on any of the m-file links on this site
and select Save Link As... to download them to
your computer for use within Matlab.
To use it, just place
isignal.m in the Matlab path and type
>> isignal(DataMatrix); or >> isignal(x,y); or >> isignal(y);
Use the cursor arrow keys to pan and zoom. The S key (or input argument "SmoothMode") cycles through four smoothing modes:
If SmoothMode=0, the signal is not smoothed.
If SmoothMode=1, rectangular (sliding-average or boxcar)
If SmoothMode=2, triangular (2 passes of sliding-average)
If SmoothMode=3, pseudo-Gaussian (3 passes of sliding-average)
The A and Z keys (or optional input argument SmoothWidth) control the SmoothWidth, w. Example shown in the figure on the right. New in Version 6: To specify a
press Shift-Q. You can specify the smooth width vector in
two ways: at the prompt you can (a) enter the number of
segments (then you'll be prompted to enter the smooth widths in the
first and last segments, and the computer will calculate integer
values of smooth widths that are evenly divided between the
specified first and last values, or (b) type in the smooth width
vector directly including the square brackets, e.g. [1 3 3
9]. In either case, subsequently adjusting the smooth width with the
A and Z keys will vary all the segments by
the same percentage factor. (To return to an ordinary single segment
smooth, enter 1 as the number of segments). Picture on right.
The X key toggles "ends" between 0 and 1. This determins how the "ends" of the signal (the first w/2 points and the last w/2 points) are handled when smoothing:
If ends=0, the ends are zero.
If ends=1, the ends are smoothed with progressively smaller smooths the closer to the end.
Generally, ends=1 is best, except in some cases using the derivative mode when ends=0 result in better vertical centering of the signal.
when smoothing peaks, you can easily measure the effect of smoothing on
peak height and width by turning on peak measure mode (press P) and then press S to cycle through the smooth modes.
Effect of smoothing on the frequency spectrum. You can also visualize the effect of smoothing on the frequency spectrum of the signal by pressing Shift-S, which displays the frequency spectrum in the lower window (temporarily replacing the full-signal display). Use the pan and zoom keys to adjust the region of the signal to be viewed. PressShift-Ato cycle through four plot modes (linear, semilog X, semilog Y, or log-log) and pressShift-Xto toggle between afrequency on the x axis andtimeon the x-axis. All signal processing functions remain active in the frequency spectrum mode (smooth, derivative, etc) so you can observe the effect of these functions on the frequency spectrum of the signal, as in the animated figure on the left.