function?[prperpulse?stdperpulse]=pulsetrain（stimpropertiespulseratepulsewidth）

%?[prperpulse?stdperpulse]=pulsetrain（stimpropertiespulseratepulsewidth）
%
%?INPUT:
%
%?stim?is?the?stimulus?intensity?in?units?of?microA?NOT?dB
%
%?properties?is?a?structure?containing?the?model?parameters:
%
%	properties.resptype?is?the?fiber?I/O?function?type:
%		‘stp‘?=?deterministic?（step-function）?model
%		‘3_p‘?=?3-piece?linear?model
%		‘erf‘?=?stochastic?（error-function）?model
%
%	properties.thrsh?is?the?matrix?of?fiber?thresholds
%
%	properties.noisestd?is?the?matrix?of?fiber?noise?standard?deviations
%
%?pulserate?is?the?rate?of?stimulation?in?pulses/s
%
%?pulsewidth?is?the?phase?duration?in?microseconds/phase
%
%?OUTPUT:
%
%?prperpulse?is?the?discharge?probability?（mean?discharge?rate）?per?pulse
%
%?stdperpulse?is?the?standard?deviation?in?discharge?rate?per?pulse
%
%?Usage?Agreement:	Any?publication?of?results?obtained?using?this?model?should?cite?both-
%
%		[1]	Bruce?I.?C.?Irlicht?L.?S.?White?M.?W.?O‘Leary?S.?J.?Dynes
%			S.?Javel?E.?Clark?G.?M.?（1999）?“A?stochastic?model?of?the
%			electrically?stimulated?auditory?nerve:?Pulse-train?response“
%			IEEE?Transactions?on?Biomedical?Engineering?46（6）:630-637.
%
%	and
%
%		[2]	Bruce?I.?C.?（1997）?Spatiotemporal?coding?of?sound?in?the?auditory
%			nerve?for?cochlear?implants?PhD?Thesis?The?University?of?Melbourne?Australia.

%?Copyright?（c）?1996-1999?Ian?Bruce
%?Created?by?Ian?Bruce?1996
%?Released?for?public?distribution?as?version?1.1?1999

Fs?=?1e5;?%?Sampling?frequency?in?Hz

resptype?=?properties.resptype;
thrsh????=?properties.thrsh;
noisestd?=?properties.noisestd;

%?Run?pulse-train?model?for?each?individual?fiber?at?each?stimulus?intensity
for?i=1:size（stim1）
???for?j=1:size（stim2）
??????disp（[num2str（pulserate）?‘Hz?-?Run?‘?num2str（size（stim2）*（i-1）+j）?‘/‘?num2str（size（stim1）*size（stim2））]）
??????[prperpulse（ij）?stdperpulse（ij）]?=?pulsetrain_perelement（stim（ij）resptypethrsh（ij）noisestd（ij）pulseratepulsewidthFs）;
???end
end

function?[prperpulse?stdperpulse?isiprs?Pinf]=pulsetrain_perelement（levelresptypethrshnoisestdpulseratepulsewidthFs）
%
%?[prperpulse?stdperpulse?isiprs?Pinf]=pulsetrain_perelement（levelresptypethrshnoisestdpulseratepulsewidthFs）

%?For?Vrefr?relative?to?Vthr
absref?=?0.7e-3;	%?Absolute?refractory?period?in?seconds
if?pulserate>1/absref
???error（[‘this?model?is?not?appropriate?for?pulserates?above?‘?num2str（round（1/absref））?‘?pulses/s‘]）
end
relref?=?20e-3;	%?Relative?refractory?period?in?seconds
magn?=?0.97;		%?Relative?magnitude?of?refactory?function
tconst?=?1.32e-3;?%?Time-constant?of?refactory?function?in?seconds
y0?=?1;				%?Asymptotic?relative?magnitude?of?refactory?function

dt?=?1/Fs;	%?bin?size?in?seconds

pulsewidth?=?pulsewidth*1e-6;?%?convert?from?microseconds/phase?to?seconds/phase

