1. Basic Signal Analysis Computations
The basic computations for analyzing signals include converting from a twosided power spectrum to a singlesided power spectrum, adjusting frequency resolution and graphing the spectrum, using the FFT, and converting power and amplitude into logarithmic units.
The power spectrum returns an array that contains the twosided power spectrum of a timedomain signal. The array values are proportional to the amplitude squared of each frequency component making up the timedomain signal. A plot of the twosided power spectrum shows negative and positive frequency components at a height
where A_{k} is the peak amplitude of the sinusoidal component at frequency k. The DC component has a height of A_{0}^{2} where A_{0} is the amplitude of the DC component in the signal.
Figure 1 shows the power spectrum result from a timedomain signal that consists of a 3 Vrms sine wave at 128 Hz, a 3 Vrms sine wave at 256 Hz, and a DC component of 2 VDC. A 3 Vrms sine wave has a peak voltage of 3.0 or about 4.2426 V. The power spectrum is computed from the basic FFT function. Refer to the Computations Using the FFT section later in this application note for an example this formula.
Figure 1. TwoSided Power Spectrum of Signal
Converting from a TwoSided Power Spectrum to a SingleSided Power Spectrum
Most realworld frequency analysis instruments display only the positive half of the frequency spectrum because the spectrum of a realworld signal is symmetrical around DC. Thus, the negative frequency information is redundant. The twosided results from the analysis functions include the positive half of the spectrum followed by the negative half of the spectrum, as shown in Figure 1.
In a twosided spectrum, half the energy is displayed at the positive frequency, and half the energy is displayed at the negative frequency. Therefore, to convert from a twosided spectrum to a singlesided spectrum, discard the second half of the array and multiply every point except for DC by two.
where S_{AA}(i) is the twosided power spectrum, G_{AA}(i) is the singlesided power spectrum, and N is the length of the twosided power spectrum. The remainder of the twosided power spectrum S_{AA}
The nonDC values in the singlesided spectrum are then at a height of
This is equivalent to
where
is the root mean square (rms) amplitude of the sinusoidal component at frequency k. Thus, the units of a power spectrum are often referred to as quantity squared rms, where quantity is the unit of the timedomain signal. For example, the singlesided power spectrum of a voltage waveform is in volts rms squared.
Figure 2 shows the singlesided spectrum of the signal whose twosided spectrum Figure 1 shows.
Figure 2. SingleSided Power Spectrum of Signal in Figure 1
As you can see, the level of the nonDC frequency components are doubled compared to those in Figure 1. In addition, the spectrum stops at half the frequency of that in Figure 1.
Refer to the Power Spectrum topic in the LabVIEW Help (linked below) for the most updated information about the power spectrum.
Adjusting Frequency Resolution and Graphing the Spectrum
Figures 1 and 2 show power versus frequency for a timedomain signal. The frequency range and resolution on the xaxis of a spectrum plot depend on the sampling rate and the number of points acquired. The number of frequency points or lines in Figure 2 equals
where N is the number of points in the acquired timedomain signal. The first frequency line is at 0 Hz, that is, DC. The last frequency line is at
where F_{s} is the frequency at which the acquired timedomain signal was sampled. The frequency lines occur at f intervals where
Frequency lines also can be referred to as frequency bins or FFT bins because you can think of an FFT as a set of parallel filters of bandwidth f centered at each frequency increment from
Alternatively you can compute f as
The computations for the frequency axis demonstrate that the sampling frequency determines the frequency range or bandwidth of the spectrum and that for a given sampling frequency, the number of points acquired in the timedomain signal record determine the resolution frequency. To increase the frequency resolution for a given frequency range, increase the number of points acquired at the same sampling frequency. For example, acquiring 2,048 points at 1.024 kHz would have yielded f = 0.5 Hz with frequency range 0 to 511.5 Hz. Alternatively, if the sampling rate had been 10.24 kHz with 1,024 points, f would have been 10 Hz with frequency range from 0 to 5.11 kHz.
Computations Using the FFT
The power spectrum shows power as the mean squared amplitude at each frequency line but includes no phase information. Because the power spectrum loses phase information, you may want to use the FFT to view both the frequency and the phase information of a signal.
The phase information the FFT yields is the phase relative to the start of the timedomain signal. For this reason, you must trigger from the same point in the signal to obtain consistent phase readings. A sine wave shows a phase of 90° at the sine wave frequency. A cosine shows a 0° phase. In many cases, your concern is the relative phases between components, or the phase difference between two signals acquired simultaneously. You can view the phase difference between two signals by using some of the advanced FFT functions. Refer to the FFTBased Network Measurement section of this application note for descriptions of these functions.
The FFT returns a twosided spectrum in complex form (real and imaginary parts), which you must scale and convert to polar form to obtain magnitude and phase. The frequency axis is identical to that of the twosided power spectrum. The amplitude of the FFT is related to the number of points in the timedomain signal. Use the following equation to compute the amplitude and phase versus frequency from the FFT.
where the arctangent function here returns values of phase between  and +, a full range of 2 radians. Using the rectangular to polar conversion function to convert the complex array
to its magnitude (r) and phase (ø) is equivalent to using the preceding formulas.
The twosided amplitude spectrum actually shows half the peak amplitude at the positive and negative frequencies. To convert to the singlesided form, multiply each frequency other than DC by two, and discard the second half of the array. The units of the singlesided amplitude spectrum are then in quantity peak and give the peak amplitude of each sinusoidal component making up the timedomain signal. For the singlesided phase spectrum, discard the second half of the array.
To view the amplitude spectrum in volts (or another quantity) rms, divide the nonDC components by the square root of two after converting the spectrum to the singlesided form. Because the nonDC components were multiplied by two to convert from twosided to singlesided form, you can calculate the rms amplitude spectrum directly from the twosided amplitude spectrum by multiplying the nonDC components by the square root of two and discarding the second half of the array. The following equations show the entire computation from a twosided FFT to a singlesided amplitude spectrum.
where i is the frequency line number (array index) of the FFT of A.
The magnitude in volts rms gives the rms voltage of each sinusoidal component of the timedomain signal.
To view the phase spectrum in degrees, use the following equation.
The amplitude spectrum is closely related to the power spectrum. You can compute the singlesided power spectrum by squaring the singlesided rms amplitude spectrum. Conversely, you can compute the amplitude spectrum by taking the square root of the power spectrum. The twosided power spectrum is actually computed from the FFT as follows.
where FFT*(A) denotes the complex conjugate of FFT(A). To form the complex conjugate, the imaginary part of FFT(A) is negated.
When using the FFT in LabVIEW and LabWindows/CVI, be aware that the speed of the power spectrum and the FFT computation depend on the number of points acquired. If N can be factored into small prime numbers, LabVIEW and LabWindows/CVI uses a highly efficient CooleyTukey mixedradix FFT algorithm. Otherwise (for large prime sizes), LabVIEW uses other algorithms to compute the discrete Fourier transform (DFT), and these methods often take considerably longer. For example, the time required to compute a 1000point and 1024point FFT are nearly the same, but a 1023point FFT may take twice as long to compute. Typical benchtop instruments use FFTs of 1,024 and 2,048 points.
So far, you have looked at display units of volts peak, volts rms, and volts rms squared, which is equivalent to meansquare volts. In some spectrum displays, the rms qualifier is dropped for Vrms, in which case V implies Vrms, and V^{2} implies Vrms^{2}, or meansquare volts.
Refer to the Computing the Amplitude and Phase Spectrums topic in the LabVIEW Help (linked below) for the most updated information about using the FFT for computations.
Converting to Logarithmic Units
Most often, amplitude or power spectra are shown in the logarithmic unit decibels (dB). Using this unit of measure, it is easy to view wide dynamic ranges; that is, it is easy to see small signal components in the presence of large ones. The decibel is a unit of ratio and is computed as follows.
where P is the measured power and P_{r} is the reference power.
Use the following equation to compute the ratio in decibels from amplitude values.
where A is the measured amplitude and A_{r} is the reference amplitude.
When using amplitude or power as the amplitudesquared of the same signal, the resulting decibel level is exactly the same. Multiplying the decibel ratio by two is equivalent to having a squared ratio. Therefore, you obtain the same decibel level and display regardless of whether you use the amplitude or power spectrum.
As shown in the preceding equations for power and amplitude, you must supply a reference for a measure in decibels. This reference then corresponds to the 0 dB level. Several conventions are used. A common convention is to use the reference 1 Vrms for amplitude or 1 Vrms squared for power, yielding a unit in dBV or dBVrms. In this case, 1 Vrms corresponds to 0 dB. Another common form of dB is dBm, which corresponds to a reference of 1 mW into a load of 50 for radio frequencies where 0 dB is 0.22 Vrms, or 600 for audio frequencies where 0 dB is 0.78 Vrms.
Refer to the Converting to Logarithmic Units topic in the LabVIEW Help (linked below) for the most updated information about converting to logarithmic units.
2. Antialiasing and Acquisition Front Ends for FFTBased Signal Analysis
FFTbased measurement requires digitization of a continuous signal. According to the Nyquist criterion, the sampling frequency, F_{s}, must be at least twice the maximum frequency component in the signal. If this criterion is violated, a phenomenon known as aliasing occurs. Figure 3 shows an adequately sampled signal and an undersampled signal. In the undersampled case, the result is an aliased signal that appears to be at a lower frequency than the actual signal.
Figure 3. Adequate and Inadequate Signal Sampling
When the Nyquist criterion is violated, frequency components above half the sampling frequency appear as frequency components below half the sampling frequency, resulting in an erroneous representation of the signal. For example, a component at frequency
appears as the frequency F_{s}  f_{0}.
Figure 4 shows the alias frequencies that appear when the signal with real components at 25, 70, 160, and 510 Hz is sampled at 100 Hz. Alias frequencies appear at 10, 30, and 40 Hz.
Figure 4. Alias Frequencies Resulting from Sampling a Signal at 100 Hz That Contains Frequency Components Greater than or Equal to 50 Hz
Before a signal is digitized, you can prevent aliasing by using antialiasing filters to attenuate the frequency components at and above half the sampling frequency to a level below the dynamic range of the analogtodigital converter (ADC). For example, if the digitizer has a fullscale range of 80 dB, frequency components at and above half the sampling frequency must be attenuated to 80 dB below full scale.
These higher frequency components, do not interfere with the measurement. If you know that the frequency bandwidth of the signal being measured is lower than half the sampling frequency, you can choose not to use an antialiasing filter. Figure 5 shows the input frequency response of the National Instruments PCI4450 Family dynamic signal acquisition boards, which have antialiasing filters. Note how an input signal at or above half the sampling frequency is severely attenuated.
Figure 5. Bandwidth of PCI4450 Family Input Versus Frequency, Normalized to Sampling Rate
Limitations of the Acquisition Front End
In addition to reducing frequency components greater than half the sampling frequency, the acquisition front end you use introduces some bandwidth limitations below half the sampling frequency. To eliminate signals at or above half of the sampling rate to less than the measurement range, antialiasing filters start to attenuate frequencies at some point below half the sampling rate. Because these filters attenuate the highest frequency portion of the spectrum, typically you want to limit the plot to the bandwidth you consider valid for the measurement.
For example, in the case of the PCI4450 Family sample shown in Figure 5, amplitude flatness is maintained to within ±0.1 dB, at up to 0.464 of the sampling frequency to 20 kHz for all gain settings, +1 dB to 95 kHz, and then the input gain starts to attenuate. The 3 dB point (or halfpower bandwidth) of the input occurs at 0.493 of the input spectrum. Therefore, instead of showing the input spectrum all the way out to half the sampling frequency, you may want to show only 0.464 of the input spectrum. To do this, multiply the number of points acquired by 0.464, respectively, to compute the number of frequency lines to display.
The characteristics of the signal acquisition front end affect the measurement. The National Instruments PCI4450 Family dynamic signal acquisition boards and the NI 4551 and NI 4552 dynamic signal analyzers are excellent acquisition front ends for performing FFTbased signal analysis measurements. These boards use deltasigma modulation technology, which yields excellent amplitude flatness, highperformance antialiasing filters, and wide dynamic range as shown in Figure 5. The input channels are also simultaneously sampled for good multichannel measurement performance.
At a sampling frequency of 51.2 kHz, these boards can perform frequency measurements in the range of DC to 23.75 kHz. Amplitude flatness is ±0.1 dB maximum from DC to 23.75 kHz. Refer to the PCI4451/4452/4453/4454 User Manual for more information about these boards.
Calculating the Measurement Bandwidth or Number of Lines for a Given Sampling Frequency
The dynamic signal acquisition boards have antialiasing filters built into the digitizing process. In addition, the cutoff filter frequency scales with the sampling rate to meet the Nyquist criterion as shown in Figure 5. The fast cutoff of the antialiasing filters on these boards means that the number of useful frequency lines in a 1,024point FFTbased spectrum is 475 lines for ±0.1 dB amplitude flatness.
To calculate the measurement bandwidth for a given sampling frequency, multiply the sampling frequency by 0.464 for the ±0.1 dB flatness. Also, the larger the FFT, the larger the number of frequency lines. A 2,048point FFT yields twice the number of lines listed above. Contrast this with typical benchtop instruments, which have 400 or 800 useful lines for a 1,024point or 2,048point FFT, respectively.
Dynamic Range Specifications
The signaltonoise ratio (SNR) of the PCI4450 Family boards is 93 dB. SNR is defined as
where V_{s} and V_{n} are the rms amplitudes of the signal and noise, respectively. A bandwidth is usually given for SNR. In this case, the bandwidth is the frequency range of the board input, which is related to the sampling rate as shown in Figure 5. The 93 dB SNR means that you can detect the frequency components of a signal that is as small as 93 dB below the fullscale range of the board. This is possible because the total input noise level caused by the acquisition front end is 93 dB below the fullscale input range of the board.
If the signal you monitor is a narrowband signal (that is, the signal energy is concentrated in a narrow band of frequencies), you are able to detect an even lower level signal than 93 dB. This is possible because the noise energy of the board is spread out over the entire input frequency range. Refer to the Computing Noise Level and Power Spectral Density section later in this application note for more information about narrowband versus broadband levels.
The spuriousfree dynamic range of the dynamic signal acquisition boards is 95 dB. Besides input noise, the acquisition front end may introduce spurious frequencies into a measured spectrum because of harmonic or intermodulation distortion, among other things. This 95 dB level indicates that any such spurious frequencies are at least 95 dB below the fullscale input range of the board.
The signaltototalharmonicdistortion (THD)plusnoise ratio, which excludes intermodulation distortion, is 90 dB from 0 to 20 kHz. THD is a measure of the amount of distortion introduced into a signal because of the nonlinear behavior of the acquisition front end. This harmonic distortion shows up as harmonic energy added to the spectrum for each of the discrete frequency components present in the input signal.
The wide dynamic range specifications of these boards is largely due to the 16bit resolution ADCs. Figure 6 shows a typical spectrum plot of the PCI4450 Family dynamic range with a fullscale 997 Hz signal applied. You can see that the harmonics of the 997 Hz input signal, the noise floor, and any other spurious frequencies are below 95 dB. In contrast, dynamic range specifications for benchtop instruments typically range from 70 dB to 80 dB using 12bit and 13bit ADC technology.
Figure 6. PCI4450 Family Spectrum Plot with 997 Hz Input at Full Scale (Full Scale = 0 dB) < /div > < /DIV >
See Also:
PCI4451/4452/4453/4454 User Manual
3. Using Windows Correctly
As mentioned in the Introduction, using windows correctly is critical to FFTbased measurement. This section describes the problem of spectral leakage, the characteristics of windows, some strategies for choosing windows, and the importance of scaling windows.
Spectral Leakage
For an accurate spectral measurement, it is not sufficient to use proper signal acquisition techniques to have a nicely scaled, singlesided spectrum. You might encounter spectral leakage. Spectral leakage is the result of an assumption in the FFT algorithm that the time record is exactly repeated throughout all time and that signals contained in a time record are thus periodic at intervals that correspond to the length of the time record. If the time record has a nonintegral number of cycles, this assumption is violated and spectral leakage occurs. Another way of looking at this case is that the nonintegral cycle frequency component of the signal does not correspond exactly to one of the spectrum frequency lines.
There are only two cases in which you can guarantee that an integral number of cycles are always acquired. One case is if you are sampling synchronously with respect to the signal you measure and can therefore deliberately take an integral number of cycles.
Another case is if you capture a transient signal that fits entirely into the time record. In most cases, however, you measure an unknown signal that is stationary; that is, the signal is present before, during, and after the acquisition. In this case, you cannot guarantee that you are sampling an integral number of cycles. Spectral leakage distorts the measurement in such a way that energy from a given frequency component is spread over adjacent frequency lines or bins. You can use windows to minimize the effects of performing an FFT over a nonintegral number of cycles.
Figure 7 shows the effects of three different windows  none (Uniform), Hanning (also commonly known as Hann), and Flat Top  when an integral number of cycles have been acquired, in this figure, 256 cycles in a 1,024point record. Notice that the windows have a main lobe around the frequency of interest. This main lobe is a frequency domain characteristic of windows. The Uniform window has the narrowest lobe, and the Hann and Flat Top windows introduce some spreading. The Flat Top window has a broader main lobe than the others. For an integral number of cycles, all windows yield the same peak amplitude reading and have excellent amplitude accuracy.
Figure 7 also shows the values at frequency lines of 254 Hz through 258 Hz for each window. The amplitude error at 256 Hz is 0 dB for each window. The graph shows the spectrum values between 240 and 272 Hz. The actual values in the resulting spectrum array for each window at 254 through 258 Hz are shown below the graph. f is 1 Hz.
Figure 7. Power Spectrum of 1 Vrms Signal at 256 Hz with Uniform, Hann, and Flat Top Windows
Figure 8 shows the leakage effects when you acquire 256.5 cycles. Notice that at a nonintegral number of cycles, the Hann and Flat Top windows introduce much less spectral leakage than the Uniform window. Also, the amplitude error is better with the Hann and Flat Top windows. The Flat Top window demonstrates very good amplitude accuracy but also has a wider spread and higher side lobes than the Hann window.
Figure 8. Power Spectrum of 1 Vrms Signal at 256.5 Hz with Uniform, Hann, and Flat Top Windows
In addition to causing amplitude accuracy errors, spectral leakage can obscure adjacent frequency peaks. Figure 9 shows the spectrum for two close frequency components when no window is used and when a Hann window is used.
Figure 9. Spectral Leakage Obscuring Adjacent Frequency Components
Window Characteristics
To understand how a given window affects the frequency spectrum, you need to understand more about the frequency characteristics of windows. The windowing of the input data is equivalent to convolving the spectrum of the original signal with the spectrum of the window as shown in Figure 10. Even if you use no window, the signal is convolved with a rectangularshaped window of uniform height, by the nature of taking a snapshot in time of the input signal. This convolution has a sine function characteristic spectrum. For this reason, no window is often called the Uniform or Rectangular window because there is still a windowing effect.
An actual plot of a window shows that the frequency characteristic of a window is a continuous spectrum with a main lobe and several side lobes. The main lobe is centered at each frequency component of the timedomain signal, and the side lobes approach zero at
Figure 10. Frequency Characteristics of a Windowed Spectrum
An FFT produces a discrete frequency spectrum. The continuous, periodic frequency spectrum is sampled by the FFT, just as the timedomain signal was sampled by the ADC. What appears in each frequency line of the FFT is the value of the continuous convolved spectrum at each FFT frequency line. This is sometimes referred to as the picketfence effect because the FFT result is analogous to viewing the continuous windowed spectrum through a picket fence with slits at intervals that corresponds to the frequency lines.
If the frequency components of the original signal match a frequency line exactly, as is the case when you acquire an integral number of cycles, you see only the main lobe of the spectrum. Side lobes do not appear because the spectrum of the window approaches zero at f intervals on either side of the main lobe. Figure 7 illustrates this case.
If a time record does not contain an integral number of cycles, the continuous spectrum of the window is shifted from the main lobe center at a fraction of f that corresponds to the difference between the frequency component and the FFT line frequencies. This shift causes the side lobes to appear in the spectrum. In addition, there is some amplitude error at the frequency peak, as shown in Figure 8, because the main lobe is sampled off center (the spectrum is smeared).
Figure 11 shows the frequency spectrum characteristics of a window in more detail. The side lobe characteristics of the window directly affect the extent to which adjacent frequency components bias (leak into) adjacent frequency bins. The side lobe response of a strong sinusoidal signal can overpower the main lobe response of a nearby weak sinusoidal signal.
Figure 11. Frequency Response of a Window
Another important characteristic of window spectra is main lobe width. The frequency resolution of the windowed signal is limited by the width of the main lobe of the window spectrum. Therefore, the ability to distinguish two closely spaced frequency components increases as the main lobe of the window narrows. As the main lobe narrows and spectral resolution improves, the window energy spreads into its side lobes, and spectral leakage worsens. In general, then, there is a trade off between leakage suppression and spectral resolution.
Refer to the Windowing Signals and Characteristics of Different Smoothing Windows topics in the LabVIEW Help (linked below) for the most updated information about window characteristics.
Defining Window Characteristics
To simplify choosing a window, you need to define various characteristics so that you can make comparisons between windows. Figure 11 shows the spectrum of a typical window. To characterize the main lobe shape, the 3 dB and 6 dB main lobe width are defined to be the width of the main lobe (in FFT bins or frequency lines) where the window response becomes 0.707 (3 dB) and 0.5 (6 dB), respectively, of the main lobe peak gain.
To characterize the side lobes of the window, the maximum side lobe level and side lobe rolloff rate are defined. The maximum side lobe level is the level in decibels relative to the main lobe peak gain, of the maximum side lobe. The side lobe rolloff rate is the asymptotic decay rate, in decibels per decade of frequency, of the peaks of the side lobes. Table 1 lists the characteristics of several window functions and their effects on spectral leakage and resolution.
Window

