Rate-Varying Poisson Process Generator 1.0

Rate-Varying Poisson-Process generator. T=PechePourPoisson(rate);T=PechePourPoisson(rate,dt);[T tvec]=PechePourPoisson(rate,dt);Example:rate=conv(randn(1,10000),ones(1,1000),'same');[T t]=PechePourPoisson(rate,.001);plot(t,rate);addline(T); addline(0,'h'); % You may not have the addline functionThis function spits out a poisson process with a rate that can vary withtime. Each sample of the "rate" vector is considered to be the beginningof a time bin of width dt. The number of events occurring in each bin issampled from a Poisson Distribution. The timing of each event within thebin is then randomized.Note: You can easily specify a constant-rate process running for acertain time by just making "rate" a scalar rate and setting dt to yourdesired time span.Inputs Description"rate" Vector representing the rate as a function of time. Negativerates are interpreted as zero rates.("dt") Optional scalar defining the time resolution of the ratevector. If neglected, dt=1, and the rate is interpreted asmeaning "events per time-sample". If dt is set, rate will bescaled by the value of dt. eg, say you want 90events/s, wherethe rate vector represents the spiking rate as a function oftime with millisecond resolution. Then dt=0.001, andrate(i)=90, where i is the particular instant with a 90event/s rate.Outputs Description"T" Vector of event times. The first sample of "rate" is considered to be time 0. If dt is specified, times will be scaled accordingly."tvec" If you're too lazy to make your own time vector for your ratevector, this function can optionally do it for you.EnjoyPeteroconnorp _at_ ethz _dot_ ch

Requirements:	No special requirements
Platforms:	Matlab
Keyword:	Defining,  Dt3d1,  Eventss,  Meaning,  Neglected,  Resolution,  Scaled,  Timesample
Users rating:	0/10