Appendix D: Exercises
Series | Geophysical References Series |
---|---|
Title | Digital Imaging and Deconvolution: The ABCs of Seismic Exploration and Processing |
Author | Enders A. Robinson and Sven Treitel |
Chapter | 4 |
DOI | http://dx.doi.org/10.1190/1.9781560801610 |
ISBN | 9781560801481 |
Store | SEG Online Store |
For the want of a nail the shoe was lost, For the want of a shoe the horse was lost, For the want of a horse the rider was lost, For the want of a rider the battle was lost, For the want of a battle the kingdom was lost, And all for the want of a horseshoe nail. —Benjamin Franklin
1. Fill in the table:
Sampling interval | Sampling rate | Nyquist frequency |
---|---|---|
0.001 s | ||
0.002 s | ||
0.003 s | ||
0.004 s | ||
0.005 s | ||
0.006 s | ||
0.007 s | ||
0.008 s |
2. For , fill in the table:
Actual frequency f | Aliased frequency resulting from sampling |
---|---|
0 | |
100 | |
200 | |
300 | |
400 | |
500 |
3. Show that Euler’s equation for gives . This equation is considered to be one of the most beautiful equations in mathematics because it relates the four special numbers to each other.
4. Discuss why the alias phenomenon represents a frequency ambiguity resulting from the sampling process.
5. Suppose there are fewer than two samples per cycle of an input sinusoidal signal of a given frequency. Show that the same sample values would occur for a sinusoidal signal of a still lower frequency. Hence, that signal looks the same as the original signal.
6. Half of the sampling frequency is called the folding or Nyquist frequency . Show that the frequency appears as the smaller apparent frequency . These two frequencies are aliases of each other.
7. To avoid aliasing, show why frequencies above the Nyquist frequency must be removed by an anti-alias filter before sampling.
8. Discuss why aliasing is an inherent property of all sampling systems. For example, aliasing applies to (1) sampling done at discrete time intervals, as with digital seismic recording, and (2) sampling done at discrete spatial intervals by separate elements of geophone and source-point arrays.
9. Suppose that . Show that f is indistinguishable from the lower (alias) frequency .
10. Describe what happens when a signal is sampled with fewer than two points per cycle. [Hint: The effect is somewhat similar to the effect observed if a wheel rotating at a high rate of speed is illuminated stroboscopically at a different rate. The wheel appears to be rotating at a speed different from its true rate.]
11. Illustrate Exercise 10 for a few specific cases, using a 0.004-s sampling interval. [Hint: Consider a 250-Hz sine wave signal sampled 250 times per second. A constant value (anywhere between -1 and +1) will result. Thus, a sampled 250-Hz signal has the appearance of a zero-frequency signal. The analogous situation for the rotating wheel analogy is that the stroboscopic rate causes the wheel to appear to be standing still.]
12. In Exercise 11, has the sampling process caused a 250-Hz signal to appear to be a zero-frequency signal? [Hint: Now consider an original 150-Hz sine-wave signal. The signal resulting from sampling at a 4-ms interval has the appearance of a 100-Hz signal, because 150 – 250 = –100 Hz, which is the same as 100 Hz. In other words, a signal above 125 Hz takes on an alias as a result of sampling. The critical frequency (125 Hz in this case) is called the Nyquist frequency, or the half-sampling frequency.]
13. At a 4-ms sampling interval, what happens to sinusoidal signals above 250 Hz in frequency? This case corresponds to the situation in the stroboscope analogy in which the wheel appears to turn counter to its true direction of rotation. Frequencies between 250 and 375 Hz appear as aliased frequencies between –125 Hz and 0 Hz.
14. In terms of frequency expressed as cycles per sampling interval, show that the Nyquist frequency always turns out to be 0.5. The corresponding Nyquist angular frequency is . For example, for a sampling interval of 0.004 s, the Nyquist frequency in Hz is given as (0.5 cycles/time-unit)/(0.004 s/time-unit) = (125 cycles/s).
Continue reading
Previous section | Next section |
---|---|
Sampling geophysical data | none |
Previous chapter | Next chapter |
Visualization | Filtering |
Also in this chapter
- Time series
- The wavelet
- Digitization
- Frequency
- Sinusoidal motion
- Aliasing
- The Nyquist frequency
- Sampling geophysical data