Code covered by the BSD License
ACdsgn(Fs)
Test_Nth_oct_filters1(res... % Test_Nth_oct_filters1: This program tests the Nth octave time filters for continuous and impulsive noise
[A2, A_str, real_digitsL,... % sd_round: Rounds an array to a specified number of Significant Digits, significant figures, digits of precision
[B, A]=Nth_octdsgn(Fs, Fc...
[B, A]=hand_arm_fil(Fs)
% hand_arm_fil: Calculates the hand arm vibrations filter coefficients
[Fsn, p, q, errors]=get_p...
[LeqA, LeqA8, LeqC, LeqC8... % Leq_all_calc: Calculates levels and peaks for A, C, and linear weighting
[SP, f, bin_size, num_ave... % pressure_spectra: Calculates an accurate estimate of the pressure spectra
[SP, f, num_averages_out]... % spectra_estimate: Is a rough estimate of the pressure spectra
[SP2, mean_array2, mean1,... % sub_mean: Removes the running average from a time record given a sampling rate and high pass cutoff frequency.
[SP2, mean_array2]=sub_me... % sub_mean: Removes the running average from a time record given a sampling rate and high pass cutoff frequency.
[SP_rms_levels_a, SP_peak... % Test_third_oct_filters: Tests Nth octave filters with pure tones and tone bursts
[Wa, Wc]=AC_weight_NB(f, ... % AC_weight_NB: Calculates the A and C frequency weights at specified frequencies
[bin_size, num_averages_o... % number_of_averages: Calculates the number of points not overlapped from the array size, bin size, and number of averages
[bz, az]=bessel_digital(F... % bessel_digital: creates a digital low pass bessel filter of order n
[f, A_atten, A_atten2, C_... % Test_ACweight: Tests the A and C weighting filters using pure tones and tone bursts
[f_sig, Wf ]=Test_hand_ar... % Test_hand_arm: Tests the accuracy of the hand-arm vibrations filters
[fc_out, SP_levels, SP_pe... % Nth_oct_time_filter: Calculates the Nth octave center frequencies, sound levels, peak levels, and time records
[fc_out, SP_levels, SP_pe... % Nth_oct_time_filter2: Calculates the Nth octave center frequencies, sound levels, peak levels, and time records
[filename_base, ext]=file... % file_extension: separates a filename and path from the file extension
[ftwcf]=window_correction... % window_correction_factor: Computes the factor for calibrating a Fourier Transform given specific processing parameters
[m2]=geospace(a, b, n, fl... % geospace: caculates a geometric sequence or psuedogeometric sequence from a to b with n elements
[ndraw, count, count2, er... % rand_int: Quckly generates n random integers from a to b integer
[prms]=rms_val(p, dim)
% rms_val: Calculates the rms value along a specific dimension
[res,raw]=fastmcd(data,op... version 22/12/2000, revised 19/01/2001,
[varargout]=convert_doubl... % This program converts the inputs into double precision arrays. Then
[w]=flat_top(N, type)
% Flat top windows are used for calibration, because the wide main lobe
[xw, xp, xb, SPLw1, SPLp1... % Test_third_oct_filters: Tests Nth octave filters with white, pink, and brown noise.
[y, t]=tone_burst(Fs, fc,... % sinusoidal_w_wave: N wave with known peak level, frequency and duration
[y, x, a]=match_height_an... % match_height_and_slopes2: creates a quartic with specifed height and slope at the end points.
[y2, num_settle_pts, sett... % filter_settling_data: Creates data to append to a time record for settling a filter
[y2]=remove_filter_settli... % remove_filter_settling_data: removes data added to time records to settle the filter
[yAC, errors]=ACweight_ti... % 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_anti... % bessel_antialias: applies an antialiasing digital Bessel filter
[y_out, t_out, b, a]=bess... % bessel_down_sample: applies an antialiasing digital Bessel filter
[y_out, x_out, y_in]=resa... % resample_interp3: resamples using interp1 with additional options
[yh, B, A, errors]=hand_a... % hand_arm_time_fil: Calculates the hand arm vibrations filter coefficients and returns the filtered time record
fastlts(x,y,options) version 22/12/2000, revised 19/01/2001, 30/01/2003
nth_freq_band(N, min_f, m... % nth_freq_band: Calculates the 1/nth octave frequency bands center, lower, and upper bandedge limits
plot_test_Nth_oct_filters... % plot_test_Nth_oct_filters: plots the third octave filter test data
progressbar(varargin)
Displays a multi leveled progressbar. This makes it easy to nest
rmean(y, db_or_lin)
% get_stats: Calculates descriptive statistics for the input variable y.
save_a_plot_reverb_time(a... % save_a_plot_reverb_time: Saves current figure to specified image type.
spatialPattern(DIM,BETA)
function x = spatialPattern(DIM, BETA),
View all files
Nth_Oct_Hand_Arm_&_AC_Filter_Tool_Box
by
Edward Zechmann
22 Dec 2008
(Updated
15 Mar 2012 )
Features Nth octave band, Hand Arm, and A and C weighting filters
[varargout]=convert_double(varargin)
function [varargout]=convert_double(varargin)
% % This program converts the inputs into double precision arrays. Then
% % outputs them. An indefinite number of inputs and outputs can be used.
% %
% Example;
% a=single([1 2]); % a is single precision vector
% [b]=convert_double(a); % varargin will have one input variable 'a'
% % varargout will have one output variable 'b'
% %
% %
% %
% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% % Program Written by Edward L. Zechmann
% % date Not certain 2007
% % modified 19 December 2007
% % added comments and an example
% %
% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% % Please feel free to modify this code.
% %
for e1=1:nargin;
ttype=class(varargin{e1});
if ~isequal(ttype, 'double')
varargout{e1}=double(varargin{e1});
else
varargout{e1}=varargin{e1};
end
end
Contact us