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 handarm 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