avgpic.m

A MATLAB tool for analysis of Transient Otoacoustic Emission signals

function avgpic(hfig, avgfile, avgtitle)

load avgfile

% get handles of  axes
haxs = get(hfig,'UserData');
haxstim  = haxs(1);
haxt = haxs(2);    
haxSrs = haxs(3);    
haxSTr = haxs(4);  
haxsmS = haxs(5);  
htab = haxs(6);  
haxfftS = haxs(7);
haxsurf = haxs(8);



%name
ht = get(haxt,'Userdata');
set(ht,'String',[int2str(ipat) '. ' name ', ' ear ', ' date])

% Stimulus
set(0,'CurrentFigure',hfig);
set(hfig,'CurrentAxes',haxstim); 
delete(get(haxstim,'Children'));
hstim=line('XData',(0:127)/25,...
           'YData',stim);
  
% FFT-Spectra
set(0,'CurrentFigure',hfig);
set(hfig,'CurrentAxes',haxfftS);
delete(get(haxfftS,'Children'));
nzd = max(size(Sr))+2;
hfftSr = patch(dupends(f),[S_ymin dbs(Sr) S_ymin],zeros(nzd,1),[.8 .8 .8]);
hfftSn = patch(dupends(f),[S_ymin dbs(Sn) S_ymin],ones(nzd,1),[.5 .5 .5]);

% Spectra
nzd = max(size(Srsc))+2;
set(0,'CurrentFigure',hfig);
set(hfig,'CurrentAxes',haxSrs);
delete(get(haxSrs,'Children'));
hSr=patch(dupends(fb),[S_ymin dbs(Srsc) S_ymin],zeros(nzd,1),[.8 .8 .8]);  
hSn=patch(dupends(fb),[S_ymin dbs(Snsc) S_ymin],ones(nzd,1),[.5 .5 .5]); 
hSs=line(fb,dbs(Sssc));                  

% Smoothed spectra
set(0,'CurrentFigure',hfig);
set(hfig,'CurrentAxes',haxsmS);
delete(get(haxsmS,'Children'));
hsmSr=patch(dupends(fb),[Ssm_ymin smooth(dbs(Srsc),gaussw(7)) Ssm_ymin],zeros(nzd,1),[.8 .8 .8]);
hsmSn=patch(dupends(fb),[Ssm_ymin smooth(dbs(Snsc),gaussw(7)) Ssm_ymin],ones(nzd,1),[.5 .5 .5]);

% T-F-Spectrum & Noise level
infb = fbtf>=1.8 & fbtf<=hz2bark(6000); 
infn = fn>=1.8 & fn<=hz2bark(6000); 
int = ttf>=2.5 & ttf<=21;
intn = tn>=2.5 & tn<=21;
rzdata = rzdata(int,infb);
rcdata = rcdata(int,infb);
Pn = Pn(intn,infn);
fbtf = fbtf(infb);
fn = fn(infn);
ttf = ttf(int);
tn = tn(intn);
rzdata(rzdata<zamin) = NaN*rzdata(rzdata<zamin);