Labview Sound and Vibration
Labview Sound and Vibration
TM
Worldwide Offices
Australia 1800 300 800, Austria 43 0 662 45 79 90 0, Belgium 32 0 2 757 00 20, Brazil 55 11 3262 3599,
Canada (Calgary) 403 274 9391, Canada (Ottawa) 613 233 5949, Canada (Québec) 450 510 3055,
Canada (Toronto) 905 785 0085, Canada (Vancouver) 514 685 7530, China 86 21 6555 7838,
Czech Republic 420 224 235 774, Denmark 45 45 76 26 00, Finland 385 0 9 725 725 11,
France 33 0 1 48 14 24 24, Germany 49 0 89 741 31 30, Greece 30 2 10 42 96 427, India 91 80 51190000,
Israel 972 0 3 6393737, Italy 39 02 413091, Japan 81 3 5472 2970, Korea 82 02 3451 3400,
Malaysia 603 9131 0918, Mexico 001 800 010 0793, Netherlands 31 0 348 433 466,
New Zealand 0800 553 322, Norway 47 0 66 90 76 60, Poland 48 22 3390150, Portugal 351 210 311 210,
Russia 7 095 783 68 51, Singapore 65 6226 5886, Slovenia 386 3 425 4200, South Africa 27 0 11 805 8197,
Spain 34 91 640 0085, Sweden 46 0 8 587 895 00, Switzerland 41 56 200 51 51, Taiwan 886 2 2528 7227,
Thailand 662 992 7519, United Kingdom 44 0 1635 523545
For further support information, refer to the Technical Support and Professional Services appendix. To comment
on the documentation, send email to [email protected].
Warranty
The media on which you receive National Instruments software are warranted not to fail to execute programming instructions, due to defects
in materials and workmanship, for a period of 90 days from date of shipment, as evidenced by receipts or other documentation. National
Instruments will, at its option, repair or replace software media that do not execute programming instructions if National Instruments receives
notice of such defects during the warranty period. National Instruments does not warrant that the operation of the software shall be
uninterrupted or error free.
A Return Material Authorization (RMA) number must be obtained from the factory and clearly marked on the outside of the package before
any equipment will be accepted for warranty work. National Instruments will pay the shipping costs of returning to the owner parts which are
covered by warranty.
National Instruments believes that the information in this document is accurate. The document has been carefully reviewed for technical
accuracy. In the event that technical or typographical errors exist, National Instruments reserves the right to make changes to subsequent
editions of this document without prior notice to holders of this edition. The reader should consult National Instruments if errors are suspected.
In no event shall National Instruments be liable for any damages arising out of or related to this document or the information contained in it.
EXCEPT AS SPECIFIED HEREIN, NATIONAL INSTRUMENTS MAKES NO WARRANTIES, EXPRESS OR IMPLIED, AND SPECIFICALLY DISCLAIMS ANY WARRANTY OF
MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. CUSTOMER’S RIGHT TO RECOVER DAMAGES CAUSED BY FAULT OR NEGLIGENCE ON THE PART OF
NATIONAL INSTRUMENTS SHALL BE LIMITED TO THE AMOUNT THERETOFORE PAID BY THE CUSTOMER. NATIONAL INSTRUMENTS WILL NOT BE LIABLE FOR
DAMAGES RESULTING FROM LOSS OF DATA, PROFITS, USE OF PRODUCTS, OR INCIDENTAL OR CONSEQUENTIAL DAMAGES, EVEN IF ADVISED OF THE POSSIBILITY
THEREOF. This limitation of the liability of National Instruments will apply regardless of the form of action, whether in contract or tort, including
negligence. Any action against National Instruments must be brought within one year after the cause of action accrues. National Instruments
shall not be liable for any delay in performance due to causes beyond its reasonable control. The warranty provided herein does not cover
damages, defects, malfunctions, or service failures caused by owner’s failure to follow the National Instruments installation, operation, or
maintenance instructions; owner’s modification of the product; owner’s abuse, misuse, or negligent acts; and power failure or surges, fire,
flood, accident, actions of third parties, or other events outside reasonable control.
Copyright
Under the copyright laws, this publication may not be reproduced or transmitted in any form, electronic or mechanical, including photocopying,
recording, storing in an information retrieval system, or translating, in whole or in part, without the prior written consent of National
Instruments Corporation.
Trademarks
LabVIEW™, National Instruments™, NI™, ni.com™, and NI-DAQ™ are trademarks of National Instruments Corporation.
Product and company names mentioned herein are trademarks or trade names of their respective companies.
Patents
For patents covering National Instruments products, refer to the appropriate location: Help»Patents in your software, the patents.txt file
on your CD, or ni.com/patents.
Chapter 1
Introduction
Sound and Vibration Toolkit .........................................................................................1-1
Toolkit Palettes ..............................................................................................................1-2
Scaling .............................................................................................................1-3
Calibration .......................................................................................................1-3
Limit Testing ...................................................................................................1-3
Weighting ........................................................................................................1-3
Integration........................................................................................................1-4
Generation .......................................................................................................1-4
Vibration Level................................................................................................1-4
Sound Level.....................................................................................................1-4
Octave Analysis...............................................................................................1-4
Frequency Analysis .........................................................................................1-4
Transient Analysis ...........................................................................................1-5
Waterfall Display.............................................................................................1-5
Swept Sine .......................................................................................................1-5
Distortion.........................................................................................................1-5
Single-Tone .....................................................................................................1-6
Front Panel Displays ......................................................................................................1-6
Examples........................................................................................................................1-6
Chapter 2
Dynamic Signals
Acquiring and Simulating Dynamic Signals .................................................................2-1
Aliasing............................................................................................................2-3
Time Continuity...............................................................................................2-4
Chapter 3
Scaling and Calibration
Scaling to EU .................................................................................................................3-1
Performing System Calibration .....................................................................................3-3
Propagation Delay Calibration ........................................................................3-3
© National Instruments Corporation v LabVIEW Sound and Vibration Toolkit User Manual
Contents
Chapter 4
Limit Testing Analysis
Limit Testing Overview ................................................................................................ 4-1
Using the SVT Limit Testing VI ................................................................................... 4-3
Chapter 5
Weighting Filters
Purpose of Weighting Filters......................................................................................... 5-1
Psophometric Weighting Filters...................................................................... 5-3
A-, B-, and C-Weighting Filters ....................................................... 5-3
Radiocommunications Weighting Filters ......................................... 5-5
Telecommunications Weighting Filters............................................ 5-6
Applying Weighting Filters........................................................................................... 5-6
Applying Weighting to Time-Domain Data ................................................... 5-8
Standards Compliance ................................................................................................... 5-9
A-, B-, and C-Weighting Filters...................................................................... 5-9
ANSI Standards ................................................................................ 5-9
ISO/IEC Standard ............................................................................. 5-10
Radiocommunications Weighting Filters........................................................ 5-10
Telecommunications Weighting Filters .......................................................... 5-10
Performing A-Weighted Sound Level Measurements...................... 5-11
Applying Weighting to an Octave Spectrum .................................................. 5-12
Errors Due to Uniform Corrections .................................................. 5-12
Applying Weighting to an FFT-Based Spectrum ........................................... 5-13
Chapter 6
Integration
Introduction to Integration............................................................................................. 6-1
Implementing Integration .............................................................................................. 6-3
Challenges When Integrating Vibration Data ................................................. 6-5
DC Component ................................................................................. 6-5
Transducers....................................................................................... 6-5
Implementing Integration using the Sound and Vibration Toolkit ................. 6-5
Time-Domain Integration.............................................................................................. 6-6
Single-Shot Acquisition and Integration......................................................... 6-6
Continuous Acquisition and Integration ......................................................... 6-7
Frequency-Domain Integration ..................................................................................... 6-11
Chapter 7
Vibration-Level Measurements
Measuring the Root Mean Square (RMS) Level ...........................................................7-2
Single-Shot Buffered Acquisition ...................................................................7-3
Continuous Signal Acquisition........................................................................7-3
Performing a Running RMS Level Measurement .........................................................7-4
Computing the Peak Level.............................................................................................7-4
Computing the Crest Factor ...........................................................................................7-5
Chapter 8
Sound-Level Measurements
Time Averaging Modes .................................................................................................8-2
Linear Averaging.............................................................................................8-3
Single-Shot Linear Averaging ..........................................................8-3
Measuring Leq Over a Longer Time Period ......................................8-4
Restart Averaging and Advanced Concepts......................................8-5
Performing a Running Leq .................................................................8-6
Exponential Averaging....................................................................................8-7
Peak Hold ........................................................................................................8-8
Considerations for Making Sound-Level Measurements ..............................................8-8
Chapter 9
Fractional-Octave Analysis
Fractional-Octave Analysis Overview...........................................................................9-2
Full-Octave Analysis in the 31.5 Hz–16 kHz Band ........................................9-3
Bandwidth and Filter Banks ..........................................................................................9-4
The Octave Filter.............................................................................................9-4
Bandedge Frequencies.....................................................................................9-5
Fractional-Octave Filters .................................................................................9-6
Filter Settling Time..........................................................................................9-7
Averaging.......................................................................................................................9-8
Linear Averaging.............................................................................................9-8
Exponential Averaging....................................................................................9-8
Equal Confidence Averaging ..........................................................................9-9
Peak-Hold Averaging ......................................................................................9-9
Resetting the Filter and Restarting the Averaging Process .............................9-9
Performing Third-Octave Analysis Outside the Audio Range ......................................9-9
ANSI and IEC Standards ...............................................................................................9-10
ANSI Standard.................................................................................................9-10
IEC Standard ...................................................................................................9-11
© National Instruments Corporation vii LabVIEW Sound and Vibration Toolkit User Manual
Contents
Chapter 10
Frequency Analysis
FFT Fundamentals......................................................................................................... 10-2
Number of Samples......................................................................................... 10-3
Frequency Resolution ..................................................................................... 10-3
Maximum Resolvable Frequency ..................................................... 10-4
Minimum Resolvable Frequency...................................................... 10-4
Number of Spectral Lines................................................................. 10-4
Relationship between Time-Domain and Frequency-Domain
Specifications and Parameters...................................................................... 10-4
Increasing Frequency Resolution .................................................................................. 10-6
Zoom FFT Analysis ........................................................................................ 10-8
Frequency Resolution of the Zoom FFT VIs.................................... 10-9
Zoom Measurement .......................................................................... 10-10
Zoom Settings................................................................................... 10-11
Subset Analysis ............................................................................................... 10-11
Using the Frequency Analysis VIs ................................................................................ 10-12
Available Measurements................................................................................. 10-12
Single-Channel Measurements ....................................................................... 10-13
Power Spectrum Measurement ......................................................... 10-14
Dual-Channel Measurements .......................................................................... 10-15
Frequency Response Function Measurement ................................... 10-15
Windowing .................................................................................................................... 10-21
Averaging Parameters ................................................................................................... 10-22
Special Considerations for Averaged Measurements ..................................... 10-23
Averaging Mode ............................................................................................. 10-23
No Averaging ................................................................................... 10-24
RMS Averaging ................................................................................ 10-24
Vector Averaging ............................................................................. 10-25
RMS versus Vector Averaging......................................................... 10-25
Peak Hold ......................................................................................... 10-27
Weighting Mode ............................................................................................. 10-27
Coherence and Coherent Output Power ........................................................................ 10-28
Extended Measurements................................................................................................ 10-28
Power in Band ................................................................................................. 10-29
Spectrum Peak Search..................................................................................... 10-29
Unit Conversion .............................................................................................. 10-30
Chapter 11
Transient Analysis
Transient Analysis with the Sound and Vibration Toolkit ............................................11-2
Performing an STFT versus Time .................................................................................11-2
Selecting the FFT Block Size ..........................................................................11-5
Overlapping .....................................................................................................11-6
Using the SVT STFT versus Time VI.............................................................11-8
Performing an STFT versus Rotational Speed ..............................................................11-9
Converting the Pulse Train to Rotational Speed .............................................11-9
STFT versus RPM ...........................................................................................11-10
Measuring a Shock Response Spectrum ........................................................................11-12
Chapter 12
Waterfall Display
Using the Display VIs....................................................................................................12-1
Initializing the Display ....................................................................................12-2
Sending Data to the Display ............................................................................12-2
Waterfall Display for Frequency Analysis........................................12-3
Waterfall Display for Transient Analysis .........................................12-3
Waterfall Display for Octave Spectra ...............................................12-5
Customizing the Waterfall Display View........................................................12-6
Closing the Waterfall Display .........................................................................12-6
Chapter 13
Swept-Sine Measurements
Swept Sine Overview.....................................................................................................13-3
Choosing Swept-Sine versus FFT Measurements .........................................................13-4
Taking a Swept Sine Measurement ...............................................................................13-6
Swept Sine Measurement Example ...............................................................................13-7
Chapter 14
Distortion Measurements
Variable Definitions.......................................................................................................14-1
Signal in Noise and Distortion (SINAD) .......................................................................14-3
Total Harmonic Distortion Plus Noise (THD+N)..........................................................14-5
Total Harmonic Distortion (THD) .................................................................................14-6
Intermodulation Distortion (IMD) .................................................................................14-9
Phase Linearity ..............................................................................................................14-12
© National Instruments Corporation ix LabVIEW Sound and Vibration Toolkit User Manual
Contents
Chapter 15
Single-Tone Measurements
Single-Tone Measurement Overview............................................................................ 15-1
Gain and Phase .............................................................................................................. 15-4
Crosstalk ........................................................................................................................ 15-5
Gain ............................................................................................................................... 15-6
Idle-Channel Noise........................................................................................................ 15-7
Dynamic Range ............................................................................................................. 15-7
Spurious Free Dynamic Range (SFDR) ........................................................................ 15-9
Appendix A
References
Appendix B
Technical Support and Professional Services
Glossary
Index
Conventions
The following conventions appear in this manual:
» The » symbol leads you through nested menu items and dialog box options
to a final action. The sequence File»Page Setup»Options directs you to
pull down the File menu, select the Page Setup item, and select Options
from the last dialog box.
bold Bold text denotes items that you must select or click on in the software,
such as menu items and dialog box options. Bold text also denotes
parameter names, controls and buttons on the front panel, dialog boxes,
menu names, and palette names.
monospace bold Bold text in this font denotes the messages and responses that the computer
automatically prints to the screen. This font also emphasizes lines of code
that are different from the other examples.
monospace Text in this font denotes text or characters that you should enter from the
keyboard, paths, directories, variables, and filenames and extensions.
© National Instruments Corporation xi LabVIEW Sound and Vibration Toolkit User Manual
About This Manual
Related Documentation
The following documents contain information that you might find helpful
as you read this manual:
• LabVIEW Help, available by selecting Help»VI, Function,
& How-To Help in LabVIEW
• Getting Started with LabVIEW
• LabVIEW User Manual
• LabVIEW Measurements Manual
You can use the Sound and Vibration Toolkit to perform measurements on
digitized or simulated data.
Note In Figure 1-1, the measurement operations shown on the Analysis line are not
necessarily performed simultaneously. The dashed boxes in Figure 1-1 indicate optional
measurement operations.
© National Instruments Corporation 1-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 1 Introduction
DSA Device
DAQ Device
Data Source WAV File
DAT Recorder
Simulation
Calibration
Calibrate
Sensor
Scale Voltage to
Scaling Engineering Units
Measure
Propagation Delay
Limit Testing
Analysis
Vibration Level Extended Transient
Weighting
Measurements Measurements Analysis
Limit Testing
Toolkit Palettes
Installing the Sound and Vibration Toolkit adds Sound & Vibration
palettes to both the LabVIEW Functions and Controls palettes. This
section briefly introduces the different palettes that compose the Sound
and Vibration Toolkit.
All the high-level VIs in the Sound and Vibration Toolkit are designed to
offer measurement capabilities. The high-level VIs perform the selected
analysis and allow you to view the results with the appropriate engineering
units in standard displays, such as magnitude/phase, real/imaginary part, or
decibels on/off.
Scaling
The SVL Scale Voltage to EU VI allows you to scale the original signal to
engineering units. The SVL Scale Voltage to EU VI is part of the Sound
and Vibration Library (SVL). The SVL is a collection of VIs shared by the
Sound and Vibration Toolkit and other National Instruments (NI) toolkits.
Refer to Chapter 3, Scaling and Calibration, and to the LabVIEW Help for
more information about the SVL Scale Voltage to EU VI.
Calibration
The Calibration VIs allow you to perform an end-to-end calibration on
a selected channel and measure the propagation delay of the measurement
device. The Calibration VIs are part of the SVL. Refer to Chapter 3,
Scaling and Calibration, for information about the calibration process.
Refer to the LabVIEW Help for information about the individual
Calibration VIs.
Limit Testing
The SVT Limit Testing VI allows you to apply limit analysis to any type of
measured result produced by the Sound and Vibration Toolkit. Refer to
Chapter 4, Limit Testing Analysis, for more information about using limit
testing to analyze measurement results. Refer to the LabVIEW Help for
information about the SVT Limit Testing VI.
Weighting
The Weighting VIs allow you to apply A-, B-, or C-weighting filters on
the time-domain signal. Additionally, ITU-R 468-4 and Dolby filters are
available for radiocommunications applications, and C-message and
CCITT filters are available for telecommunications applications. Refer to
Chapter 5, Weighting Filters, for more information about applying
weighting to a signal. Refer to the LabVIEW Help for information about the
individual Weighting VIs.
© National Instruments Corporation 1-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 1 Introduction
Integration
The SVT Integration VI allows you to perform single or double integration
on the time-domain signal. Refer to Chapter 6, Integration, for information
about the integration process. Refer to the LabVIEW Help for information
about the SVT Integration VI.
Generation
The SVT Pink Noise Waveform VI allows you to generate a continuous
pink noise waveform. The Generation palette also contains a subpalette
linked to the standard Waveform Generation palette in LabVIEW. Refer
to the LabVIEW Help for information about the Generation VIs.
Vibration Level
The Vibration Level VIs offer level measurements typically used for
vibration measurements, including measuring the crest factor. Averaging
modes include RMS averaging, exponential averaging, and peak hold.
Refer to Chapter 7, Vibration-Level Measurements, for more information
about performing vibration level measurements. Refer to the LabVIEW
Help for information about the individual Vibration Level VIs.
Sound Level
The Sound Level VIs offer typical sound-level measurements, including
equivalent continuous averaging (Leq), exponential averaging, and peak
hold. Refer to Chapter 8, Sound-Level Measurements, for information
about performing sound-level measurements. Refer to the LabVIEW Help
for information about the individual Sound Level VIs.
Octave Analysis
The Octave Analysis VIs offer a set of tools to perform fractional-octave
analysis, including 1/1, 1/3, 1/6, 1/12, and 1/24 octave-band analysis. The
Octave Analysis VIs can accommodate any sampling frequency and any
number of fractional-octave bands. Refer to Chapter 9, Fractional-Octave
Analysis, for information about performing octave analyses. Refer to
the LabVIEW Help for information about the individual Octave
Analysis VIs.
Frequency Analysis
The Frequency Analysis VIs are a collection of frequency-analysis tools
based on the discrete Fourier transform (DFT) and the fast Fourier
transform (FFT). The Frequency Analysis VIs also provide zoom FFT
frequency measurements and extended measurements. Refer
to Chapter 10, Frequency Analysis, for information about performing
frequency analyses. Refer to the LabVIEW Help for information about the
individual Frequency Analysis VIs.
Transient Analysis
The Transient Analysis VIs offer two techniques for obtaining information
about transient signals. Use the short-time Fourier transform (STFT) to
extract frequency information as a function of time or rotational speed. Use
the shock response spectrum (SRS) to evaluate the severity of a shock
signal. Refer to Chapter 11, Transient Analysis, for information about
performing transient analyses. Refer to the LabVIEW Help for information
about the individual Transient Analysis VIs.
Waterfall Display
The Waterfall Display VIs allow you to display the results of frequency
analyses and octave analyses as waterfall graphs. The Waterfall Display
VIs generate and manage the external window of the waterfall graph. Refer
to Chapter 12, Waterfall Display, for information about displaying analysis
results in a waterfall graph. Refer to the LabVIEW Help for information
about the individual Waterfall Display VIs.
Swept Sine
The Swept-Sine VIs allow you to characterize the frequency response of a
device under test (DUT). The swept-sine measurements include dynamic
measurements for stimulus level, response level, frequency response (gain
and phase), total harmonic distortion (THD), and individual harmonic
distortion. Refer to Chapter 13, Swept-Sine Measurements, for information
about performing swept-sine measurements. Refer to the LabVIEW Help
for information about the individual Swept-Sine VIs.
Distortion
The Distortion VIs allow you to measure the harmonic, intermodulation,
and broadband noise components due to nonlinearities in the DUT. Refer
to Chapter 14, Distortion Measurements, for information about performing
distortion analyses. Refer to the LabVIEW Help for information about the
individual Distortion VIs.
© National Instruments Corporation 1-5 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 1 Introduction
Single-Tone
The Single-Tone VIs allow you to perform single-tone measurements,
defined as a group of measurements where the excitation is a single tone.
These measurements are often used to measure the linear response,
nonlinear distortion, and noise of audio devices. Refer to Chapter 15,
Single-Tone Measurements, for information about performing single-tone
measurements. Refer to the LabVIEW Help for information about the
individual Single-Tone VIs.
Examples
The Sound and Vibration Toolkit includes examples to help you get started
using the toolkit. Select Help»Find Examples in LabVIEW to launch the
NI Example Finder. Select Toolkits and Modules»Sound and Vibration
in the Browse tab to view all of the available examples, or use the Search
tab to locate a specific example. The examples demonstrate the following
Sound and Vibration Toolkit capabilities:
• Display
• Frequency analysis
• Integration
• Level measurements
• Octave analysis
• Scaling
• Transient analysis
• Weighting filters
• Swept-Sine measurements
• Audio measurements
• Limit testing
Refer to the LabVIEW Help for more information about acquiring and
simulating data.
Figure 2-1 illustrates how the data source, either acquired or simulated, fits
into the sound and vibration measurement process.
© National Instruments Corporation 2-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 2 Dynamic Signals
DSA Device
DAQ Device
Data Source WAV File
DAT Recorder
Simulation
Calibration
Calibrate
Sensor
Scale Voltage to
Scaling Engineering Units
Measure
Propagation
Delay
Limit Testing
Limit Testing
Whether you are obtaining the data from a DAQ system, reading the data
from a file, or simulating the data, aliasing and time continuity are common
issues which you should consider in your measurement analysis.
Aliasing
When a dynamic signal is discretely sampled, aliasing is the phenomenon
in which frequency components greater than the Nyquist frequency are
erroneously shifted to lower frequencies. The Nyquist frequency is
calculated with the following formula:
Simulated data also can exhibit aliasing. The signals often are generated
according to a time-domain expression and, therefore, have high-frequency
components that are aliased in the discretely sampled data. Figure 2-2
shows an example of this aliasing for a simulated square wave.
© National Instruments Corporation 2-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 2 Dynamic Signals
The only way to protect data from aliasing is to apply appropriate aliasing
protection before the data are generated or acquired. Aliasing occurs when
the data are generated or sampled, and it is not possible to remove aliased
components from the data without detailed knowledge of the original
signal. In general, it is not possible to distinguish between true frequency
components and aliased frequency components. Therefore, accurate
frequency measurements require adequate alias protection.
Time Continuity
When you acquire data in a continuous acquisition, you can use the t0
parameter in the waveform datatype to ensure there are no gaps between
successive blocks of waveforms returned by sequential calls to the DAQmx
Read VI or AI Read VI. When signals are generated with one of the
waveform generation VIs in the Generation palette or the Waveform
Generation palette, the t0 of the current waveform is one sample period
later than the timestamp of the last sample in the previous waveform.
Continuity is enforced in this way until the generation is reset.
The waveform datatype is integral for testing time continuity in the Sound
and Vibration Toolkit. If you read data from a file or simulate a signal using
one of the VIs in the Signal Generation palette, wire a t0 that meets the
continuous timestamp condition to the waveform datatype connected to the
measurement analysis VIs. This action prevents unexpected resets of the
measurement analysis due to detected discontinuities in the input signal.
© National Instruments Corporation 2-5 LabVIEW Sound and Vibration Toolkit User Manual
Scaling and Calibration
3
This chapter discusses using the SVL Scale Voltage to EU VI located on
the Scaling palette to scale a signal to engineering units (EU) and using the
Calibration VIs located on the Calibration palette.
Refer to the LabVIEW Help for more information about the SVL Scale
Voltage to EU VI and the Calibration VIs.
Scaling to EU
This section discusses scaling data to the appropriate EU so you can
perform measurement analysis.
Figure 3-1 illustrates how scaling and calibration fit into the sound and
vibration measurement process.
© National Instruments Corporation 3-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 3 Scaling and Calibration
DSA Device
DAQ Device
Data Source WAV File
DAT Recorder
Simulation
Calibration
Calibrate
Sensor
Scale Voltage to
Scaling Engineering Units
Measure
Propagation
Delay
Limit Testing
Limit Testing
Figure 3-1. Relationship of Scaling to the Sound and Vibration Measurement Process
All measurement VIs in the Sound and Vibration Toolkit expect input
signals and return results with the appropriate units, such as time-domain
signals in the correct EU, frequency spectra in decibels with the proper
reference, phase information in degrees or radians, and so on. To handle
units properly, the high-level VIs need the signal to be scaled to the
appropriate EU.
Note If you use any method outside of the Sound and Vibration Toolkit to apply scaling
to a waveform, do not use the SVL Scale Voltage to EU VI. NI provides several tools and
methods to apply scaling to a waveform. These include, but are not limited to, NI-DAQmx
tasks or global channels created with Measurement & Automation Explorer (MAX), the
DAQ Assistant, or the DAQmx Create Virtual Channel VI.
There are two ways to determine the propagation delay of the DSA device.
You can refer to the documentation for the DSA device to find the
propagation delay specifications, also referred to as group delay. You also
can measure the propagation delay in samples with the SVL Measure
Propagation Delay VIs. The SVL Measure Propagation Delay VIs allow
you to measure the delay introduced in the input and output circuitry for a
© National Instruments Corporation 3-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 3 Scaling and Calibration
specific device at the desired sample rate. Connect the DSA device output
channel directly to the input channel, as displayed in Figure 3-2, to measure
the device propagation delay.
PFI0
AI0
AI1
AO0
AO1
Note Do not put a DUT in the signal path when measuring the propagation delay for the
DAQ device.
For an E or S Series DAQ device from NI, you should expect to measure a
one-sample propagation delay due to the time required for the signal to
traverse the signal path between the D/A converter (DAC) on the analog
output channel and the A/D converter (ADC) on the analog input channel.
Figure 3-3 shows the time domain data for the propagation delay
measurement of an NI PCI-6052E.
For DSA devices, or any other device which has onboard filtering on either
the input, output, or both channels, you should expect to measure a
propagation delay consistent with the sum of the delays specified for the
onboard filters on the input and output channels. Figure 3-4 shows the
delay of a smooth pulse generated and acquired by an NI PXI-4461 with a
204.8 kHz sample rate.
Figure 3-4. NI PXI-4461 Propagation Delay with a 204.8 kHz Sample Rate
Not all DSA devices have a constant propagation delay across the entire
range of supported sample rates. For example, the NI PXI-4461
propagation delay is dependent on the output update rate. Figure 3-5 shows
the total propagation delay versus sample rate relationship for the
NI PXI-4461 from output to input as a function of the sample rate.
© National Instruments Corporation 3-5 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 3 Scaling and Calibration
As illustrated by Figures 3-3, 3-4, and 3-5, the propagation delay can vary
significantly with different sample rates and devices. To ensure
measurement accuracy in your I/O applications, determine and account for
the propagation delay of the DAQ device at the same sample rate used in
your application.
It is important to remove the effects of the delay due to the data acquisition
system for two reasons. First, there is always a delay between the generated
output signal and the acquired input on the device even when the output and
input channels are hardware synchronized. Second, the anti-imaging and
anti-aliasing filters of the device introduce additional delays. You must
account for this delay to perform accurate dynamic measurements. Use the
device propagation delay [samples] input on the examples found in the
LabVIEW program directory under \examples\Sound and
Vibration\Audio Measurements\ to remove the delay due to the DAQ
device.
The anti-imaging and anti-aliasing filters have a lowpass filter effect on the
data. This effect results in a transient response at sharp transitions in the
data. These transitions are common at the start and stop of a generation,
at a change in frequency (swept sine), and when the amplitude changes
(amplitude sweep). The swept-sine analysis and audio measurements
examples in the Sound and Vibration Toolkit account for this transient
behavior in the device response to achieve the highest degree of accuracy.
PFI0
AI0
DUT
AI1
In Out
AO0
AO1
The DUT propagation delay is the delay of the entire system minus the
device delay. Remember to measure the device delay without the DUT
connected.
© National Instruments Corporation 3-7 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 3 Scaling and Calibration
The propagation delay for an analog DUT is a constant time delay rather
than a delay of samples. Use the following equation to convert the
measured delay in samples to the equivalent delay in seconds:
Note The swept sine VIs expect the DUT propagation delay measurement in seconds and
use the equation to convert the delay in seconds to samples.
Use the SVT Limit Testing VI to perform analysis on any type of measured
result produced by the Sound and Vibration Toolkit, including the
following measurements:
• Waveform
• Spectrum
• Peak
• Octave
• Swept sine
• Scalar
Refer to the LabVIEW Help for information about the individual SVT Limit
Testing VI instances.
© National Instruments Corporation 4-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 4 Limit Testing Analysis
DSA Device
DAQ Device
Data Source WAV File
DAT Recorder
Simulation
Calibration
Calibrate
Sensor
Scale Voltage to
Scaling Engineering Units
Measure
Propagation Delay
Limit Testing
Limit Testing
Waveform Waveform
Visualization Chart Graph
XY Graph
You can use the SVT Limit Testing VI to analyze almost any measured
result produced by the Sound and Vibration Toolkit. Refer to Table 4-1 for
examples of datatypes supported by the SVT Limit Testing VI and VIs that
generate supported datatypes.
© National Instruments Corporation 4-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 4 Limit Testing Analysis
both to the SVT Limit Testing VI to define a pass range that varies in
shape and level based on acceptable results at any given point in the
measurement. You also can create a discontinuous mask which allows you
to perform limit testing on only a part of the results while ignoring the rest.
You must enter at least one limit, or the SVT Limit Testing VI returns an
error. You can visually display the input signal, failures, upper limit, and
lower limit by creating an indicator from the output values terminal.
The upper limit and lower limit inputs to the SVT Limit Testing VI must be
compatible with the input signal. Table 4-2 lists the criteria that must be
met for each input signal type that is compatible with the SVT Limit
Testing VI.
Limit testing covers a broad range of data testing from range detection to
discontinuous mask testing of a swept-sine frequency response spectrum.
Figures 4-2, 4-4, 4-6, and 4-8 illustrate some, but not all, of the different
ways you can use the SVT Limit Testing VI in your application.
© National Instruments Corporation 4-5 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 4 Limit Testing Analysis
Figure 4-4 shows a pass/fail test on the measured THD. This test only
checks the upper limit of the measurement, therefore, only the upper limit
is wired to the VI. The upper limit should have the same units as the input
measurement. In this case both THD and the upper limit are expressed as
percentages. Figure 4-5 shows the THD test output results.
© National Instruments Corporation 4-7 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 4 Limit Testing Analysis
© National Instruments Corporation 4-9 LabVIEW Sound and Vibration Toolkit User Manual
Weighting Filters
5
This chapter discusses using weighting filters in sound and vibration
analysis, including describing the purpose of weighting filters, the types of
weighting filters, and applying weighting to time-domain data, FFT-based
spectra, and octave spectra.
© National Instruments Corporation 5-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 5 Weighting Filters
120
100
80
Intensity (dB)
60
40
20
0
10 100 1000 10,000
Frequency (Hz)
© National Instruments Corporation 5-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 5 Weighting Filters
© National Instruments Corporation 5-5 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 5 Weighting Filters
The frequency response of the CCITT and C-message weighting filters are
specified in the ITU-T O.41 standard and Bell System Technical Reference
41009, respectively. Figure 5-5 shows the relative attenuation defined for
the CCITT and C-message weighting filters.
DSA Device
DAQ Device
Data Source WAV File
DAT Recorder
Simulation
Calibration
Calibrate
Sensor
Scale Voltage to
Scaling Engineering Units
Measure
Progagation Delay
Weighting Weighting
Limit Testing
© National Instruments Corporation 5-7 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 5 Weighting Filters
Note The weighting filter parameter in the SVL Scale Voltage to EU VI channel info
control assigns the correct units to the waveform, but it does not cause the VI to perform
any filtering.
For A, B, or C-weighting, you can use the VIs on the Weighting (Arbitrary
Rate) palette to apply weighting for sample rates not listed in Table 5-1.
Supported Filters
Sample Rates A, B, C-weighting ITU-R 468-4/Dolby CCITT/C-message
4 kHz to 20 kHz Yes — Yes
4 kHz, 8 kHz, 10 kHz,
11.025 kHz, 12.8 kHz,
20 kHz to 1 MHz Yes Yes Yes
20 kHz, 22.05 kHz,
25.6 kHz, 40 kHz,
44.1 kHz, 48 kHz, 50 kHz,
51.2 kHz, 80 kHz, 96 kHz,
100 kHz, 102.4 kHz,
192 kHz, 200 kHz,
204.8 kHz, 500 kHz, 1 MHz
NI recommends using the fixed-rate weighting filter VIs if the VIs support
the desired sample rate. These VIs offer two advantages over the arbitrary
rate VIs: compliance with the appropriate standards over the entire
frequency range and slightly faster execution due to precomputed filter
coefficients.
Note The filter design algorithms used by the fixed and arbitrary rate weighting
approaches are different. Using a fixed-rate weighting filter with a supported frequency or
using the equivalent arbitrary rate filter at the same sample rate achieve different results.
Each implementation offers compliance with the appropriate standard over the frequency
range specified in the Standards Compliance section.
Standards Compliance
This section discusses the standards to which the various weighting filter
VIs comply.
ANSI Standards
When combined with any DSA device, the weighting filter used by the
SVT A, B, C Weighting Filter (Fixed Rates) VI or designed by the SVT
Weighting Filter VI complies with the following standards:
• ANSI S1.4-1983
• ANSI S1.42-2001
The SVT Weighting Filter VI accommodates any sample rate greater than
4 kHz and designs the filter coefficients to target the attenuation curves
defined by the ANSI standards. Given the selected sampling frequency,
compliance with a particular filter type, either Type 1 or Type 0, is ensured
up to a specific frequency. This frequency is the maximum frequency
within tolerances. Use the SVT Max Frequency Within Tolerances [ANSI]
VI located on the Weighting (Arbitrary Rate) palette to determine the
maximum frequency within tolerances.
© National Instruments Corporation 5-9 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 5 Weighting Filters
ISO/IEC Standard
Use the SVT A, B, or C Weighting Filter (Fixed Rates) VI or the SVT
Weighting Filter VI to apply an A-, B-, or C-weighting filter to
time-domain signals. When combined with any DSA device, the weighting
filter used by the SVT A, B, or C Weighting Filter (Fixed Rates) VI or
designed by the SVT Weighting Filter VI complies with the IEC
61672-1:2002 standard.
The SVT Weighting Filter VI accommodates any sample rate greater than
4 kHz and designs the filter coefficients to target the attenuation curves
defined by the IEC standards. Given the selected sampling frequency,
compliance with a particular filter type, either Class 2 or Class 1, is ensured
up to a specific frequency. This frequency is the maximum frequency
within tolerances. Use the SVT Max Frequency Within Tolerances [IEC]
VI to determine the maximum frequency within tolerances.
Time-domain data is simulated and scaled before being sent to the SVT A,
B, C Weighting Filter (Fixed Rates) VI. The weighted signal is sent to the
Sound Level Measurement VI.
Figure 5-8 shows the time-domain input and output waveforms when
a 250 Hz sine wave is sent to the SVT A, B, C Weighting Filter (Fixed
Rates) VI using the A-weighting filter.
© National Instruments Corporation 5-11 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 5 Weighting Filters
There is a phase difference between the input and output signals because a
filter applies the time-domain weighting. The transient behavior at the
beginning of the filtered waveform corresponds to the filter settling time.
0
Potential
Measurement
Error A-Weighted
Filter
–20 Response
Proportional
Bandwidth
Filter
Correction
–40
Note The same type of measurement error in Figure 5-10 can occur when applying
weighting to FFT-based spectra. However, the error is almost always negligible as long as
the frequency resolution of the spectrum is reasonable. For example, the error is negligible
with a frequency resolution of 10 Hz or finer.
© National Instruments Corporation 5-13 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 5 Weighting Filters
The Sound and Vibration Toolkit contains the following integration VIs:
• SVT Integration VI located on the Integration palette for time-domain
integration
• SVT Integration (frequency) VI located on the Frequency
Analysis»Extended Measurements palette for frequency-domain
integration
Refer to the LabVIEW Help for information about the integration VIs.
Introduction to Integration
The conversion between acceleration, velocity, and displacement is based
on one of the fundamental laws in Newtonian physics, represented by the
following equations:
. d
x = ----- ( x )
dt
2
.. d . d
x = ----- x = ------2- ( x )
dt dt
© National Instruments Corporation 6-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 6 Integration
a = A sin (ω t )
A A π
v = – ---- cos (ω t ) = ---- sin ω t – --- (6-1)
ω ω 2
A A
d = – ------ sin (ω t ) = -----2- sin (ω t – π ) (6-2)
2
ω ω
Note Initial condition is arbitrarily set to zero in Equations 6-1 and 6-2.
Implementing Integration
If you need to perform measurements on velocity or displacement data
when you have only acquired acceleration or velocity data, respectively,
integrate the measured signal to yield the desired data. You can perform
integration either in the time domain as a form of signal conditioning or
in the frequency domain as a stage of analysis. When performed in the
frequency domain, integration is one of the extended measurements for
frequency analysis.
Figure 6-2 illustrates how integration fits into the sound and vibration
measurement process.
© National Instruments Corporation 6-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 6 Integration
DSA Device
DAQ Device
Data Source WAV File
DAT Recorder
Simulation
Calibration
Calibrate
Sensor
Scale Voltage to
Scaling Engineering Units
Measure
Propagation Delay
Waveform
Integration
Conditioning
Integration
Extended
Measurements
Limit Testing
DC Component
Even though a DC component in the measured signal might be valid, the
presence of a DC component indicates that the DUT has a net acceleration
along the axis of the transducer. For a typical vibration measurement, the
DUT is mounted or suspended in the test setup. The net acceleration of the
DUT is zero. Therefore, any DC component in the measured acceleration
is an artifact and should be ignored.
Transducers
Most acceleration and velocity transducers are not designed to accurately
measure frequency components close to DC. Closeness to DC is relative
and depends on the specific transducer. A typical accelerometer can
accurately measure components down to about 10 Hz. A typical velocity
probe can accurately measure components down to 2–3 Hz. Inaccurately
measured low-frequency vibrations can dominate the response when the
signal is integrated because integration attenuates low-frequency
components less than high-frequency components.
© National Instruments Corporation 6-5 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 6 Integration
Time-Domain Integration
This section presents examples of and discussion about time-domain
integration.
© National Instruments Corporation 6-7 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 6 Integration
In this example, the highpass cut-off frequency used for the integration is
10 Hz. Additionally, the integration is explicitly reset in the first iteration
of the VI and performed continuously thereafter. In this example, this
additional wiring is optional because the SVT Integration VI automatically
resets the first time it is called and runs continuously thereafter.
If you use the block diagram in Figure 6-5 in a larger application that
requires starting and stopping the data acquisition process more than once,
NI suggests setting the reset filter control to TRUE for the first iteration of
the While Loop. Setting the reset filter control to TRUE causes the filter to
reset every time the data acquisition process starts. Set the reset filter
control to FALSE for subsequent iterations of the While Loop.
Figure 6-6 shows the results of the continuous acquisition and integration
of the same 38 Hz sinusoid used in the single-shot acquisition and
integration example.
© National Instruments Corporation 6-9 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 6 Integration
In Figure 6-7, you can see the characteristic 20 dB per decade roll-off of the
magnitude response of the single integration. In Figure 6-8, you can see
the characteristic 40 dB per decade roll-off of the magnitude response of
the double integration.
Upper and lower frequency limits exist for which you can obtain a specified
degree of accuracy in the magnitude response. For example, sampling at
a rate of 51.2 kHz, the magnitude response of the integrator is accurate
to within 1 dB from 1.17 Hz to 9.2 kHz for single integration and from
1.14 Hz to 6.6 kHz for double integration. The accuracy ranges change
with the sampling frequency and the highpass cut-off frequency.
The attenuation of the single integration filter at 9.2 kHz is –95 dB.
The attenuation of the double integration filter at 6.6 kHz is –185 dB.
Accuracy at high frequencies usually is not an issue.
Frequency-Domain Integration
You can use the following strategies to obtain the spectrum of an integrated
signal:
• Perform the integration in the time domain before computing the
spectrum.
• Compute the spectrum before performing the integration in the
frequency domain.
Figure 6-9. Integration in the Time Domain and in the Frequency Domain
© National Instruments Corporation 6-11 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 6 Integration
Figure 6-10 shows the results of integrating in the time and frequency
domains.
Figure 7-1 illustrates how vibration-level measurement fits into the sound
and vibration measurement process.
© National Instruments Corporation 7-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 7 Vibration-Level Measurements
DSA Device
DAQ Device
Data Source WAV File
DAT Recorder
Simulation
Calibration
Calibrate
Sensor
Scale Voltage to
Scaling Engineering Units
Measure
Propagation Delay
Waveform
Integration
Conditioning
Vibration Level
Analysis Measurements
Limit Testing
Waveform Waveform
Visualization Chart Graph
© National Instruments Corporation 7-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 7 Vibration-Level Measurements
Note Set the restart averaging control on the SVT RMS Level VI to TRUE. Otherwise,
the SVT RMS Level VI accumulates intermediate results to compute the RMS vibration
level over the entire data acquisition instead of just over the last block of data.
V pk
F c = ---------
-
V rms
where
Fc is the crest factor.
Vpk is the peak value of the signal.
Vrms is the RMS value of the signal.
© National Instruments Corporation 7-5 LabVIEW Sound and Vibration Toolkit User Manual
Sound-Level Measurements
8
This chapter discusses some of the analysis concepts associated with
performing sound-level measurements and how you can use the Sound
Level VIs located on the Sound Level palette to perform sound-level
measurements. You can combine different sound-level measurements and
use them simultaneously to provide flexibility with acoustic measurements.
Figure 8-1 illustrates how sound-level measurement fits into the sound and
vibration measurement process.
© National Instruments Corporation 8-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 8 Sound-Level Measurements
DSA Device
DAQ Device
Data Source WAV File
DAT Recorder
Simulation
Calibration
Calibrate
Sensor
Scale Voltage to
Scaling Engineering Units
Measure
Propagation Delay
Waveform
Weighting Filter
Conditioning
Sound Level
Analysis Measurements
Limit Testing
Waveform Waveform
Visualization Chart Graph
Linear Averaging
You compute the Leq by integrating the square of the signal over a
fixed-time interval and dividing by the time interval. When you select
linear averaging, the Sound Level VIs return a single value. The value
returned represents the continuous decibel level that would have produced
the same sound energy in the same time T as the actual noise history. To
obtain intermediate results, you must split a long time record into several
smaller records. Linear averaging is represented by the following equation.
T
1 P rms
2
∫
L eq = 10log10 --- --------------
T P2
0 0
- dt
where
P0 is the reference pressure of 20 µPa for acoustics.
© National Instruments Corporation 8-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 8 Sound-Level Measurements
Figure 8-3 displays the resulting Leq measurement and the instantaneous
sound pressure level.
Figure 8-3. Leq and Instantaneous Sound Pressure Level versus Time
For one hour of data with a sampling frequency of 51.2 kS/s, you need to
acquire more than 184 million samples, as calculated by the following
equation:
60 minutes 60
--------------------------
seconds 51,200
--------------------------
samples-
------------------------------------ 184,320,000 samples
= --------------------------------------------------
hour minute second hour
The VI in Figure 8-4 performs an Leq over one second and repeats the
operation 3,600 times using a For Loop. The last result returned by the
SVT Leq Sound Level VI is the Leq over the one-hour period.
In order for the SVT Leq Sound Level VI to accumulate the intermediate
results, set the restart averaging control to FALSE or leave the control
unwired. You can make the intermediate results available by using the
auto-indexing capability of the For Loop. Refer to the LabVIEW Help for
information about auto-indexing.
Instead of performing an Leq over one second and repeating the operation
3,600 times, you can perform the measurement over two seconds and repeat
it 1,800 times, or perform the measurement over four seconds and repeat it
900 times, and so forth.
© National Instruments Corporation 8-5 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 8 Sound-Level Measurements
Exponential Averaging
Exponential averaging is a continuous averaging process that weights
current and past data differently. The amount of weight given to past data
as compared to current data depends on the exponential time constant.
In exponential averaging, the averaging process continues indefinitely.
Refer to the LabVIEW Help for information about the SVT Exp Avg Sound
Level VI and the SVT Decimated Exp Avg Sound Level VI.
© National Instruments Corporation 8-7 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 8 Sound-Level Measurements
Peak Hold
In peak-hold averaging, the largest measured sound pressure level value of
all previous values is computed and returned until a new value exceeds the
current maximum. The new value becomes the new maximum value and
is the value returned until a new value exceeds it. Peak hold actually is
not a true form of averaging because successive measurements are not
mathematically averaged. However, as with other averaging processes,
peak-hold averaging combines the results of several measurements into one
final measurement. As with exponential averaging, the averaging process
continues indefinitely. The formula for peak averaging is defined by the
following equation.
y [ k ] = max ( y [ k – 1 ], x [ k ] )
where
x[k] is the new measurement.
y[k] is the new average.
y[k – 1] is the previous average.
Use the Octave Analysis VIs located on the Octave Analysis palette to
perform the following analyses:
• Full-octave
• Third-octave
• Fractional-octave, including 1/1, 1/3, 1/6, 1/12, and 1/24 octave
The Octave Analysis VIs can accommodate any sampling frequency and
any number of fractional-octave bands.
Refer to the LabVIEW Help for information about the individual Octave
Analysis VIs.
© National Instruments Corporation 9-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 9 Fractional-Octave Analysis
DSA Device
DAQ Device
Data Source WAV File
DAT Recorder
Simulation
Calibration
Calibrate
Sensor
Scale Voltage to
Scaling Engineering Units
Measure
Propagation Delay
Waveform
Weighting Filter
Conditioning
Octave
Analysis Analysis
Weighting
Limit Testing
Waterfall Colormap /
Visualization XY Graph
Display Intensity Graph
© National Instruments Corporation 9-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 9 Fractional-Octave Analysis
The power in each band is computed and displayed in a bar graph with a
log scale for the x-axis. Figure 9-3 shows this bar graph.
A 1 Octave
F
f/4 f/2 f 2f 4f
1 Octave
Bandedge Frequencies
The quality constant Q is defined as the ratio of the bandwidth over the
center frequency of the filter. Q remains constant across all octave bands
for octave filters. For example, an octave filter with a center frequency of
1,000 Hz leads to the following bandedge frequencies:
f 1 = 1,000
------------- = 707 Hz
2
f 2 = ( 1,000 ) ( 2 ) = 1,414 Hz
BW = f 2 – f 1 = 707 Hz
707
Q = ------------- = 0.707
1,000
where
f1 and f2 are bandedge frequencies.
BW is the bandwidth.
Q is the quality constant.
© National Instruments Corporation 9-5 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 9 Fractional-Octave Analysis
8,000
f 1 = ------------- = 5,657 Hz
2
f 2 = ( 8,000 ) ( 2 ) = 11,314 Hz
BW = f 2 – f 1 = 5,657 Hz
5,657
Q = ------------- = 0.707
8,000
where
f1 and f2 are bandedge frequencies.
BW is the bandwidth.
Q is the quality constant.
The results obtained from calculating the bandedge frequencies indicate the
following bandwidth characteristics:
• The bandwidth of the octave filters is narrow if the center frequency
is low.
• The bandwidth of the octave filters is wider when the center frequency
is higher.
Fractional-Octave Filters
Octave filter resolution is limited, being that there are only 11 octaves in
the 16 Hz–16 kHz range. To overcome the limited resolution of octave
filters, you can use other filters known as fractional-octave filters. Rather
than covering one octave with a single filter, N filters are applied per octave
in order to improve resolution. Of the fractional-octave filters, the
third-octave (1/3) filter is widely used for fractional-octave analysis.
Figure 9-5 shows the 1/3 octave response at frequencies of 500 Hz, 630 Hz,
and 800 Hz.
© National Instruments Corporation 9-7 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 9 Fractional-Octave Analysis
Averaging
The Octave Analysis VIs support the following averaging types:
• Linear averaging
• Exponential averaging
• Equal confidence averaging
• Peak-hold averaging
Linear Averaging
Linear averaging is computed by integrating the square of the filtered
signal over a fixed time interval and dividing by the time interval. Refer to
Chapter 8, Sound-Level Measurements, for more information about linear
averaging.
Exponential Averaging
Exponential averaging is a continuous averaging process that weights
current and past data differently. The amount of weight given to past data
as compared to current data depends on the exponential time constant.
Figure 9-6 illustrates the block diagram for a VI performing 1/3 octave
analysis in the 20Hz–20kHz range using fast exponential averaging.
Note In order to use exponential averaging, you must connect the exp avg settings
parameter to a control or constant. Refer to the LabVIEW Help for information about the
exp avg settings control.
Peak-Hold Averaging
In peak-hold averaging, the largest measured level value of all previous
values is returned for each band until a new value exceeds the current
maximum. The new value becomes the new maximum value and is the
value returned until a new value exceeds it. Refer to Chapter 8,
Sound-Level Measurements, for more information about peak-hold
averaging.
You can use the restart averaging control to restart the averaging process
without resetting the filter. Restarting the averaging without resetting the
filter avoids the transient phase associated with the settling of the filter.
Refer to the Filter Settling Time section for information about filter settling
time.
© National Instruments Corporation 9-9 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 9 Fractional-Octave Analysis
Note The ANSI and IEC Octave Analysis VIs use the base 2 system for midband
frequencies instead of the base 10 system. Using the base 2 system means that the ratio of
two midband frequencies is a fractional power of 2 instead of a fractional power of 10.
ANSI Standard
According to the ANSI S1.11-1986 standard, the midband, or center,
frequency of the bandpass filter is defined by the following equation.
ib
f i = 1000 ⋅ 2
where
f i is the center frequency of the ith band-pass filter expressed in hertz.
i is an integer when i = 0 , f 0 = 1 kHz, which is the reference
frequency for the audio range.
b is the bandwidth designator and equals 1 for octave, 1/3 for
1/3 octave, 1/6 for 1/6 octave, 1/12 for 1/12 octave, and 1/24
for 1/24 octave.
IEC Standard
According to the IEC 1260:1995 standard, the midband frequency,
or center, frequency of the bandpass filter is defined by the following
equations.
ib
f i = 1000 ⋅ 2 for 1/N octave filters when N is odd
(------------------
i + 1 )b-
2
f i = 1000 ⋅ 2 for 1/N octave filters when N is even
where
f i is the center frequency of the ith band-pass filter expressed in hertz.
i is an integer when i = 0 , f 0 = 1 kHz, which is the reference
frequency for the audio range.
b is the bandwidth designator and equals 1 for octave, 1/3 for
1/3 octave, 1/6 for 1/6 octave, 1/12 for 1/12 octave, and 1/24
for 1/24 octave.
Nominal Frequencies
The exact midband frequencies are used to design the filters for
fractional-octave analysis. However, all the Octave Analysis VIs return
the nominal midband frequencies, also called the preferred frequencies.
In the case of octave and 1/3 octave analyses, the nominal frequencies are
calculated in accordance with the ANSI and IEC standards. In the case of
1/6, 1/12, and 1/24 octave analyses, the nominal frequencies are calculated
in accordance with the Annex A (informative) of the IEC 1260:1995
standard.
In the case of the ANSI standard, the default order of the bandpass filter
is 3, which leads to a Type 1-D filter.
© National Instruments Corporation 9-11 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 9 Fractional-Octave Analysis
For octave filters, choose a sampling frequency at least three times the exact
center frequency of the highest frequency band. For fractional-octave
filters, choose a sampling frequency at least 2.5 times the exact center
frequency of the highest frequency band.
Displaying Results
The Octave Analysis VIs produce a cluster containing center frequencies
and band power. The Sound & Vibration controls palette contains two
graphs designed to display octave results, the Octave Graph and the
Multiplot Octave Graph. Refer to the Sound and Vibration Toolkit Help for
information about the Octave Graph and the Multiplot Octave Graph.
You can connect the cluster containing the center frequencies and band
power directly to an XY graph. However, you must modify the XY graph
in order to display octave results.
Weighting Filters
Use the SVT Weighting Filter (octave) VI to apply A-, B-, or C-weighting
filters to a previously computed octave spectrum or spectra. Refer to the
LabVIEW Help for more information about the SVT Weighting Filter
(octave) VI. Refer to Chapter 5, Weighting Filters, for more information
about weighting filters.
Refer to the Sound and Vibration Toolkit Help for information about
individual Frequency Analysis VIs.
Note For simplicity, the remainder of this document uses the term FFT to denote both the
FFT and the DFT.
© National Instruments Corporation 10-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
DSA Device
DAQ Device
Data Source WAV File
DAT Recorder
Simulation
Calibration
Calibrate
Sensor
Scale Voltage to
Scaling Engineering Units
Measure
Propagation Delay
Frequency
Analysis Analysis
Extended
Measurements
Limit Testing
FFT Fundamentals
The FFT resolves the time waveform into its sinusoidal components. The
FFT takes a block of time-domain data and returns the frequency spectrum
of that data. The FFT is a digital implementation of the Fourier transform.
Thus, the FFT does not yield a continuous spectrum. Instead, the FFT
returns a discrete spectrum where the frequency content of the waveform
is resolved into a finite number of frequency lines, or bins.
Number of Samples
The computed spectrum is completely determined by the sampled time
waveform input to the FFT. If an arbitrary signal is sampled at a rate equal
to fs over an acquisition time T, N samples are acquired. Compute T with
the following equation:
N
T = ----
fs
where
T is the acquisition time.
N is the number of samples acquired.
fs is the sampling frequency.
N = T ⋅ fs
where
N is the number of samples acquired.
T is the acquisition time.
fs is the sampling frequency.
Frequency Resolution
Because of the properties of the FFT, the spectrum computed from the
sampled signal has a frequency resolution df. Calculate the frequency
resolution with the following equation:
f
df = --1- = ----s
T N
where
df is the frequency resolution.
T is the acquisition time.
fs is the sampling frequency.
N is the number of samples.
© National Instruments Corporation 10-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
f
f max = f Nyquist = ---s
2
where
fmax is the maximum resolvable frequency.
fNyquist is the Nyquist frequency.
fs is the sampling frequency.
Directly specify the number of lines in the spectrum for the Zoom FFT VIs
in the Sound and Vibration Toolkit. Specify the number of data samples to
control the number of spectral lines for the Baseband FFT and
the Baseband Subset VIs.
Frequency Domain
fNyquist fmax # lines df
fs f 1 f
Time N ---s fs ⋅ Eb Eb ⋅ N --- = ----s
2 T N
Domain T
where
fs is the sampling frequency.
Eb is the effective bandwidth.
N is the number of samples acquired.
T is the acquisition time.
Use the information in Table 10-2 if you prefer to specify the spectrum
parameters and determine the required data-acquisition parameters from
these specifications.
Time Domain
fs N T
fmax f max # lines 1- = #---------------
lines
Frequency # lines --------- --------------- ----
Eb Eb df f max
Domain df
where
fmax is the maximum resolvable frequency.
Eb is the effective bandwidth.
# lines is the number of lines in the spectrum.
df is the frequency resolution.
© National Instruments Corporation 10-5 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
# lines = E b ⋅ N
1
# lines = ---------- ⋅ 1,024 = 400 lines
2.56
Note Table 10-1 shows that the sampling frequency and the block size acquired during
each cycle of a continuous acquisition completely determine the frequency-domain
parameters in baseband FFT analysis. However, many stand-alone instruments are
operated by specifying the frequency range of interest and the number of lines in the FFT.
Table 10-2 shows how a stand-alone instrument uses the frequency range of interest and
the number of lines in the FFT to determine an appropriate sampling frequency and block
size.
© National Instruments Corporation 10-7 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
Implement the decreased fs strategy with zoom FFT analysis. Use baseband
FFT and FFT-subset analyses to implement the increased N strategy.
Baseband FFT analysis and FFT-subset analysis both achieve the same
frequency resolution. However, FFT subset analysis only computes a
narrow subset of the spectrum.
Refer to the Zoom FFT Analysis section and the Subset Analysis section
for examples that demonstrate the importance of frequency resolution in
frequency analysis. The examples illustrate how to achieve a finer
frequency resolution with the frequency analysis tools in the Sound and
Vibration Toolkit.
Zoom FFT analysis achieves a finer frequency resolution than the baseband
FFT. The Zoom FFT VI acquires multiple blocks of data and downsamples
to simulate a lower sampling frequency. The block size is decoupled from
the achievable frequency resolution because the Zoom FFT VI accumulates
the decimated data until you acquire the required number of points.
Because the transform operates on a decimated set of data, you only need
to compute a relatively small spectrum. The data is accumulated, so do not
think of the acquisition time as the time required to acquire one block of
samples. Instead, the acquisition time is the time required to accumulate the
required set of decimated samples.
The Zoom FFT VIs complete the following steps to process the
sampled data:
1. Modulate the acquired data to center the analysis band at 0 frequency.
2. Filter the modulated data in the time domain to isolate the analysis
band and prevent aliasing when the data is resampled at a lower
sampling frequency.
3. Decimate the filtered data to reduce the effective sampling frequency.
4. Accumulate the decimated data until sufficient samples are available to
compute the spectrum.
5. Use the Discrete Zak Transform (DZT) to efficiently compute the
desired spectral lines.
6. Demodulate, or shift, the computed spectrum.
Note The exact frequency resolution is returned as df in the spectrum computed by the
Zoom FFT VIs.
© National Instruments Corporation 10-9 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
Zoom Measurement
The following example demonstrates a zoom measurement of the power
spectrum.
Acquire the signal at 51.2 kHz. The VI reads the data in blocks of
2,048 samples. Compute the frequency resolution of this measurement
using baseband analysis with the following equation:
51,200 Hz- = 25 Hz
------------------------
2,048
Use the SVT Zoom Power Spectrum VI located on the Zoom FFT palette
to analyze a narrow band with a much finer frequency resolution.
Figure 10-5 shows the result of limiting the measurement to the frequency
band between 1 kHz and 2 kHz and computing 400 lines. Derive the
frequency resolution of the computed spectrum with the following
equation:
Zoom Settings
Figure 10-6 shows the zoom settings control used to acquire the zoom
measurement results displayed in Figure 10-5. Use this control to specify
the frequency range, window, number of lines, and percent overlap used in
the zoom analysis. Refer to the LabVIEW Help for information about the
zoom settings control.
Subset Analysis
The Baseband Subset VIs located on the Baseband Subset palette allow you
to compute a subset of the baseband FFT measurement. Subset analysis
uses the DZT to compute a subset of the baseband FFT. The frequency
resolution for spectral measurements computed with the Baseband Subset
VIs equals the frequency resolution for measurements made with the
© National Instruments Corporation 10-11 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
The Baseband Subset VIs algorithm computes only the desired spectral
lines. The only programming difference between the Baseband Subset VIs
and the Baseband FFT VIs is the additional parameter frequency range.
The frequency range parameter specifies which spectral lines the
Baseband Subset VI computes. The computed spectral lines are always
inclusive of the start frequency and the stop frequency.
Note Setting the start frequency to 0 Hz and the stop frequency to fmax yields the same
spectrum as the corresponding Baseband FFT VI. If you set the stop frequency to –1, the
baseband subset VIs return the Nyquist frequency as the highest frequency in the computed
spectrum.
The following consideration can help you decide when to use the Baseband
Subset VIs instead of the Baseband FFT VIs:
• The required block size yields an acceptable frequency resolution.
• The analysis of a narrow subset of the baseband span requires better
processing performance than the Baseband FFT VI can provide.
Available Measurements
The following Frequency Analysis palettes offer the same basic
measurements but designed for specific measurement needs:
• Baseband FFT
• Baseband Subset
• Zoom FFT
For example, each of the three palettes contains a VI for measuring the
power spectrum. The SVT Power Spectrum VI located on the Baseband
FFT palette computes the power spectrum of the input signal. The SVT
Power Spectrum Subset VI located on the Baseband Subset palette
computes a subset of the power spectrum of the input signal. The SVT
Zoom Power Spectrum VI located on the Zoom FFT palette computes a
zoom power spectrum of the input signal. Refer to the LabVIEW Help for
a complete listing of the Frequency Analysis VIs and information about
each VI.
The Baseband FFT, Baseband Subset, and Zoom FFT VIs all share the
same basic relationships between the input signal and the computed
spectrum. For the baseband and subset analyses, you can obtain a tighter
frequency resolution only by increasing the block size. Increasing the block
size results in an FFT computed with more lines. Zoom analysis internally
reduces the sampling frequency by decimating the data. In baseband FFT,
baseband subset, and zoom FFT analyses, the frequency resolution is the
inverse of the measurement duration.
Single-Channel Measurements
You can perform the following single-channel measurements with the
Frequency Analysis VIs:
• Power spectrum computes the power present within each spectral bin.
All phase information is lost in the computation. This measurement is
a useful tool for examining the various frequency components in a
signal.
• Power spectral density computes the power present within each bin
normalized by the bin width. All phase information is lost in the
computation. This measurement is a useful tool for examining the
noise floor in a signal or the power in a specific frequency range.
Normalizing the power spectrum by the bin width decouples the result
of this measurement from the block size N.
• FFT spectrum computes either the magnitude and phase or the real and
imaginary parts of the spectrum of the input signal. Phase information
is retained depending on the selected averaging mode. This
measurement is most often used by more advanced measurements that
require magnitude and phase information.
© National Instruments Corporation 10-13 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
Note In this example, the output of the read VI is an array of waveforms. You can wire the
output array from the read VI directly to the SVL Scale Voltage to EU VI. If the read VI
output array contains only a single waveform, use the instance of the polymorphic read VI
that returns only a single channel before wiring the waveform to the SVL Scale Voltage to
EU VI.
Dual-Channel Measurements
Dual-channel measurements differ from single-channel measurements
because the output spectrum of the dual-channel measurements is
dependent on the relationship between two input channels. Typically,
the input signals are a stimulus and a response. Some form of broadband
excitation usually is required to obtain accurate results. Broadband signals
include noise, chirps, multi-tone signals, impulses, as well as others.
Refer to the LabVIEW Help for more information about the view parameter.
© National Instruments Corporation 10-15 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
PFI0
AI0
AI1
DUT
AO0
In Out
AO1
Figure 10-9. Connection Scheme for FRF Measurement with an NI 4461 Device
The DSA device converts the stimulus signal from digital to analog. Analog
output channel 0, AO0, sends the stimulus signal to the DUT. Analog input
channel 0, AI0, receives the stimulus signal. Analog input channel 1, AI1,
receives the response of the DUT.
Figure 10-10 shows the block diagram of the VI used to perform the FRF
measurement.
The While Loop in Figure 10-10 controls both the generation and the
acquisition of the signal. For each iteration of the While Loop, the Uniform
White Noise Waveform VI generates a white-noise signal. The white-noise
signal is sent to the DUT on analog output channel 0, AO0 in Figure 10-9.
Analog input channel 0, AI0 in Figure 10-9, acquires the same number of
samples as the buffer that is generated. Simultaneously, analog input
channel 1, AI1 Figure 10-9, acquires the response signal from the DUT.
Figure 10-11 shows the measured time-domain stimulus and response
signals.
© National Instruments Corporation 10-17 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
Note The SVT Frequency Response [Mag-Phase] VI uses the stimulus and response
signals from the DUT to compute the FRF. In the following examples, only the magnitude
of the frequency response function is displayed.
© National Instruments Corporation 10-19 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
Windowing
Periodicity is one of the basic assumptions made in FFT-based frequency
analysis. The FFT algorithm implicitly assumes that every block of
acquired data indefinitely repeats in both positive and negative time.
Windowing is one method of ensuring periodicity.
© National Instruments Corporation 10-21 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
The windows supported by the Sound and Vibration Toolkit and their
equivalent noise bandwidths (ENBW) are listed in Table 10-3. ENBW is
a property of the window applied to the signal.
Window ENBW
None 1
Hanning 1.50
Hamming 1.36
Blackman-Harris 1.71
Exact Blackman 1.69
Blackman 1.73
Flat Top 3.77
4 Term Blackman-Harris 2.00
7 Term Blackman-Harris 2.63
Low Sidelobe 2.22
Force-Exponential N/A
Averaging Parameters
Each of the Frequency Analysis VIs supports averaging. The averaging
parameters control in Figure 10-16 defines how the averaged spectrum is
computed. averaging mode, weighting mode, number of averages, and
linear mode each control a particular feature of the averaging process.
Refer to the LabVIEW Help for more information about the averaging
parameters control and for information about number of averages and
linear mode.
Averaging Mode
You can choose from the following averaging modes when performing
frequency analysis with the Frequency Analysis VIs:
• No averaging
• RMS averaging
• Vector averaging
• Peak hold
© National Instruments Corporation 10-23 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
Note Not all of the Frequency Analysis VIs support all of the averaging modes listed
above. Refer to the LabVIEW Help for information about the averaging modes supported
by a specific VI.
No Averaging
No averaging is the default setting and does not apply any averaging to
the measurement. You can use the No averaging setting for quick
computations or when the signal-to-noise ratio is high.
RMS Averaging
RMS averaging reduces signal fluctuations but not the noise floor. The
noise floor is not reduced because RMS averaging averages the power
of the signal. Because RMS averaging averages the power of the signal,
averaged RMS quantities of single-channel measurements have zero phase.
RMS averaging for dual-channel measurements preserves important phase
information.
FFT spectrum 〈 X∗ ⋅ X 〉
power spectrum 〈 X∗ ⋅ X 〉
cross spectrum 〈 X∗ ⋅ Y 〉
frequency response 〈 X∗ ⋅ Y 〉-
H1 = --------------------
〈 X∗ ⋅ X 〉
〈 Y∗ ⋅ Y 〉
H2 = ---------------------
〈 Y∗ ⋅ X 〉
( H1 + H2 )
H3 = --------------------------
2
where
X is the complex FFT of the stimulus signal x.
Y is the complex FFT of the response signal y.
X∗ is the complex conjugate of X.
Y∗ is the complex conjugate of Y.
〈 X 〉 is the average of X, real and imaginary parts being averaged
separately.
Vector Averaging
Vector averaging, also called coherent averaging or time synchronous
averaging, reduces the amount of noise in synchronous signals. Vector
averaging computes the average of complex quantities directly. The real
and imaginary parts are averaged separately, which preserves phase
information. However, for single-channel measurements, using vector
averaging without a triggered acquisition can cause strong spectral
components to be eliminated in the averaged spectrum.
FFT spectrum 〈X 〉
power spectrum 〈 X∗ 〉 ⋅ 〈 X 〉
cross spectrum 〈 X∗ 〉 ⋅ 〈 Y 〉
〈 X 〉-
frequency response ---------- ( H1 = H2 = H3 )
〈Y 〉
where
X is the complex FFT of the stimulus signal x.
Y is the complex FFT of the response signal y.
〈 X 〉 is the average of X, real and imaginary parts being averaged
separately.
〈 Y 〉 is the average of Y, real and imaginary parts being averaged
separately.
© National Instruments Corporation 10-25 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
RMS Averaging does not reduce the noise floor. However, RMS averaging
does smooth the noise out enough to unmask the tone at 15 kHz.
Peak Hold
Peak hold is performed at each individual frequency line and retains the
RMS peak levels of the averaged quantities from one FFT record to the next
record. Peak-hold averaging is most useful when configuring a
measurement system or when applying limit or upper limit testing to a
frequency spectrum.
where
X is the complex FFT of the stimulus signal x.
X∗ is the complex conjugate of X.
Weighting Mode
Linear weighting weights each individual spectrum by the same amount in
the averaged spectrum. Linear weighting is most often used for analysis.
N–1 1
Y i = ------------- Y i – 1 + ---- X i ,
N N
where
X i is the result of the analysis performed on the ith block.
Y i is the result of the averaging process from X 1 to X i .
N = i for linear weighting.
N is a constant for exponential weighting with N = 1 for i = 1.
© National Instruments Corporation 10-27 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
2
Coherence
2 〈 X∗ ⋅ Y 〉
γ = ------------------------------------------
-
〈 X ⋅ X 〉 〈 Y∗ ⋅ Y 〉
∗
2
Coherent output COP = γ 〈 Y∗ ⋅ Y 〉
power
where
X is the complex FFT of the stimulus signal x.
Y is the complex FFT of the response signal y.
X∗ is the complex conjugate of X.
Y∗ is the complex conjugate of Y.
〈 X 〉 is the average of X, real and imaginary parts being averaged
separately.
Extended Measurements
You can use the VIs located on the Frequency Analysis»Extended
Measurements palette to perform the following measurements:
• Limit testing
• Frequency-domain weighting
• Spectrum peak search
• Power in band
• Units conversion
• Frequency-domain integration
Refer to the LabVIEW Help for information about the individual Extended
Measurements VIs.
Power in Band
The SVT Power in Band VI measures the total power within the specified
frequency range. Table 10-4 shows the equations for computing power in
band based on the type of input spectrum.
∑ PS
start frequency
--------------------------------
ENBW
∑
start frequency
( MS )
2
-----------------------------------------
ENBW
where
ENBW is the equivalent noise bandwidth.
df is the frequency resolution of the measurement.
© National Instruments Corporation 10-29 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 10 Frequency Analysis
The SVT Spectrum Peak Search VI uses an algorithm to estimate the actual
frequency and amplitude of each tone that meets the search criteria. The
algorithm also estimates the frequency and amplitude of tones that are not
exactly periodic within the measurement window. The ability to identify
nonperiodic tones that are not exactly on-bin allows you to use a smaller
block size, therefore a smaller FFT, and still accurately identify tones
present within the spectrum.
Unit Conversion
The SVT Unit Conversion VI operates on complex, magnitude, and power
spectra. You can use the VI to convert the complex spectra produced by
distortion and single-tone measurements to a magnitude or power spectrum
more suitable for display on a waveform graph. You also can use the
conversion utility to switch between magnitude and power spectra, switch
between linear and logarithmic scaling, change the decibel reference for the
spectrum, and change the peak units.
Note You cannot use the SVT Unit Conversion VI to scale the data to engineering units
or change the spectrum physical units. Use MAX, the DAQ Assistant, or the SVL Scale
Voltage to EU VI to scale your data to the proper engineering units.
When using the SVT Unit Conversion VI, input the desired scaling
parameters as they should apply to the converted spectrum parameter.
You do not need to perform units conversions on the decibel reference. The
decibel reference always has the same units as the converted spectrum. For
example, if you have a magnitude spectrum in Vrms, you can convert this to
a power spectrum in decibels referenced to 10 Vpk by specifying the
desired scaling parameters shown in Table 10-5.
desired scaling
linear/decibels dB
decibel reference (0: no change) 10
power/magnitude power
view peak
DSA Device
DAQ Device
Data Source WAV File
DAT Recorder
Simulation
Calibration
Calibrate
Sensor
Scale Voltage to
Scaling Engineering Units
Measure
Propagation Delay
Waveform
Integration Weighting Filter
Conditioning
Transient
Analysis Analysis
Waterfall Colormap /
Visualization Display Intensity Graph
XY Graph
© National Instruments Corporation 11-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 11 Transient Analysis
You can use the STFT VIs to extract frequency information as a function of
time directly from the signal of interest. Additionally, in the case of a
rotating machine where a tachometer signal is simultaneously acquired
with the signal of interest, the STFT VIs can extract frequency information
as a function of the rotational speed.
You can use the SVT Shock Response Spectrum VI to evaluate the severity
of a shock signal. The results generated by the SRS are typically displayed
on an XY graph.
Note Other LabVIEW toolkits are available that provide additional transient analysis
capabilities. The Order Analysis Toolkit is designed for rotating machinery analysis and
monitoring. The Signal Processing Toolkit has tools, such as wavelets and joint
time-frequency analysis (JTFA), for the analysis of fast transients.
Figure 11-2 shows the signal corresponding to the first 200 ms of the
waveform.
Figure 11-3 shows the result of applying a baseband FFT on the entire
waveform.
© National Instruments Corporation 11-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 11 Transient Analysis
Instead of computing a single FFT on the whole data set, you can divide the
data set into smaller blocks and compute FFTs on these smaller data blocks.
For example, divide the signal into 100 ms blocks and perform an FFT on
each of the blocks with the SVT STFT versus Time VI.
Leave from [s] and to [s] each equal to –1.00 to ensure that all of the signal
is used in the STFT computation. In this particular example, the –1.00
setting in both from [s] and to [s] is equivalent to setting from [s] to 0 and
to [s] to 10.
Because the time increment is 100 ms and a 1,024 sample FFT only
requires a 20 ms block, only one block out of five is used for computation.
Figure 11-5 shows the result obtained with a 1,024 sample FFT.
If you select an FFT Block size of 4,096 samples, or 1,600 alias-free lines,
the resolution improves, as illustrated in Figure 11-6. However, the
increased resolution comes with the expense of extra processing.
© National Instruments Corporation 11-5 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 11 Transient Analysis
Overlapping
Overlapping is a method that uses a percentage of the previous data block
to compute the FFT of the current data block. When combined with
windowing, overlapping maximizes the use of the entire data set. If no
overlapping is used, the part of the signal close to the window edges
becomes greatly attenuated. The attenuation of the signal near the window
edges could result in the loss of information in the region near the
window edges.
Note Set the desired overlapping rate by specifying % in the time increment units (%)
in the time increment control. Refer to Figure 11-4 for the location of this control. No
overlapping, or 0%, corresponds to a time increment of 100%. An overlapping of 75%
corresponds to a time increment of 25%. An overlapping of 50% corresponds to a time
increment of 50%, and so forth. The advantage of using the time increment control is that
you can specify values greater than 100%. For example, a time increment of 200%
corresponds to computing an FFT on every other block of data.
Figure 11-7 and Figure 11-8 illustrate the overlapping process. Figure 11-7
shows a 50% overlap.
Figure 11-8 shows the resulting subdivisions when you use a 50% overlap
and a Hamming window.
© National Instruments Corporation 11-7 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 11 Transient Analysis
Note Use the X Scale and Y Scale offset and multiplier properties to properly scale the
axes of the intensity graph. In this example, the X Scale range is 0 to 10 s. The Y Scale
range is 0 to 25,600 Hz. 25,600 Hz is the Nyquist frequency. You can adjust the Z Scale so
that only the relevant part of the signal is displayed. In other words, you can hide noise in
the displayed signal by increasing the minimum limit of the z-axis. Refer to the LabVIEW
Help for information about the offset and multiplier properties for graph controls.
Note For simplicity, the remainder of this chapter uses the term tachometer to denote both
a tachometer and an encoder.
© National Instruments Corporation 11-9 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 11 Transient Analysis
You can use the SVT Convert to RPM (analog) VI to measure the rotational
speed in RPM as a function of time. Figure 11-13 shows the result obtained
with the SVT Convert to RPM (analog) VI and a simulated tachometer
signal.
The signal from the tachometer is also acquired. The measured tachometer
signal is converted to RPM with the SVT Convert to RPM (analog) VI.
Figure 11-15 shows the rotational speed as a function of time as computed
by the SVT Convert to RPM (analog) VI.
© National Instruments Corporation 11-11 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 11 Transient Analysis
Using the SVT STFT versus RPM (analog) VI allows you to measure the
frequency content of the signal as a function of the rotational speed of the
engine. Figure 11-16 displays the results obtained with the SVT STFT
versus RPM (analog) VI on an intensity graph.
Figure 11-16. Intensity Graph of Sound Pressure Level for an Engine Run-Up
The resonance frequency, fN, and the critical damping factor, ζ, characterize
a SDOF system, where:
1 k
f N = ------ ---
-
2π m
c
ζ = --------------
2 km
For light damping where ζ is less than or equal to 0.05, the peak value of
the frequency response occurs in the immediate vicinity of fN and is given
by the following equation, where Q is the resonant gain:
1
Q = ------
2ζ
© National Instruments Corporation 11-13 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 11 Transient Analysis
Use the signals shown in Figure 11-17 to construct the SRS. For example,
the maximax, the absolute maximum response of the calculated shock
response signal over the entire signal duration, uses the absolute maximum
system response as a function of the system natural frequency.
Figure 11-18 illustrates the maximax SRS for the same half-sine pulse.
Note Each computed SRS is specific to the pulse used to perform the measurement.
You can use other types of shock spectra depending on the application.
These spectra include the initial shock response from the system response
over the pulse duration or from the residual shock spectrum from the
system response after the pulse. You can use the positive maximum, the
negative maximum, or the absolute maximum response signal value.
The Sound and Vibration Toolkit uses the Smallwood algorithm to compute
the SRS. The SVT Shock Response Spectrum VI also offers the ability to
preprocess the time-domain signal to improve SRS results. You can remove
any DC component or apply a lowpass filter with a selectable cut-off
frequency.
© National Instruments Corporation 11-15 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 11 Transient Analysis
The SVT Shock Response Spectrum VI can compute the SRS from the
absolute acceleration response or from the relative displacement response.
Use the model control on the SVT Shock Response Spectrum VI to select
the appropriate response.
Figure 11-19 shows how to use the SVT Shock Response Spectrum VI.
The example in Figure 11-19 acquires 1,000 samples of data from an
accelerometer during a shock. The shock signal triggers the acquisition.
The program stores 100 samples before the trigger to properly capture the
entire shock signal.
Figure 11-20 displays the acquired time-domain signal and the computed
SRS.
© National Instruments Corporation 11-17 LabVIEW Sound and Vibration Toolkit User Manual
Waterfall Display
12
This chapter discusses using the Waterfall Display VIs located on the
Waterfall Display palette.
Refer to the LabVIEW Help for information about the individual Waterfall
Display VIs.
Use the Waterfall Display VIs to display FFT spectra from frequency
analysis or transient analysis and octave spectra from octave analysis in
waterfall graphs. Refer to Front Panel Displays in the LabVIEW Help for
more information about displaying octave results.
Specific Waterfall Display VIs display the results of frequency analysis and
octave analysis in a waterfall graph. The waterfall display opens in an
external window called the Waterfall window.
© National Instruments Corporation 12-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 12 Waterfall Display
Twenty data blocks of 1,024 samples are acquired. The power spectrum is
computed on each block. The auto-indexing capability of the For Loop is
used to build an array of 20 spectra. The array or spectra is sent to the
waterfall display. Refer to the LabVIEW Help for information about
auto-indexing.
© National Instruments Corporation 12-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 12 Waterfall Display
The data are scaled and sent to the SVT STFT versus Time VI. The SVT
STFT versus Time VI returns a 2D array. You can use the results in the 2D
array in an intensity graph or connect the 2D array directly to the SVL Send
Data to Waterfall VI. Figure 12-3 shows the 2D array connected directly to
the Waterfall VI. The While Loop keeps the waterfall display open until the
Stop control is set to TRUE.
Note Connect f0 and delta f and y0 and delta y on the SVL Send Data to Waterfall VI to
ensure the graph shows the proper scales.
Figure 12-4 shows the result obtained with the STFT VI illustrated in
Figure 12-3.
Figure 12-6 shows the waterfall display created by the VI in Figure 12-5.
© National Instruments Corporation 12-5 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 12 Waterfall Display
3 4 8
1 9
2 5 6 7
You can use the Change vertical perspective and Change horizontal
perspective sliders to change the perspective of the display. You can store
and recall any view at a later time and turn transparency on or off. You also
can autoscale each axis independently or all three axes simultaneously.
© National Instruments Corporation 13-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 13 Swept-Sine Measurements
Initialize
Calibrate
Sensor
Configure
Measure
Propagation Delay
Start
Generate
and Acquire
Process Data
Graph
Measurements
No
Measurements
Done?
Yes
Close
© National Instruments Corporation 13-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 13 Swept-Sine Measurements
Figure 13-4 shows the simulated frequency response function for a four
DOF system. The peak at 17.6 Hz has a magnitude roughly 1,000 times
larger than the peak at 5.8 Hz. To use an FFT-based technique, use
broadband excitation to excite the entire frequency range of interest to
measure the frequency response. This situation forces you to set the input
range so that the overall response does not overload the DUT or the
acquisition device. Therefore, when you measure the response at 5.8 Hz,
you lose 60 dB of measurement dynamic range. The swept-sine technique
allows you to tailor the excitation amplitude to the specific test frequency,
preserving the full measurement dynamic range.
© National Instruments Corporation 13-5 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 13 Swept-Sine Measurements
Use the Configure Swept Sine VIs in the Configure Swept Sine palette to
configure the scaling, test frequencies, averaging, delays, and other
measurement settings. These configuration VIs allow control over basic
and advanced measurement parameters. The order in which you place the
configuration VIs is important, as it allows you to customize a swept-sine
measurement. For example, you can easily generate 100 logarithmically
spaced test frequencies in the audio range, then apply inverse A-weighted
scaling to the excitation level by adding code similar to that in Figure 13-5
into your swept-sine application.
You can use the swept-sine configuration VIs to customize your swept-sine
application. For example, to speed up a swept-sine measurement, reduce
the settling or integration time specified by the SVT Set Swept Sine
Averaging VI. You also can configure the device integrated electronic
piezoelectric (IEPE) excitation with the SVT Set Swept Sine Coupling and
IEPE Excitation (DAQmx) VI. You also can reduce the block duration
input to SVT Set Swept Sine Block Duration VI.
Note The minimum block duration is limited by the capabilities of the computer
processing the measurement. A very small block duration can result in a loss of continuous
processing, causing the swept-sine measurement to stop and return an error.
Use the SVT Start Swept Sine VI to begin the generation and acquisition.
Initially, the VI fills the device output buffer with zeros before writing the
first test frequency excitation.
The SVT Swept Sine Engine VI continually acquires data and processes it
to remove samples acquired during delays, transitions, and settling periods.
The SVT Swept Sine Engine VI performs measurement analysis on
samples acquired during integration periods. The SVT Swept Sine Engine
VI updates the excitation to excite the DUT at the next test frequency after
it integrates sufficient data at the current test frequency.
Note The transition to the next excitation tone, both frequency and amplitude, always
occurs at a zero crossing to minimize transients introduced to the DUT.
Use the Read Swept Sine Measurements VIs in the Read Swept Sine
Measurements palette to read the raw measurements, scale the
measurements, and perform additional conversions to display and report
the data in the desired format.
Use the SVT Close Swept Sine VI to stop the generation and acquisition,
and clear the swept-sine task.
© National Instruments Corporation 13-7 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 13 Swept-Sine Measurements
PFI0
AI0
AI1
DUT
AO0
In Out
AO1
The NI PXI-4461 converts the desired stimulus signal from digital data to
an analog signal and outputs that signal on AO0. The excitation signal is
connected to both the stimulus input channel AI0 and the input terminal of
the DUT. The response signal is connected from the output terminal of the
DUT to the response input channel AI1.
© National Instruments Corporation 13-9 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 13 Swept-Sine Measurements
Table 13-2 documents the actions performed by the VIs in Figure 13-7.
Some steps are required and must be done for the VI to function correctly.
The optional steps allow you to customize your measurement.
The While Loop in Figure 13-7 controls the synchronized generation and
acquisition. Display controls and measurement indicators are updated
inside the While Loop. This loop allows for the monitoring of intermediate
results.
Many of the steps in Table 13-2 are configuration steps. Through the Sound
and Vibration Toolkit swept-sine configuration VIs, you can specify
numerous configuration parameters to achieve fine control of the
swept-sine measurement parameters. For many applications two or three
configuration VIs are sufficient.
© National Instruments Corporation 13-11 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 13 Swept-Sine Measurements
test frequency. From the time-domain data, you can see that the notch filter
has attenuated the signal and introduced a phase shift.
Figure 13-9 shows the magnitude and phase responses of the notch filter at
all the test frequencies in the magnitude and phase spectra in the Bode plot.
You expect to see a peak in the THD at the notch frequency. The peak
occurs because the fundamental frequency is attenuated at the notch
frequency. However, the graph indicates that this measurement has failed to
accurately identify the power in the harmonic distortion components. For
the example in Figure 13-10, the number of integration cycles is two.
More integration cycles must be specified to perform accurate harmonic
distortion measurements. If you change the number of integration cycles to
10 and rerun the example, you obtain the THD versus frequency results
displayed in Figure 13-11.
© National Instruments Corporation 13-13 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 13 Swept-Sine Measurements
Now, with a sufficient number of integration cycles specified, you can see
the characteristic peak in the THD at the center frequency of the notch filter.
Variable Definitions
In the equations in this chapter, the variables are defined as:
∑P Harmonics = ( H 2 ) 2 + ( H 3 ) 2 + ... + ( H K ) 2
∑P ∑P
2
Harmonics + Noise = ( H 2 ) 2 + ( H 3 ) 2 + ... + ( H K ) 2 + ( N )
∑P Total = ( F ) 2 + ( H 2 ) 2 + ... + ( H K ) 2 + ( N ) 2
© National Instruments Corporation 14-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 14 Distortion Measurements
Figure 14-1 shows where distortion analysis fits into the sound and
vibration measurement process.
DSA Device
DAQ Device
Data Source WAV File
DAT Recorder
Simulation
Calibration
Calibrate
Sensor
Scale Voltage to
Scaling Engineering Units
Measure
Propagation Delay
Distortion
Analysis Measurements
Extended
Measurements
Limit Testing
Waveform Waveform
Visualization Chart
XY Graph
Graph
∑
P Total
SINAD = -----------------------------------------------------------
∑
P Harmonics + ∑ P Noise
You can input a frequency range to refine the measurement. The measured
power in the total signal is distributed across the entire measurement
bandwidth as a result of the noise in the measured signal. You can limit the
measurement to a particular frequency range, exclude low frequency noise,
or exclude high frequency noise by selecting an appropriate frequency
range. For example, if you want to specifically exclude DC energy from the
measurement, you should set the start frequency to at least 8(df), where df
is the frequency resolution of the FFT used by the SVT SINAD VI. If there
are frequency components of interest below 8(df) that are not DC, reduce
df by increasing the measurement duration. Use the following equation to
determine the relationship between df and the measurement duration,
where T is the measurement duration:
1
df = ---
T
The example illustrated in Figures 14-2 and 14-3 shows how you can
effectively remove DC energy from the measurement and still measure
SINAD with low-frequency test signals.
© National Instruments Corporation 14-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 14 Distortion Measurements
You can use the VI in Figure 14-2 to measure the SINAD for the left
channel of a device like a two-channel audio equalizer. In this example, the
left channel level settings are set to 0 dB for all available octave bands. The
measurement result displayed in Figure 14-3 has a 1 s integration time that
results in a frequency resolution of 1 Hz. The start and stop frequencies are
set to 20 Hz and 20 kHz, respectively. The tone frequency used in the test
is 400 Hz. With these settings, the SINAD measurement includes the entire
audio band and prevents DC energy from giving an inaccurately low result.
The reported SINAD for the left channel of the audio equalizer is
SINAD=73.4 dB, 20 Hz–20 kHz, test amplitude 2 Vpk, test frequency
400 Hz, and unity gain.
∑
P Harmonics + ∑ P Noise
THD+N = ----------------------------------------------------------- (14-1)
∑
P Total
∑
P Harmonics + ∑ P Noise
THD+N = ----------------------------------------------------------- (14-2)
F
© National Instruments Corporation 14-5 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 14 Distortion Measurements
When the fundamental tone dominates the signal, these two definitions
are equivalent. The Sound and Vibration Toolkit uses the definition in
Equation 14-1. In this equation the power in the harmonic components
and noise is divided by the total signal power to compute the THD+N.
Note that when we use the definition in Equation 14-1, the THD+N is the
inverse of the SINAD. When these measurements are expressed in decibels,
the THD+N is the negative of the SINAD.
∑
P Harmonics
THD = ------------------------------
F
THD only includes the energy in the harmonics and does not include the
energy in the broadband noise. Therefore, the measured THD is less than
the measured THD+N for the same input signal.
The SVT THD and Harmonic Components VI allows you to specify the
maximum harmonic used to compute the THD. You can wire –1 to this
input to include all the harmonics up to the Nyquist frequency. If the
maximum harmonic input to the VI is greater than the Nyquist frequency,
the function uses all the harmonics up to the Nyquist frequency and returns
a warning with code 1947. This warning indicates a coerced maximum
harmonic.
You can view any of the measured harmonics by specifying the desired
harmonics as an array to the harmonics to visualize terminal. If you wire
a single-element array with a –1 in the first element to this terminal, all the
measured harmonics will be output to the harmonic components terminal.
Figure 14-5 and Table 14-1 show the magnitude spectrum, measured gain,
THD+N, and THD of the left channel of the same two-channel audio
equalizer used for the example in Figure 14-3.
The VI in Figure 14-4 uses the SVXMPL_One Shot Gain and Distortion
example VI modified by the addition of the SVT THD+N VI. In the
measurement analysis section of this VI, the audio measurements are
applied sequentially to the same channel. The Scale to EU VI scales the
time-domain data to engineering units. The Gain VI computes the complex
spectrum and measures the gain. The SVT THD+N VI then uses the
spectrum to compute the THD+N. The SVT THD and Harmonic
Components VI uses the same spectrum to compute the THD, and the Unit
Conversion VI converts the complex spectrum to a magnitude spectrum for
display.
© National Instruments Corporation 14-7 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 14 Distortion Measurements
Variable Result
Fundamental Frequency 1,000.00 Hz
Fundamental Amplitude 1.82 Vpk
gain –0.80 dB
THD+N –71.69 dB
THD –76.16 dB
h2 0.000021 Vpk
h3 0.000120 Vpk
h4 0.000040 Vpk
h5 0.000059 Vpk
Figure 14-5 illustrates that the measured THD is lower than the THD+N
over the same bandwidth of 20 Hz to 20 kHz. This relationship always
holds true because the THD+N measurement includes the energy in the
harmonics as well as the energy in the noise.
IMD is often used to measure the distortion of the DUT near the
high-frequency limit of the DUT or the measurement system. You also
can arrange the test so that many IMD components occur within the
measurement bandwidth.
The IMD measurement is dictated by the IMD standard used for the
measurement. There are several standard configurations for IMD
measurements, and these configurations use one of two types of IMD
test signals.
© National Instruments Corporation 14-9 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 14 Distortion Measurements
In the first type of IMD test signal, with a low-frequency modulation tone
and a high-frequency carrier tone, the intermodulation components appear
as side bands around the high-frequency carrier tone.
The Figure 14-6 shows the side bands around the 8 kHz carrier tone, as
computed by the following equation:
f2 ± n • f1
Figure 14-6. IMD Test Signal with High-Frequency Carrier Tone Sidebands
In the second type of IMD test signal, with two closely spaced tones
near the high-frequency limit of the measurement bandwidth, the
intermodulation components appear at multiples of the difference
frequency calculated by the following equation:
m • ( f2 – f1 )
Test Signal
IMD Standard (f1, f2, ratio) Typical Applications
SMPTE/DIN (60, 7000, 4) Excite low-frequency
(250, 8000, 4) distortion mechanisms
others (such as thermal distortion
in power amplifiers)
Disk recording and film
ITU-R (CCIF) (11000, 12000, 1) ADC and DAC
(14000, 15000, 1) Slope-induced distortion
(19000, 20000, 1)
others
© National Instruments Corporation 14-11 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 14 Distortion Measurements
The frequencies f1 and f2 are the frequencies, in hertz, of the low- and
high-frequency tones, respectively. The ratio is the ratio of the amplitude
of the low-frequency tone to the amplitude of the high-frequency tone.
Do not expect to get the same measured IMD when performing the
measurement with different standards. It is important that you specify
the test signal used to perform the measurement when you report IMD
measurement results.
Phase Linearity
Phase linearity is one way to measure the amount of group delay distortion
introduced by the DUT. The phase linearity measurement uses a multitone
signal to measure the relative delay of each component through the DUT.
The multitone test signal should be constructed with zero relative phase
between the individual frequency components. In addition, the signal
should be constructed with evenly spaced frequency components given by
the following equation:
fi = i • f
This test signal has a high crest factor. It is important that you construct the
multitone signal with frequency components appropriate for the DUT. In
general, f should be low enough to measure the group delay introduced by
the DUT, as calculated by the following equation:
1
group delay ≤ ---
f
The SVT Phase Linearity VI returns the group delay. The group delay is the
slope of the best linear fit to the plot of phase versus frequency. Next, the
VI computes the difference between the measured phase and the best linear
fit. You can plot this difference against frequency to show the deviation of
the measured phase from the best fit line. Figure 14-9 shows this plot.
© National Instruments Corporation 14-13 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 14 Distortion Measurements
The maximum absolute value of the phase deviation is the value returned
as the phase linearity. In this case the maximum phase deviation value is
0.215°. Phase linearity measurements are always dependent on the
bandwidth of the measurement. To completely report the result, specify that
the phase linearity is 0.215° with a 20 kHz bandwidth, 1 Vpk test amplitude,
and unity gain set for the DUT.
© National Instruments Corporation 15-1 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 15 Single-Tone Measurements
PFI0
AI0
AI1
DUT
AO0
In Out
AO1
Typically, you know the frequency and amplitude of the excitation signal.
The real advantage of connecting the excitation signal to one of the input
channels is the ability to measure the relative phase between the stimulus
and response signals. In a two-channel measurement configuration the
input channel connected to the excitation signal is the stimulus channel,
and the input channel connected to the output of the DUT is the response
channel.
PFI0
AI0
DUT
AI1
In Out
AO0
AO1
Figure 15-3 illustrates where single-tone measurements fit into the sound
and vibration measurement process.
© National Instruments Corporation 15-3 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 15 Single-Tone Measurements
DSA Device
DAQ Device
Data Source WAV File
DAT Recorder
Simulation
Calibration
Calibrate
Sensor
Scale Voltage to
Scaling Engineering Units
Measure
Propagation Delay
Single-Tone
Analysis Measurements
Extended
Measurements
Limit Testing
Waveform Waveform
Visualization Chart
XY Graph
Graph
searches for a response channel tone at the same frequency as the stimulus
tone. The SVT Gain and Phase VI is a dual-channel measurement which
requires you to wire both stimulus and response signals. A single-channel
measurement reads only the response signal.
You have the option of generating the excitation signal from the acquisition
device, a separate output device, or an external source since the stimulus
and response measurements are acquired. If you are using a DSA device,
NI recommends you perform a dual-channel versus a single-channel
measurement when possible. Performing a dual-channel measurement
takes full advantage of the excellent interchannel gain mismatch and
simultaneous sampling of the analog input channels. Refer to the DSA
device specifications for more information. For example, in single-channel
mode the measurement uncertainty when using an NI PXI-4461 is 0.1 dB,
in dual-channel mode the measurement uncertainty is typically 0.01 dB.
Crosstalk
The SVT Crosstalk VI measures the amplification and phase lag between
the reference channel and the idle response channel at the test frequency.
The SVT Crosstalk VI identifies the detected tone amplitude and phase on
the stimulus channel. The stimulus tone amplitude and phase are treated as
references when specifying the idle response channel tone relative
amplitude and relative phase. The SVT Crosstalk VI searches for an idle
response tone at the same frequency as the stimulus tone.
© National Instruments Corporation 15-5 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 15 Single-Tone Measurements
Gain
The SVT Gain VI measures the amplification between the amplitude of the
identified tone in the input signal and the reference amplitude. The VI
automatically identifies the tone with the highest amplitude as the
fundamental tone, unless you wire a value to the expected fundamental
frequency [Hz] terminal. You must connect the amplitude of the excitation
tone to the amplitude input terminal regardless of whether you use the
output of the acquisition device or an external source. Figure 15-5 shows
the required connection.
1 Amplitude Connection
Figure 15-5. Connecting the Amplitude of the Excitation Tone to the Amplitude Input
Idle-Channel Noise
The SVT Idle Channel Noise VI measures the total power within the
specified frequency range. For practical purposes, the SVT Idle Channel
Noise VI is a power-in-band measurement that can operate on waveform
data or the complex spectrum output by the other audio measurements. You
can measure idle-channel noise with the input channel to the DUT floating,
terminated, or excited with a zero amplitude signal.
Dynamic Range
A dynamic range measurement quantifies the distortion and noise found in
the DUT output when a low-level signal is provided. Typically, the
dynamic range is measured with the amplitude of the excitation signal set
to –60 dB FS, where FS is the DUT input channel full-scale amplitude.
© National Instruments Corporation 15-7 LabVIEW Sound and Vibration Toolkit User Manual
Chapter 15 Single-Tone Measurements
Figure 15-6 shows the final measured spectrum for the dynamic range
measurement.
© National Instruments Corporation 15-9 LabVIEW Sound and Vibration Toolkit User Manual
References
A
This appendix lists the reference material used for the Sound and Vibration
Toolkit for LabVIEW. Refer to the following documents for more
information about the theory implemented in this toolkit.
Randall, R.B. 1987. Frequency Analysis. Nærum, Denmark: Brüel & Kjær.
© National Instruments Corporation A-1 LabVIEW Sound and Vibration Toolkit User Manual
Appendix A References
If you searched ni.com and could not find the answers you need, contact
your local office or NI corporate headquarters. Phone numbers for our
worldwide offices are listed at the front of this manual. You also can visit
the Worldwide Offices section of ni.com/niglobal to access the branch
office Web sites, which provide up-to-date contact information, support
phone numbers, email addresses, and current events.
© National Instruments Corporation B-1 LabVIEW Sound and Vibration Toolkit User Manual
Glossary
Numbers/Symbols
2D Two-dimensional.
3D Three-dimensional.
A
acceleration The rate of change of velocity.
anti-aliasing filter To avoid aliasing, analog lowpass filters are used before A/D conversion to
filter out the frequencies greater than half the sampling frequency. Since
they are used to prevent aliasing, these analog lowpass filters are known as
anti-aliasing filters. See also Nyquist frequency.
B
bandedge frequency The upper and lower cutoff frequencies of an ideal bandpass filter.
bandpass filter A filter with a single transmission band extending from a lower bandedge
frequency greater than zero to a finite upper bandedge frequency.
© National Instruments Corporation G-1 LabVIEW Sound and Vibration Toolkit User Manual
Glossary
C
calibrator A controlled source generating a known level of excitation used to calibrate
a sensor.
coherence Gives a measure of the degree of linear dependence between two signals,
as a function of frequency.
coherent output The coherent output power spectrum gives a measure of what part of the
power spectrum (output) power spectrum is fully coherent with the input signal.
crest factor The ratio of the peak value of a signal to its RMS value. For a sine wave,
the crest factor is 1.414. For a square wave, the crest factor is 1.
cross power The cross power spectrum of two signals has an amplitude that is the
spectrum product of the two amplitudes, and a phase that is the difference of the two
phases.
cut-off frequency The frequency at which the filter attenuates the input 3 dB, or half of its
original power.
D
DAQ Data acquisition.
dB Decibels. A logarithmic unit for measuring ratios of levels. If the levels are
specified in terms of power, then 1 dB = 10*log10 (P/Pr) where P is the
measured power and Pr is the reference power. If the levels are specified in
terms of amplitude, then 1 dB = 20*log10(A/Ar) where A is the measured
amplitude and Ar is the reference amplitude.
Distortion The production of signal components not in the original signal due to
non-linearities in the system or transmission path.
Dynamic range The ratio of the largest signal level a circuit can handle to the smallest
signal level it can handle (usually taken to be the noise level), normally
expressed in decibels.
E
ENBW Equivalent noise bandwidth of the time-domain window applied to the
signal in order to perform frequency analysis.
equal confidence Special exponential averaging mode used for fractional-octave analysis.
For equal confidence the time constant for each band is set individually
so that the relative confidence in the measurement is equal across all the
bands. There is a 68% probability that the results will be within ± the
specified confidence level of the true mean value.
equivalent continuous The energy average level of a signal over a given time interval.
level (Leq)
exponential Time-averaging technique that gives recent data more importance than
averaging older data.
F
Fast Exponential averaging using a time constant of 125 ms.
FFT Fast Fourier Transform—an efficient and fast method for calculating the
Discrete Fourier Transform. The Fast Fourier Transform determines the
amplitude and phase of frequency components present in a time domain
digital signal.
© National Instruments Corporation G-3 LabVIEW Sound and Vibration Toolkit User Manual
Glossary
FFT lines The number of FFT lines is related to the FFT block size. Theoretically the
number of lines is half of the block size, but it is practically reduced to 80%
of that value due to the anti-aliasing filter. For example, a 400 line FFT is
based on a block size of 1,024 points.
filtering A type of signal conditioning that allows you to modify the frequency
content of a signal.
fractional-octave The interval between two frequencies, one of which is a fractional power of
two of the other.
frequency response Represents the ratio of output-to-input in the frequency domain, and fully
function characterizes linear, time-invariant systems.
fundamental The frequency of the dominant (highest amplitude) tone in the signal. For
frequency harmonic analysis the fundamental frequency is the greatest common
divisor of the harmonic frequencies.
G
g Unit for measuring acceleration. One g = 9.81 m/s2, the acceleration due to
gravity at the surface of the Earth.
H
H1 The frequency response function computed as the ratio of the cross
spectrum to the input autospectrum: Gxy/Gxx. This technique gives the
best performance in the presence of noise for measuring anti-resonances,
where the signal to noise ratio tends to be poor. For measurement of
resonances, the frequency response function H2, gives a better estimate.
In a noise free environment, both techniques give the same result. Since
both measurements are based on the same data set, the choice of technique
can be made after the data acquisition is completed.
Harmonic A signal whose frequencies are integer multiples of the input signal.
I
IEC International Electrotechnical Commission.
intensity graph A method of displaying three dimensions of data on a 2D plot with the use
of color to indicate the value in the third dimension.
intermodulation The distortion that arises as the result of the modulation between two or
distortion more signals.
L
Leq See equivalent continuous level (Leq).
linear averaging Time-averaging technique that weights all data in the average equally.
© National Instruments Corporation G-5 LabVIEW Sound and Vibration Toolkit User Manual
Glossary
M
MAX Measurement & Automation Explorer.
maximax The absolute maximum of the calculated shock response signal over the
entire signal duration.
microphone Sensor used to convert sound pressure variations into an electrical signal,
usually when the acoustic medium is air.
midband frequency The center frequency of a bandpass filter, defined as the geometric mean of
the bandedge frequencies.
ms Millisecond.
N
noise Any unwanted signal. Noise can be generated by internal sources such as
semiconductors, resistors, and capacitors, or from external sources such
as the AC power line, motors, generators, thunderstorms, and radio
transmitters.
nonstationary signal Signal whose frequency content changes within a captured frame.
Nyquist frequency Half the sampling frequency. Any analog frequency component above
the Nyquist frequency will, after sampling, be converted (aliased) to a
frequency below the Nyquist frequency. See also aliasing and anti-aliasing
filter.
O
octave Refers to the interval between two frequencies, one of which is twice the
other. For example, frequencies of 250 Hz and 500 Hz are one octave apart,
as are frequencies of 1 kHz and 2 kHz.
overlapping A method that uses a portion of the previous data block to compute the FFT
of the current data block.
P
Pa Pascal. International unit of pressure.
peak hold Peak detection process retaining the maximum value of a signal.
phon The unit of loudness on a scale corresponding to the decibel scale of sound
pressure level with the number of phons of a given sound being equal to the
decibels of a pure 1 kHz sine tone judged by the average listener to be equal
in loudness to the given sound.
pink noise Noise for which the spectral energy per octave or any fractional-octave
band is independent of the band. The spectrum looks flat on an octave or
fractional-octave band display.
© National Instruments Corporation G-7 LabVIEW Sound and Vibration Toolkit User Manual
Glossary
R
reference sound A reference pressure of 20E-6 Pa. This reference pressure was
pressure conventionally chosen to correspond to the quietest sound at 1,000 Hz
that the human ear can detect.
reverberation time T60[s]. At a point in an enclosure and for a stated frequency or frequency
band, the time required for the pressure level to decrease by 60 dB after the
source has been stopped.
RMS averaging RMS averaging is used to average the power of a signal. RMS averaging
reduces fluctuations.
S
s Seconds.
S/s Samples/second.
sensor A device that converts a physical stimulus (such as force, sound, pressure,
motion) into a corresponding electrical signal.
settling time The amount of time required for a signal to reach its final value within
specified limits.
Shannon Sampling States that to properly sample a signal, the signal must not contain
Theorem frequencies above the Nyquist frequency.
shock response spectrum A processing method which evaluates the severity of a shock signal.
signal in noise and The ratio of the input signal to the sum of noise and harmonics.
distortion
sound level meter A device used to measure sound pressure levels. Sound level meters usually
consist of a microphone, a preamplifier, a set of standardized frequency
weighting filters, standardized exponential time weighting circuits, a
logarithmic amplitude detector, and a display in decibels.
sound pressure level In decibels, 20 times the base 10 logarithm of the ratio of the sound
pressure, in a stated frequency band, to the reference sound pressure.
spectral leakage A phenomenon whereby the measured spectral energy appears to leak from
one frequency into other frequencies. It occurs when a sampled waveform
does not contain an integral number of cycles over the time period during
which it was sampled. The technique used to reduce spectral leakage is to
multiply the time-domain waveform by a window function. See also
windowing.
spurious free The dynamic range from full-scale deflection to the highest spurious signal
dynamic range in the frequency domain.
SVL Sound and Vibration Library. SVL is a collection of VIs shared by the
Sound and Vibration Toolkit and other NI Toolkits.
T
tach See tachometer.
total harmonic The ratio, in decibels, of the sum of noise and harmonics to the input signal.
distortion plus noise
© National Instruments Corporation G-9 LabVIEW Sound and Vibration Toolkit User Manual
Glossary
time constant A standardized time constant used in exponential time weighting for
acoustical analysis. The standard time constants are Slow = 1,000 ms,
Fast = 125 ms, and Impulse = 35 ms while the signal level is increasing
or 1500 ms while the signal level is decreasing.
total harmonic The ratio of the sum of harmonics to the fundamental tone.
distortion
V
V Volts.
vector averaging Computes the average of complex quantities directly, that is, the real and
imaginary parts are averaged separately. Vector averaging eliminates noise
from synchronous signals and usually requires a trigger.
velocity A vector quantity whose magnitude is a body's speed and whose direction
is the body's direction of motion.
W
waterfall A 3-dimensional plot displaying the amplitude of spectral components as a
function of both time and frequency. The frequency spectrum is displayed
as a curve for each specified time instant. Several such curves (for different
time instants) are displayed simultaneously.
weighting filter Filter used to reproduce the varying sensitivity of the human ear to sound
at different frequencies. Originally, A-weighting was intended to represent
the varying sensitivity of the ear to sound pressure levels ranging between
40 and 60 dB ref 20 µPa. Subsequently, B-weighting and C-weighting were
developed to represent the varying sensitivity of the ear over higher sound
pressure level ranges.
white noise Noise that has the same power spectral density at all frequencies. As an
example, the average power of white noise in a 100 Hz bandwidth between
300 Hz and 400 Hz, is the same as the average power of white noise in the
100 Hz bandwidth between 10,000 Hz and 10,100 Hz.
window function A smooth waveform that generally has zero value at the edges.
See also windowing.
© National Instruments Corporation G-11 LabVIEW Sound and Vibration Toolkit User Manual
Index
© National Instruments Corporation I-1 LabVIEW Sound and Vibration Toolkit User Manual
Index
© National Instruments Corporation I-3 LabVIEW Sound and Vibration Toolkit User Manual
Index
© National Instruments Corporation I-5 LabVIEW Sound and Vibration Toolkit User Manual
Index
O linearity, 14-12
velocity, 6-2
Octave Analysis palette, 1-4, 9-1
phons, 5-3
Extended Measurements palette, 5-1
power in band, 10-29
octave analysis. See fractional-octave analysis
equations for, 10-29
Octave Graph, 1-6
power-in-band measurement, 15-7
overlapping, 11-6
programming examples (NI resources), B-1
time increment control, use in, 11-6
propagation delay, 3-3
P
R
palettes
related documentation, xii
Calibration, 1-3
reset filter control, 9-9
Distortion, 1-5
restart averaging control
Frequency Analysis, 1-4, 10-1
fractional-octave analysis, use in, 9-9
Extended Measurements, 6-1
linear averaging, use in, 8-5, 8-6
front panel display, 1-6
vibration level measurements, use in, 7-4
Generation, 1-4
RMS averaging, 10-24
Integration, 1-4, 6-1
RMS level, 7-2
Limit Testing, 1-3, 4-1
continuous signal acquisition, 7-3
Octave Analysis, 1-4, 9-1
running, 7-4
Scaling, 1-3
single-shot buffered acquisition, 7-3
Single-Tone, 1-6
root mean square. See RMS level
Single-tone Measurements, 15-1
Sound & Vibration, 1-2
Sound Level, 1-4, 8-1 S
Swept Sine, 1-5, 13-1
sampling frequency
Transient Analysis, 1-5, 11-1
for octave filters, 9-12
Vibration Level, 1-4, 7-1
selection of, making sound level
Waterfall Display, 1-5, 12-1
measurements, 8-8
Weighting, 1-3
scaling, 3-1
peak averaging. See peak hold averaging
signals to engineering units, 3-1
peak detection, 7-4
Scaling palette, 1-3, 3-1
peak hold averaging, 8-8
separation. See crosstalk
equation for, 8-8
settling time, 9-7
fractional-octave analysis, use in, 9-9
SFDR, 15-9
frequency analysis, use in, 10-27
Shannon Sampling Theorem, 10-4
periodicity, 10-21
shock response spectrum, 11-2
phase
signal in noise and distortion, 14-3
displacement, 6-2
lag, 15-4
© National Instruments Corporation I-7 LabVIEW Sound and Vibration Toolkit User Manual
Index
SVT STFT versus RPM (analog), 11-9 applying to time-domain data, 5-8
SVT STFT versus Time, 11-8 errors due to uniform corrections, 5-12
SVT Weighting Filter (frequency), 5-13 fractional-octave analysis, use in, 9-12
SVT Weighting Filter (octave), 5-12 purpose of, 5-1
SVT Zoom Power Spectrum, 10-10 weighting filter control, 5-7
weighting filters
See filters
W See weighting
waterfall display, 12-1 weighting mode control, 10-22
autoscaling z-axis, 12-6 Weighting palette, 1-3, 5-1
closing, 12-6 windowing, 10-21
customizing view, 12-6 periodicity, 10-21
defining graph properties, 12-2 supported windows (table), 10-22
definition of, 12-1
example
frequency analysis, 12-3 Z
octave spectra, 12-5 zoom FFT analysis, 10-8
transient analysis, 12-3 controls, 10-11
initializing, 12-2 example, 10-10
procedure for creating display, 12-1 frequency resolution, equation for, 10-9
sending data to, 12-2 settings, 10-11
waterfall window, 12-1 steps in, 10-9
Waterfall Display palette, 1-5, 12-1 VIs, 10-9
Web resources, B-1
weighting
applying to an octave spectrum, 5-12
applying to FFT-based spectrum, 5-13
© National Instruments Corporation I-9 LabVIEW Sound and Vibration Toolkit User Manual