matlab时间序列分析工具程序

包含了用分析时间序列是常用的程序源代码，及程序的请详细说明

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