%?Need?to?calculate?number?of?bins?per?pulse?to?calculate?mean?and?variance?of?renewal?time
%?depending?on?which?bin?the?prev

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????122712??2000-07-12?11:16??matlab_download\xg_timeseries\tstool\poster.pdf
?????目錄???????????0??2000-07-12?11:18??matlab_download\xg_timeseries\tstool\TSTOOL?Demo?2.files\
?????文件???????10585??2000-07-12?11:17??matlab_download\xg_timeseries\tstool\TSTOOL?Demo?2.files\demo2a.gif
?????文件????????1962??2000-07-12?11:17??matlab_download\xg_timeseries\tstool\TSTOOL?Demo?2.files\demo2b.gif
?????文件????????2372??2000-07-12?11:17??matlab_download\xg_timeseries\tstool\TSTOOL?Demo?2.files\demo2c.gif
?????文件????????4399??2000-07-12?11:17??matlab_download\xg_timeseries\tstool\TSTOOL?Demo?2.files\demo2d.gif
?????文件????????2110??2000-07-12?11:17??matlab_download\xg_timeseries\tstool\TSTOOL?Demo?2.files\demo2e.gif
?????文件????????2096??2000-07-12?11:17??matlab_download\xg_timeseries\tstool\TSTOOL?Demo?2.files\demo2f.gif
?????文件???????73912??2000-07-12?11:17??matlab_download\xg_timeseries\tstool\TSTOOL?Demo?2.files\logo.gif
?????文件????????3340??2000-07-12?11:18??matlab_download\xg_timeseries\tstool\TSTOOL?Demo?2.htm
?????文件?????4571598??2000-07-05?20:43??matlab_download\xg_timeseries\tstool\OpenTSTOOL1_02.tgz
?????文件??????959280??2000-07-05?20:27??matlab_download\xg_timeseries\tstool\pentium_mexdll1_02.tgz
?????文件????????8225??2000-07-05?20:14??matlab_download\xg_timeseries\tstool\TSTOOL?Installation?Page.htm
?????文件?????????310??2000-07-15?11:29??matlab_download\xg_timeseries\tstool\readme.txt
?????文件?????2320325??2000-07-15?11:14??matlab_download\xg_timeseries\tftb\tftb.zip
?????文件???????93211??2000-07-15?11:17??matlab_download\xg_timeseries\tftb\newfile.zip
?????文件?????????277??2000-07-15?11:16??matlab_download\xg_timeseries\tftb\readme.txt
?????文件???????54219??2000-07-05?20:15??matlab_download\xg_timeseries\tsa\tsa244.zip
?????文件????????4492??2000-07-12?11:12??matlab_download\xg_timeseries\tsa\Time?Series?Analysis?with?Matlab.htm
?????文件?????????217??2000-07-15?11:28??matlab_download\xg_timeseries\tsa\readme.txt
?????文件????????3825??2000-07-15?11:20??matlab_download\xg_timeseries\cochlear\Ian?Bruce‘s?Cochlear?Implant?Models.htm
?????目錄???????????0??2000-07-15?11:20??matlab_download\xg_timeseries\cochlear\Ian?Bruce‘s?Cochlear?Implant?Models.files\
?????文件????????1132??2000-07-15?11:18??matlab_download\xg_timeseries\cochlear\Ian?Bruce‘s?Cochlear?Implant?Models.files\Home.gif
?????文件????????9016??2000-07-15?11:18??matlab_download\xg_timeseries\cochlear\Ian?Bruce‘s?Cochlear?Implant?Models.files\Pink_and_White.gif
?????文件????????1782??2000-07-15?11:20??matlab_download\xg_timeseries\cochlear\singlepulse.m
?????文件????????7611??2000-07-15?11:20??matlab_download\xg_timeseries\cochlear\totalan.m
?????文件????????7229??2000-07-15?11:21??matlab_download\xg_timeseries\cochlear\pulsetrain.m
?????文件?????????188??2000-07-15?11:21??matlab_download\xg_timeseries\cochlear\readme.txt
?????文件????????3447??2000-07-12?11:09??matlab_download\xg_timeseries\Subject?Re?ARMA?matlab?code.htm