Download All

Code covered by the BSD License

Highlights from
Continuous Sound and Vibration Analysis

choosebox(varargin) CHOOSEBOX Two-listed item selection dialog box.
selectdlg2( cellItems, guiTit... selectdlg2 Generate a scrolled matrix of choices for user selection.
tableGUI(varargin) TABLEGUI - Spreadsheet like display and edition of a generic 2D array. By generic it is
ACdsgn(Fs)
ArgStruct=parseArgs(args,ArgS... Helper function for parsing varargin.
[ ha_chan2, ha_axes2, ha_acce... % config_ha_accels: Configures the axis directions and channels for accelerometers measuring hand arm vibrations
[ wb_chan2, wb_axes2, wb_acce... % config_wb_accels: Configures the axis directions and channels for accelerometers measuring whole body vibrations
[A2, A_str, real_digitsL, rea... % m_round: Rounds an array to a specified number of significant digits or a specified digits place: significant figures, sigfigs
[A2, A_str, real_digitsL, rea... % pow10_round: Round a numeric array to a Specified Digits Place
[A2, A_str, real_digitsL, rea... % sd_round: Rounds an array to a specified number of Significant Digits, significant figures, digits of precision
[B, A]=Nth_octdsgn(Fs, Fc, N,...
[B, A]=hand_arm_fil(Fs) % hand_arm_fil: Calculates the hand arm vibrations filter coefficients
[Fsn, p, q, errors]=get_p_q2(...
[LeqA, LeqA8, LeqC, LeqC8, Le... % Leq_all_calc: Calculates levels and peaks for A, C, and linear weighting
[SP, vibs, Fs_SP, Fs_vibs, de... % data_loader2: Loads data specified as sound or vibrations and further specified as data, time increment or sampling rate
[SP2, mean_array2, mean1, arr... % sub_mean: Removes the running average from a time record given a sampling rate and high pass cutoff frequency.
[SP2, mean_array2]=sub_mean(S... % sub_mean: Removes the running average from a time record given a sampling rate and high pass cutoff frequency.
[alpha]=t_alpha(t, nu) % t_alpha: cumulative probability function of the t-distribution,
[arms_w, arms8h_w, Dy_w, Dy_5... % Vibs_calc_hand_arm: Calculates the metrics for the hand arm vibrations
[arms_wb, arms8h_wb, VDV_wb, ... % Vibs_calc_whole_body: Calculates the metrics for the whole body vibrations
[buf2, num_files_not_empty]=s... % splat_cell: Returns a row vector of numbers containng all of the numbers in a cell array
[bz, az]=bessel_digital(Fs, F... % bessel_digital: creates a digital low pass bessel filter of order n
[comb_data_buf2]=combine_acce... % combine_accel_directions_ha: Calculates overall hand arm vibrations acceleration metrics from multiple axis accelerometers.
[comb_data_buf2]=combine_acce... % combine_accel_directions_wb: Calculates overall whole body vibrations acceleration metrics from multiple axis accelerometers.
[fc_out, SP_levels, SP_peak_l... % Nth_oct_time_filter2: Calculates the Nth octave center frequencies, sound levels, peak levels, and time records
[fid, mean_vals, max_num_chan... % print_channel_stats: Print the ststistics for each channel for the impulsive sound table
[fid]=print_overall_stats(fid... % print_overall_stats: Print the overall statistics for the impulsive sound table
[filename_base, ext]=file_ext... % file_extension: separates a filename and path from the file extension
[filtercoeffs, filter_types, ... % whole_Body_Filter: Calculates the filter coefficients for whole body vibrations.
[h, h2, wb_th_array]=plot_snd... % plot_snd_vibs(: plots sound and vibrations data in the time domain
[interval, error]=t_confidenc... % t_confidence_interval: One or two sided confidence interval of standard error with t distribution
[m2]=geospace(a, b, n, flag) % geospace: caculates a geometric sequence or psuedogeometric sequence from a to b with n elements
[num_accels, accel_num_chan, ... % accel_config: Configures the accelerometer channels given the cell array axes_config
[num_accels, accel_num_chan, ... % config_accels: Configures the accelerometer channels using an input dialog box.
[num_channels_a, num_diff_cha... % table_append_channels: caculates the number of channels of data for making a table
[num_samples_ca, s]=num_impul... % num_impulsive_samples: Counts the number of impulsive noise samples
[nums, matches]=find_nums(str... % find_nums: Finds floating point complex, real, integer, and currency (dollars) numbers in a string
[out]=print_data_loader_confi... % print_data_loader_configuration_table: Prints the variable configuration to a table
[pfq]=genHyper(a,b,z,lnpfq,ix... function [pfq]=genHyper(a,b,z,lnpfq,ix,nsigfig)
[prms]=rms_val(p, dim) % rms_val: Calculates the rms value along a specific dimension
[ptsa, nptsa, rt_stats, other... % data_outliers: Determine the outliers and return descriptive statistics
[res,raw]=fastmcd(data,option... version 22/12/2000, revised 19/01/2001,
[s]=calc_diff_metrics_cont(s,... % calc_diff_metrics_cont: Calculates the differences in metrics between pairs of channels: Impulsive Noise Meter
[snd, vibras_ha, vibras_wb, v... % Continuous_Sound_and_Vibrations_Analysis: Loads time records and timeseries, Calculates metrics for sound and vibrations (hand arm and whole body).
[t_SP_rs, SP_rs]=resample_plo... % resample_plot: Resamples a plot to 10000 data points using the max and%min from each bin
[t_out, T_out, error1_out, er... % inverse cumulative probability function, t-distribution
[varargout]=convert_double(va... % This program converts the inputs into double precision arrays. Then
[y, x, a]=match_height_and_sl... % match_height_and_slopes2: creates a quartic with specifed height and slope at the end points.
[y2, num_settle_pts, settling... % filter_settling_data: Creates data to append to a time record for settling a filter
[y2]=remove_filter_settling_d... % remove_filter_settling_data: removes data added to time records to settle the filter
[yAC, errors]=ACweight_time_f... % ACweight_time_filter: Applies an A or C weighting time filter to a sound recording
[y]=moving(x,m,fun) MOVING will compute moving averages of order n (best taken as odd)
[y_out, b, a]=bessel_antialia... % bessel_antialias: applies an antialiasing digital Bessel filter
[y_out, t_out, b, a]=bessel_d... % bessel_down_sample: applies an antialiasing digital Bessel filter
[y_out, x_out, y_in]=resample... % resample_interp3: resamples using interp1 with additional options
[yh, B, A, errors]=hand_arm_t... % hand_arm_time_fil: Calculates the hand arm vibrations filter coefficients and returns the filtered time record
[yhw, filtercoeffs, filter_ty... % whole_body_time_filter: Calculates the whole body vibrations filter coefficients and returns the filtered time record
[ytick_m, YTickLabel1, ytick_... % fix_YTick: Sets appropriate Y-Tick values for small plots
channel_data_type_selection(n... % channel_data_type_selection
fastlts(x,y,options)
h=subaxis(varargin) SUBAXIS Create axes in tiled positions. (just like subplot)
kurtosis2(x, dimension) % This program calculates the kurtosis.
loc=LMSloc(X)
main_sound_and_vibs % Main_snd_and_vibs: Main Program for continuous sound and vibrations analysis
make_summary_cont_stats_table... % make_summary_cont_stats_table: Makes a table of the output of the
nth_freq_band(N, min_f, max_f... % nth_freq_band: Calculates the 1/nth octave frequency bands center, lower, and upper bandedge limits
psuedo_box(h_array)
rmean(y, db_or_lin) % get_stats: Calculates descriptive statistics for the input variable y.
save_a_plot2_audiological(a, ... % save_a_plot2_audiological: Saves current figure to specified image type.
variable_data_type_selection(... % channel_data_type_selection
View all files

5.0

5.0 | 3 ratings Rate this file 73 Downloads (last 30 days) File Size: 1.48 MB File ID: #21384

Continuous Sound and Vibration Analysis

by Edward Zechmann

09 Sep 2008 (Updated 11 Jan 2012)

This program analyzes sound and vibrations data using metrics for continuous noise and vibrations.

| Watch this File

File Information

Description

The main program is Main_sound_and_vibs.m. Up to 100 metrics are for each of sound, hand arm vibrations, whole body vibrations, and motion sickness. Main_sound_and_vibs calls the primary depedent function
"Continuous_Sound_and_Vibrations_Analysis" which contains the switch statement controlling the data flow.

This program calculates metrics for continuous sound and vibrations. Sound metrics include: peaks, Leq, LeqA, LeqC, kurtosis, third octave band peaks and levels, and more.

Vibration metrics for hand-arm include: arms, armq, Dy, peak, crest factor, kurtosis, third octave band levels and peaks, and more.

Vibration metrics for whole-body include: arms, armq, VDV, MSDV, crest factor, kurtosis, third octave band levels and peaks, and more.

The vibration metrics except for the third octave bands are calculated using both the weighted and unweighted filters.

This program prompts the user for all of the inputs to the Continuous_Sound_and_Vibrations_Analysis.m program.

The user selects either matlab or wav files to analyze for sound and vibrations. There are a series of prompts for output filenames, formating of output images of the figures, and the Name of the device under test (Tool Name).

Required Products

Signal Processing Toolbox

MATLAB release

MATLAB 7.5 (R2007b)

Tags for This File
Everyone's Tags	dsp, hand arm, leq, level, noise(2), peak, sound(2), vibrations, whole body
Tags I've Applied
Add New Tags	Please login to tag files.

Comments and Ratings (3)
09 Oct 2008	Seeker Tang	good job!
27 Oct 2009	Florentino Ruiz	Estoy empezando a investigar sobre el analisis de vibraciones en maquinas rotativas, me interesa las cosas que haces, voy a analizar tus programas y si puedes apoyarme en algo, ahi te comento, gracias
23 Aug 2011	Benjamin	Seems to do the job. Good effort.

Please login to add a comment or rating.

Updates
09 Oct 2008	Added Comments and modified the subprogram find_nums. Fixed a bug with reporting the scaling factors for whole body vibration. Fixed bug where values of the mtrics were not reported to the metrics field.
12 Nov 2008	Fixed a bug in the running average.
21 Jan 2009	Fixed a bug in loading wave files.
27 Apr 2010	Fixed a bug in calculating VDV for whole body vibrations.
18 Jan 2011	Fixed abug in the selection of the posture type. This affects posture types 2-8 including the filter coefficients.
21 Feb 2011	Fixed bugs in outputting vibrations metrics.
11 Jan 2012	Fixed a bug in the calculation of VDV. Thanks to Longqi Wang for finding the bug.

Tag Activity for this File
Tag	Applied By	Date/Time
dsp	Edward Zechmann	22 Oct 2008 10:18:31
sound	Edward Zechmann	22 Oct 2008 10:18:31
vibrations	Edward Zechmann	22 Oct 2008 10:18:31
hand arm	Edward Zechmann	22 Oct 2008 10:18:31
whole body	Edward Zechmann	22 Oct 2008 10:18:31
level	Edward Zechmann	22 Oct 2008 10:18:31
leq	Edward Zechmann	22 Oct 2008 10:18:31
peak	Edward Zechmann	22 Oct 2008 10:18:31
noise	Edward Zechmann	22 Oct 2008 10:18:31
noise	serkan	22 Feb 2011 04:36:10
sound	owais hamza	11 Jun 2011 02:37:18

Highlights from Continuous Sound and Vibration Analysis

Continuous Sound and Vibration Analysis

Highlights from
Continuous Sound and Vibration Analysis