%
%    the dynamic time constant of the OHC control path of the AN model by
%    Zhang et al. 2001
%
DC = 0.37; R0 = 0.05; R1=1-R0;
TN = 1.0; TW = 3.0; % TN and TW chosen arbitrarily
time = 0:0.0001:0.005;
V_LP = DC * exp(-800.*time); % the LP filter output runs between [1,0]

Tau = TN*(R0+ R1*( ((TW/TN)-R0)/R1).^(V_LP/DC));

figure(1); clf;
plot(time, Tau,'r-','linewidth',2);
hold on
plot([0 0.005],[TN TN],'k--')
plot([0 0.005],[TW TW],'k--')
axis([0 0.005 0.5*TN 1.2*TW])
text(0.001, TN+0.1, 'TN','fontsize',12,'fontweight','bold');
text(0.001, TW-0.1, 'TW','fontsize',12,'fontweight','bold');


xlabel('TIME [s]','fontsize',13,'fontweight','bold');
ylabel('OHC TIME CONSTANT [ms]','fontsize',13,'fontweight','bold');
set(gca, 'fontsize',12,'fontweight','bold');
print -djpeg exc6_5