simoae.m

A MATLAB Simulation of Transient Otoacoustic Emission signals

function simoae(command_str, radbutt)
%
%			            	      PROGRAM SIMOAE( ) 
%
%	Graphic User Interface (GUI) for the simulation  of transient evoked otoacoustic emissions.
%   
%   Input signal files have extension   *.sig. The following three signal files
%   are available: click.sig (regular click), cl_filt.sig (click filtered by a middle ear with
%   irregular transfer function), noise.sig (noise-like input signal).
%   
%   Auxiliary files (see comments in files): 
%
%   ALLDATA.M  loads  all  program data. 
%   BMDATA.M loads data generated by BMOPS.M and prepare data for backward Euler method 
%   BMOPS.M builds the distributed parameters for the basilar membrane (BM) motion equation;
%   TMOPS.M  builds the distributed parameters for the tectorial membrane (TM) motion equation;  
%   GAUSSGRID.M creates a set of BM coordinates that covers the axis 0-1 with Gaussian density 
%   INTERPOL.M routine for data interpolation 
%   GREENF.M, G0.M, GPEAKS.M, PSIGMA.M routines for the generation of Green's function and
%                   stapes-BM coupling coefficients.
%   MASS.M    Organ of Corti  mass per unit length.
%   SIGMOID.M  The profile of  the mechano-mechanic saturation as a function of stereocilia diplacement  
%   BMWDATA.MAT, CDDATA.MAT Data on the human cochlea geometry. 
%   LAMBDATA.MAT  Coefficients for OHC force control as determined by a recursive procedure
%                               This file contains "lam" (smooth CA gain) and "lam1" (damaged CA gain)   
% 
%   


%   ------------------------------  NOTE ON THIS STRUCTURE OF THIS ROUTINE ---------------
%   All variables shared by the different calls of  SIMOAE.M routine must be declared as global variables or
%   retrived by "get(gcf, 'Userdata')".
%
%	CommandString: string specifying the GUI mode to be called. The calling string
%   can be one of the following:
%
%   'Initialize' - This is the default call. It generates the graphic environment and stores all
%   GUI's item handles (axes, lines, menus, buttons etc.) to the figure's Userdata buffer.
%   Any command-string different from this causes the retrieval of  such handles. 
%
%	'Load' - This command-string is passed when button LOAD is pushed.
%   It activates the procedure for loading input signal files.
%
%   'InputLevel' - Command string passed when the vertical slider at the right of the top diagram is used.
%
%   'Help' - Command string passed when button Help is pressed. Displays a help dialog box. 
%
%	'Start' - This command-string is passed when button (RE)START is pushed.
%   It sets START_FLAG =1, resets all dynamical  variables to zero initial values and calls
%   the routine  passing the command-string 'Compute'.
%
%	'Continue' - This command-string is passed when button CONTINUE, is pushed.
%   It sets START_FLAG = 1 without resetting the dynamical variables to initial values, thus
%   the computation continues from the time at which it was left by the command 'Stop'.	
%
%	'Stop' - This command-string is passed when button STOP, is pushed.
%    It sets START_FLAG = 0, thus interrupting computations.	
%
%   'BMScale' - Command string passed when the vertical slider at the right of middle diagram is used.
%
%   'BMPopup' - Command string passed when the popup menu on the top of the bottom diagram is used.
%
%   'VPScale' - Command string passed when the vertical slider at the right of bottom diagram is used.
%
%   'RadioCtl' - Command string passed when one of the three radio buttons at bottom-right is chosen .
%
%	'Compute' - This command-string is passed by other commands. It activates the computation cycle
%   and updates graphics via functions "set(h_Lo,'Ydata', ydata1)"; "set(h_L1,'Ydata', ydata1)",
%   "set(h_L2,'Ydata', ydata2)", "set(h_txt,'String',num2str(t))", etc. 
%
%	 Second parameter RadioButton = auxiliary parameter to handle radio button mutual exclusion
%
%    〤opyright by Renato Nobili, Padova University, Italy, 2003 - Not for commercial use.
%   The data file LAMBDAS.MAT, which assures a regular cochlear amplifier gain, was computed by a recursive
%   procedure created by Ales Vetesnik - Palack