3 dB Main Lobe Width (bins)

6 dB Main Lobe Width (bins)

Maximum Side Lobe Level (dB)

Side Lobe RollOff Rate (dB/decade)

Uniform (None) 
0.89

1.21

13

20

Hanning (Hann) 
1.44

2.00

31

60

Hamming 
1.30

1.82

43

20

BlackmanHarris 
1.62

2.27

71

20

Exact Blackman 
1.61

2.25

68

20

Blackman 
1.64

2.30

58

60

Flat Top 
3.72

4.58

93

20

Strategies for Choosing Windows
Each window has its own characteristics, and different windows are used for different applications. To choose a spectral window, you must guess the signal frequency content. If the signal contains strong interfering frequency components distant from the frequency of interest, choose a window with a high side lobe rolloff rate. If there are strong interfering signals near the frequency of interest, choose a window with a low maximum side lobe level.
If the frequency of interest contains two or more signals very near to each other, spectral resolution is important. In this case, it is best to choose a window with a very narrow main lobe. If the amplitude accuracy of a single frequency component is more important than the exact location of the component in a given frequency bin, choose a window with a wide main lobe. If the signal spectrum is rather flat or broadband in frequency content, use the Uniform window (no window). In general, the Hann window is satisfactory in 95% of cases. It has good frequency resolution and reduced spectral leakage.
The Flat Top window has good amplitude accuracy, but because it has a wide main lobe, it has poor frequency resolution and more spectral leakage. The Flat Top window has a lower maximum side lobe level than the Hann window, but the Hann window has a faster roll offrate. If you do not know the nature of the signal but you want to apply a window, start with the Hann window. Figures 7 and 8 contrast the characteristics of the Uniform, Hann, and Flat Top windows.
If you are analyzing transient signals such as impact and response signals, it is better not to use the spectral windows because these windows attenuate important information at the beginning of the sample block. Instead, use the Force and Exponential windows. A Force window is useful in analyzing shock stimuli because it removes stray signals at the end of the signal. The Exponential window is useful for analyzing transient response signals because it damps the end of the signal, ensuring that the signal fully decays by the end of the sample block.
Selecting a window function is not a simple task. In fact, there is no universal approach for doing so. However, Table 2 can help you in your initial choice. Always compare the performance of different window functions to find the best one for the application. Refer to the references at the end of this application note for more information about windows.
Table 2. Initial Window Choice Based on Signal Content
Signal Content

