Discrete Fourier series
The discrete Fourier series (DFS), also referred to as the discrete-time Fourier series (DTFS), is a mathematical tool in signal processing and mathematics that represents a periodic discrete-time signal as a finite linear combination of complex exponential functions with frequencies that are integer multiples of a fundamental frequency \omega_0 = 2\pi / N, where N is the period of the signal.[1][2] This representation decomposes the signal into its frequency components, analogous to the continuous-time Fourier series but adapted for sequences defined on integer indices rather than continuous time.[3] The core of the DFS consists of two equations: the analysis equation, which computes the Fourier coefficients X for k = 0, 1, \dots, N-1 as X = \sum_{n=0}^{N-1} x e^{-j 2\pi k n / N}, and the synthesis equation, which reconstructs the original periodic sequence x as x = \frac{1}{N} \sum_{k=0}^{N-1} X e^{j 2\pi k n / N}.[3][2] For real-valued signals, the coefficients satisfy X[N-k] = X^*, where ^* denotes the complex conjugate, ensuring symmetry in the frequency domain.[2] Unlike the continuous Fourier series, the DFS involves no convergence issues for finite-period signals, as the representation is exact due to the finite orthonormal basis of N harmonic sinusoids.[1] Key properties of the DFS include linearity, allowing superposition of signals; time-shifting, which shifts the coefficients by a phase factor; frequency-shifting, modulating the signal to alter its frequency content; and Parseval's theorem, stating that \sum_{n=0}^{N-1} |x|^2 = \frac{1}{N} \sum_{k=0}^{N-1} |X|^2, preserving energy between time and frequency domains.[2] The DFS relates closely to the discrete Fourier transform (DFT), which extends it to aperiodic finite-length sequences, and forms the basis for efficient algorithms like the fast Fourier transform (FFT) used in practical computations.[3] In engineering applications, it is fundamental for analyzing periodic signals in digital systems, enabling tasks such as spectral estimation, filtering, and system identification in fields like communications and audio processing.[3][1]Introduction and Background
Definition and Motivation
The discrete Fourier series (DFS), also known as the discrete-time Fourier series (DTFS), provides a representation for periodic discrete-time signals, expressing a sequence x with period N as a finite sum of harmonically related complex exponentials. Specifically, the synthesis equation is given by x = \sum_{k=0}^{N-1} c_k e^{j 2\pi k n / N}, where c_k are the DFS coefficients, and the summation is over one period, with the exponentials having fundamental frequency \omega_0 = 2\pi / N.[4][5] This formulation assumes familiarity with basic complex numbers, infinite sequences, and the concept of discrete-time signals, which are sequences of values indexed by integers, often arising from sampled continuous-time data.[4] The primary motivation for the DFS stems from the need to analyze periodic signals in the digital domain, particularly those obtained by sampling continuous-time periodic functions, such as in audio recordings or communication waveforms. By decomposing these signals into their frequency components, the DFS enables efficient frequency-domain processing, where operations like filtering become straightforward multiplications in the frequency space, leveraging the fact that complex exponentials are eigenfunctions of linear time-invariant systems.[4][5] In practical applications, this is crucial for digital signal processing tasks in audio analysis, where it helps identify tonal frequencies in music or speech, and in communications, where it supports modulation and spectrum estimation for efficient data transmission.[6][7] A distinguishing feature of the DFS is its inherent periodicity: both the time-domain sequence x and the frequency-domain coefficients c_k repeat every N samples, which contrasts with aperiodic discrete transforms that handle non-repeating signals without such repetition. This periodicity aligns naturally with the sampled nature of digital signals, ensuring that the representation captures the full harmonic content within one period while avoiding redundancy beyond N terms.[4][5]Historical Context
The discrete Fourier series emerged in the mid-20th century as digital computing advanced, building on the foundational continuous Fourier series developed by Joseph Fourier in his 1822 treatise Théorie analytique de la chaleur, which established the decomposition of periodic functions into trigonometric series.[8] This continuous framework provided the theoretical basis for extending Fourier analysis to discrete domains, particularly for periodic sequences in sampled signals. Early theoretical explorations of discrete versions predated widespread computation, serving as analytical tools in nascent signal processing.[9] Key milestones trace back to the 1940s during World War II, when research on sampled-data systems for radar and control applications spurred interest in discrete-time analysis. Pioneers like John R. Ragazzini and Lotfi A. Zadeh advanced sampled-data theory in a 1952 paper, introducing methods like the z-transform to model discrete systems, motivated by postwar radar needs and the limitations of analog processing.[10] These efforts laid groundwork for handling periodic discrete signals, though full formalization of the discrete Fourier series awaited the digital era. By the 1960s, recognition grew of its equivalence to the discrete Fourier transform (DFT) under periodic extension, as detailed in emerging digital signal processing literature, enabling spectral analysis of finite-length sequences as approximations of infinite periodic ones.[11] The 1965 publication of the Cooley-Tukey algorithm marked a computational breakthrough, providing an efficient method to calculate the DFT—and by extension, discrete Fourier series coefficients—reducing complexity from O(N²) to O(N log N) and facilitating practical implementation on early computers. This innovation, rooted in earlier 20th-century ideas but popularized post-1965, accelerated the shift from analog to digital spectral analysis.[9] Formalization in signal processing texts solidified the discrete Fourier series in the 1970s, notably through Alan V. Oppenheim and Ronald W. Schafer's 1975 book Digital Signal Processing, which systematically presented its synthesis and analysis for periodic discrete-time signals, integrating it with DFT computations and FFT efficiency.[12] This work influenced the broader transition to the digital era, enabling real-time applications in fields like communications and enabling scalable spectral processing on emerging hardware.[9]Mathematical Foundations
Synthesis and Analysis Equations
The synthesis and analysis equations provide the mathematical framework for representing a discrete-time periodic signal x with period N in the frequency domain using its Fourier coefficients X, and vice versa. These equations leverage the orthogonality of the complex exponential basis functions over one period to ensure exact reconstruction and decomposition. The synthesis equation reconstructs the time-domain signal from its frequency-domain coefficients as x = \frac{1}{N} \sum_{k=0}^{N-1} X \, e^{j 2\pi k n / N}, \quad n = 0, 1, \dots, N-1. The factor of $1/N in the synthesis equation normalizes the contribution of each frequency component, arising directly from the orthogonality relation of the basis functions \sum_{n=0}^{N-1} e^{j 2\pi (k - l) n / N} = N \delta_{k,l \mod N}.[11] The corresponding analysis equation extracts the frequency-domain coefficients from the time-domain signal: X = \sum_{n=0}^{N-1} x \, e^{-j 2\pi k n / N}, \quad k = 0, 1, \dots, N-1. To derive this, substitute the synthesis form into the left-hand side: multiplying both sides by e^{-j 2\pi l n / N} and summing over n yields \sum_n x e^{-j 2\pi l n / N} = \sum_k X \delta_{k,l \mod N} = X, confirming the analysis form via the same orthogonality.[11] Here, X denotes the frequency-domain representation, capturing the amplitude and phase of the k-th harmonic. The coefficients X are inherently periodic with period N, satisfying X[k + mN] = X for any integer m, due to the repeating nature of the exponential basis functions modulo N.[11] The finite period N in the discrete Fourier series formulation ensures exact reconstruction of the periodic signal via the synthesis equation without time-domain aliasing, as the summation over precisely one period aligns with the signal's periodicity, avoiding overlap from periodic replicas.[13]Coefficient Computation
The computation of the discrete Fourier series (DFS) coefficients from a given periodic sequence x with period N is performed using the analysis equation, which extracts the frequency-domain representation X for k = 0, 1, \dots, N-1: X = \sum_{n=0}^{N-1} x e^{-j 2\pi k n / N}. This formula arises from the orthogonality of the complex exponential basis functions, analogous to an inner product in the discrete space, where multiplying the synthesis equation by e^{-j 2\pi m n / N} and summing over one period yields the coefficient via the Kronecker delta property: the sum equals N when m = k and zero otherwise.[14] To compute the coefficients, one evaluates the finite summation over a single period of N samples, weighting each time-domain sample x by the complex exponential kernel; this discrete process directly avoids the infinite integrals required in the continuous-time Fourier series, where coefficients involve integration over the period instead of summation.[14][15] The resulting coefficients X are themselves periodic with period N, satisfying X[k + N] = X due to the periodicity of the exponential term in the summation.[14] For real-valued sequences x, the coefficients exhibit conjugate symmetry, where X[N - k] = X^* (with ^* denoting the complex conjugate), halving the independent values needed for representation.[16] Parseval's theorem relates the energy in the time and frequency domains, stating that \sum_{n=0}^{N-1} |x|^2 = \frac{1}{N} \sum_{k=0}^{N-1} |X|^2, preserving total energy across the transform pair.[14]Properties and Characteristics
Periodicity and Symmetry
The discrete Fourier series (DFS) is fundamentally tied to the periodicity of the underlying signal in the time domain. Specifically, the DFS represents a discrete-time signal \tilde{x} that is periodic with period N, satisfying \tilde{x}[n + mN] = \tilde{x} for any integer m. This periodicity is inherently enforced by the synthesis equation of the DFS, which reconstructs the signal as a sum of N complex exponentials over one fundamental period: \tilde{x} = \frac{1}{N} \sum_{k=0}^{N-1} \tilde{X} e^{j 2\pi k n / N}, \quad 0 \leq n < N. Any extension beyond one period repeats the signal periodically, ensuring the representation captures the repeating nature of the sequence without loss of information within the periodic framework.[17] In the frequency domain, the DFS coefficients \tilde{X} exhibit analogous periodicity with period N, such that \tilde{X}[k + mN] = \tilde{X} for integer m. This property arises directly from the analysis equation: \tilde{X} = \sum_{n=0}^{N-1} \tilde{x} e^{-j 2\pi k n / N}, \quad 0 \leq k < N, where shifting k by multiples of N yields identical results due to the exponential's periodicity. Consequently, only N unique coefficients are required to fully describe the spectrum, limiting the representation to discrete harmonics spaced by $2\pi / N radians per sample. This frequency-domain periodicity reflects the finite resolution imposed by the signal's period N.[17][11] Symmetry properties further simplify the DFS representation, particularly for real-valued signals. For a real periodic sequence \tilde{x}, the coefficients satisfy Hermitian symmetry: \tilde{X} = \tilde{X}^*[N - k], where * denotes the complex conjugate. This relation implies that the real part of \tilde{X} is even (symmetric around k = 0) and the imaginary part is odd (antisymmetric), reducing the number of independent coefficients by approximately half. In the context of even and odd extensions, a real even sequence—defined discretely as \tilde{x} = \tilde{x}[-n \mod N]—yields purely real coefficients, while a real odd sequence \tilde{x} = -\tilde{x}[-n \mod N] produces purely imaginary coefficients. These symmetries exploit the signal's structure to enhance computational efficiency and interpretability.[17][11][18] The discrete sampling inherent to the time-domain signal introduces potential aliasing in the frequency domain, where higher-frequency components could fold into the principal range [0, N-1]. However, within the periodic framework of the DFS, this aliasing is resolved by the analysis equation's summation over exactly one period, which aggregates all contributions from aliased replicas into the discrete harmonic coefficients without distortion for truly periodic signals. This encapsulation ensures the DFS provides an exact representation, distinguishing it from non-periodic transforms where aliasing may require additional mitigation.[11]Linearity and Convolution Theorem
The discrete Fourier series (DFS) exhibits linearity, meaning that if two periodic sequences x_1 and x_2, each with period N, have DFS coefficients X_1 and X_2, then the linear combination z = a x_1 + b x_2 (where a and b are constants) has DFS coefficients Z = a X_1 + b X_2. This property follows directly from the superposition principle inherent in the DFS representation, as the synthesis equation for z is the sum of scaled versions of the individual exponential components, and the analysis equation preserves this additivity due to the linearity of summation.[19][20] A related property is the time-shift theorem, which states that shifting a periodic sequence in time by m samples, yielding z = x[n - m], results in DFS coefficients multiplied by a phase factor: Z = X e^{-j 2\pi k m / N}. This arises because the shift modifies the phase of each harmonic exponential in the synthesis equation, introducing a linear phase term that factors out in the analysis sum; the periodicity ensures the shift remains within one period. This property highlights how temporal displacements manifest as frequency-dependent phase rotations in the DFS domain.[19][21] The convolution theorem for DFS asserts that the circular (or periodic) convolution of two periodic sequences x and y, defined as z = \sum_{l=0}^{N-1} x y[(n - l) \mod N], has DFS coefficients given by the pointwise product Z = X Y. This equivalence stems from the orthogonality of the exponential basis functions, where the convolution sum in time aligns with multiplication in the frequency domain upon applying the analysis equation; the inverse also holds, with time-domain multiplication corresponding to frequency-domain circular convolution. Unlike linear convolution for aperiodic signals, all DFS operations are inherently circular due to the enforced periodicity of the sequences, which wraps contributions around the period N and distinguishes DFS from non-periodic transforms.[19][20]Relations to Other Transforms
Comparison with Continuous Fourier Series
The discrete Fourier series (DFS) and the continuous Fourier series (CFS) share fundamental conceptual similarities in their approach to signal decomposition. Both represent periodic signals as sums of orthogonal complex exponential basis functions, enabling the analysis of frequency content through harmonic components.[18] In the CFS, a continuous-time periodic signal x(t) with period T is expressed as an infinite sum: x(t) = \sum_{k=-\infty}^{\infty} c_k e^{j 2\pi k t / T}, where the coefficients c_k are computed via integration: c_k = \frac{1}{T} \int_0^T x(t) e^{-j 2\pi k t / T} \, dt. [18] Similarly, the DFS decomposes a discrete-time periodic sequence x with period N into a finite sum of orthogonal discrete exponentials: x = \frac{1}{N} \sum_{k=0}^{N-1} X e^{j 2\pi k n / N}, with coefficients obtained through summation: X = \sum_{n=0}^{N-1} x e^{-j 2\pi k n / N}. [18] This parallelism underscores their roles as frequency-domain tools for periodic signals, differing primarily in the domain (continuous versus discrete) and the nature of the basis functions (infinite versus finite harmonics).[22] Key differences arise from the discrete nature of the DFS, which eliminates certain convergence challenges inherent to the CFS. The CFS synthesis involves an infinite series that may not converge uniformly for discontinuous signals, potentially leading to issues like the Gibbs phenomenon—persistent oscillations near discontinuities even as more terms are added. In contrast, the DFS uses a finite sum limited to N harmonics, inherently bandlimiting the representation to frequencies up to the Nyquist rate relative to the sampling, and thus avoids the Gibbs phenomenon entirely due to the absence of infinite partial sums. Additionally, the integral-based coefficient computation in the CFS requires exact continuous integration, while the DFS employs discrete summation, simplifying numerical evaluation but introducing quantization effects.[18] The DFS can be viewed as arising from sampling a continuous-time periodic signal, where its coefficients correspond to sampled versions of the CFS coefficients. Specifically, if a bandlimited continuous signal is sampled at N points over one period, the DFS coefficients X equal N times the CFS coefficients c_k, provided the signal's bandwidth does not exceed the sampling rate to prevent aliasing.[22] Undersampling leads to aliasing, where higher-frequency components fold into lower frequencies, distorting the discrete spectrum—a phenomenon absent in the purely continuous CFS but critical in bridging the two domains.[22] This sampling relationship highlights the DFS as an approximation of the CFS for digital processing, preserving essential frequency information under proper conditions.[22]Connection to Discrete Fourier Transform
The discrete Fourier series (DFS) and the discrete Fourier transform (DFT) are closely related mathematical tools for analyzing discrete signals, with the DFS specifically tailored for periodic sequences and the DFT providing a generalization for finite-length sequences. For a periodic signal x with period N, the DFS coefficients X are computed over one full period, and this computation is identical to applying the N-point DFT to that single period of the signal. This equivalence arises because both operations project the signal onto the same set of orthogonal basis functions, differing only by a scaling factor of N in the inverse transform.[11][17] In practice, the DFT of a finite sequence x of length M (where M \leq N) implicitly treats the sequence as one period of a periodic signal by periodically extending it with repetitions every N samples, effectively converting the finite data into a periodic form suitable for DFS analysis. When the sequence length exactly matches the transform size (M = N), the DFT yields the precise DFS coefficients without approximation. The defining formula for this case is the N-point DFT: X = \sum_{n=0}^{N-1} x e^{-j 2\pi k n / N}, \quad k = 0, 1, \dots, N-1, with the inverse given by x = \frac{1}{N} \sum_{k=0}^{N-1} X e^{j 2\pi k n / N}, \quad n = 0, 1, \dots, N-1. This exact match holds because the finite sequence aligns perfectly with the periodicity assumption of the DFS.[11][17] For aperiodic finite sequences, the DFT can approximate a DFS representation by zero-padding the signal to length N (where N > M), which interpolates the frequency-domain coefficients but introduces periodicity artifacts unless N is sufficiently large. This extension allows the DFT to handle non-periodic data by enforcing artificial periodicity, bridging the gap between finite-duration analysis and periodic decomposition. Overall, the DFS emphasizes the periodic nature of signals for exact harmonic representation, whereas the DFT offers flexibility for processing arbitrary finite data, often serving as the computational framework for DFS in digital signal processing. Both discrete tools originate from sampling the continuous Fourier series, adapting its principles to discrete domains.[11][17][23]Computation and Implementation
Direct Calculation Methods
Direct calculation methods for discrete Fourier series coefficients involve straightforward numerical evaluation of the analysis equation through summation, without exploiting any symmetries or recursive structures. This approach computes each coefficient X by directly summing the contributions from all N time-domain samples, typically implemented as nested loops over the indices k (frequencies) and n (time samples). For a periodic discrete-time signal x with period N, the computation requires performing N complex multiplications and additions for each of the N coefficients, resulting in an overall time complexity of O(N^2).[18] This method can be viewed as a matrix-vector multiplication, where the discrete Fourier series coefficients vector \mathbf{X} is obtained by \mathbf{X} = \mathbf{W} \mathbf{x}, with \mathbf{x} being the input signal vector and \mathbf{W} the N \times N Vandermonde matrix whose entries are W_{k,n} = e^{-j 2 \pi k n / N}. The direct evaluation thus corresponds to explicit multiplication of this matrix by the input vector, which is computationally intensive but conceptually simple.[24] The following pseudocode illustrates a basic implementation in a programming language like Python or MATLAB:This brute-force procedure is suitable for small values of N, such as N < 100, where the quadratic complexity remains manageable on modern hardware and the simplicity aids in educational or prototyping contexts.[18] However, even in direct computation, numerical errors arise due to round-off in floating-point arithmetic, particularly from the repeated complex multiplications involving exponentials, which can accumulate and degrade accuracy as N increases. While the method yields exact results in arbitrary-precision or exact arithmetic environments, practical implementations in standard floating-point formats (e.g., IEEE 754 double precision) introduce these errors, making the direct approach less stable compared to optimized alternatives for moderate to large N.[25]for k in 0 to N-1: X[k] = 0 for n in 0 to N-1: X[k] += x[n] * [exp](/page/Exp)(-1j * 2 * pi * n * k / N)for k in 0 to N-1: X[k] = 0 for n in 0 to N-1: X[k] += x[n] * [exp](/page/Exp)(-1j * 2 * pi * n * k / N)
Efficient Algorithms Using FFT
The Fast Fourier Transform (FFT) provides an efficient method for computing the Discrete Fourier Series (DFS) coefficients, significantly reducing the computational burden compared to direct summation methods. The seminal Cooley-Tukey radix-2 algorithm decomposes the N-point discrete Fourier transform into smaller sub-transforms by exploiting the symmetry and periodicity of the exponential basis functions, achieving a time complexity of O(N log N) operations for sequences where N is a power of 2, as opposed to the O(N²) complexity of naive evaluation.[26][27] In the context of DFS, which represents a periodic discrete-time signal over one period of N samples, the coefficients can be computed directly by applying the FFT to these N samples, yielding the frequency-domain representation equivalent to the DFS analysis formula.[28] This approach leverages the fact that the DFS and discrete Fourier transform are mathematically equivalent for finite-length periodic sequences.[29] The radix-2 Cooley-Tukey FFT requires the sequence length N to be a power of 2 for optimal efficiency, but variants such as Bluestein's algorithm extend fast computation to arbitrary N, including prime lengths, by reformulating the transform as a convolution computable via larger power-of-2 FFTs, maintaining O(N log N) complexity.[30] This enables practical DFS computation without unnecessary zero-padding, which would otherwise inflate the transform size. For implementation, widely-used libraries such as FFTW and NumPy's fft module incorporate optimized versions of these algorithms, including Cooley-Tukey and Bluestein variants, to perform real-time DFS computations on large datasets across various hardware platforms.[31][32] These libraries provide substantial speedups—often orders of magnitude—over direct calculation methods, making FFT-based DFS indispensable for applications requiring high performance.[28]Applications and Examples
Signal Processing Uses
In digital signal processing, the discrete Fourier series (DFS) plays a fundamental role in spectral analysis by decomposing periodic discrete-time signals into their constituent frequency components. For signals such as audio tones or mechanical vibrations, the DFS coefficients provide the amplitudes and phases of harmonic frequencies, allowing engineers to identify dominant spectral lines and harmonic content essential for tasks like noise detection or vibration monitoring. This frequency-domain representation enables precise characterization of periodic phenomena, where the magnitude spectrum highlights energy distribution across discrete frequencies.[11] The convolution theorem for DFS further extends its utility to filtering applications, where multiplication of DFS coefficients in the frequency domain implements circular convolution in the time domain. This property supports the design of frequency-selective filters for periodic signals, such as low-pass or band-stop filters, by attenuating unwanted spectral components while preserving others, thereby enabling efficient removal of interference in applications like audio equalization or seismic data processing.[11] A notable application rooted in DFS principles is audio compression via the modified discrete cosine transform (MDCT) in the MP3 standard. The MDCT, a real-valued transform derived from the discrete Fourier transform through shifted variants, represents audio blocks in a critically sampled frequency domain, facilitating perceptual coding by quantizing less audible spectral components and achieving significant bitrate reduction without substantial quality loss.[33] In wireless communications, DFS basis functions underpin orthogonal frequency-division multiplexing (OFDM), where multiple orthogonal subcarriers—essentially the complex exponentials of the DFS—carry parallel data streams. This modulation scheme, implemented via the inverse discrete Fourier transform, combats multipath fading by converting frequency-selective channels into flat-fading subchannels, enhancing data rates and reliability in standards like Wi-Fi and 4G LTE.[34]Numerical Examples and Illustrations
To illustrate the computation of discrete Fourier series (DFS) coefficients, consider a simple periodic discrete-time signal with period N=4, defined as x = \cos\left(\frac{2\pi n}{4}\right) for n = 0, 1, 2, 3. This yields the sequence values x{{grok:render&&&type=render_inline_citation&&&citation_id=0&&&citation_type=wikipedia}} = 1, x{{grok:render&&&type=render_inline_citation&&&citation_id=1&&&citation_type=wikipedia}} = 0, x{{grok:render&&&type=render_inline_citation&&&citation_id=2&&&citation_type=wikipedia}} = -1, x{{grok:render&&&type=render_inline_citation&&&citation_id=3&&&citation_type=wikipedia}} = 0.[35] The DFS coefficients X are obtained via the analysis equation: X = \sum_{n=0}^{N-1} x e^{-j \frac{2\pi k n}{N}}. For this signal:- X{{grok:render&&&type=render_inline_citation&&&citation_id=0&&&citation_type=wikipedia}} = 1 \cdot 1 + 0 \cdot 1 + (-1) \cdot 1 + 0 \cdot 1 = 0
- X{{grok:render&&&type=render_inline_citation&&&citation_id=1&&&citation_type=wikipedia}} = 1 \cdot 1 + 0 \cdot e^{-j \pi/2} + (-1) \cdot e^{-j \pi} + 0 \cdot e^{-j 3\pi/2} = 1 + 0 + 1 + 0 = 2
- X{{grok:render&&&type=render_inline_citation&&&citation_id=2&&&citation_type=wikipedia}} = 1 \cdot 1 + 0 \cdot e^{-j \pi} + (-1) \cdot e^{-j 2\pi} + 0 \cdot e^{-j 3\pi} = 1 + 0 - 1 + 0 = 0
- X{{grok:render&&&type=render_inline_citation&&&citation_id=3&&&citation_type=wikipedia}} = 1 \cdot 1 + 0 \cdot e^{-j 3\pi/2} + (-1) \cdot e^{-j 3\pi} + 0 \cdot e^{-j 9\pi/2} = 1 + 0 + 1 + 0 = 2
| k | X |
|---|---|
| 0 | 0 |
| 1 | 2 |
| 2 | 0 |
| 3 | 2 |
- X{{grok:render&&&type=render_inline_citation&&&citation_id=0&&&citation_type=wikipedia}} = 1 + 1 + 0 + 0 = 2
- X{{grok:render&&&type=render_inline_citation&&&citation_id=1&&&citation_type=wikipedia}} = 1 \cdot 1 + 1 \cdot e^{-j \pi/2} + 0 + 0 = 1 - j
- X{{grok:render&&&type=render_inline_citation&&&citation_id=2&&&citation_type=wikipedia}} = 1 \cdot 1 + 1 \cdot e^{-j \pi} + 0 + 0 = 1 - 1 = 0
- X{{grok:render&&&type=render_inline_citation&&&citation_id=3&&&citation_type=wikipedia}} = 1 \cdot 1 + 1 \cdot e^{-j 3\pi/2} + 0 + 0 = 1 + j
| k | X |
|---|---|
| 0 | 2 |
| 1 | $1 - j |
| 2 | 0 |
| 3 | $1 + j |