| [SP, vibs, Fs_SP, Fs_vibs, default_wav_settings_out, default_mat_config_out]=data_loader2(filename, default_wav_settings_in, default_mat_config_in)
|
function [SP, vibs, Fs_SP, Fs_vibs, default_wav_settings_out, default_mat_config_out]=data_loader2(filename, default_wav_settings_in, default_mat_config_in)
% %
% % data_loader2.m loads data stored in a matlab file or a wave file
% %
% % If the file is a wave file it is loaded
% % Then the user is prompted to specify which tracks are sound and
% % which are vibrations.
% %
% % Then the user is prompted for the calibration sensitivies for each
% % channel listed by accelerometer channels and microphone channels.
% %
% %
% % This program finds the numeric variables stored in the
% % specified data file then prompts the user to select the data variables
% %
% % Then the user is prompted to specify whether each variable is by
% % typing a single letter for each variable in the dialog box input
% %
% % s for sound
% % v for vibrations
% % b for both
% % n for neither
% %
% % In teh case of both soudn and vibrations the user is prompted to
% % specify which channels are sound and which channels are vibrations
% %
% % Then the program prompts the user to select the time or sampling rate
% % variables. These variables are used to specify the sampling rate for
% % each data variable.
% %
% % The user is prompted to relate each data variable to the corresponding
% % time or sampling rate variable.
% %
% % If there are no time or sampling rate variables then the program
% % prompts the user to enter the sampling rates for each data variable.
% %
% %
% %
% %
% Example=''; % for matlab data
% filename='data.mat'; % filename of the matlab file to open
% default_wav_settings_in={[1 0 ], [40.1], [10.1]};
%
% % used to automatically select sensor types
% % and specify the calibration sensitivities
% % default is empty cell array.
% % default is to prompt the user with
% % dialog boxes.
% %
% default_mat_config_in={{'SP'}, {'Fs_SP'}; {'vibs'}, {'Fs_vibs'}; {}, {}; {1}, {};};
% %
% % used to automatically select sensor types
% % and specify the calibration sensitivities
% % default is empty cell array.
% % default is to prompt the user with
% % dialog boxes.
% %
% %
% % Output Variables
% %
% % SP={}; % Pa Sound pressure
% % vibs={}; % m/s^2 acceleration
% % Fs_SP={}; % Hz smpling rate of sound pressure
% % Fs_vibs={}; % Hz smpling rate of acceleration
% %
% % default_wav_settings_out={};
% %
% % % used to automatically select sensor types
% % % and specify the calibration sensitivities
% % % default is to set the values to
% % % the user entered data from the dialog
% % % boxes then use the code
% %
% % default_wav_settings_in=default_wav_settings_out;
% %
% % default_mat_config_out=cell(4, 5);
% %
% % % used to automatically match the data variables
% % % with the sensor types
% % % specify the time or sampling rate
% % % variables
% % % match the data variables with the
% % % time or sampling rate
% % % The best way to set the valeus is to
% % % run the file through the program then
% % % run the code
% %
% % default_mat_config_in=default_mat_config_out;
% %
% % List of Sub Programs
% %
% % file_extension
% % convert_double
% %
% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% %
% % Program Written by Edward L. Zechmann
% %
% % Date 18 December 2007
% % modified 21 December 2007 updated comments
% % added use default settings
% % for reading wave files
% %
% % modified 2 January 2008 added option to select the data variable
% % and time or frequency variable
% %
% % modified 4 January 2008 added more code to find and select
% % the data variables and time or sampling
% % variables.
% %
% % modified 6 January 2008 finished the variable selection code
% % and updated comments
% %
% % modified 7 January 2008 fixed variaous bugs when configuring
% % data for both sound and vibrations
% % fixed bugs for specifying sampling rates
% %
% % modified 13 January 2008
% %
% % modified 13 January 2008 added the uigetfile program to take care of
% % the problem of a lack of files to open
% %
% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %
% % Please Feel Free to Modify This Program
% %
if logical(nargin < 1) || logical(length(filename ) < 1)
[filename, pathname] = uigetfile( { '*.mat','*.wav'}, 'Select the files To Process', 'MultiSelect', 'on');
cd(pathname);
end
% Make sure that the configuration variables exist for the wav file
if nargin < 2 || isempty(default_wav_settings_in)
default_wav_settings_in=cell(4,1);
end
default_wav_settings_out=default_wav_settings_in;
% configuring the matlab variables is quite invovled and requires
% checking initializiing the configuration variables
% Make sure that the configuration variables exist
if nargin < 3 || isempty(default_mat_config_in)
default_mat_config_in=cell(4,5);
end
[cfm1 cfn1]=size(default_mat_config_in);
% Make sure input automated configuration meets the minimum size
% requirements
for e1=1:4;
for e2=1:5;
if logical(e1 > cfm1) || logical(e2 > cfn1)
default_mat_config_in{e1, e2}={};
end
end
end
default_mat_config_out=default_mat_config_in;
% Make sure that the output variables exist
SP=[];
vibs=[];
Fs_SP=[];
Fs_vibs=[];
[filename_base, ext]=file_extension(filename);
fexist=exist(filename);
if isequal(fexist, 2) && ~isempty(ext)
switch ext
case 'mat'
keep_config=2;
if ~isequal(default_mat_config_in{4,1}, 1)
while keep_config == 2
% listing of some variables in the matlab file
SP=[];
vibs=[];
Fs_SP=[];
Fs_vibs=[];
t_SP=[];
t_vibs=[];
F_SP=[];
F_vibs=[];
% determine the variable names in tbe data file
bb=whos( '-file', filename);
data_vars=[];
cl_n={'single', 'double', 'int8', 'uint8', 'int16', 'uint16', 'int32', 'uint32', 'int64', 'uint64'};
% make a list of the data variable names containing
% numeric data
for e1=1:length(bb);
data_class=bb(e1).class;
% make sure the data variables are numeric
bb3=strmatch(data_class, cl_n, 'exact');
% make sure the data variables are not empty
data_not_empty=(max(bb(e1).size) > 0);
if ~isempty(bb3) && data_not_empty
data_vars=[data_vars e1];
end
end
% get the data variable names
if ~isempty(data_vars)
num_data_vars=length(data_vars);
var_names=cell(num_data_vars, 1);
% var_names is a cell array of the data variable
% names.
for e1=1:num_data_vars;
var_names{e1}=bb(data_vars(e1)).name;
end
[data_var_ix,ok] = listdlg('Name', 'Choose the data variables do not choose the time or Sampling Rate variables', 'PromptString','Select all the data variables not the Time or sampling rate variables:','SelectionMode','multiple', 'ListSize', [500, 500],'ListString', var_names);
if isequal(ok, 1)
select_data_var=data_vars(data_var_ix);
% update the data_vars2 array
ix_set=setdiff(1:num_data_vars, data_var_ix);
data_vars2=data_vars(ix_set);
num_data_vars2=length(data_vars2);
var_names2=cell(num_data_vars2, 1);
% update the var_names2
for e1=1:num_data_vars2;
var_names2{e1}=bb(data_vars2(e1)).name;
end
else
select_data_var=[];
var_names2=var_names;
data_vars2=data_vars;
end
num_select_data_var=length(select_data_var);
prompt=cell(num_select_data_var,1);
defAns=cell(num_select_data_var,1);
% determine if each data variable contains
% sound data
% vibrations data
% both sound and vibrations
% neither
%
for e1=1:num_select_data_var;
prompt{e1, 1}=['Variable ', bb(select_data_var(e1)).name, ', s for sound, v for vibrations, b for both, n for none'];
switch e1
case 1
defAns{e1, 1}='s';
case 2
defAns{e1, 1}='v';
otherwise
defAns{e1, 1}='n';
end
end
dlg_title='Data Variable Type, s, v, n';
num_lines=1;
% The inputdlg prompts the user to specify the sensors for each
% channel
datatype_str=inputdlg(prompt,dlg_title,num_lines,defAns);
% convert the data_type string into integers
data_type=zeros(num_select_data_var, 1);
for e1=1:num_select_data_var;
switch datatype_str{e1}
case 's'
data_type(e1)=1;
case 'v'
data_type(e1)=2;
case 'b'
data_type(e1)=3;
otherwise
data_type(e1)=4;
end
end
snd_var=select_data_var(find(data_type == 1));
vibs_var=select_data_var(find(data_type == 2));
both_var=select_data_var(find(data_type == 3));
num_snd_vars=length(snd_var);
num_vibs_vars=length(vibs_var);
num_both_vars=length(both_var);
% store the sound variables in the configuration
% cell array
snd_var_name=cell(num_snd_vars, 1);
for e1=1:num_snd_vars;
snd_var_name{e1, 1}=bb(snd_var(e1)).name;
end
% cell array of sound variable names
default_mat_config_out{1,1}=snd_var_name;
vibs_var_name=cell(num_vibs_vars, 1);
for e1=1:num_vibs_vars;
vibs_var_name{e1, 1}=bb(vibs_var(e1)).name;
end
% cell array of vibs variable names
default_mat_config_out{2,1}=vibs_var_name;
both_var_name=cell(num_both_vars, 1);
for e1=1:num_both_vars;
both_var_name{e1, 1}=bb(both_var(e1)).name;
end
% cell array of both variable names
default_mat_config_out{3,1}=both_var_name;
end
% for the data variables that have both microphone and
% accelerometer data, determine which channels
% are microphones and accelerometers
%
both_snd_ch=cell(num_both_vars, 1);
both_vibs_ch=cell(num_both_vars, 1);
if num_both_vars > 0
for e1=1:num_both_vars;
buf=bb(both_var(e1)).size;
bvm1=buf(1);
bvn1=buf(2);
num_channels=min(bvm1, bvn1);
% prompt the user for which channels are microphones and
% accelerometers
prompt=cell(num_channels,1);
defAns=cell(num_channels,1);
for e2=1:num_channels;
prompt{e2, 1}=['Channel ', num2str(e2), ', s for sound, v for vibrations'];
defAns{e2, 1}='s';
end
dlg_title=['For Variable ', both_var_name{e1, 1}, 'Input Sensor Type, s or v'];
num_lines=1;
% The inputdlg prompts the user to specify the sensors for each
% channel
datatype_str=inputdlg(prompt,dlg_title,num_lines,defAns);
data_type=ones(num_channels, 1);
for e2=1:num_channels;
k = strfind(datatype_str{e2}, 's');
if ~isempty(k)
data_type(e2)=1;
else
data_type(e2)=0;
end
end
% Calculate the number of microphones and accelerometers
both_snd_ch{e1, 1}=find(data_type == 1);
both_vibs_ch{e1, 1}=find(data_type == 0);
end
end
default_mat_config_out{3,4}=both_snd_ch;
default_mat_config_out{3,5}=both_vibs_ch;
% get the Time or Sampling Rate variable names
% tosr is an acronym for Time Or Sampling Rate
num_data_vars2=length(data_vars2);
var_names2=cell(num_data_vars2, 1);
for e1=1:num_data_vars2;
var_names2{e1}=bb(data_vars2(e1)).name;
end
% tosr acronym for (time or sampling rate)
[tosr_var, ok] = listdlg('Name', 'Choose the Time or Rate variables', 'PromptString','Select the Time and Sampling Rate variables:','SelectionMode','multiple', 'ListSize', [500, 500],'ListString', var_names2);
% If there are not any time or frequency variables prompt user
% to input the sampling rate for each data variable
if isempty(tosr_var) || isequal(ok, 0)
snd_tosr_bool={};
Fs_SP={};
vibs_tosr_bool={};
Fs_vibs={};
both_tosr_bool={};
Fs_both={};
% Enter the sampling rates for the sound variables
if num_snd_vars > 0
prompt=cell(num_snd_vars,1);
defAns=cell(num_snd_vars,1);
snd_tosr_bool=cell(num_snd_vars,1);
for e1=1:num_snd_vars;
snd_tosr_bool{e1}=0;
prompt{e1, 1}=['Enter Sampling Rate (Hz) for Variable ', default_mat_config_out{1,1}{e1,1}];
defAns{e1, 1}='50000';
end
dlg_title='Enter Sampling Rates for the Sound Variables';
num_lines=1;
Fs_SP_var=inputdlg(prompt,dlg_title,num_lines,defAns);
num_Fs_SP=length(Fs_SP_var);
Fs_SP=cell(num_Fs_SP, 1);
for e1=1:num_Fs_SP;
Fs_SP{e1}=str2double(Fs_SP_var{e1});
end
end
% Enter the sampling rates for the vibrations variables
if num_vibs_vars > 0
prompt=cell(num_vibs_vars,1);
defAns=cell(num_vibs_vars,1);
vibs_tosr_bool=cell(num_vibs_vars,1);
for e1=1:num_vibs_vars;
vibs_tosr_bool{e1}=0;
prompt{e1, 1}=['Enter Sampling Rate (Hz) for Variable ', default_mat_config_out{2,1}{e1,1}];
defAns{e1, 1}='5000';
end
dlg_title='Enter Sampling Rates for the Vibrations Variables';
num_lines=1;
Fs_vibs_var=inputdlg(prompt,dlg_title,num_lines,defAns);
num_Fs_vibs=length(Fs_vibs_var);
Fs_vibs=cell(num_Fs_vibs, 1);
for e1=1:num_Fs_vibs;
Fs_vibs{e1}=str2double(Fs_vibs_var{e1});
end
end
% Enter the sampling rates for the both sound and vibrations variables
if num_both_vars > 0
prompt=cell(num_both_vars,1);
defAns=cell(num_both_vars,1);
both_tosr_bool=cell(num_both_vars,1);
for e1=1:num_both_vars;
both_tosr_bool{e1}=0;
prompt{e1, 1}=['Enter Sampling Rate (Hz) for Variable ', default_mat_config_out{3,1}{e1,1}];
defAns{e1, 1}='50000';
end
dlg_title='Enter Sampling Rates for the Variables with Both Sound and Vibrations';
num_lines=1;
Fs_both_var=inputdlg(prompt,dlg_title,num_lines,defAns);
num_Fs_both=length(Fs_both_var);
Fs_both=cell(num_Fs_both, 1);
for e1=1:num_Fs_both;
Fs_both{e1}=str2double(Fs_both_var{e1});
end
end
% Store sampling rate arrays to the
% configuration cell array
default_mat_config_out{1,2}=Fs_SP;
default_mat_config_out{2,2}=Fs_vibs;
default_mat_config_out{3,2}=Fs_both;
default_mat_config_out{1,3}=snd_tosr_bool;
default_mat_config_out{2,3}=vibs_tosr_bool;
default_mat_config_out{3,3}=both_tosr_bool;
else
num_tosr_var=length(tosr_var);
tosr_var_names=cell(num_tosr_var, 1);
for e1=1:num_tosr_var;
tosr_var_names{e1}=var_names2{tosr_var(e1)};
end
% determine if the variables are time or sampling rate variables.
tosr_bool1=ones(num_tosr_var, 1);
for e1=1:num_tosr_var;
tosr_bool1(e1) = menu(['Is the variable ', tosr_var_names{e1} ,' Time or Sampling Rate'], 'Time', 'Sampling Rate');
end
% determine the correspondence between data variables and
% Time or Sampling Rate variables
snd_tosr_var=cell(num_snd_vars, 1);
snd_tosr_bool=cell(num_snd_vars, 1);
for e1=1:num_snd_vars;
[buf,ok] = listdlg('Name', 'Choose the corresponsing Time or Sampling Rate Variable', 'PromptString',['For Sound Variable ',default_mat_config_out{1,1}{e1,1} ],'SelectionMode','single', 'ListSize', [500, 500],'ListString', tosr_var_names);
if ok
snd_tosr_var{e1}=tosr_var_names{buf};
snd_tosr_bool{e1}=tosr_bool1(buf);
else
snd_tosr_var{e1}={};
snd_tosr_bool{e1}=[];
end
end
vibs_tosr_var=cell(num_vibs_vars, 1);
vibs_tosr_bool=cell(num_vibs_vars, 1);
for e1=1:num_vibs_vars;
[buf,ok] = listdlg('Name', 'Choose the corresponsing Time or Sampling Rate Variable', 'PromptString',['For Vibrations Variable ',default_mat_config_out{2,1}{e1,1} ],'SelectionMode','single', 'ListSize', [500, 500],'ListString', tosr_var_names);
if ok
vibs_tosr_var{e1}=tosr_var_names{buf};
vibs_tosr_bool{e1}=tosr_bool1(buf);
else
vibs_tosr_var{e1}={};
vibs_tosr_bool{e1}=[];
end
end
both_tosr_var=cell(num_both_vars, 1);
both_tosr_bool=cell(num_both_vars, 1);
for e1=1:num_both_vars;
[buf,ok] = listdlg('Name', 'Choose the corresponsing Time or Sampling Rate Variable', 'PromptString',['For Variable for both Sound and Vibrations ',default_mat_config_out{3,1}{e1,1} ],'SelectionMode','single', 'ListSize', [500, 500],'ListString', tosr_var_names);
if ok
both_tosr_var{e1}=tosr_var_names{buf};
both_tosr_bool{e1}=tosr_bool1(buf);
else
both_tosr_var{e1}={};
both_tosr_bool{e1}=[];
end
end
default_mat_config_out{1,2}=snd_tosr_var;
default_mat_config_out{2,2}=vibs_tosr_var;
default_mat_config_out{3,2}=both_tosr_var;
default_mat_config_out{1,3}=snd_tosr_bool;
default_mat_config_out{2,3}=vibs_tosr_bool;
default_mat_config_out{3,3}=both_tosr_bool;
end
keep_config=menu('Keep Variable Configuration', 'Yes', 'No');
end
automatic_cal=menu('Use Sensor Configuration for all files', 'Yes', 'No');
default_mat_config_out{4,1}=automatic_cal;
end
% load matlab file
load(filename);
% check the configuration of the data
bool_config=ones(3,2);
bool_config2=ones(3, 1);
for e1=1:3;
for e2=1:3;
if iscell(default_mat_config_out{e1,e2})
if length(default_mat_config_out{e1,e2}) < 1
bool_config(e1, e2)=0;
else
if iscell(default_mat_config_out{e1,e2}{1,1})
if length(default_mat_config_out{e1,e2}{1,1}{1,1}) < 1
bool_config(e1, e2)=0;
else
if iscell(default_mat_config_out{e1,e2}{1,1}{1,1})
bool_config(e1, e2)=0;
else
if ~isempty(default_mat_config_out{e1,e2}{1,1}{1,1})
bool_config(e1, e2)=1;
else
bool_config(e1, e2)=0;
end
end
end
else
if ~isempty(default_mat_config_out{e1,e2}{1,1})
bool_config(e1, e2)=1;
else
bool_config(e1, e2)=0;
end
end
end
else
if ~isempty(default_mat_config_out{e1,e2})
bool_config(e1, e2)=1;
else
bool_config(e1, e2)=0;
end
end
end
bool_config2(e1, 1)=prod(bool_config(e1, :));
end
% Load the data variables
% Determine the number of sound variables that can be loaded
% correctly
% Append the sound data to a cell array
% Append the sound sampling rates to a cell array
SP_var={};
Fs_SP_var={};
if isequal(bool_config2(1,1), 1)
num_snd_vars=length(default_mat_config_out{1,1});
num_snd_tosr_vars=length(default_mat_config_out{1,2});
num_snd_tosr_bools=length(default_mat_config_out{1,3});
num_snd_vars=min([num_snd_vars, num_snd_tosr_vars, num_snd_tosr_bools]);
SP_var=cell(num_snd_vars, 1);
Fs_SP_var=cell(num_snd_vars, 1);
for e1=1:num_snd_vars;
SP_var_buf=eval(default_mat_config_out{1,1}{e1,1});
[svb_m1, svbn1]=size(SP_var_buf);
if svb_m1> svbn1
SP_var_buf=SP_var_buf';
[svb_m1, svbn1]=size(SP_var_buf);
end
SP_var{e1}=SP_var_buf;
if ischar(default_mat_config_out{1,2}{e1,1})
Fs_buf=eval(default_mat_config_out{1,2}{e1,1});
else
Fs_buf=default_mat_config_out{1,2}{e1,1};
end
if max(size(Fs_buf)) > 1
Fs_buf=Fs_buf(2)-Fs_buf(1);
else
Fs_buf=Fs_buf(1);
end
if isequal(default_mat_config_out{1,3}{e1, 1}, 1)
Fs_SP_var{e1}=1./Fs_buf;
else
Fs_SP_var{e1}=Fs_buf;
end
end
end
% Append the vibrations data to a cell array
% Append the vibrations sampling rates to a cell array
vibs_var={};
Fs_vibs_var={};
if isequal(bool_config2(2,1), 1)
num_vibs_vars=length(default_mat_config_out{2,1});
num_vibs_tosr_vars=length(default_mat_config_out{2,2});
num_vibs_tosr_bools=length(default_mat_config_out{2,3});
num_vibs_vars=min([num_vibs_vars, num_vibs_tosr_vars, num_vibs_tosr_bools]);
vibs_var=cell(num_vibs_vars, 1);
Fs_vibs_var=cell(num_vibs_vars, 1);
for e1=1:num_vibs_vars;
vibs_var_buf=eval(default_mat_config_out{2,1}{e1,1});
[vvb_m1, vvbn1]=size(vibs_var_buf);
if vvb_m1> vvbn1
vibs_var_buf=vibs_var_buf';
[vvb_m1, vvbn1]=size(vibs_var_buf);
end
vibs_var{e1}=vibs_var_buf;
if ischar(default_mat_config_out{2,2}{e1,1})
Fs_buf=eval(default_mat_config_out{2,2}{e1,1});
else
Fs_buf=default_mat_config_out{2,2}{e1,1};
end
if max(size(Fs_buf)) > 1
Fs_buf=Fs_buf(2)-Fs_buf(1);
else
Fs_buf=Fs_buf(1);
end
if isequal(default_mat_config_out{2,3}{e1, 1}, 1)
Fs_vibs_var{e1}=1./Fs_buf;
else
Fs_vibs_var{e1}=Fs_buf;
end
end
end
% Append the both data to a cell array
% Append the both sampling rates to a cell array
if isequal(bool_config2(3,1), 1)
num_both_vars=length(default_mat_config_out{3,1});
num_both_tosr_vars=length(default_mat_config_out{3,2});
num_both_tosr_bools=length(default_mat_config_out{3,3});
num_both_vars=min([num_both_vars, num_both_tosr_vars, num_both_tosr_bools]);
for e1=1:num_both_vars;
both_var=eval(default_mat_config_out{3,1}{e1,1});
[bv_m1, bvn1]=size(both_var);
if bv_m1> bvn1
both_var=both_var';
[bv_m1, bvn1]=size(both_var);
end
if ischar(default_mat_config_out{3,2}{e1,1})
Fs_buf=eval(default_mat_config_out{3,2}{e1,1});
else
Fs_buf=default_mat_config_out{3,2}{e1,1};
end
if max(size(Fs_buf)) > 1
Fs_buf=Fs_buf(2)-Fs_buf(1);
else
Fs_buf=Fs_buf(1);
end
if isequal(default_mat_config_out{3,3}{e1, 1}, 1)
Fs_both_var=1./Fs_buf;
else
Fs_both_var=Fs_buf;
end
both_snd_ch=default_mat_config_out{3,4}{e1,1};
if length(both_snd_ch) > 1
buf=(both_snd_ch > 0);
if isequal(all(buf), 1)
SP_m1=length(SP_var);
SP_var{SP_m1+1}=both_var(both_snd_ch, :);
Fs_SP_var{SP_m1+1}=Fs_both_var;
end
end
both_vibs_ch=default_mat_config_out{3,5}{e1,1};
if length(both_vibs_ch) > 1
buf=(both_vibs_ch > 0);
if isequal(all(buf), 1)
vibs_m1=length(vibs_var);
vibs_var{vibs_m1+1}=both_var(both_vibs_ch, :);
Fs_vibs_var{vibs_m1+1}=Fs_both_var;
end
end
end
end
% data may be stored in single precision or other formats
% and must be changed to double precision.
%
SP={};
for e1=1:length(SP_var);
[buf]=convert_double(SP_var{e1});
SP_var{e1}=[];
SP{e1}=buf;
end
vibs={};
for e1=1:length(vibs_var);
[buf]=convert_double(vibs_var{e1});
vibs_var{e1}=[];
vibs{e1}=buf;
end
Fs_SP={};
for e1=1:length(Fs_SP_var);
[buf]=convert_double(Fs_SP_var{e1});
Fs_SP_var{e1}=[];
Fs_SP{e1}=buf;
end
Fs_vibs={};
for e1=1:length(Fs_vibs_var);
[buf]=convert_double(Fs_vibs_var{e1});
Fs_vibs_var{e1}=[];
Fs_vibs{e1}=buf;
end
case 'wav'
% Read the wave file
% y is the data
% fs is the sampling rate
[y, fs, nbits, opts] = wavread(filename);
[num_data, num_channels]=size(y);
cal_snd=[];
cal_vibs=[];
Fs_SP=fs;
Fs_vibs=fs;
keep_config=2;
% Determine whether to use the automated configuration
if ~isequal(default_wav_settings_in{4,1}, 1)
while keep_config == 2
% prompt the user for which channels are microphones and
% accelerometers
prompt=cell(num_channels,1);
defAns=cell(num_channels,1);
for e1=1:num_channels;
prompt{e1, 1}=['Channel ', num2str(e1), ', s for sound, v for vibrations'];
defAns{e1, 1}='s';
end
dlg_title='Input Sensor Type, s or v';
num_lines=1;
% The inputdlg prompts the user to specify the sensors for each
% channel
datatype_str=inputdlg(prompt,dlg_title,num_lines,defAns);
data_type=ones(num_channels, 1);
for e1=1:num_channels;
k = strfind(datatype_str{e1}, 's');
if ~isempty(k)
data_type(e1)=1;
else
data_type(e1)=0;
end
end
% Calculate the number of microphones and accelerometers
num_snd=length(data_type(data_type==1));
num_vibs=length(data_type(data_type==0));
if num_snd+num_vibs < num_channels
data_type=ones(num_channels, 1);
num_vibs=0;
num_snd=num_channels;
end
% get the calibration sensitivities for each sound channel
if num_snd > 0
prompt=cell(num_snd, 1);
defAns=cell(num_snd, 1);
for e1=1:num_snd;
prompt{e1, 1}=['Mic ', num2str(e1), ', mV/Pa'];
defAns{e1, 1}=num2str(40.00, 1);
end
dlg_title='Input Calibration Sensitivity mV/Pa';
num_lines=1;
cal_snd_str=inputdlg(prompt,dlg_title,num_lines,defAns);
cal_snd=ones(num_snd, 1);
for e1=1:num_snd;
cal_snd(e1)=str2double(cal_snd_str{e1});
end
end
% get the calibration sensitivities for each accelerometer
% channel
if num_vibs > 0
prompt=cell(num_vibs, 1);
defAns=cell(num_vibs, 1);
for e1=1:num_vibs;
prompt{e1, 1}=['Accel ', num2str(e1), ', mV/g'];
defAns{e1, 1}=num2str(10.00, 1);
end
dlg_title='Input Calibration Sensitivity mV/g Program converts from g''s to m/s^2 ';
num_lines=1;
% The inputdlg prompts the user to specify the calibration
% sensitivities
cal_vibs_str=inputdlg(prompt,dlg_title,num_lines,defAns);
cal_vibs=ones(num_vibs, 1);
for e1=1:num_vibs;
cal_vibs(e1)=str2double(cal_vibs_str{e1});
end
end
keep_config=menu('Keep Sensor Configuration', 'Yes', 'No');
end
automatic_cal=menu('Use Sensor Configuration for all files', 'Yes', 'No');
else
% Apply the input settings
data_type=default_wav_settings_in{1, 1};
cal_snd=default_wav_settings_in{2, 1};
cal_vibs=default_wav_settings_in{3, 1};
automatic_cal=default_wav_settings_in{4, 1};
% Calculate the number of microphones and accelerometers
num_snd=length(data_type(data_type==1));
num_vibs=length(data_type(data_type==0));
end
% Set the output settings
default_wav_settings_out{1,1}=data_type;
default_wav_settings_out{2,1}=cal_snd;
default_wav_settings_out{3,1}=cal_vibs;
default_wav_settings_out{4,1}=automatic_cal;
[cal_snd_a]=meshgrid(cal_snd, 1:num_data );
[cal_vibs_a]=meshgrid(cal_vibs, 1:num_data );
% Calibrate the wave files and convert to engineering units
if num_snd > 0
SP=1./cal_snd_a.*y(:, data_type==1); % Pa
else
SP=[]; % Pa
end
if num_vibs > 0
vibs=9.806./cal_vibs_a.*y(:, data_type==0); % m/s^2
else
vibs=[]; % m/s^2
end
otherwise
end
end
[m1 n1]=size(SP);
if m1 > n1
SP=SP';
[m1 n1]=size(SP);
end
[m1 n1]=size(vibs);
if m1 > n1
vibs=vibs';
[m1 n1]=size(vibs);
end
|
|