Window

Sine wave or combination of sine waves  Hann 
Sine wave (amplitude accuracy is important)  Flat Top 
Narrowband random signal (vibration data)  Hann 
Broadband random (white noise)  Uniform 
Closely spaced sine waves  Uniform, Hamming 
Excitation signals (Hammer blow)  Force 
Response signals  Exponential 
Unknown content  Hann 
Windows are useful in reducing spectral leakage when using the FFT for spectral analysis. However, because windows are multiplied with the acquired timedomain signal, they introduce distortion effects of their own. The windows change the overall amplitude of the signal. The windows used to produce the plots in Figures 7 and 8 were scaled by dividing the windowed array by the coherent gain of the window. As a result, each window yields the same spectrum amplitude result within its accuracy constraints.
You can think of an FFT as a set of parallel filters, each f in bandwidth. Because of the spreading effect of a window, each window increases the effective bandwidth of an FFT bin by an amount known as the equivalent noisepower bandwidth of the window. The power of a given frequency peak is computed by adding the adjacent frequency bins around a peak and is inflated by the bandwidth of the window. You must take this inflation into account when you perform computations based on the spectrum. Refer to the Computations on the Spectrum section for sample computations.
Table 3 lists the scaling factor (or coherent gain), the noise power bandwidth, and the worstcase peak amplitude accuracy caused by offcenter components for several popular windows.
Table 3. Correction Factors and WorstCase Amplitude Errors for Windows
Window

