First-order exponential smoother
The simplest model how the voltage $x$ at a synapse responds to input $u$ is a first-order filter:
$$\tau \dot x = -x + u.$$
This corresponds to convolving signal $u(t)$ with exponential filter $\operatorname H(t) \exp(-t/\tau)$, where $\operatorname H(\cdot)$ is the Heaviside step function:
$$\begin{aligned}x(t) &= h(t) * u(t)\\h(t)&=\operatorname H(t) \exp(-t/\tau).\end{aligned}$$
The alpha function
A first-order filter has a discontinuous jump in response to an abrupt inputs (like spikes). A more realistic response is the "alpha function" $t\cdot \exp(-t)$. The alpha function can be obtained by convolving two first decay functions (i.e. chaining together two first-order filters):
$$\begin{aligned}\tau \dot x_1 &= -x_1 + u\\\tau \dot x_2 &= -x_2 + x_1.\end{aligned}$$
This is sometimes written in the compact notation
$$ \left(\tau \frac{d }{dt} + 1\right)^2 x = u.$$
Higher orders
You can repeat this operation many times, obtaining responses with increasing smoothness. The family $t^n\cdot \exp(-t)$ reflects $n+1$ feed-forward variables coupled by exponential decay $\dot x_n=x_{n-1}-x_n$. The integral of $t^n \exp(-t)$ grows with $n$. To normalize, divide by $n! = \Gamma(n+1)$:
$$h(t) = \operatorname H(t) \frac {t^n}{\Gamma(n+1)}e^{-t}.$$
The response of $t^n\exp(-t)$ peaks at time $t=n$. Rescale time with $t\gets nt$ to get a peak response at $t=1$. To keep the integral of the response normalized when rescaling, multiply by $n$.
$$h(t) = \operatorname H(t) n \frac {(nt)^n}{\Gamma(n+1)}e^{-nt}$$
This is equivalent to choosing a time constant $\tau=1/n$ for each of the $n$ filtering stages. To place the peak response at time $t_0$, set $\tau = t_0/n$.
To stabilize the variance instead of time-to-peak, rescale time by $1/\sqrt{n+1}$. This corresponds to a gamma distribution with $k=n+1$ and $\theta=1/\sqrt{k}$. The time-to-peak in this case diverges as $n\to\infty$.
No comments:
Post a Comment