Home / Articles / Structural Health Monitoring SHM-01
Sampling, Nyquist and Aliasing: The Error You Cannot Fix Afterwards
The first link in the measurement chain, and the only decision in it that cannot be undone. A badly chosen sampling rate embeds a falsehood in the data that can never be separated out again — and no amount of software will help.
1. The measurement chain
In structural health monitoring, data passes through these links before it becomes a mode:
That distinction is the whole subject of this article. If you chose the wrong window, you recompute the spectrum. If you chose the wrong sampling rate, you have no option but to repeat the measurement.
2. What sampling does
A continuous acceleration signal is read at equal intervals and converted to numbers. The number of readings per second is the sampling rate fs. What happens between readings is unknown to us — and the data does not know it is unknown either. The digital record carries no trace of the gap.
Intuitively: to recognise a cycle of vibration you need at least two points per cycle. One peak, one trough. With fewer than that, there is no mathematical way to tell a fast oscillation from a slow one.
3. The Nyquist limit and folding
That intuition is made precise by the Nyquist–Shannon sampling theorem. In a signal sampled at fs, the highest frequency that can be represented correctly is the Nyquist limit:
The apparent frequency follows from the folding relation:
4. Worked case: folded components
The table below is computed for fs = 50 Hz; the figures are produced by the engine when the page loads. The Nyquist limit is —.
| True frequency | Appears in spectrum at | Status |
|---|---|---|
| computing… | ||
The row to watch is the 45 Hz one. A high-frequency noise source — a machine, a pump, a mains harmonic — appears in the spectrum as though it were a low-frequency structural mode. The engineer looking at the spectrum has no information with which to tell them apart.
5. Why it cannot be fixed afterwards
This is why the remedy is not numerical but analogue: before sampling, while the signal is still continuous, you suppress content above Nyquist with a low-pass filter. This is the anti-aliasing filter, and it sits in hardware between the accelerometer and the ADC.
In practice the filter is not ideal; suppression after the cut-off is gradual. The sampling rate is therefore usually chosen appreciably higher than twice the highest frequency of interest — to leave room for the filter's transition band.
6. Collision: the worst case
The most insidious form of aliasing is when several components fold to the same apparent frequency.
In that case you see a single peak in the spectrum, and that peak is the sum of three separate physical events with nothing to do with one another. Interpreting its amplitude is meaningless; even attributing the existence of the peak to a mode is wrong.
7. Record length is a separate decision
The sampling rate determines how high you can look. The record length determines how finely you can discriminate. These are two independent decisions and they are often confused.
A worked figure: separating two modes 0.02 Hz apart requires a record of at least 50 seconds; raising the sampling rate from 100 to 1000 Hz changes nothing about that. This distinction is treated in detail in SHM-02.
8. Setting up a measurement
- Identify the highest mode you care about. In structural health monitoring the first three to five modes are usually enough; for building periods that is a region of a few Hz.
- Choose a sampling rate at least twice that, and in practice appreciably above. Leave room for the filter's transition band.
- Confirm that an anti-aliasing filter exists. This is not a software setting; it is looked up in the instrument's datasheet. Without one, no sampling rate will save you.
- Set the record length from the resolution you need, not from the sampling rate.
- Note the high-frequency sources in the environment. Machines, generators, traffic. Work out in advance where they will land when folded — which is exactly what the table in this article is for.
9. Test yourself
- You are sampling at fs = 200 Hz. At what frequency does a 130 Hz component appear in the spectrum?
- In that situation, does raising the sampling rate to 400 Hz solve the problem? What about doubling the record length?
- You see an unexpected peak at 5 Hz in a spectrum. How would you test whether it is folded noise? (Hint: change the sampling rate and measure again — a genuine mode stays put, a folded component moves.)
- Why must the anti-aliasing filter be analogue rather than digital?
- What is the minimum record length needed to separate two modes 0.01 Hz apart?
References
- Nyquist–Shannon sampling theorem. The derivation is universal and belongs to no national code.
- The folding relation fapparent = |f − fs·round(f/fs)| follows directly from the periodicity of the spectrum.
All figures in this article are produced by engine/dsp.js. Verification is not by comparison with third-party output but by closed loop: we generate a known signal ourselves, pass it through the chain and measure whether the engine recovers the known values.
Aliasing cannot be corrected by any numerical operation after the recording. The remedy is in hardware.