Scaling Factor
(Coherent Gain) 
Noise Power
Bandwidth 
WorstCase Amplitude
Error (dB) 
Uniform (none) 
1.00

1.00

3.92

Hann 
0.50

1.50

1.42

Hamming 
0.54

1.36

1.75

BlackmanHarris 
0.42

1.71

1.13

Exact Blackman 
0.43

1.69

1.15

Blackman 
0.42

1.73

1.10

Flat Top 
0.22

3.77

< 0.01

4. Computations on the Spectrum
When you have the amplitude or power spectrum, you can compute several useful characteristics of the input signal, such as power and frequency, noise level, and power spectral density.
Estimating Power and Frequency
The preceding windowing examples demonstrate that if you have a frequency component in between two frequency lines, it appears as energy spread among adjacent frequency lines with reduced amplitude. The actual peak is between the two frequency lines. In Figure 8, the amplitude error at 256.5 Hz is due to the fact that the window is sampled at ±0.5 Hz around the center of its main lobe rather than at the center where the amplitude error would be 0. This is the picketfence effect explained in the Window Characteristics section of this application note.
You can estimate the actual frequency of a discrete frequency component to a greater resolution than the f given by the FFT by performing a weighted average of the frequencies around a detected peak in the power spectrum.
Similarly, you can estimate the power in Vrms2 of a given peak discrete frequency component by summing the power in the bins around the peak (computing the area under the peak)
Notice that this method is valid only for a spectrum made up of discrete frequency components. It is not valid for a continuous spectrum. Also, if two or more frequency peaks are within six lines of each other, they contribute to inflating the estimated powers and skewing the actual frequencies. You can reduce this effect by decreasing the number of lines spanned by the preceding computations. If two peaks are that close, they are probably already interfering with one another because of spectral leakage.
Similarly, if you want the total power in a given frequency range, sum the power in each bin included in the frequency range and divide by the noise power bandwidth of the windows.
Refer to the Computations on the Spectrum topic in the LabVIEW Help (linked below) for the most updated information about estimating power and frequency.
Computing Noise Level and Power Spectral Density
The measurement of noise levels depends on the bandwidth of the measurement. When looking at the noise floor of a power spectrum, you are looking at the narrowband noise level in each FFT bin. Thus, the noise floor of a given power spectrum depends on the f of the spectrum, which is in turn controlled by the sampling rate and number of points. In other words, the noise level at each frequency line reads as if it were measured through a f Hz filter centered at that frequency line. Therefore, for a given sampling rate, doubling the number of points acquired reduces the noise power that appears in each bin by 3 dB. Discrete frequency components theoretically have zero bandwidth and therefore do not scale with the number of points or frequency range of the FFT.
To compute the SNR, compare the peak power in the frequencies of interest to the broadband noise level. Compute the broadband noise level in Vrms^{2} by summing all the power spectrum bins, excluding any peaks and the DC component, and dividing the sum by the equivalent noise bandwidth of the window. For example, in Figure 6 the noise floor appears to be more than 120 dB below full scale, even though the PCI4450 Family dynamic range is only 93 dB. If you were to sum all the bins, excluding DC, and any harmonic or other peak components and divide by the noise power bandwidth of the window you used, the noise power level compared to full scale would be around 93 dB from full scale.
Because of noiselevel scaling with f, spectra for noise measurement are often displayed in a normalized format called power or amplitude spectral density. This normalizes the power or amplitude spectrum to the spectrum that would be measured by a 1 Hzwide square filter, a convention for noiselevel measurements. The level at each frequency line then reads as if it were measured through a 1 Hz filter centered at that frequency line.
Power spectral density is computed as
The spectral density format is appropriate for random or noise signals but inappropriate for discrete frequency components because the latter theoretically have zero bandwidth.
Refer to the Computations on the Spectrum topic in the LabVIEW Help (linked below) for the most updated information about computing noise levels and the power spectral density.
5. FFTBased Network Measurement
When you understand how to handle computations with the FFT and power spectra, and you understand the influence of windows on the spectrum, you can compute several FFTbased functions that are extremely useful for network analysis. These include the frequency response, impulse response, and coherence functions. Refer to the Frequency Response and Network Analysis section of this application note for more information about these functions. Refer to the Signal Sources for Frequency Response Measurement section for more information about Chirp signals and broadband noise signals.
Cross Power Spectrum
One additional building block is the cross power spectrum. The cross power spectrum is not typically used as a direct measurement but is an important building block for other measurements.
The twosided cross power spectrum of two timedomain signals A and B is computed as
The cross power spectrum is in twosided complex form. To convert to magnitude and phase, use the RectangularToPolar conversion function. To convert to a singlesided form, use the same method described in the Converting from a TwoSided Power Spectrum to a SingleSided Power Spectrum section of this application note. The units of the singlesided form are in volts (or other quantity) rms squared.
The power spectrum is equivalent to the cross power spectrum when signals A and B are the same signal. Therefore, the power spectrum is often referred to as the auto power spectrum or the auto spectrum. The singlesided cross power spectrum yields the product of the rms amplitudes of the two signals, A and B, and the phase difference between the two signals.
When you know how to use these basic blocks, you can compute other useful functions, such as the Frequency Response function.
Refer to the Cross Power Spectrum conceptual topic in the LabVIEW Help (linked below) for the most updated information about the cross power spectrum.
Frequency Response and Network Analysis
Three useful functions for characterizing the frequency response of a network are the frequency response, impulse response, and coherence functions.
The frequency response of a network is measured by applying a stimulus to the network as shown in Figure 12 and computing the frequency response from the stimulus and response signals.
Figure 12. Configuration for Network Analysis
Frequency Response Function
The frequency response function (FRF) gives the gain and phase versus frequency of a network and is typically computed as
where A is the stimulus signal and B is the response signal.
The frequency response function is in twosided complex form. To convert to the frequency response gain (magnitude) and the frequency response phase, use the RectangularToPolar conversion function. To convert to singlesided form, simply discard the second half of the array.
You may want to take several frequency response function readings and then average them. To do so, average the cross power spectrum, S_{AB}(f), by summing it in the complex form then dividing by the number of averages, before converting it to magnitude and phase, and so forth. The power spectrum, S_{AA}(f), is already in real form and is averaged normally.
Refer to the Frequency Response and Network Analysis topic in the LabVIEW Help (linked below) for the most updated information about the frequency response function.
Impulse Response Function
The impulse response function of a network is the timedomain representation of the frequency response function of the network. It is the output timedomain signal generated when an impulse is applied to the input at time t = 0.
To compute the impulse response of the network, take the inverse FFT of the twosided complex frequency response function as described in the Frequency Response Function section of this application note.
Refer to the Frequency Response and Network Analysis topic in the LabVIEW Help (linked below) for the most updated information about the impulse response function.
Coherence Function
The coherence function is often used in conjunction with the frequency response function as an indication of the quality of the frequency response function measurement and indicates how much of the response energy is correlated to the stimulus energy. If there is another signal present in the response, either from excessive noise or from another signal, the quality of the network response measurement is poor. You can use the coherence function to identify both excessive noise and causality, that is, identify which of the multiple signal sources are contributing to the response signal. The coherence function is computed as
The result is a value between zero and one versus frequency. A zero for a given frequency line indicates no correlation between the response and the stimulus signal. A one for a given frequency line indicates that the response energy is 100 percent due to the stimulus signal; in other words, there is no interference at that frequency.
For a valid result, the coherence function requires an average of two or more readings of the stimulus and response signals. For only one reading, it registers unity at all frequencies. To average the cross power spectrum, S_{AB}(f), average it in the complex form then convert to magnitude and phase as described in the Frequency Response Function section of this application note. The auto power spectra, S_{AA}(f) and S_{BB}(f), are already in real form, and you average them normally.
Refer to the Frequency Response and Network Analysis topic in the LabVIEW Help (linked below) for the most updated information about the coherence function.
Signal Sources for Frequency Response Measurements
To achieve a good frequency response measurement, significant stimulus energy must be present in the frequency range of interest. Two common signals used are the chirp signal and a broadband noise signal. The chirp signal is a sinusoid swept from a start frequency to a stop frequency, thus generating energy across a given frequency range. White and pseudorandom noise have flat broadband frequency spectra; that is, energy is present at all frequencies.
It is best not to use windows when analyzing frequency response signals. If you generate a chirp stimulus signal at the same rate you acquire the response, you can match the acquisition frame size to match the length of the chirp. No window is generally the best choice for a broadband signal source. Because some stimulus signals are not constant in frequency across the time record, applying a window may obscure important portions of the transient response.
6. Conclusion
There are many issues to consider when analyzing and measuring signals from plugin DAQ devices. Unfortunately, it is easy to make incorrect spectral measurements. Understanding the basic computations involved in FFTbased measurement, knowing how to prevent antialiasing, properly scaling and converting to different units, choosing and using windows correctly, and learning how to use FFTbased functions for network measurement are all critical to the success of analysis and measurement tasks. Being equipped with this knowledge and using the tools discussed in this application note can bring you more success with your individual application.
7. Related Links
LabVIEW Help: Power Spectrum
LabVIEW Help: Computing the Amplitude and Phase Spectrums
LabVIEW Help: Converting to Logarithmic Units
LabVIEW Help: Spectral Leakage
LabVIEW Help: Windowing Signals
LabVIEW Help: Characteristics of Different Smoothing Windows
LabVIEW Help: Computations on the Spectrum
LabVIEW Help: Cross Power Spectrum
LabVIEW Help: Frequency Response and Network Analysis
8. References
Harris, Fredric J. "On the Use of Windows for Harmonic Analysis with the Discrete Fourier Transform" in Proceedings of the IEEE Vol. 66, No. 1, January 1978.
Horowitz, Paul, and Hill, Winfield, The Art of Electronics, 2nd Edition, Cambridge University Press, 1989.
Nuttall, Albert H. "Some Windows with Very Good Sidelobe Behavior," IEEE Transactions on Acoustics, Speech, and Signal Processing Vol. 29, No. 1, February 1981.
Randall, R.B., and Tech, B. Frequency Analysis, 3rd Edition, Bruël and Kjær, September 1979.
The Fundamentals of Signal Analysis, Application Note 243, HewlettPackard, 1985.