Sampling of a continuous function: Kronecker's or Dirac's delta? I've been reading some papers in signal proccesing and I'm very confused about the issue in the title of my question. Consider a continuous function of time $t$, $f(t)$, that I sample at uneven times $t_k$, where $k=1,2,...,N$. To me, it makes sense that the sampled function is:
$$f_s(t)=\sum_{k=1}^N\delta_{t,t_k}f(t),\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1)$$
where $\delta_{t,t_k}$ is Kronecker's delta (equals $1$ when $t=t_k$, zero elsewhere). However, in this paper, the author defines the sampled signal as:
$$f_s(t)=\frac{1}{N}\sum_{k=1}^{N}f(t)\delta(t-t_k),\ \ \ (2)$$
where $\delta(t-t_k)$ is Dirac's delta function and I really don't get why the $1/N$ appears here (the author claims that the sampling function is actually a weighted sum of delta functions
$$s(t)=C\frac{\sum_{k=1}^{N}w_k\delta(t-t_k)}{\sum_{k=1}^{N}w_k},$$
and here he choose $C=w_k=1$. I really didn't understand why). This last statement doesn't make much sense to me: the sampled signal would have infinite amplitude at $t=t_k$! 

Despite of all this, it is much easier to define the Fourier Transform of $f_s(t)$ in the second case (equation $(2)$), because it is just the convolution of the window function (the FT of the Dirac comb) and the FT of the continuous signal $f(t)$, while on equation $(1)$ the FT is a bit more complicated because we have an integer function (Kronecker's delta) multiplied by a continuous function ($f(t)$). Any highlights on this?
Modeling of the sampling process via multiplication of a continuous-time signal by a train of Dirac impulses is the most common interpretation in my experience. If you dig into it deeply enough, you'll find some disagreement about the mathematical precision of this approach*, but I wouldn't worry about it; it's just a convenient model for the process. There are no impulse generators inside your cell phone's ADC generating periodic bolts of lightning that multiply their analog inputs.

As you noted, you can't calculate the continuous-time Fourier transform of the Kronecker delta function, as its domain isn't continuous (it is limited to the integers). The Dirac delta function, in contrast, has a simple Fourier transform, and the effect of multiplying a signal by a train of Dirac impulses is easy to show due to its sifting property.

*: As an example, if you're going to be mathematically precise, you would say that the Dirac delta isn't a function at all, but a distribution instead. But at an engineering level, these issues are really just semantics.

Edit: I'll address the comment below. You gave your mental model of the sampling process as:

$$
f_s(t) = \sum_{k=1}^N \int_{t_k-\epsilon_k}^{t_k+\epsilon_k} f(t) \delta(t-t_k) dt.
$$

The problem with this interpretation is that the typical ideal sampling model does not have that integration built in. Instead, it's a pure multiplication of the input signal by a Dirac impulse train. If you look more closely at the equation you showed for $f_s(t)$, you'll see that the right hand side actually has no independent variable; $t$ is the dummy variable of integration. For any $\epsilon_k &gt; 0$ in the above, according to the Dirac impulse's sifting property, you would get: 

$$
f_s(t)=\sum_{k=1}^N f(t_k),
$$ 

which isn't correct. Instead, the model for the sampled signal is:

$$
f_s(t)=\sum_{k=-\infty}^{\infty} f(t) \delta(t-kT)
$$

Which is very similar to the above, except generalizing for an infinitely long impulse train along the time axis and assuming that the data is uniformly sampled at time instants $t_k = kT$. The Fourier transform of the resulting signal is:

$$
\begin{align}
F_s(\omega) &amp;= \int_{-\infty}^{\infty}f_s(t) e^{-j \omega t} dt \\
&amp;= \int_{-\infty}^{\infty} \sum_{k=-\infty}^{\infty} f(t) \delta(t-kT) e^{-j \omega t} dt \\
&amp;= \sum_{k=-\infty}^{\infty} \int_{-\infty}^{\infty} f(t) \delta(t-kT) e^{-j \omega t} dt \\
&amp;= \sum_{k=-\infty}^{\infty} f(kT) e^{-j \omega kT}
\end{align}
$$

If we define the discrete, sampled version of the signal $f(t)$ to be $x[n] = f(nT)$, then you are left with:

$$
F_s(\omega) = \sum_{n=-\infty}^{\infty} x[n] e^{-j \omega n}
$$

which is exactly the definition of the discrete-time Fourier transform.