LabVIEW Sound and Vibration Analysis User Manual (2007)
LabVIEW Sound and Vibration Analysis User Manual (2007)
TM
Support Worldwide Technical Support and Product Information ni.com National Instruments Corporate Headquarters 11500 North Mopac Expressway Worldwide Offices Australia 1800 300 800, Austria 43 662 457990-0, Belgium 32 (0) 2 757 0020, Brazil 55 11 3262 3599, Canada 800 433 3488, China 86 21 5050 9800, Czech Republic 420 224 235 774, Denmark 45 45 76 26 00, Finland 385 (0) 9 725 72511, France 33 (0) 1 48 14 24 24, Germany 49 89 7413130, India 91 80 41190000, Israel 972 3 6393737, Italy 39 02 413091, Japan 81 3 5472 2970, Korea 82 02 3451 3400, Lebanon 961 (0) 1 33 28 28, Malaysia 1800 887710, Mexico 01 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 495 783 6851, Singapore 1800 226 5886, Slovenia 386 3 425 42 00, South Africa 27 0 11 805 8197, Spain 34 91 640 0085, Sweden 46 (0) 8 587 895 00, Switzerland 41 56 2005151, Taiwan 886 02 2377 2222, Thailand 662 278 6777, Turkey 90 212 279 3031, United Kingdom 44 (0) 1635 523545 For further support information, refer to the Technical Support and Professional Services appendix. To comment on National Instruments documentation, refer to the National Instruments Web site at ni.com/info and enter the info code feedback. 19992007 National Instruments Corporation. All rights reserved. Austin, Texas 78759-3504 USA Tel: 512 683 0100
Important Information
Warranty
The media is warranted against defects in materials and workmanship for a period of 90 days from the date of shipment, as evidenced by receipts or other documentation. National Instruments will, at its option, repair or replace equipment that proves to be defective during the warranty period. This warranty includes parts and labor. 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 . C USTOMERS 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 owners failure to follow the National Instruments installation, operation, or maintenance instructions; owners modification of the product; owners 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. National Instruments respects the intellectual property of others, and we ask our users to do the same. NI software is protected by copyright and other intellectual property laws. Where NI software may be used to reproduce software or other materials belonging to others, you may use NI software only to reproduce materials that you may reproduce in accordance with the terms of any applicable license or other legal restriction.
Trademarks
National Instruments, NI, ni.com, and LabVIEW are trademarks of National Instruments Corporation. Refer to the Terms of Use section on ni.com/legal for more information about National Instruments trademarks. Other product and company names mentioned herein are trademarks or trade names of their respective companies. Members of the National Instruments Alliance Partner Program are business entities independent from National Instruments and have no agency, partnership, or joint-venture relationship with National Instruments.
Patents
For patents covering National Instruments products, refer to the appropriate location: HelpPatents in your software, the patents.txt file on your CD, or ni.com/patents.
THE SUITABILITY OF NATIONAL INSTRUMENTS PRODUCTS WHENEVER NATIONAL INSTRUMENTS PRODUCTS ARE INCORPORATED IN A SYSTEM OR APPLICATION, INCLUDING, WITHOUT LIMITATION, THE APPROPRIATE DESIGN, PROCESS AND SAFETY LEVEL OF SUCH SYSTEM OR APPLICATION.
Contents
About This Manual
Conventions ...................................................................................................................xi Related Documentation..................................................................................................xii
Chapter 1 Introduction
Sound and Vibration Toolkit .........................................................................................1-1 Toolkit Palettes ..............................................................................................................1-3 S&V Express Measurements...........................................................................1-3 UFF58 File I/O ................................................................................................1-3 Scaling .............................................................................................................1-3 Calibration .......................................................................................................1-3 Limit Testing ...................................................................................................1-4 Weighting ........................................................................................................1-4 Integration........................................................................................................1-4 Generation .......................................................................................................1-4 Vibration Level................................................................................................1-4 Sound Level.....................................................................................................1-5 Octave Analysis...............................................................................................1-5 Frequency Analysis .........................................................................................1-5 Transient Analysis ...........................................................................................1-5 Waterfall Display.............................................................................................1-5 Swept Sine .......................................................................................................1-6 Distortion.........................................................................................................1-6 Single-Tone .....................................................................................................1-6 Front Panel Displays ......................................................................................................1-6 Examples........................................................................................................................1-7
Contents
vi
ni.com
Contents
Chapter 7 Integration
Introduction to Integration .............................................................................................7-1 Implementing Integration ..............................................................................................7-3 Challenges When Integrating Vibration Data .................................................7-5 DC Component .................................................................................7-5 Transducers .......................................................................................7-5 Implementing Integration using the Sound and Vibration Toolkit .................7-5 Time-Domain Integration ..............................................................................................7-6 Single-Shot Acquisition and Integration .........................................................7-6 Continuous Acquisition and Integration..........................................................7-7 Frequency-Domain Integration......................................................................................7-11
vii
Contents
Bandedge Frequencies .................................................................................... 10-5 Fractional-Octave Filters................................................................................. 10-6 Filter Settling Time ......................................................................................... 10-7 Averaging ...................................................................................................................... 10-8 Linear Averaging ............................................................................................ 10-8 Exponential Averaging ................................................................................... 10-8 Equal Confidence Averaging .......................................................................... 10-9 Peak-Hold Averaging...................................................................................... 10-9 Resetting the Filter and Restarting the Averaging Process............................. 10-9 Performing Third-Octave Analysis Outside the Audio Range...................................... 10-9 ANSI and IEC Standards............................................................................................... 10-10 Midband Frequencies ...................................................................................... 10-10 Nominal Frequencies ...................................................................................... 10-11 Compliance with ANSI and IEC Standards .................................................... 10-11 Displaying Results......................................................................................................... 10-11 Octave Extended Measurements ................................................................................... 10-12 Weighting Filters............................................................................................. 10-12 Octave Band Power Value .............................................................................. 10-12 Octave Subset.................................................................................................. 10-13
viii
ni.com
Contents
Vector Averaging ..............................................................................11-17 RMS versus Vector Averaging .........................................................11-18 Peak Hold ..........................................................................................11-19 Weighting Mode..............................................................................................11-20 Using the Frequency Analysis VIs ................................................................................11-20 Available Measurements .................................................................................11-20 Single-Channel Measurements........................................................................11-21 Power Spectrum Measurement .........................................................11-22 Dual-Channel Measurements ..........................................................................11-23 Frequency Response Function Measurement ...................................11-24 Frequency Response Function Measurement Example ....................11-25 Cross Spectrum .................................................................................11-30 Coherence..........................................................................................11-30 Coherent Output Power.....................................................................11-31 Spectrum Extended Measurements................................................................................11-31 Get Spectrum Value ........................................................................................11-32 Get Spectrum Subset .......................................................................................11-32 Power in Band .................................................................................................11-32 Spectrum Peak Search .....................................................................................11-33 Unit Conversion...............................................................................................11-33
ix
Contents
ni.com
Conventions
The following conventions appear in this manual: [] Square brackets enclose the units of measure specified or returned by a parameterfor example, [dB]. The symbol leads you through nested menu items and dialog box options to a final action. The sequence FilePage SetupOptions directs you to pull down the File menu, select the Page Setup item, and select Options from the last dialog box. This icon denotes a note, which alerts you to important information. 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. Italic text denotes variables and cross-references. Text in this font denotes text or characters that you should enter from the keyboard, paths, directories, variables, and filenames and extensions. 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.
italic
monospace
monospace bold
xi
Related Documentation
The following documents contain information that you might find helpful as you read this manual: LabVIEW Help, available by selecting HelpSearch the LabVIEW Help in LabVIEW Getting Started with LabVIEW LabVIEW Fundamentals
xii
ni.com
Introduction
This chapter introduces the Sound and Vibration Toolkit, the toolkit functions and controls palettes, and where to find examples to help you get started.
You can use the Sound and Vibration Toolkit to perform measurements on digitized or simulated data.
1-1
Chapter 1
Introduction
necessarily performed simultaneously. The dashed boxes in Figure 1-1 indicate optional measurement operations.
Data Source
Calibration
Calibrate Sensor
Scaling
Waveform Conditioning
Integration
Weighting Filter
Limit Testing
Octave Analysis
Frequency Analysis
Single-Tone Measurements
Transient Analysis
Analysis
Vibration Level Measurements Weighting Extended Measurements Distortion Measurements
Limit Testing
Visualization
Waveform Chart
Waveform Graph
XY Graph
Waterfall Display
1-2
ni.com
Chapter 1
Introduction
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 4, 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 4, Scaling and Calibration, for information about the calibration process.
1-3
Chapter 1
Introduction
Refer to the LabVIEW Help for information about the individual Calibration VIs.
Limit Testing
The Limit Test Express VI and the SVL Limit Testing VI enable you to apply limit analysis to any type of measured result produced by the Sound and Vibration Toolkit. Refer to Chapter 5, Limit Testing Analysis, for more information about using limit testing to analyze measurement results. Refer to the LabVIEW Help for information about the SVL 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 6, Weighting Filters, for more information about applying weighting to a signal. Refer to the LabVIEW Help for information about the individual Weighting VIs.
Integration
The SVL Integration VI allows you to perform single or double integration on the time-domain signal. Refer to Chapter 7, Integration, for information about the integration process. Refer to the LabVIEW Help for information about the SVL Integration VI.
Generation
The SVT Pink Noise Waveform VI allows you to generate a continuous pink noise waveform.
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 8, Vibration-Level Measurements, for more information about performing vibration level measurements. Refer to the LabVIEW Help for information about the individual Vibration Level VIs.
1-4
ni.com
Chapter 1
Introduction
Sound Level
The Sound Level VIs offer typical sound-level measurements, including equivalent continuous averaging (Leq), exponential averaging, and peak hold. Refer to Chapter 9, 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 10, 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 11, 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 12, 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 buffer the results of frequency analyses before displaying them on a Waterfall Graph, convert Waterfall Graphs to Colormaps, and convert Colormaps to Waterfall Graphs. Refer to the LabVIEW Help for information about the individual Waterfall Display VIs.
1-5
Chapter 1
Introduction
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 14, Distortion 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 15, Single-Tone Measurements, for information about performing distortion analyses. Refer to the LabVIEW Help for information about the individual Distortion VIs.
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 used to measure the linear response and noise characteristics 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.
1-6
ni.com
Chapter 1
Introduction
Examples
The Sound and Vibration Toolkit includes examples to help you get started using the toolkit. Select HelpFind Examples in LabVIEW to launch the NI Example Finder. Select Toolkits and ModulesSound 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 Single-tone measurements Dual-tone measurements Multitone measurements Audio test Limit test Vibration test File input and output Filtering and conditioning
1-7
Dynamic Signals
This chapter discusses how to obtain data to analyze with the LabVIEW Sound and Vibration Toolkit, as well as issues that can affect the quality of your data. You can simulate data with the Generation VIs and with other LabVIEW VIs. 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.
2-1
Chapter 2
Dynamic Signals
Data Source
Calibration
Calibrate Sensor
Scaling
Waveform Conditioning
Integration
Weighting Filter
Limit Testing
Octave Analysis
Frequency Analysis
Single-Tone Measurements
Distortion Measurements
Analysis
Vibration Level Measurements Weighting Extended Measurements Transient Analysis
Limit Testing
Visualization
Waveform Chart
Waveform Graph
XY Graph
Waterfall Display
Figure 2-1. Relationship of the Data Source to the Sound and Vibration Measurement Process
It is important that you keep certain considerations in mind when you obtain your data. The measurement and analysis VIs in the Sound and Vibration Toolkit do not compensate for inaccurate data. Therefore, the test equipment and test procedure should be calibrated to ensure accurate results. Generally, the test equipment should have specifications at least
2-2
ni.com
Chapter 2
Dynamic Signals
10 times better than those of the DUT. Use a verifiable and repeatable test procedure to get accurate results. 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: fNyquist = sample rate/2 When acquiring data with an NI Dynamic Signal Acquisition (DSA) device, aliasing protection is automatic in any acquisition. The sharp anti-aliasing filters on DSA devices track the sample rate and filter out (attenuate) all frequencies above the Nyquist frequency. When performing frequency measurements with an NI E Series DAQ device, you must take steps to eliminate aliasing. These anti-aliasing steps can include the following actions: Increasing the sample rate Applying an external lowpass filter Using an inherently bandlimited DUT
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.
2-3
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 Generation VIs, 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.
2-4
ni.com
Chapter 2
Dynamic Signals
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 Signal Generation VIs, 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.
2-5
The S&V Express Measurements Express VIs are configurable Express VIs designed to significantly reduce the amount of time required to implement sound and vibration measurements. Express VIs appear on the block diagram as expandable nodes with icons surrounded by a blue field. The S&V Express Measurements Express VIs enable you to configure each measurement through an interactive dialog box, the configuration view, which provides visual feedback on the specified measurement parameters.
SignalExpress. Refer to the LabVIEW Help for information about the SignalExpress Express VIs. The Sound and Vibration Toolkit includes Express VIs for specific sound and vibration operations, such as the following: Load from UFF58 Save to UFF58 Sound Level Vibration Level Octave Analysis Power Spectrum Zoom Power Spectrum Peak Search Power in Band Tone Measurements Noise Measurements
3-1
Chapter 3
Refer to the LabVIEW Help for information about the S&V Express Measurements Express VIs.
run the VI. If you view the configuration dialog box a second time without running the VI again, the configuration dialog box displays simulated data.
Simulated Data 100 Hz tone with harmonics and noise in units of g 1 kHz tone with harmonics and noise in units of Pa
3-2
ni.com
Chapter 3
Express VI Power Spectrum Zoom Power Spectrum Power in Band Peak Search Tone Measurements Noise Measurements Frequency Response
Stimuluschirp from 0 Hz to Nyquist in units of V Responsefiltered chirp signal in units of V. The filter is an IIR notch filter centered at 1 kHz.
You can configure most Express VIs to perform measurements on one channel or multiple (N) channels. You can configure this setting on the Input tab of the configuration dialog box. The Express VIs report any errors or warnings at the bottom of the configuration dialog box. The configuration dialog box for these Express VIs also contains a context help window with two portions. The upper portion of the context help window displays basic information about the Express VI. The lower portion of the context help window displays information about specific controls in the configuration dialog box. Click the Show Help button that appears at the top of the configuration dialog box to open the context help window. Click the Hide Help button to open the context help window. Click the OK button to set the configuration of S&V Express Measurements VI. The selected measurement parameters become the default configuration for that instance of the Express VI on the block diagram. Click the Cancel button to revert any changes you made to the configuration of the S&V Express Measurements VI. You can have multiple instances of the same S&V Express Measurements Express VI and each instance can have different configuration settings.
3-3
Chapter 3
3-4
ni.com
This chapter discusses using the SVL Scale Voltage to EU VI to scale a signal to engineering units (EU) and using the Calibration VIs. 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 4-1 illustrates how scaling and calibration fit into the sound and vibration measurement process.
4-1
Chapter 4
Data Source
Calibration
Calibrate Sensor
Scaling
Waveform Conditioning
Integration
Weighting Filter
Limit Testing
Octave Analysis
Frequency Analysis
Single-Tone Measurements
Distortion Measurements
Analysis
Vibration Level Measurements Weighting Extended Measurements Transient Analysis
Limit Testing
Visualization
Waveform Chart
Waveform Graph
XY Graph
Waterfall Display
Figure 4-1. Relationship of Scaling to the Sound and Vibration Measurement Process
Typically, scaling a signal to the appropriate EU occurs before any analysis is performed. Use the SVL Scale Voltage to EU VI to scale the signal to the appropriate EU. All measurement VIs in the Sound and Vibration Toolkit expect input signals and return results with the appropriate units, such as time-domain
4-2
ni.com
Chapter 4
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. The Sound and Vibration Toolkit software maintains engineering units applied by NI-DAQmx through subsequent measurement analysis. In particular, the Sound and Vibration Toolkit handles data scaled to Pascals in NI-DAQmx Sound Pressure measurements and g in NI-DAQmx Acceleration measurements in the same manner as if you used the SVL Scale Voltage to EU VI to scale the data.
4-3
Chapter 4
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 specific device at the desired sample rate. Connect the DSA device output channel directly to the input channel, as displayed in Figure 4-2, to measure the device propagation delay.
PFI0
AI0
AI1
AO0
AO1
Figure 4-2. Measuring the Device Propagation Delay 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 4-3 shows the time domain data for the propagation delay measurement of an NI PCI-6052E.
LabVIEW Sound and Vibration Analysis User Manual 4-4 ni.com
Chapter 4
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 4-4 shows the delay of a smooth pulse generated and acquired by an NI PXI-4461 with a 204.8 kHz sample rate.
Figure 4-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 4-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.
4-5
Chapter 4
As illustrated by Figures 4-3, 4-4, and 4-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
4-6
ni.com
Chapter 4
examples in the Sound and Vibration Toolkit account for this transient behavior in the device response to achieve the highest degree of accuracy. The propagation delay of the DUT is also an important specification in some applications. For example, the propagation delay for the DUT is a required input when performing audio measurements and when measuring the frequency response using swept sine. If the DUT and the propagation medium can successfully pass the pulse signal used by the SVL Measure Propagation Delay VIs without excessive attenuation, then this measurement also applies when measuring the propagation delay of the DUT and the propagation medium. Figure 4-6 shows the wiring diagram for this configuration.
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.
4-7
Chapter 4
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: delay[s] = delay[samples]/sample rate[Hz]
Note The swept sine VIs expect the DUT propagation delay measurement in seconds and
4-8
ni.com
This chapter discusses using the Limit Test Express VI, and the polymorphic SVL Limit Testing VI. Use the Limit Test Express VI to interactively configure a limit test on the following measurement results: Waveform Spectrum Scalar
Refer to Chapter 3, S&V Express Measurements, and to the LabVIEW Help for more information about the Limit Test Express VI. Use the SVL 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 SVL Limit Testing VI instances.
5-1
Chapter 5
Data Source
Calibration
Calibrate Sensor
Scaling
Waveform Conditioning
Integration
Weighting Filter
Limit Testing
Octave Analysis
Frequency Analysis
Single-Tone Measurements
Distortion Measurements
Analysis
Vibration Level Measurements Weighting Extended Measurements
Limit Testing
Visualization
Waveform Chart
Waveform Graph
XY Graph
Figure 5-1. Relationship of Limit Testing to the Sound and Vibration Measurement Process
5-2
ni.com
Chapter 5
You can use the SVL Limit Testing VI to analyze almost any measured result produced by the Sound and Vibration Toolkit. Refer to Table 5-1 for examples of datatypes supported by the SVL Limit Testing VI and VIs that generate supported datatypes.
Table 5-1. Compatible Data Types for SVL Limit Testing VI
Output VIs AI Read, DAQmx Read, Waveform Generation, Weighting, Integration, Vibration Level, Sound Level
Baseband FFT, Baseband Subset FFT, Zoom FFT, Spectrum Extended Measurements
XY Data
Peak Measurement
Distortion, Single-Tone, Spectrum Extended Measurements Calibration, Vibration Level, Sound Level, Octave, Distortion, Single-Tone, Spectrum and Octave Extended Measurements
Scalar Measurement
5-3
Chapter 5
is the signal is the upper mask limit is the lower mask limit
5-4
ni.com
Chapter 5
Criteria on Input Limit Masks dt > 0 dtS = dtU = dtL NS = NU = NL f0S = f0U = f0L dfS = dfU = dfL NS = NU = NL [X]S = [X]U = [X]L NS = NU = NL f(i)S = f(i)U = f(i)L NS = NU = NL
Octave Spectrum, Swept Sine Spectrum, XY data ([X], [Y]) Identified Peaks, Harmonic Components, Multitone Phases [(frequency, amplitude)]
Limit testing covers a broad range of data testing from range detection to discontinuous mask testing of a swept-sine frequency response spectrum. Figures 5-2, 5-4, 5-6, and 5-8 illustrate some, but not all, of the different ways you can use the SVL Limit Testing VI in your application.
Figure 5-2 illustrates a range-detection test. Scaled waveform data and upper and lower limits are input to the SVL Limit Testing VI. The VI checks that the data falls within the envelope specified by the upper and lower limits. Figure 5-2 shows the output results for the range detection test.
5-5
Chapter 5
Figure 5-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 5-5 shows the THD test output results.
5-6
ni.com
Chapter 5
Figure 5-6 shows a continuous mask test on a power spectrum. Formula nodes define both the upper and lower limits in this VI, making this a more complex test than the one in Figure 5-4. Figure 5-7 shows the output graph for the power spectrum continuous mask test.
5-7
Chapter 5
5-8
ni.com
Chapter 5
Figure 5-8 shows a discontinuous mask test on a swept-sine frequency response. A discontinuous mask test can track and test the results at different magnitudes and ranges, as well as stop testing at defined intervals. For example, you might use the envelope defined by the upper and lower limit masks in this example for a DUT such as a notch filter. Figure 5-9 shows the output graph for the discontinuous mask test.
5-9
Weighting Filters
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. Use the Sound Level Express VI to apply A-, B-, or C-weighting to a sound level measurement. Refer to Chapter 3, S&V Express Measurements, and to the LabVIEW Help for more information about the Sound Level Express VI. The following Sound and Vibration Toolkit palettes contain weighting filter VIs: Weighting Octave Extended Measurements Spectrum Extended Measurements S&V Express Measurements
Refer to the LabVIEW Help for information about individual weighting filter VIs.
6-1
Chapter 6
Weighting Filters
The frequency response of this microphone is designed to be as flat as possible in the frequency range of 10 Hz to 10 kHz. Compare this frequency response with the equal loudness curves shown in Figure 6-2.
100
Intensity (dB)
80
60
40
20
10
100
10,000
6-2
ni.com
Chapter 6
Weighting Filters
Loudness is a subjective indicator of the perceived noise level expressed in phons. The loudness level in phons is the sound pressure level in decibels of a 1 kHz tone having the same perceived loudness as the tone being evaluated. Thus, a 1 kHz tone with a loudness level of 30 phons is equally as loud as a 1 kHz tone with a sound pressure level of 30 dB referenced to 20 Pa. However, a signal frequency of 100 Hz requires a sound pressure level of 44 dB referenced to 20 Pa to provide the same loudness level.
6-3
Chapter 6
Weighting Filters
The frequency responses of these filters are standardized according to the ANSI S1.4 standard. Figure 6-3 shows the relative attenuation defined for A-, B-, and C-weighting filters.
Figure 6-3. Relative Attenuation of A-, B-, and C-Weighting Filters Note Each of these filters has a relative attenuation of 0 dB at 1,000 Hz.
6-4
ni.com
Chapter 6
Weighting Filters
6-5
Chapter 6
Weighting Filters
6-6
ni.com
Chapter 6
Weighting Filters
Data Source
Calibration
Calibrate Sensor
Scaling
Waveform Conditioning
Weighting Filter
Analysis
Single-Tone Measurements
Frequency Analysis
Octave Analysis
Distortion Measurements
Transient Analysis
Weighting
Weighting
Limit Testing
Visualization
Waveform Chart
Waveform Graph
XY Graph
Waterfall Display
Figure 6-6. Relationship of Weighting to the Sound and Vibration Measurement Process
Traditionally, weighting filters are built using analog components. If you use an external weighting filter, use the weighting filter parameter in the channel info control of the SVL Scale Voltage to EU VI to ensure proper display of the selected units. Also, the weighting VIs use the channel info parameter to report an error if the application attempts to apply additional weighting to a previously weighted signal.
6-7
Chapter 6
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. Refer to Chapter 4, Scaling and Calibration, for information about scaling a signal.
Supported Filters Sample Rates 4 kHz to 20 kHz 4 kHz, 8 kHz, 10 kHz, 11.025 kHz, 12.8 kHz, 20 kHz to 1 MHz 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 A, B, C-weighting Yes ITU-R 468-4/Dolby CCITT/C-message Yes
Yes
Yes
Yes
6-8
ni.com
Chapter 6
Weighting Filters
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 to determine the maximum frequency within tolerances. The SVT A, B, C Weighting Filter (Fixed Rates) VI supports the sample rates listed in Table 6-1. For all supported sample rates the VI achieves Type 0 compliance to the Nyquist frequency.
6-9
Chapter 6
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. 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. The SVT A, B, C Weighting Filter (Fixed Rates) VI supports the sample rates listed in Table 6-1. For all supported sample rates the VI achieves Class 1 compliance to the Nyquist frequency.
The SVT Telecommunications Weighting Filter (Fixed Rates) VI supports the sample rates listed in Table 6-1. For all supported sample rates, the VI achieves compliance to the Nyquist frequency.
6-10
ni.com
Chapter 6
Weighting Filters
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 6-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.
6-11
Chapter 6
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.
6-12
ni.com
Chapter 6
Weighting Filters
A-weighting filter is used at frequencies below 500 Hz. At frequencies below 500 Hz, the slope of the A-weighting curve is steep. Figure 6-10 shows how the slope of the A-weighting curve can contribute to potential measurement errors at low frequencies.
40
2k
4k
8k
Figure 6-10. Potential Measurement Error for A-Weighting at Low Frequency Note The same type of measurement error in Figure 6-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.
6-13
Chapter 6
Weighting Filters
Figure 6-12 shows a different implementation based on applying the weighting filter on the time-domain signal and then computing the power spectrum.
In Figure 6-11, a single power spectrum is computed. In Figure 6-12, the power spectrum is computed twice, leading to more CPU usage and increased processing time. By applying the weighting filter in the frequency domain, as in Figure 6-11, you can decrease CPU usage and processing time.
6-14
ni.com
Integration
7
This chapter discusses the integration process, including basic theory and implementation in the time and frequency domains. The Sound and Vibration Toolkit contains the following VIs for integration: SVL Integration VI for time-domain integration SVL Integration (frequency) VI for frequency-domain integration Vibration Level Express VI for interactively converting acceleration data to velocity and displacement data, and configuring vibration level measurements. Refer to Chapter 3, S&V Express Measurements, for more information about the Vibration Level Express VI.
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: . dx = ---- ( x ) dt
2 .. d- . dx = ---- x = ------ ( x ) 2 dt dt
Velocity is the first derivative of displacement with respect to time. Acceleration is the first derivative of velocity and the second derivative of displacement with respect to time. Therefore, given acceleration, perform a single integration with respect to time to compute the velocity or perform a double integration with respect to time to compute the displacement.
7-1
Chapter 7
Integration
When representing the acceleration of a point by a simple sinusoid, the velocity and the displacement of the point are well known and represented by the following equations: a = A sin ( t ) A A v = --- cos ( t ) = --- sin t -- 2 AAd = ----- sin ( t ) = ----- sin ( t ) 2 2
Note Initial condition is arbitrarily set to zero in Equations 7-1 and 7-2.
(7-1)
(7-2)
The amplitude of the velocity is inversely proportional to the frequency of vibration. The amplitude of the displacement is inversely proportional to the square of the frequency of vibration. Furthermore, the phase of the velocity lags the acceleration by 90. The phase of the displacement lags the acceleration by 180. Figure 7-1 illustrates the relationship between acceleration, velocity, and displacement.
7-2
ni.com
Chapter 7
Integration
The integration of a sinusoid is known in closed form. Integration of an arbitrary waveform typically requires a numerical approach. You can use several numerical integration schemes to evaluate an integral in the time domain. In the frequency domain, you can define any arbitrary band-limited waveform as a sum of sinusoids. Because the amplitude and phase relationships are known for sinusoids, you can carry out the integration in the frequency domain, as in Figures 7-9 and 7-10.
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 7-2 illustrates how integration fits into the sound and vibration measurement process.
7-3
Chapter 7
Integration
Data Source
Calibration
Calibrate Sensor
Scaling
Waveform Conditioning
Integration
Analysis
Frequency Analysis
Transient Analysis
Integration
Extended Measurements
Limit Testing
Visualization
Waveform Chart
Waveform Graph
Waterfall Display
Figure 7-2. Relationship of Integration to the Sound and Vibration Measurement Process
7-4
ni.com
Chapter 7
Integration
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 23 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.
7-5
Chapter 7
Integration
Time-Domain Integration
This section presents examples of and discussion about time-domain integration.
Because the integration is implemented with filters, there is a transient response associated with integration while the filters settle. You should take care to avoid the transient region when making further measurements. Figure 7-4 shows the results of a single-shot acquisition and integration of a 38 Hz sine wave. You can see the transient response in the first 200 ms of the integrated signal.
7-6
ni.com
Chapter 7
Integration
7-7
Chapter 7
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 SVL Integration VI automatically resets the first time it is called and runs continuously thereafter. If you use the block diagram in Figure 7-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 7-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.
As in single-shot acquisition and integration, continuous acquisition and integration has an initial transient response. Take care to avoid making additional measurements until the response of the filters settles. Once the filters settle, you can use the integrated signals for additional analysis.
7-8
ni.com
Chapter 7
Integration
Figure 7-7 shows the frequency response for time-domain single integration. Figure 7-8 shows the frequency response for time-domain double integration.
7-9
Chapter 7
Integration
In Figure 7-7, you can see the characteristic 20 dB per decade roll-off of the magnitude response of the single integration. In Figure 7-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.
7-10
ni.com
Chapter 7
Integration
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.
The following example demonstrates the implementation of the strategies used to obtain the spectrum of an integrated signal. Figure 7-9 shows the block diagram for the example VI.
Figure 7-9. Integration in the Time Domain and in the Frequency Domain
The highpass cutoff frequency parameter of the SVL Integration VI is wired with a constant of 10 Hz. The SVL Integration (frequency) VI does not have a highpass cutoff frequency parameter. Instead, the SVL Integration (frequency) VI sets the DC component of the integrated signal to zero if the spectrum scale is linear or to negative infinity (Inf) if the spectrum scale is in decibels.
7-11
Chapter 7
Integration
Figure 7-10 shows the results of integrating in the time and frequency domains.
The power spectrum is computed after the time-domain integration filters settle. The frequency-domain integration scales the spectrum at each frequency line. No settling time is necessary for the frequency-domain integration because integration filters are not involved in the frequency-domain integration. Perform frequency-domain integration in the following situations to maximize performance: When the integrated signal is not needed in the time domain When spectral measurements are made
7-12
ni.com
Vibration-Level Measurements
This chapter briefly discusses the analysis concepts associated with performing vibration-level measurements and how you can use the Vibration Level Express VI and the Vibration Level VIs to perform vibration-level measurements. Use the Vibration Level Express VI to quickly develop and interactively configure the following vibration level measurements: RMS Running RMS Peak Max-Min Exponential average
You can also use the Vibration Level Express VI to export a single-integrated or double-integrated signal. Refer to Chapter 3, S&V Express Measurements, and to the LabVIEW Help for more information about the individual Vibration Level VIs. Figure 8-1 illustrates how vibration-level measurement fits into the sound and vibration measurement process.
8-1
Chapter 8
Vibration-Level Measurements
Data Source
Calibration
Calibrate Sensor
Scaling
Waveform Conditioning
Integration
Analysis
Limit Testing
Visualization
Waveform Chart
Waveform Graph
Figure 8-1. Relationship of Vibration-Level Measurements to the Sound and Vibration Measurement Process
8-2
ni.com
Chapter 8
Vibration-Level Measurements
The sampling frequency is 10 kS/s. A buffer containing 1 s of data is returned by the read VI.
8-3
Chapter 8
Vibration-Level Measurements
Note Set the restart averaging control on the SVL RMS Level VI to TRUE. Otherwise,
the SVL 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.
8-4
ni.com
Chapter 8
Vibration-Level Measurements
8-5
Sound-Level Measurements
This chapter discusses some of the analysis concepts associated with performing sound-level measurements and how you can use the Sound Level VIs to perform sound-level measurements. You can combine different sound-level measurements and use them simultaneously to provide flexibility with acoustic measurements. Use the Sound Level Express VI to develop and interactively configure the following sound level measurements: Leq Running Leq Peak Exponential average
You also can use the Sound Level Express VI to export the A-, B-, or C-weighted signal. Refer to Chapter 3, S&V Express Measurements, and to the LabVIEW Help for more information about the individual Sound Level VIs. When combined with any National Instruments DSA device, the appropriate microphone, and proper signal conditioning, the Sound Level VIs are compliant with the IS0/IEC 61672-1:2002 standard.
9-1
Chapter 9
Sound-Level Measurements
Figure 9-1 illustrates how sound-level measurement fits into the sound and vibration measurement process.
DSA Device DAQ Device WAV File DAT Recorder Simulation
Data Source
Calibration
Calibrate Sensor
Scaling
Waveform Conditioning
Weighting Filter
Analysis
Limit Testing
Visualization
Waveform Chart
Waveform Graph
Figure 9-1. Relationship of Sound-level Measurements to the Sound and Vibration Measurement Process
9-2
ni.com
Chapter 9
Sound-Level Measurements
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.
2 1 P rms L eq = 10log10 -- -------------- dt T P2 0 0 T
The single-shot acquisition is configured with the following parameters: Sampling frequency fs = 51,200 samples per second (S/s). Buffer size = 51,200 samples, which is one second of data at the specified sampling frequency.
9-3
Chapter 9
Sound-Level Measurements
Figure 9-3 displays the resulting Leq measurement and the instantaneous sound pressure level.
Figure 9-3. Leq and Instantaneous Sound Pressure Level versus Time
9-4
ni.com
Chapter 9
Sound-Level Measurements
The memory required to accumulate such a large number of samples might prohibit you from accumulating the number of samples required for your analysis. An alternate method is to process small chunks of data, keep track of the intermediate results, and integrate the measurement over time. You can use the Sound and Vibration Toolkit to perform this alternate method of Leq measurement. Figure 9-4 shows the block diagram for a VI designed to compute the Leq for a one-hour period.
The VI in Figure 9-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 SVL 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.
9-5
Chapter 9
Sound-Level Measurements
Because you divide the 10 s measurement period into 50 ms blocks, you need 200 iterations (10 s/50 ms) of the For Loop. To continuously measure the Leq, set the restart averaging control to TRUE. When the restart averaging control is set to TRUE, the SVL Leq Sound Level VI does not accumulate intermediate results but restarts the averaging process with each iteration of the For Loop. Figure 9-6 shows the results of performing a reverberation time measurement for 10 s.
9-6
ni.com
Chapter 9
Sound-Level Measurements
Exponential Averaging
Exponential averaging is a continuous averaging process that weighs 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. The exponential averaging mode supports the following time constants: Slow uses a time constant of 1,000 ms. Slow averaging is useful for tracking the sound pressure levels of signals with sound pressure levels that vary slowly. Fast uses a time constant of 125 ms. Fast averaging is useful for tracking the sound pressure of signals with sound pressure levels that vary quickly. Impulse uses a time constant of 35 ms if the signal is rising and 1,500 ms if the signal is falling. Impulse averaging is useful for tracking sudden increases in the sound pressure level and recording the increases so that you have a record of the changes. Custom allows you to specify a time constant suitable for your particular application.
Use the following sound-level measurement VIs for exponential time averaging: SVL Exp Avg Sound Level VI SVL Decimated Exp Avg Sound Level VI
Refer to the LabVIEW Help for information about the SVL Exp Avg Sound Level VI and the SVL Decimated Exp Avg Sound Level VI. Figure 9-7 illustrates a block diagram for a VI that performs an exponential averaged measurement using the SVL Exp Avg Sound Level VI with a Slow time constant. Use the restart averaging control to reset the averaging.
9-7
Chapter 9
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.
9-8
ni.com
Fractional-Octave Analysis
10
This chapter discusses fractional-octave analysis, including fractional-octave analysis theory, averaging modes supported by the Octave Analysis VIs, settling time, and ANSI and IEC standards. Use the Octave Express VI to develop and interactively configure 1/1, 1/3, 1/6, 1/12, and 1/24 fractional-octave analyses. Refer to Chapter 3, S&V Express Measurements, and to the LabVIEW Help for more information about the Octave Express VI. Use the Octave Analysis VIs 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. These VIs offer compliance with both the ANSI S1.11-2004, class 1 and the IEC 1260:1995, class 1 standards. Refer to the LabVIEW Help for information about the individual Octave Analysis VIs.
10-1
Chapter 10
Fractional-Octave Analysis
Figure 10-1 illustrates the relationship of fractional-octave analysis to the sound and vibration measurement process.
DSA Device DAQ Device WAV File DAT Recorder Simulation
Data Source
Calibration
Calibrate Sensor
Scaling
Waveform Conditioning
Weighting Filter
Analysis
Octave Analysis
Weighting
Limit Testing
Visualization
XY Graph
Waterfall Display
Figure 10-1. Relationship of Fractional-Octave Analysis to the Sound and Vibration Process
10-2
ni.com
Chapter 10
Fractional-Octave Analysis
The number of filters to apply, as well as their passband, is determined by the type of fractional-octave measurement you are making, 1/1 through 1/24 octave. The filters used are known as constant Q filters. Refer to the Bandedge Frequencies section of this chapter for more information about constant Q filters.
10-3
Chapter 10
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 10-3 shows this bar graph.
10-4
ni.com
Chapter 10
Fractional-Octave Analysis
A 0
1 Octave
Just as a tuning fork provides a 440 Hz reference frequency for musicians, a reference frequency is needed to fully define octave filters. Instrumentation equipment typically uses a 1 kHz reference frequency.
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: 1,000 f 1 = ------------ = 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.
10-5
Chapter 10
Fractional-Octave Analysis
An octave filter with a center frequency of 8,000 Hz leads to the following bandedge frequencies: 8,000 f 1 = ------------ = 5,657 Hz 2 f 2 = ( 8,000 ) ( 2 ) = 11,314 Hz BW = f 2 f 1 = 5,657 Hz Q = 5,657 = 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 when the center frequency is low. The bandwidth of the octave filters is wider when the center frequency is higher.
Because of the bandwidth characteristics, fractional-octave analysis uses a logarithmic frequency scale to compute and display octave spectra.
Fractional-Octave Filters
Octave filter resolution is limited, being that there are only 11 octaves in the 16 Hz16 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 10-5 shows the 1/3 octave response at frequencies of 500 Hz, 630 Hz, and 800 Hz.
10-6
ni.com
Chapter 10
Fractional-Octave Analysis
Fractional-octave analysis is a CPU-intensive operation. Increasing the number of filters applied to a signal increases the demands placed on the CPU and can result in increased computation time.
10-7
Chapter 10
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 9, 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 10-6 illustrates the block diagram for a VI performing 1/3 octave analysis in the 20Hz20kHz range using fast exponential averaging.
Figure 10-6. 1/3 Octave Analysis with Fast Exponential Averaging VI 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.
10-8
ni.com
Chapter 10
Fractional-Octave Analysis
Refer to Chapter 9, Sound-Level Measurements, for more information about exponential averaging.
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 9, Sound-Level Measurements, for more information about peak-hold averaging.
10-9
Chapter 10
Fractional-Octave Analysis
The fractional-octave analysis VIs accept numeric values for the low-frequency and high-frequency limits and enable the measurement of signals outside the 20 Hz20 kHz range.
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.
Midband Frequencies
According to the IEC 1260:1995 and the ANSI S1.11:2004 standards, the midband frequency, or center, frequency of the bandpass filter is defined by the following equations. f i = 1000 2 f i = 1000 2 where f i is the center frequency of the i th 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.
ib
( i + 1 )b -----------------2
10-10
ni.com
Chapter 10
Fractional-Octave Analysis
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 and the ANSI S1.11:2004 standards.
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 LabVIEW 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.
10-11
Chapter 10
Fractional-Octave Analysis
To modify an XY graph to display octave results, complete the following steps: 1. 2. 3. 4. 5. 6. 7. 8. Right-click the XY graph. Select X ScaleMappingLogarithmic from the shortcut menu. Right-click the Plot Legend. Select Bar Plots and the style for the bar graph from the shortcut menu. Right-click the Plot Legend. Select Fill Base Line-Infinity from the shortcut menu. Right-click the Plot Legend. Select Interpolation and the style with no line connecting the points.
Refer to Chapter 5, Limit Testing Analysis, for information about limit testing. Refer to Chapter 6, Weighting Filters, for information about octave weighting.
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 6, Weighting Filters, for more information about weighting filters.
10-12
ni.com
Chapter 10
Fractional-Octave Analysis
returns the band power value of the octave spectrum for the specified input frequency band. If the input frequency does not exactly coincide with a nominal frequency, the VI finds the octave band that includes the input frequency. Using the results of this calculation, you can use the band power value to verify or perform limit testing on the power of that specified frequency band.
Octave Subset
If you have a range of frequency bands of interest, and you want to analyze a subset of the octave spectrum, you can use the SVT Get Octave Subset VI. The SVT Get Octave Subset VI extracts the subset of the octave spectrum specified by a frequency range and limits subsequent spectrum operations to the frequency range of interest. This VI also computes the total band power of the returned frequency bands. You can use this information to compare how the total band power varies at different frequency ranges and to see the effect of a frequency range on the whole signal. You also can apply different weighting factors to different frequency ranges. The SVT Get Octave Subset VI is useful in applications where you need to combine fractional octave bands, for example, in loudness measurements. The following example shows an application where you divide the full octave spectrum into three different frequency ranges and determine the total band power for each range of frequencies. Figure 10-8 shows the full octave spectrum and the three different frequency ranges.
10-13
Chapter 10
Fractional-Octave Analysis
Using the SVT Get Octave Subset VI, you can divide the octave spectrum into three rangeslow, middle, and high frequency ranges. The block diagram code in Figure 10-9 defines the low range from 20 Hz to 200 Hz, the middle range from 200 Hz to 2000 Hz, and the high range from 2000 Hz to 20,000 HZ.
The total band powers indicator returns the band power value for each range of frequencies. The band power values for the low, middle, and high frequency ranges are 50.6 dB, 88.3 dB, and 49.1 dB, respectively. If you add the band power value for each range together, the result is the total band power value, 88.3 dB, for the full octave spectrum.
10-14
ni.com
Frequency Analysis
11
This chapter discusses methods used by the Frequency Analysis VIs for windowing, averaging, and performing frequency-domain measurements. Use the Express VIs on the Frequency Analysis palette to develop and interactively configure the following measurements: Power spectrum Zoom power spectrum Frequency response Peak search Power in band
Refer to the LabVIEW Help for information about individual Frequency Analysis VIs. The Frequency Analysis VIs offer various frequency measurements based on the discrete Fourier transform (DFT). Refer to the LabVIEW Help for information about the DFT and Fast Fourier Transform (FFT).
Note For simplicity, the remainder of this document uses the term FFT to denote both the
11-1
Chapter 11
Frequency Analysis
Figure 11-1 illustrates the relationship of frequency analysis to the sound and vibration measurement process.
DSA Device DAQ Device WAV File DAT Recorder Simulation
Data Source
Calibration
Calibrate Sensor
Scaling
Waveform Conditioning
Integration
Weighting Filter
Analysis
Limit Testing
Visualization
Waveform Graph
Waterfall Display
Figure 11-1. Relationship of Frequency Analysis to the Sound and Vibration Measurement Process
11-2
ni.com
Chapter 11
Frequency Analysis
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. Compute N with the following equation: 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: fs 1 df = -- = --T N
11-3
Chapter 11
Frequency Analysis
where df is the frequency resolution. T is the acquisition time. fs is the sampling frequency. N is the number of samples.
Note The frequency resolution is determined solely by the acquisition time.
11-4
ni.com
Chapter 11
Frequency Analysis
Frequency Domain fNyquist Time Domain 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 11-2 if you prefer to specify the spectrum parameters and determine the required data-acquisition parameters from these specifications.
Table 11-2. Frequency-Domain Specifications to Time-Domain Parameters
fmax fs Eb
# lines Eb N
df fs 1 -- = --N T
fs N T
fs -2
Time Domain fs Frequency Domain 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. fmax # lines df f max -------Eb N # lines --------------Eb T 1 # lines ---- = --------------df f max
11-5
Chapter 11
Frequency Analysis
In Tables 11-1 and 11-2, Eb is the ratio of the alias-free bandwidth to 1 the sampling frequency. In traditional instruments, Eb is typically ---------- . 2.56 However, the alias-free bandwidth depends on the hardware used to acquire the dynamic signal. Eb can have a maximum value of 0.5. This maximum value corresponds to a perfect anti-aliasing filter. For example, if 1,024 samples are input to the FFT algorithm, the computed spectrum has 512 non-DC spectral lines. The computed spectrum has a total of 513 lines including the DC component. Acquire these same 1,024 samples on an 1 instrument with a standard ---------- effective bandwidth, then use the equation 2.56 in Table 11-1 to find the expected number of alias-free lines in the computed spectrum. Complete the necessary calculations with the following equations: # lines = E b N 1# lines = --------- 1,024 = 400 lines 2.56 Eb is entirely a hardware property. However, mathematically, you can use the FFT to compute the frequency spectrum up to the Nyquist frequency. Remember to account for the presence or absence of an anti-aliasing filter when performing frequency analysis. The Frequency Analysis VIs compute every spectral line, alias-free or not. Use the frequency range to limit the analysis to the alias-free region of the spectrum with FFT subset and zoom FFT measurements. Use the SVFA Get Spectrum Subset VI to limit the analysis to the alias-free region of the spectrum with baseband FFT measurements.
Note Table 11-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 11-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.
11-6
ni.com
Chapter 11
Frequency Analysis
Increase the acquisition time to 1 s to achieve a frequency resolution of 1 Hz. Figure 11-3 shows the results obtained with a 1 s acquisition time. You can distinguish the individual tones with the increased acquisition time.
11-7
Chapter 11
Frequency Analysis
The following strategies achieve a finer frequency resolution: Decrease the sampling frequency, fs. Decreasing fs usually is not practical because decreasing fs reduces the frequency range. Increase the number of samples, N. Increasing N yields an increased number of lines over the original frequency range.
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 Windowing 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.
11-8
ni.com
Chapter 11
Frequency Analysis
The Zoom FFT VIs complete the following steps to process the sampled data: 1. 2. Modulate the acquired data to center the analysis band at 0 frequency. 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. Decimate the filtered data to reduce the effective sampling frequency. Accumulate the decimated data until sufficient samples are available to compute the spectrum. Use the Discrete Zak Transform (DZT) to efficiently compute the desired spectral lines. Demodulate, or shift, the computed spectrum.
3. 4. 5. 6.
11-9
Chapter 11
Frequency Analysis
Zoom Measurement
The following example demonstrates a zoom power spectrum measurement. Acquire a sine wave at 1,390 Hz with an National Instruments DSA device and the VI displayed in Figure 11-4.
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 SVFA Zoom Power Spectrum VI to analyze a narrow band with a much finer frequency resolution. Figure 11-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: 2,000 1,000 Hz ----------------------------------------- = 2.5 Hz 400 lines
11-10
ni.com
Chapter 11
Frequency Analysis
Use the Zoom Power Spectrum Express VI to interactively configure zoom measurements. Refer to Chapter 3, S&V Express Measurements, and to the LabVIEW Help for more information about the Zoom Power Spectrum Express VI.
Zoom Settings
Figure 11-6 shows the zoom settings control used to acquire the zoom measurement results displayed in Figure 11-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.
11-11
Chapter 11
Frequency Analysis
Subset Analysis
The Baseband Subset VIs 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 Baseband FFT VIs. The acquisition time determines the frequency resolution. The only way to achieve a finer frequency resolution is to increase the length of the time record. In the case of baseband or subset analysis, a longer time record implies a larger block size. 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 offers better processing performance than the Baseband FFT VI can provide.
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. Windowing multiplies the time-domain data by a window function before the FFT is performed. Window functions typically have a value of zero at the start and end of the measurement period. Figure 11-7 shows how a signal that is not the same at the start and end of the measurement period
11-12
ni.com
Chapter 11
Frequency Analysis
appears not to be periodic and how this same signal becomes periodic when multiplied by a window function.
The windows supported by the Sound and Vibration Toolkit and their equivalent noise bandwidths (ENBW) are listed in Table 11-3. ENBW is a property of the window applied to the signal.
Table 11-3. Supported Time-Domain Windows
11-13
Chapter 11
Frequency Analysis
Window Flat Top 4 Term Blackman-Harris 7 Term Blackman-Harris Low Sidelobe Blackman Nuttall Triangle Barlett-Hanning Bohman Parzen Welch Kaiser Dolph-Chebyshev Gaussian Force-Exponential
ENBW 3.77 2.00 2.63 2.22 1.98 1.33 1.46 1.79 1.92 1.20 3.85 1.40 1 N/A
Note You can find a complete list of supported time-domain window in the window
11-14
ni.com
Chapter 11
Frequency Analysis
Averaging Parameters
Each of the Frequency Analysis VIs supports averaging. The averaging parameters control in Figure 11-8 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.
11-15
Chapter 11
Frequency Analysis
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
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. RMS averaged measurements are computed according to the following equations: FFT spectrum power spectrum cross spectrum frequency response X X X X X Y X Y H1 = ------------------- X X Y Y H2 = ------------------- Y X
11-16
ni.com
Chapter 11
Frequency Analysis
( H1 + H2 ) H3 = ------------------------2 coherence coherent output 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.
2 X Y = ----------------------------------------- X Y Y X 2
COP = Y Y
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. Vector averaged measurements are computed according to the following equations: FFT spectrum power spectrum cross spectrum frequency response where X is the complex FFT of the stimulus signal x. Y is the complex FFT of the response signal y. X X X X Y X ---------( H1 = H2 = H3 ) Y
11-17
Chapter 11
Frequency Analysis
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.
11-18
ni.com
Chapter 11
Frequency Analysis
No Averaging only identifies the dominant tone. RMS Averaging does not reduce the noise floor. However, RMS averaging does smooth the noise out enough to unmask the tone at 15 kHz. Vector Averaging, Untriggered Acquisition underestimates the energy present at 10 kHz. Also, the tone at 15 kHz is indistinguishable from the noise. Vector averaging, Triggered Acquisition with triggered acquisition accurately computes the energy of the tones, reduces the noise floor by 20 dB, and reveals the tone at 15 kHz. The 20 dB reduction in the noise floor corresponds to a factor of 10, or 100 , where 100 is the number of averages completed.
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. Peak-hold averaged measurements are computed according to the following equations: FFT spectrum power spectrum where X is the complex FFT of the stimulus signal x. X is the complex conjugate of X.
Note Dual-channel, stimulus-response measurements do not support peak-hold
MAX X X MAX ( X X )
averaging.
11-19
Chapter 11
Frequency Analysis
Weighting Mode
Linear weighting weights each individual spectrum by the same amount in the averaged spectrum. Linear weighting is most often used for analysis. Exponential weighting weights the most recent spectrum more than previous spectra. Weighting the most recent spectrum more than previous spectra makes the averaged spectrum more responsive to changes in the input signal. This responsiveness makes exponential weighting ideal for the configuration phase of a measurement. Exponential weighting also is useful for monitoring applications because the averaged spectrum responds to a singular event. A linearly averaged spectrum might not respond noticeably to a singular event, especially with a large number of averages. Weighting is applied according to the following equation: N1 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.
Available Measurements
The following Frequency Analysis palettes offer the same basic measurements but designed for specific measurement needs: Baseband FFT Zoom FFT Baseband Subset
For example, each of the three palettes contains a VI for measuring the power spectrum. The SVFA Power Spectrum VI computes the power spectrum of the input signal. The SVFA Power Spectrum Subset VI computes a subset of the power spectrum of the input signal. The SVFA
11-20
ni.com
Chapter 11
Frequency Analysis
Zoom Power Spectrum VI 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, the real and imaginary parts, or the complex form 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.
11-21
Chapter 11
Frequency Analysis
Figure 11-11. Continuous Acquisition 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.
11-22
ni.com
Chapter 11
Frequency Analysis
Use the Power Spectrum Express VI to interactively configure power spectrum measurements. Refer to Chapter 3, S&V Express Measurements, and to the LabVIEW Help for more information about the Power Spectrum Express 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. Use the Frequency Analysis VIs to make the following dual-channel measurements: Frequency response computes the transfer function of the response to the stimulus. Use the frequency response as a general tool to characterize the dynamic response of a system. The coherence often is used to validate the frequency response results. The coherence quantifies the portion of the response that is linearly dependent on, or coherent with, the stimulus. The coherence ranges from 0 to 1. Cross spectrum computes the cross spectrum of the two inputs. Use the cross spectrum in some advanced analyses. Dynamic data often is stored in terms of cross spectra. Use the cross spectra to compute other useful measurements, such as frequency-response functions. Coherent output power computes the portion of the response power that is coherent with the stimulus.
You can view most dual-channel measurement results as magnitude and phase or as real and imaginary parts. Use the view parameter of the Frequency Analysis VIs returning magnitude and phase to specify the following viewing options: Whether the magnitude is expressed in decibels. Convert the magnitudes to a logarithmic scale where: magnitude decibels = 20 log ------------------------ reference Converting the magnitudes to a logarithmic scale helps to visualize small and large spectrum components at the same time.
11-23
Chapter 11
Frequency Analysis
Whether the phase is unwrapped. Unwrapping the phase removes discontinuities greater than (180). The phase is wrapped between and by default. Whether the phase is returned in degrees or radians
Refer to the LabVIEW Help for more information about the view parameter.
11-24
ni.com
Chapter 11
Frequency Analysis
PFI0
AI0
AI1 DUT
AO0
In
Out
AO1
Figure 11-12. 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.
11-25
Chapter 11
Frequency Analysis
Figure 11-13 shows the block diagram of the VI used to perform the FRF measurement.
The AO physical channel control is set to Dev1/ao0. The AI physical channels control is set to Dev1/ai0:1. The sampling frequency (fs) is 51,200 Hz and is the same for the output and input channels. This sampling frequency means the measurement is performed in the audio domain. The buffers are automatically configured by NI-DAQmx. The While Loop in Figure 11-13 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 11-12. Analog input channel 0, AI0 in Figure 11-12, acquires the same number of samples as the buffer that is generated. Simultaneously, analog input channel 1, AI1 Figure 11-12, acquires the response signal from the DUT. Figure 11-14 shows the measured time-domain stimulus and response signals.
11-26
ni.com
Chapter 11
Frequency Analysis
It is difficult to interpret the time-domain signals in Figure 11-14. Looking at the results in the frequency domain provides more insight into the DUT dynamic response.
signals from the DUT to compute the FRF. In the following examples, only the magnitude of the frequency response function is displayed.
11-27
Chapter 11
Frequency Analysis
11-28
ni.com
Chapter 11
Frequency Analysis
Using averaging improves the results as compared to the results obtained with no averaging. Refer to the Averaging Parameters section for more information about averaging and frequency analysis.
Increasing the frequency resolution results in a more accurate measurement of the attenuation of the filter around its notch frequency.
11-29
Chapter 11
Frequency Analysis
Test Type Single Input, Single Output (SISO) Multiple Input, Multiple Output (MIMO)
Common Use Cases Validating single-channel devices such as audio filters, amplifiers, equalizers, speakers, and cell phones Performing structural tests, or any type of test with multiple excitation locations and multiple response locations, such as a ground vibration test (GVT) for airplanes, bridge response, and modal analysis Validating multiple single-channel devices Performing structural tests with a single excitation location and multiple response locations Performing structural tests with multiple excitation locations and a single response location
NxM
N pairs Nx1
SISO with n number of DUTs Single Input, Multiple Output (SIMO) Multiple Input, Single Output (MISO)
1xM
Cross Spectrum
A cross spectrum measurement is an important building block for other measurements. Typically, you do not use the cross power spectrum as a direct measurement. The magnitude of the cross spectrum is equivalent to the power spectrum when the two input signals to the cross spectrum are the same signal, and for the general case, the magnitude of the cross spectrum is the product of the RMS amplitudes of the input signals, X and Y. The phase of the cross spectrum is equal to the phase of the frequency response.
Coherence
The coherence function indicates the portion of the response energy correlated to the stimulus energy. You can use the coherence function to identify excess noise and identify which of the multiple stimulus inputs is contributing to the response signal. For more information about calculating coherence, including the equation, refer to the RMS Averaging section. The coherence equation always yields a value for coherence between zero and
11-30
ni.com
Chapter 11
Frequency Analysis
one. A value of one at a given frequency indicates that all of the response energy is due to the stimulus signal and that no interference is occurring at that frequency. A value of zero indicates that no correlation exists between the response signal and the stimulus signal at that frequency. Use coherence with averaged measurements. The coherence is unity at all frequencies for applications with only one average.
Refer to Chapter 5, Limit Testing Analysis, for information about limit testing. Refer to Chapter 6, Weighting Filters, for information about weighting. Refer to Chapter 7, Integration, for information about frequency-domain integration. Refer to the LabVIEW Help for information about the individual Spectrum Extended Measurements VIs.
11-31
Chapter 11
Frequency Analysis
Power in Band
The SVFA Power in Band VI measures the total power within the specified frequency range. Table 11-5 shows the equations for computing power in band based on the type of input spectrum.
Table 11-5. Power in Band Equations Based on Input Type
PS
------------------------------ENBW Power Spectral Density (PSD) stop frequency PSD start frequency df
( MS )
---------------------------------------ENBW
where ENBW is the equivalent noise bandwidth. df is the frequency resolution of the measurement.
11-32
ni.com
Chapter 11
Frequency Analysis
Use the Power in Band Express VI to interactively configure power in band computations. Refer to Chapter 3, S&V Express Measurements, and to the LabVIEW Help for more information about the Power in Band Express VI.
Unit Conversion
The SVL 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, switch between peak units of RMS and peak-to-peak, and switch between spectrum and spectral density.
Note You cannot use the SVL 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.
11-33
Chapter 11
Frequency Analysis
When using the SVL 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 11-6.
Table 11-6. Desired Scaling Parameters
desired scaling linear/decibels dB reference (0: no change) power/magnitude view dB 10 power peak
11-34
ni.com
Transient Analysis
12
DSA Device DAQ Device WAV File DAT Recorder Simulation
This chapter discusses performing transient analysis with the Transient Analysis VIs. Refer to the LabVIEW Help for information about individual Transient Analysis VIs. Figure 12-1 illustrates the relationship of transient analysis to the sound and vibration measurement process.
Data Source
Calibration
Calibrate Sensor
Scaling
Waveform Conditioning
Integration
Weighting Filter
Analysis
Transient Analysis
Visualization
Waterfall Display
XY Graph
Figure 12-1. Relationship of Transient Analysis to the Sound and Vibration Measurement Process
12-1
Chapter 12
Transient Analysis
Figure 12-2 shows the signal corresponding to the first 200 ms of the waveform.
12-2
ni.com
Chapter 12
Transient Analysis
Figure 12-3 shows the result of applying a baseband FFT on the entire waveform.
Figure 12-3. Baseband FFT on a Chirp Signal Note No window is applied on the signal.
12-3
Chapter 12
Transient Analysis
The spectrum is flat from 10 Hz to 10 kHz. Only noise is measured at frequencies above 10 kHz. Unfortunately, this measurement does not provide any information about how the frequency content of the signal changes with time. However, the STFT can reveal useful information about the time dependence of the frequency content. 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. Subdivide the time-domain signal by configuring the time segment control displayed in Figure 12-4.
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. Create a 100 ms time increment by setting time increment to 100.00 and time increment units (%) to ms. The 100 ms time increment causes the SVT STFT versus Time VI to compute one FFT every 100 ms. Setting time increment is independent from selecting the FFT block size.
12-4
ni.com
Chapter 12
Transient Analysis
The measurement is performed using the following settings: Acquisition time = 10 s Sampling frequency = 51.2 kS/s FFT block size = 1,024 samples, or 512 lines (400 alias-free lines) Time increment= 100 ms
Based on the sampling frequency of 51,200 Hz, a 1,024 sample FFT requires a 20 ms block of data, leading to a frequency resolution of 50 Hz. 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 12-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 12-6. However, the increased resolution comes with the expense of extra processing.
12-5
Chapter 12
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 12-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.
12-6
ni.com
Chapter 12
Transient Analysis
Figure 12-7 and Figure 12-8 illustrate the overlapping process. Figure 12-7 shows a 50% overlap.
Figure 12-8 shows the resulting subdivisions when you use a 50% overlap and a Hamming window.
12-7
Chapter 12
Transient Analysis
This example uses the instance of the SVT STFT versus Time VI that produces a two-dimensional array that can be displayed directly on an SV Intensity Graph. Use the Colormap or Waterfall instance if you want to display results directly on the Colormap or Waterfall Graph, respectively. The example in Figure 12-9 acquires 10 s of data at a sample rate of 51.2 kHz. After scaling, the signal is sent to the SVT STFT versus Time VI. The result is displayed on the SV Intensity Graph in Figure 12-10.
12-8
ni.com
Chapter 12
Transient Analysis
Figure 12-10. STFT versus Time Graph Note Use the X Scale and Y Scale offset and multiplier properties to properly scale the axes of the SV 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.
12-9
Chapter 12
Transient Analysis
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 Figure 12-11 illustrates the response of a SDOF system to a half-sine pulse with a 10 g acceleration amplitude and 10 ms duration. The top graph shows the time-domain acceleration. The middle graph is the SDOF system response with a 50 Hz resonance frequency. The bottom graph is the SDOF system response with a 150 Hz resonance frequency. In both cases, is 0.05.
12-10
ni.com
Chapter 12
Transient Analysis
Use the signals shown in Figure 12-11 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 12-12 illustrates the maximax SRS for the same half-sine pulse.
Note Each computed SRS is specific to the pulse used to perform the measurement.
12-11
Chapter 12
Transient Analysis
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.
12-12
ni.com
Chapter 12
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 12-13 shows how to use the SVT Shock Response Spectrum VI. The example in Figure 12-13 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.
12-13
Chapter 12
Transient Analysis
Figure 12-14 displays the acquired time-domain signal and the computed SRS.
12-14
ni.com
Swept-Sine Measurements
13
This chapter discusses using the Swept Sine VIs. Refer to the LabVIEW Help for information about the individual Swept-Sine VIs. Use the Swept Sine Express VI on the Swept Sine (DAQmx) palette to develop and interactively configure a complete swept sine measurement. The swept-sine measurements include dynamic measurements for stimulus level, response level, frequency response (gain and phase), THD, and individual harmonic distortion. Figure 13-1 illustrates the application flow used when performing measurements with swept sine.
13-1
Chapter 13
Swept-Sine Measurements
Initialize
Start
Process Data
Graph Measurements
Close
13-2
ni.com
Chapter 13
Swept-Sine Measurements
13-3
Chapter 13
Swept-Sine Measurements
Swept-sine frequency-response measurements compare a response signal to the stimulus tone in order to compute the FRF of the DUT. The magnitude of the FRF is equivalent to gain and represents the ratio of the output level to the input level for each test frequency. The phase of the FRF is equivalent to the phase lag introduced by the DUT for each test frequency. Swept-sine measurements require a signal source. The stimulus signal is always a single tone that excites the DUT at the test frequency. Since the stimulus is a single tone, swept-sine analysis can measure the harmonic distortion while simultaneously measuring the linear response.
Swept-Sine Frequency Response Single-tone excitation Can measure harmonics Arbitrary test frequencies Longer test time for many test frequencies Better dynamic range possibility
FFT-based Frequency Response Broadband excitation Cannot measure harmonics Linearly spaced frequency resolution
Swept-sine measurements offer superior dynamic range over FFT-based measurements because you can optimize the signal level and input ranges at each test frequency. FFT-based techniques must specify a signal level and input ranges appropriate for the maximum broadband response.
13-4
ni.com
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.
FFT-based measurements are limited to a linear-spaced frequency resolution determined by the sample rate and the block size. Refer to Chapter 11, Frequency Analysis, for more information on FFT-based measurements. When the response changes rapidly, this frequency resolution may not yield enough information about the dynamic response. Also, a linear resolution may yield an excessive amount of information in frequency regions where the dynamic response is relatively constant. Swept-sine analysis has the ability to test arbitrary frequency resolutions that are linear, logarithmic, or adapted to the dynamic response of the DUT. When the frequency resolution is adapted to the DUT dynamic response, you can test more frequencies in regions where the dynamic response is of interest to your application and less where it is not. The main benefit of swept-sine analysis is the ability to measure harmonic distortion simultaneously with linear response. FFT-based analysis offers a speed advantage for broadband measurements with many test frequencies.
13-5
Chapter 13
Swept-Sine Measurements
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.
13-6
ni.com
Chapter 13
Swept-Sine Measurements
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 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.
13-7
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. The DUT for this example is a notch filter centered at 1 kHz. Figure 13-7 shows the block diagram of the example SVXMPL_Swept Sine FRF DAQmx VI, which ships with the Sound and Vibration Toolkit.
13-8
ni.com
Chapter 13
Swept-Sine Measurements
13-9
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.
Table 13-2. Swept Sine Measurement Steps
Step Number 1 2 3
Description Initialize a swept-sine measurement by specifying the hardware device and channel settings. Specify the scaling that will be applied to the acquired stimulus and response data. Configure the source by specifying the test frequencies, amplitude, and whether or not the sweep automatically restarts after completion. Set the settling and integration parameters to allow sufficient time for the DUT to settle before the measurement is performed at the new test frequency and sufficient integration time to achieve the desired level of accuracy. Set the block duration input terminal for the measurement to be small enough to give a reasonable test time and large enough so it does not put the test computer at risk of being unable to continuously generate and read the signals. The smaller the block size, the faster the swept sine can transition from one test frequency to the next. Explicitly set the sample rate for the measurement. The rate is automatically selected if this VI is not used. The same rate is used for input and output channels. Specify the propagation time terminal input specific to the DAQ device being used for the measurement. You can measure the device propagation time using the SVL Measure Propagation Delay VI. Refer to Chapter 4, Scaling and Calibration, for more information.
Required
Optional
Optional
Optional
13-10
ni.com
Chapter 13
Swept-Sine Measurements
Step Number 8
Description Configure the harmonic distortion measurement by specifying the maximum harmonic to use in the computation of the THD. Only those harmonics specified in the harmonics to visualize array return individual harmonic components. Start the swept sine to perform the hardware configuration and start the output and input tasks. Channel synchronization is performed internally in this VI. Generate the excitation and acquire the stimulus and response data at each test frequency. Convert the raw data to the specified format in order to display and report measurement results. Stop the swept-sine measurement and clear the output and input tasks to release the device.
Required
10 11 12
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. It is important to allow for the propagation delay of the DAQ or DSA device. This delay is specific to the device used to perform the measurement. To determine the device propagation delay, refer to the device documentation or measure the delay with the SVL Measure Propagation Delay VI. Refer to Chapter 4, Scaling and Calibration, for more information about the SVL Measure Propagation Delay VI. Figures 13-8 to 13-11 display measurement results obtained with the SVXMPL_Swept Sine FRF DAQmx VI example program. Figure 13-8 shows the time-domain stimulus and response signals for the 138.49 Hz
13-11
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.
13-12
ni.com
Chapter 13
Swept-Sine Measurements
In addition to measuring the frequency response, this example simultaneously measures the harmonic distortion at each test frequency. Figure 13-10 shows the graph of THD versus frequency.
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.
13-13
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.
13-14
ni.com
Distortion Measurements
14
This chapter discusses using the Distortion VIs. Refer to the LabVIEW Help for information about the individual Distortion VIs. Use the Tone Measurements Express VI to develop and interactively configure distortion and channel evaluation measurements.
Variable Definitions
In the equations in this chapter, the variables are defined as: F is the fundamental tone root-mean-square (RMS) level K is the maximum harmonic in the measurement bandwidth Hi is the ith harmonic RMS level, where i is an integer <2..K> N is the noise RMS level
P P
Harmonics
Harmonics
( H 2 ) 2 + ( H 3 ) 2 + ... + ( H K ) 2 ( H 2 ) 2 + ( H 3 ) 2 + ... + ( H K ) 2 + ( N )
2
P
=
Noise
Total
( F ) 2 + ( H 2 ) 2 + ... + ( H K ) 2 + ( N ) 2
14-1
Chapter 14
Distortion Measurements
Figure 14-1 shows where distortion analysis fits into the sound and vibration measurement process.
DSA Device DAQ Device WAV File DAT Recorder Simulation
Data Source
Calibration
Calibrate Sensor
Scaling
Waveform Conditioning
Integration
Weighting Filter
Analysis
Limit Testing
Visualization
Waveform Chart
XY Graph
Waveform Graph
Figure 14-1. Relationship of Distortion Analysis to the Sound and Vibration Measurement Process
14-2
ni.com
Chapter 14
Distortion Measurements
The SVT SINAD VI identifies the fundamental tone and effectively removes the fundamental tone from the input signal to measure the power in the harmonics and noise. 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 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.
14-3
Chapter 14
Distortion Measurements
The VI in Figure 14-2 is a modified version of the SVXMPL_One Shot Gain and Distortion (Traditional DAQ) audio example VI in the Sound and Vibration Toolkit. Figure 14-2 uses the same single-tone excitation signal and finite acquisition used in the example VI. However, the SVT SINAD VI replaces SVT Gain and SVT THD and Harmonic Components VIs the used in the example. The added frequency range control allows you to control the bandwidth of the SINAD measurement. 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.
14-4
ni.com
Chapter 14
Distortion Measurements
The reported SINAD for the left channel of the audio equalizer is SINAD=73.4 dB, 20 Hz20 kHz, test amplitude 2 Vpk, test frequency 400 Hz, and unity gain.
(14-1)
Alternatively, other sources define THD+N with the following equation: P Harmonics + P Noise THD+N = ---------------------------------------------------------F
(14-2)
14-5
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. The Sound and Vibration Toolkit THD+N function identifies the fundamental tone and effectively removes the fundamental tone from the input signal to measure the power in the harmonics and noise. 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. 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 8(df). If there are frequency components of interest below 8(df) that are not DC, reduce df by increasing the measurement duration.
14-6
ni.com
Chapter 14
Distortion Measurements
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. In this way you can perform several audio measurements sequentially to take advantage of the polymorphic Distortion VIs and Single Tone Measurements VIs. These VIs accept a waveform or a complex spectrum, and their polymorphic features allow multiple measurement analysis VIs to act sequentially on the data without recomputing any additional FFTs.
14-7
Chapter 14
Distortion Measurements
Figure 14-5 displays the distortion test measured magnitude spectrum. The fundamental frequency of 1,000 Hz is clearly visible followed by harmonics in multiples of 1,000. Refer to Table 14-1 for more information.
14-8
ni.com
Chapter 14
Distortion Measurements
Result 1,000.00 Hz 1.82 Vpk 0.80 dB 71.69 dB 76.16 dB 0.000021 Vpk 0.000120 Vpk 0.000040 Vpk 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.
14-9
Chapter 14
Distortion Measurements
You can approximate the SNR for signals dominated by a fundamental tone. For these signals the SNR is high and approximately the inverse of the difference between the THD+N and the THD.
2 2 1 -2
Use the SINAD result returned by the SVT SINAD VI if you require that the fundamental tone harmonics be included in the SNR computation. This measurement is nearly equivalent to the dynamic range measurement. The differences are that the dynamic range is a property of the channel and usually measured with a fundamental that is 60 dB FS, and the SNR is a property of the signal usually measured with a fundamental that is 0 dB FS. Some specifications call for measuring SNR in a two-step process whereby you measure the full-scale level with the signal applied, and then measure the noise with no signal present. Measure the full-scale level and then use the SVT Idle Channel Noise VI to compute the SNR with no signal present.
14-10
ni.com
Chapter 14
Distortion Measurements
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 )
14-11
Chapter 14
Distortion Measurements
Figure 14-7 shows a typical IMD test spectrum. The intermodulation components are found at multiples of 1,000 Hz, the difference frequency between the twin 11 kHz and 12 kHz tones.
IMD measurements are defined by several standards. Each of these standards have particular advantages and disadvantages when trying to reveal distortion. The standards, associated test signals, and typical use cases are listed in Table 14-2.
Table 14-2. IMD Standards and Typical Applications
Test Signal (f1, f2, ratio) (60, 7000, 4) (250, 8000, 4) others
Typical Applications Excite low-frequency distortion mechanisms (such as thermal distortion in power amplifiers) Disk recording and film ADC and DAC Slope-induced distortion
ITU-R (CCIF)
14-12
ni.com
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. You also might use a square wave to measure phase linearity. Square waves are generated easily and have a nominal crest factor of one. 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 1 Aim to use a fundamental frequency, f, such that f -------------------------------------- . 10 group delay The phase linearity is defined as the maximum deviation from a straight line of the plot of the phase versus frequency. This measurement is performed in steps. First, the multitone signal is generated, passed through the DUT, and acquired on the input channel. The phase of each tone in the multitone signal is measured and plotted versus frequency to yield a graph similar to the one in Figure 14-8.
14-13
Chapter 14
Distortion Measurements
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.
14-14
ni.com
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.
14-15
Single-Tone Measurements
15
This chapter discusses using the Single Tone Measurement VIs. Refer to the LabVIEW Help for information about the individual Single Tone Measurement VIs. Use the Tone Measurements and Noise Measurements Express VIs to develop and interactively configure single tone and noise measurements.
15-1
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.
15-2
ni.com
Chapter 15
Single-Tone Measurements
Examples of single-channel measurements are measurements of gain, idle-channel noise, dynamic range, and spurious-free dynamic range. Single-channel measurements do not require that the excitation signal be connected to an input channel of the acquisition device. Figure 15-2 illustrates the connection scheme for single-channel measurements.
PFI0
AI0 DUT
AI1 In Out
AO0
AO1
15-3
Chapter 15
Single-Tone Measurements
Figure 15-3 illustrates where single-tone measurements fit into the sound and vibration measurement process.
Data Source
Calibration
Calibrate Sensor
Scaling
Waveform Conditioning
Integration
Weighting Filter
Analysis
Limit Testing
Visualization
Waveform Chart
XY Graph
Waveform Graph
Figure 15-3. Relationship of Single-Tone Measurements to the Sound and Vibration Measurement Process
15-4
ni.com
Chapter 15
Single-Tone Measurements
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.
15-5
Chapter 15
Single-Tone Measurements
The measurement analysis is equivalent to that performed for a gain and phase measurement. The difference between these two measurements is that the gain and phase measurement is applied on the active response channel, while the crosstalk measurement is applied to the idle response channel. Figure 15-4 illustrates this difference.
Figure 15-4. Difference between Gain and Phase Measurement and Crosstalk Measurement
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.
15-6
ni.com
Chapter 15
Single-Tone Measurements
Amplitude Connection
Figure 15-5. Connecting the Amplitude of the Excitation Tone to the Amplitude Input
The accuracy of this measurement is dependent on the absolute amplitude accuracies of the output channel and the input channel. When generating and acquiring the signal with an NI PXI-4461, the expected measurement uncertainty is less than 0.3 dB.
Idle-Channel Noise
The SVT Idle Channel Noise VI measures the total power, the median noise density level, the maximum spur, and the SNR within the specified frequency range. This VI can operate on waveform data or the complex spectrum output by the other audio measurements. For practical purposes, the returned idle-channel noise is a power-in-band measurement that returns the aggregate noise energy in the specified frequency range. The SVT Idle Channel Noise VI computes the SNR as the full-scale level divided by the idle-channel noise. The median noise level is normalized by the bin width, df, and is one measure of the noise floor, or noise density level. You might find noise floor knowledge useful in advanced measurements. For example, you can specify the threshold input to the SVFA Spectrum Peak Search VI in order to identify all the spectral components that exceed the noise floor plus 12 dB, or four times the noise floor.
15-7
Chapter 15
Single-Tone Measurements
Unlike gain, phase, and crosstalk measurements, idle-channel noise is a broadband measurement in which the amplitude of the fundamental frequency is expected to be zero. Use the frequency range control to limit the bandwidth of the measurement. To exclude DC power, set the start frequency to 8(df), where df is the frequency resolution of the spectrum. Increase the start frequency in the frequency range control to exclude additional low-frequency components from the measurement. You can enter the default values of -1 for the start frequency and 1 for the stop frequency to perform the measurement over the entire acquisition bandwidth.
Note You can measure the idle-channel noise with the input channel to the DUT floating,
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. The connection scheme for a dynamic range measurement follows that of the single-channel measurement in Figure 15-2. Configure the measurement by specifying the measurement settings in the SVXMPL_Settled Dynamic Range (DAQmx) VI. It is important to allow for the propagation delay of the DAQ or DSA device. This delay is specific to the device used to perform the measurement. To determine the device propagation delay, refer to the device documentation or measure the delay with the SVL Measure Propagation Delay VI. Refer to Chapter 4, Scaling and Calibration, for more information about the SVL Measure Propagation Delay VI. The SVXMPL_Settled Dynamic Range (DAQmx) VI performs a continuous generation and acquisition. The measurement continues until the measured dynamic range is settled according to the tolerances specified under the Settling Parameters tab. For the left channel of this DUT, a two-channel audio equalizer, a typical measurement configuration results in a dynamic range of 83 dB FS, 20 Hz20 kHz, test amplitude 60 dB FS, test frequency 1.0 kHz, and unity gain.
15-8
ni.com
Chapter 15
Single-Tone Measurements
Figure 15-6 shows the final measured spectrum for the dynamic range measurement.
15-9
Chapter 15
Single-Tone Measurements
15-10
ni.com
References
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.
Note The following resources offer useful background information on the general
concepts discussed in this documentation. These resources are provided for general informational purposes only and are not affiliated, sponsored, or endorsed by National Instruments. The content of these resources is not a representation of, may not correspond to, and does not imply current or future functionality in the Sound and Vibration Toolkit or any other National Instruments product. American National Standards Institute. 1983. Specification for Sound Level Meters. ANSI S1.4-1983. Washington: American National Standards Institute. . 1986a. Design Response of Weighting Networks for Acoustical Measurements. ANSI S1.42-2001. Washington: American National Standards Institute. . 2004. Specification for Octave-Band and Fractional Octave-Band Analog and Digital Filters. ANSI S1.11-2004. New York: Acoustical Society of America. Crocker, Malcolm J. 1998. Handbook of Acoustics. New York: John Wiley & Sons, Inc. Hassall, J. R., and K. Zaveri. 1988. Acoustic Noise Measurements. Nrum, Denmark: Brel & Kjr. International Electrotechnical Commission. 1975. Preferred Frequencies for Measurements. International Standard IEC 266, 1st ed. 1975-07-15. Geneva, Switzerland: International Electrotechnical Commission. . 2002. Electroacoustics - Sound Level Meters. International Standard IEC 61672, 1st ed, 2002-05. Geneva, Switzerland: International Electrotechnical Commission.
A-1
Appendix A
References
. 1995. Octave-Band and Fractional Octave-Band Filters. International Standard IEC 1260, 1st ed. 1995-07. Geneva, Switzerland: International Electrotechnical Commission. Randall, R.B. 1987. Frequency Analysis. Nrum, Denmark: Brel & Kjr. . 1986. Measurement of Audio-Frequency Noise Voltage Level in Sound Broadcasting. ITU-R Recommendation 468-4. . 1995. Psophometer for Use on Telephone-Type Circuits. ITU-T Recommendation O.41, Revised, 1993-1996. Telecommunication Standardization Sector of the International Telecommunication Union. Smallwood, D. 1980. An Improved Recursive Formula for Calculating Shock Response Spectra. 51st Shock and Vibration Bulletin.
A-2
ni.com
Visit the following sections of the National Instruments Web site at ni.com for technical support and professional services: SupportOnline technical support resources at ni.com/support include the following: Self-Help ResourcesFor answers and solutions, visit the award-winning National Instruments Web site for software drivers and updates, a searchable KnowledgeBase, product manuals, step-by-step troubleshooting wizards, thousands of example programs, tutorials, application notes, instrument drivers, and so on. Free Technical SupportAll registered users receive free Basic Service, which includes access to hundreds of Application Engineers worldwide in the NI Discussion Forums at ni.com/forums. National Instruments Application Engineers make sure every question receives an answer. For information about other technical support options in your area, visit ni.com/services or contact your local office at ni.com/contact. Training and CertificationVisit ni.com/training for self-paced training, eLearning virtual classrooms, interactive CDs, and Certification program information. You also can register for instructor-led, hands-on courses at locations around the world. System IntegrationIf you have time constraints, limited in-house technical resources, or other project challenges, National Instruments Alliance Partner members can help. To learn more, call your local NI office or visit ni.com/alliance.
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.
B-1
Glossary
Numbers/Symbols
2D 3D Pa Two-dimensional. Three-dimensional. Micro-Pascal. See also Pa.
A
acceleration accelerometer aliasing The rate of change of velocity. A sensor mounted on a structure to measure the acceleration at a particular location in one or multiple directions. Aliasing is a phenomenon whereby an analog signal of frequency greater than half the sampling frequency (Nyquist frequency) appears, after sampling, at a frequency less than half the sampling frequency. See also anti-aliasing filter. American National Standards Institute. 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 bandpass filter baseband analysis The upper and lower cutoff frequencies of an ideal bandpass filter. A filter with a single transmission band extending from a lower bandedge frequency greater than zero to a finite upper bandedge frequency. The range of frequencies extending from 0 Hz (DC) to a maximum frequency, often the Nyquist frequency.
G-1
Glossary
C
calibrator coherence A controlled source generating a known level of excitation used to calibrate a sensor. Gives a measure of the degree of linear dependence between two signals, as a function of frequency. The coherent output power spectrum gives a measure of what part of the (output) power spectrum is fully coherent with the input signal. A method of displaying 3-dimensional data on a 2-dimensional plot using color to indicate the value in the third dimension. The type of filters used in fractional-octave measurements. Coherent output power. 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. The cross power spectrum of two signals has an amplitude that is the product of the two amplitudes, and a phase that is the difference of the two phases. The frequency at which the filter attenuates the input 3 dB, or half of its original power.
D
DAQ dB Data acquisition. 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. Discrete Fourier Transformdetermines the amplitude and phase of frequency components present in a time-domain digital signal. The production of signal components not in the original signal due to non-linearities in the system or transmission path.
DFT Distortion
G-2
ni.com
Glossary
Dynamic signal acquisition. Device under test. 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. Discrete Zak Transform.
DZT
E
ENBW engineering unit (EU) equal confidence Equivalent noise bandwidth of the time-domain window applied to the signal in order to perform frequency analysis. Term of data measurement, such as Pa, m/s2, g, and so on.
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. The energy average level of a signal over a given time interval. See engineering unit (EU). Time-averaging technique that gives recent data more importance than older data.
F
Fast FFT Exponential averaging using a time constant of 125 ms. Fast Fourier Transforman 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. The number of samples used to compute an FFT.
G-3
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. A group of filters. A type of signal conditioning that allows you to modify the frequency content of a signal. The interval between two frequencies, one of which is a fractional power of two of the other. Represents the ratio of output-to-input in the frequency domain, and fully characterizes linear, time-invariant systems. See frequency response function. The frequency of the dominant (highest amplitude) tone in the signal. For harmonic analysis the fundamental frequency is the greatest common divisor of the harmonic frequencies.
G
g gain Unit for measuring acceleration. One g = 9.81 m/s2, the acceleration due to gravity at the surface of the Earth. The amplification or attenuation of a signal.
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.
G-4
ni.com
Glossary
H2
The frequency response function computed as the ratio of the output autospectrum to the backwards cross spectrum: Gyy/Gyx. This technique gives the best performance in the presence of noise for measuring resonances, where the signal to noise ratio tends to be best. For measurement of anti-resonances, the frequency response function H1, 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. The frequency response function computed as an average of H1 and H2. A signal whose frequencies are integer multiples of the input signal. A form of distortion in analog circuits that generates harmonics. It is calculated as the ratio of the level of a single harmonic to the level of the original signal. Hertz. Cycles per second.
Hz
I
IEC IEPE IMD Impulse International Electrotechnical Commission. Integrated Electronic Piezoelectric. See intermodulation distortion. Exponential averaging using a time constant of 35 ms if the signal level is rising and 1,500 ms if the signal level is falling. 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. The distortion that arises as the result of the modulation between two or more signals.
L
Leq linear averaging See equivalent continuous level (Leq). Time-averaging technique that weights all data in the average equally.
G-5
Glossary
M
MAX maximax Measurement & Automation Explorer. The absolute maximum of the calculated shock response signal over the entire signal duration. Sensor used to convert sound pressure variations into an electrical signal, usually when the acoustic medium is air. The center frequency of a bandpass filter, defined as the geometric mean of the bandedge frequencies. Maximum Length SequenceA pseudo-random sequence of ones and zeroes you can use as a test signal to perform delay and frequency response measurements. Millisecond.
MLS
ms
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. Rounded midband frequency for the designation of a particular fractional-octave filter. This term is used by the IEC standards, but nominal frequencies are identical to the preferred frequencies defined in the ANSI standards. Signal whose frequency content changes within a captured frame. 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.
nominal frequency
G-6
ni.com
Glossary
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. Harmonic of the rotational speed of rotating machinery. As an example, in the case of the shaft rotating at 6000 rpm, the first order component occurs at a frequency of 100 Hz (6000/60), whereas the third order component would occur at a frequency of 300 Hz. A method that uses a portion of the previous data block to compute the FFT of the current data block.
order
overlapping
P
Pa peak hold periodicity phon Pascal. International unit of pressure. Peak detection process retaining the maximum value of a signal. One of the basic assumptions made in FFT-based frequency analysis. 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. 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. Microphone calibrator generating a known sound pressure level, typically at a certain reference frequency. Rounded midband frequency for the designation of a particular fractional-octave filter. This term is used by the ANSI standards, but preferred frequencies are identical to the nominal frequencies defined in the IEC standards. Any gain applied to a signal by an external device (amplifier, preamplifier, signal conditioning, and so on) before the data acquisition device. Power spectral density.
pink noise
pregain PSD
G-7
Glossary
R
reference sound pressure reverberation time A reference pressure of 20E-6 Pa. This reference pressure was conventionally chosen to correspond to the quietest sound at 1,000 Hz that the human ear can detect. 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. Root Mean Square. RMS averaging is used to average the power of a signal. RMS averaging reduces fluctuations. Revolutions per minute.
rpm
S
s S/s sampling frequency SDOF sensor Seconds. Samples/second. The rate at which a continuous waveform is digitized. Single degree of freedom. A device that converts a physical stimulus (such as force, sound, pressure, motion) into a corresponding electrical signal. The amount of time required for a signal to reach its final value within specified limits. See spurious free dynamic range. States that to properly sample a signal, the signal must not contain frequencies above the Nyquist frequency.
shock response spectrum A processing method which evaluates the severity of a shock signal. signal in noise and distortion SINAD The ratio of the input signal to the sum of noise and harmonics. See signal in noise and distortion.
G-8
ni.com
Glossary
An acquisition that runs only one time. Exponential averaging using a time constant of 1,000 ms. 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. 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. 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. The dynamic range from full-scale deflection to the highest spurious signal in the frequency domain. See shock response spectrum. Short-Time Fourier Transform. Sound and Vibration Frequency Analysis. SVFA is a collection of VIs shared by the Sound and Vibration Toolkit and other NI products. Sound and Vibration Library. SVL is a collection of VIs shared by the Sound and Vibration Toolkit and other NI toolkits. Sound and Vibration Toolkit.
SVL SVT
T
tach tachometer THD THD+N See tachometer. Device used to measure the rotational speed of a rotating part. See total harmonic distortion. See total harmonic distortion plus noise.
G-9
Glossary
Ratio between two frequencies, equal to 21/3. 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. The ratio of the sum of harmonics to the fundamental tone.
total harmonic distortion total harmonic distortion plus noise transducer transient
The ratio, in decibels, of the sum of noise and harmonics to the input signal. See sensor. A very short-duration signal. Normally only occurs once, or very infrequently. Transition-transistor logica typical medium speed digital technology. Nominal TTL logic levels are 0 and 5 V.
TTL
V
V vector averaging Volts. 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. A vector quantity whose magnitude is a bodys speed and whose direction is the bodys direction of motion.
velocity
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.
G-10
ni.com
Glossary
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. 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. A smooth waveform that generally has zero value at the edges. See also windowing. Technique used to reduce spectral leakage by multiplying the time-domain waveform by a window function. The process of windowing reduces the amplitudes of discontinuities at the edges of a waveform, thereby reducing spectral leakage. If the waveform contains an integral number of cycles, there is no spectral leakage. See also spectral leakage.
white noise
window function
windowing
Z
Z-weighting Zero-weighting. Weighting with unity gain at all frequencies; equivalent to linear or no weighting.
G-11