Alternative expression pulse train
In the module sampling and reconstruction we have seen that the continuous time signal x(t)x(t) can be represented by samples x[n]x[n] when these samples are obtained by using an ideal C-to-D convertor which runs at a sampling frequency fsfs which is larger than 2 times the maximum frequency fmaxfmax of the continuous time signal x(t)x(t). Here we will assume that this is the case, so there is no aliasing.
In the development of the FTD we will use the model of an ideal C-to-D convertor as depicted in Fig. 1. In this model we first represent the discrete-time samples x[n]x[n] by the continuous time signal xs(t)xs(t), which is obtained by multiplying the continuous time signal x(t)x(t) by a continuous-time signal s(t)s(t) which is an infinite train of continuous-time delta pulses with an inter-pulse distance of Ts=1/fsTs=1/fs seconds.
Although this intermediate continuous-time signal xs(t)xs(t) can never be realized in practice, it is very useful in the following conceptual derivation of the FTD.
In the module Fourier Transform for Continuous-time signals we have derived the following alternative expression for this infinite train of continuous-time delta pulses s(t)s(t) as an infinite sum of phasors: s(t)=∞∑n=−∞δ(t−n⋅Ts)ˆ=1Ts∞∑k=−∞ej2πTskts(t)=∞∑n=−∞δ(t−n⋅Ts)^=1Ts∞∑k=−∞ej2πTskt The reasoning behind this equality can be shown by Fig. 2. The figure at the left hand side indicates an example of the phasor terms of the summation for any continuous time variable tt which is not equal to n×Tsn×Ts.
In such cases, the phasors terms in the summation, which are unit length vectors, can have any angle. There will be an infinite amount of these unit length vectors with positive imaginary part, colored in blue, and with negative imaginary part, colored in red. As a result of this infinite amount of vectors, for each blue vector there will be a red vector in opposite direction which cancel each other and the result is that the summation of all these phasors will add up to zero. Only for values of the continuous-time variable t≠n×Tst≠n×Ts, all phasor terms of the summation are equal to one and are all on top of each other. The result of the summation for this case is infinite, which is indicated at the right hand side of the figure. Thus the infinite summation of phasors can indeed be represented by an infinite train of delta pulses with an inter distance of TsTs seconds, as depicted in the equation above.
In the following subsection we will use this equality in the derivation of the FTD.
FTD concept
In this subsection we will derive the FTD by using frequency domain descriptions of signals x(t),xs(t)x(t),xs(t) and x[n]x[n], as indicated in Fig. 1 of the mathematical model of the ideal C-to-D convertor.
Assume we are given a continuous time signal x(t)x(t) and its FTC X(ω)X(ω). The maximum radian frequency of x(t)x(t) is ωmωm, as depicted in the top figure Fig. 3.
As a result of the definition of an impulse train we can write the intermediate continuous-time signal xs(t)xs(t) as follows: xs(t)=x(t)⋅1Ts∞∑k=−∞ej2πTskt=1Ts∞∑k=−∞x(t)⋅ejωsktwithωs=2πTsxs(t)=x(t)⋅1Ts∞∑k=−∞ej2πTskt=1Ts∞∑k=−∞x(t)⋅ejωsktwithωs=2πTs This equation allows us to apply the frequency shift property to obtain the FTC of xs(t)xs(t): Xs(ω)=1Ts⋅∞∑k=−∞X(ω−k⋅ωs),Xs(ω)=1Ts⋅∞∑k=−∞X(ω−k⋅ωs), which is a periodic function of ωω with period ωs=2πTsωs=2πTs. As depicted in the middle figure of Fig. 3, Xs(ω)Xs(ω) (in green) gives the spectral representation of xs(t)xs(t) in terms of shifted copies of X(ω)X(ω), the FTC of the continuous-time signal x(t)x(t). On the other hand if we express xs(t)xs(t) as an infinite train of continuous-time delta pulses xs(t)=x(t)⋅∑∞n=−∞δ(t−n⋅Ts)xs(t)=x(t)⋅∑∞n=−∞δ(t−n⋅Ts) and using the delay property of the FTC, it follows that an equivalent relation for Xs(ω)Xs(ω) can be given by the following equation: Xs(ω)=∞∑n=−∞x(t)|t=nTse−jωnTsXs(ω)=∞∑n=−∞x(t)∣∣∣t=nTse−jωnTs which expresses Xs(ω)Xs(ω) in terms of the sampled values x(t)∣t=nTsx(t)∣t=nTs of the continuous-time signal x(t)x(t). Because of the exponential functions, as well as the description of Xs(ω)Xs(ω) this form of Xs(ω)Xs(ω) is periodic in ωω with period ωs=2πTsωs=2πTs. Furthermore by representing the sampled values x(t)∣t=nTsx(t)∣t=nTs by the samples x[n]x[n] and using the relative frequency θ=ω⋅Tsθ=ω⋅Ts we can rewrite the above equation as: X(ejθ)=∞∑n=−∞x[n]e−jnθX(ejθ)=∞∑n=−∞x[n]e−jnθ in which the periodicity of this function is denoted by expressing the continuous function variable θθ as an exponent. This equation transforms the discrete-time samples x[n]x[n] into a continuous periodic relative frequency domain function X(ejθ)X(ejθ) and as such it represents the Fourier Transform of Discrete-time signals (FTD). From the equivalence of the spectra Xs(ω)Xs(ω) and X(ejθ)X(ejθ) it follows that the FTD can be viewed in terms of shifted copies of X(ω)X(ω). Finally, as depicted in the lower figure of Fig. 3, it follows from θ=ω⋅Tsθ=ω⋅Ts that the period of X(ejθ)X(ejθ) is equal to 2π2π. Only one period is necessary to represent X(ejθ)X(ejθ), which is called the Fundamental Interval (FI). In most cases we will use for the FI the period in between −π−π and +π+π.
Inverse FTD
By using the fact that X(ejθ)X(ejθ) is a periodic function we can use the Fourier Series (FS) representation to find inverse of the FTD expression. For this reason we recall the FS definition which implies that any periodic function xp(t)xp(t), with period T0=1/F0T0=1/F0, can be represented as an infinite sum of weighted phasors while in the inverse case the weights αkαk can be evaluated by an integral over one period, as denoted in the following equations: xp(t)=∞∑k=−∞αkej2πkF0t∘−∘αk=1T0∫T0/2−T0/2xp(t)e−j2πkF0tdtxp(t)=∞∑k=−∞αkej2πkF0t∘−∘αk=1T0∫T0/2−T0/2xp(t)e−j2πkF0tdt
In our case we change time and frequency dimensions of the FS equations. Furthermore the periodic function is X(ejθ)X(ejθ) with period 2π2π and the ‘weights’ are represented by the samples x[n]x[n]. Now by using the FS equation for the ‘weights’ x[n]x[n] this results in the following expression for the Inverse FTD: x[n]=12π∫π−πX(ejθ)ejnθdθx[n]=12π∫π−πX(ejθ)ejnθdθ Given the periodic expression X(ejθ)X(ejθ) as a function of the continuous relative frequency θθ this equation defines the way to compute the discrete-time samples x[n]x[n].
Differences FTC and FTD
Both FTC and FTD have similarities but are different. That is why the most important aspects of both transformations are summarized here.
Based on absolute frequency ff in Hertz, the FTC and IFTC equations are: FTC:X(f)=∫∞−∞x(t)e−j2πftdt∘−∘x(t)=∫∞−∞X(f)ej2πftdfFTC:X(f)=∫∞−∞x(t)e−j2πftdt∘−∘x(t)=∫∞−∞X(f)ej2πftdf which alternatively can be expressed based on the absolute radian frequency ωω in radians per second as follows: FTC:X(ω)=∫∞−∞x(t)e−jωtdt∘−∘x(t)=12π∫∞−∞X(ω)ejωtdωFTC:X(ω)=∫∞−∞x(t)e−jωtdt∘−∘x(t)=12π∫∞−∞X(ω)ejωtdω
Important aspects of these transforms are:
- The continuous time signal x(t)x(t) is not periodic,
- The frequency representation X(ω)X(ω) is a continuous function and not periodic,
- Both FTC and IFTC are evaluated by an integral,
- Finally we will use these notations for the FTC and IFTC: X(ω)=F{x(t)} and x(t)=F−1{X(ω)}
The FTD and IFTD equations are as follows: FTD:X(ejθ)=∞∑n=−∞x[n]e−jnθ∘−∘x[n]=12π∫π−πX(ejθ)ejnθdθ Important aspects of these transforms, which differ from the FTC and IFTC equations, are:
- The signal samples x[n] are discrete in time and not periodic,
- The normalized frequency θ=ω⋅Ts,
- The Fundamental Interval is usually chosen as: −π≤θ<π,
- The frequency representation X(ejθ) is a continuous function and periodic,
- The FTD is evaluated over an infinite summation of samples x[n],
- The IFTD is evaluated by an integral over one period of the frequency representation X(ejθ).
- Finally, to make a clear distinction we will use the following notations for the FTD and IFTD: X(ejθ)=F{x[n]} and x[n]=F−1{X(ejθ)}