Skip to main content

Chapter 2 · 9 hours

Fourier Series

IOE past exam questions

Past questions and answers

46 questions set from this chapter, 13 of them more than once. Most asked first.

  • Asked 9 times
  • 2082 Chaitra · 5 marks
  • 2082 Kartik · 6 marks
  • 2083 Bhadra (new course) · 4 marks
  • 2082 Bhadra (new course) · 4 marks
  • 2075 Baisakh · 5 marks
  • 2074 Bhadra · 5 marks
  • 2072 Asoj · 5 marks
  • 2071 Bhadra · 7 marks
  • 2069 Bhadra · 8 marks

Derive the exponential form of the Fourier series of a continuous-time periodic signal x(t) (synthesis and analysis equations).

Answer

A continuous-time periodic signal x(t)x(t) with fundamental period TT and fundamental frequency ω0=2π/T\omega_0 = 2\pi/T can be written as a weighted sum of harmonically related complex exponentials. This is the exponential Fourier series.

Synthesis equation

The set ϕk(t)=ejkω0t\phi_k(t) = e^{jk\omega_0 t}, k=0,±1,±2,…k = 0, \pm1, \pm2, \dots, contains signals that are all periodic with period TT (the kkth one has frequency kω0k\omega_0). A linear combination of them is also periodic with period TT:

x(t)=∑k=−∞∞ak ejkω0t=∑k=−∞∞ak ejk(2π/T)tx(t) = \sum_{k=-\infty}^{\infty} a_k\, e^{jk\omega_0 t} = \sum_{k=-\infty}^{\infty} a_k\, e^{jk(2\pi/T)t}

This is the synthesis equation. a0a_0 is the dc term, k=±1k = \pm1 are the fundamental components, and k=±Nk = \pm N are the NNth harmonics.

Analysis equation

To find aka_k, multiply both sides by e−jnω0te^{-jn\omega_0 t} and integrate over one period TT:

∫0Tx(t)e−jnω0t dt=∑k=−∞∞ak∫0Tej(k−n)ω0t dt\int_0^T x(t)e^{-jn\omega_0 t}\,dt = \sum_{k=-\infty}^{\infty} a_k \int_0^T e^{j(k-n)\omega_0 t}\,dt

Evaluate the integral on the right (orthogonality):

∫0Tej(k−n)ω0tdt=∫0Tcos⁡((k−n)ω0t) dt+j∫0Tsin⁡((k−n)ω0t) dt\int_0^T e^{j(k-n)\omega_0 t}dt = \int_0^T \cos((k-n)\omega_0 t)\,dt + j\int_0^T \sin((k-n)\omega_0 t)\,dt
  • For k≠nk \ne n, both integrals cover a whole number (k−n)(k-n) of periods, so they are 0.
  • For k=nk = n, the integrand is 1 and the integral is TT.
∫0Tej(k−n)ω0tdt={T,k=n0,k≠n\int_0^T e^{j(k-n)\omega_0 t}dt = \begin{cases} T, & k = n \\ 0, & k \ne n \end{cases}

So only the k=nk = n term survives:

∫0Tx(t)e−jnω0tdt=anT⇒an=1T∫0Tx(t)e−jnω0tdt\int_0^T x(t)e^{-jn\omega_0 t}dt = a_n T \quad\Rightarrow\quad a_n = \frac{1}{T}\int_0^T x(t)e^{-jn\omega_0 t}dt

Since any interval of length TT gives the same result, the analysis equation is

ak=1T∫Tx(t) e−jkω0t dta_k = \frac{1}{T}\int_{T} x(t)\,e^{-jk\omega_0 t}\,dt

The Fourier series pair

x(t)⟷FSakx(t) \overset{FS}{\longleftrightarrow} a_k
EquationFormula
Synthesisx(t)=∑kakejkω0tx(t) = \sum_k a_k e^{jk\omega_0 t}
Analysisak=1T∫Tx(t)e−jkω0tdta_k = \frac{1}{T}\int_T x(t)e^{-jk\omega_0 t}dt
DC valuea0=1T∫Tx(t) dta_0 = \frac{1}{T}\int_T x(t)\,dt (average value)

Notes:

  • The aka_k are the Fourier series coefficients or spectral coefficients; they are complex in general, ak=∣ak∣ej∠aka_k = |a_k|e^{j\angle a_k}, giving the magnitude and phase spectra.
  • For real x(t)x(t), a−k=ak∗a_{-k} = a_k^*, so ∣ak∣|a_k| is even and ∠ak\angle a_k is odd in kk.

Example: x(t)=cos⁡ω0t=12ejω0t+12e−jω0tx(t) = \cos\omega_0 t = \frac{1}{2}e^{j\omega_0 t} + \frac{1}{2}e^{-j\omega_0 t} gives a1=a−1=1/2a_1 = a_{-1} = 1/2 and all other ak=0a_k = 0.

  • Asked 8 times
  • 2082 Kartik · 3 marks
  • 2079 Jestha · 6 marks
  • 2078 Baisakh · 6 marks
  • 2076 Bhadra · 5 marks
  • 2073 Bhadra · 5 marks
  • 2072 Magh · 2+4 marks
  • 2072 Asoj · 5 marks
  • 2083 Bhadra (new course) · 4 marks

State and prove Parseval's relation for Discrete Time Fourier Series (discrete time periodic signals).

Answer

Parseval's relation for the DTFS: the average power of a discrete-time periodic signal over one period equals the sum of the squared magnitudes of its Fourier series coefficients:

1N∑n=⟨N⟩∣x[n]∣2=∑k=⟨N⟩∣ak∣2\frac{1}{N}\sum_{n=\langle N\rangle}|x[n]|^2 = \sum_{k=\langle N\rangle}|a_k|^2

Here x[n]x[n] has period NN, ω0=2π/N\omega_0 = 2\pi/N, and

x[n]=∑k=⟨N⟩akejkω0n,ak=1N∑n=⟨N⟩x[n]e−jkω0nx[n] = \sum_{k=\langle N\rangle} a_k e^{jk\omega_0 n}, \qquad a_k = \frac{1}{N}\sum_{n=\langle N\rangle}x[n]e^{-jk\omega_0 n}

Proof

Start from the left side and write ∣x[n]∣2=x[n] x∗[n]|x[n]|^2 = x[n]\,x^*[n]:

1N∑n=⟨N⟩∣x[n]∣2=1N∑n=⟨N⟩x[n] x∗[n]\frac{1}{N}\sum_{n=\langle N\rangle}|x[n]|^2 = \frac{1}{N}\sum_{n=\langle N\rangle}x[n]\,x^*[n]

Replace x∗[n]x^*[n] by the conjugate of the synthesis equation, x∗[n]=∑k=⟨N⟩ak∗e−jkω0nx^*[n] = \sum_{k=\langle N\rangle}a_k^* e^{-jk\omega_0 n}:

1N∑n∣x[n]∣2=1N∑n=⟨N⟩x[n]∑k=⟨N⟩ak∗e−jkω0n=∑k=⟨N⟩ak∗[1N∑n=⟨N⟩x[n]e−jkω0n]=∑k=⟨N⟩ak∗ ak=∑k=⟨N⟩∣ak∣2\begin{aligned} \frac{1}{N}\sum_{n}|x[n]|^2 &= \frac{1}{N}\sum_{n=\langle N\rangle}x[n]\sum_{k=\langle N\rangle}a_k^* e^{-jk\omega_0 n} \\ &= \sum_{k=\langle N\rangle}a_k^*\left[\frac{1}{N}\sum_{n=\langle N\rangle}x[n]e^{-jk\omega_0 n}\right] \\ &= \sum_{k=\langle N\rangle}a_k^*\, a_k \\ &= \sum_{k=\langle N\rangle}|a_k|^2 \end{aligned}

The bracket is exactly the analysis equation for aka_k, which completes the proof.

Meaning

  • ∣ak∣2|a_k|^2 is the average power in the kkth harmonic akejkω0na_k e^{jk\omega_0 n}, because ∣ejkω0n∣=1|e^{jk\omega_0 n}| = 1.
  • So total average power = sum of powers of the NN harmonic components; power can be computed in either the time or the frequency domain.
  • Only NN distinct coefficients exist (ak+N=aka_{k+N} = a_k), so the sum is finite.

Example: x[n]=cos⁡(πn/2)x[n] = \cos(\pi n/2), N=4N = 4: a1=a−1=a3=1/2a_1 = a_{-1} = a_3 = 1/2. Time domain: samples 1,0,−1,01, 0, -1, 0 give 14(1+0+1+0)=0.5\frac{1}{4}(1 + 0 + 1 + 0) = 0.5. Frequency domain: ∣a1∣2+∣a−1∣2=0.25+0.25=0.5|a_1|^2 + |a_{-1}|^2 = 0.25 + 0.25 = 0.5. Both agree.

  • Asked 6 times
  • 2077 Chaitra · 4 marks
  • 2076 Baisakh · 6 marks
  • 2075 Bhadra · 6+4 marks
  • 2072 Magh · 5 marks
  • 2071 Magh · 7 marks
  • 2083 Baisakh (new course) · 4 marks

Derive the exponential form of the Fourier series (synthesis and analysis equations) for a discrete time periodic signal x[n] with period N.

Answer

A discrete-time signal x[n]x[n] is periodic with period NN if x[n+N]=x[n]x[n+N] = x[n]; its fundamental frequency is ω0=2π/N\omega_0 = 2\pi/N. The DTFS writes it as a sum of harmonically related complex exponentials ϕk[n]=ejkω0n\phi_k[n] = e^{jk\omega_0 n}.

Only NN distinct exponentials

ϕk+N[n]=ej(k+N)(2π/N)n=ejk(2π/N)n ej2πn=ϕk[n]\phi_{k+N}[n] = e^{j(k+N)(2\pi/N)n} = e^{jk(2\pi/N)n}\, e^{j2\pi n} = \phi_k[n]

since ej2πn=1e^{j2\pi n} = 1 for integer nn. So only NN different exponentials exist, e.g. k=0,1,…,N−1k = 0, 1, \dots, N-1. The sum therefore runs over any NN consecutive values of kk, written k=⟨N⟩k = \langle N\rangle.

Synthesis equation

x[n]=∑k=⟨N⟩ak ejkω0n=∑k=⟨N⟩ak ejk(2π/N)nx[n] = \sum_{k=\langle N\rangle} a_k\, e^{jk\omega_0 n} = \sum_{k=\langle N\rangle} a_k\, e^{jk(2\pi/N)n}

Orthogonality result

For integer rr, using the finite geometric series:

∑n=0N−1ejr(2π/N)n={N,r=0,±N,±2N,…1−ej2πr1−ej2πr/N=0,otherwise\sum_{n=0}^{N-1} e^{jr(2\pi/N)n} = \begin{cases} N, & r = 0, \pm N, \pm 2N, \dots \\ \dfrac{1 - e^{j2\pi r}}{1 - e^{j2\pi r/N}} = 0, & \text{otherwise} \end{cases}

Analysis equation

Multiply the synthesis equation by e−jr(2π/N)ne^{-jr(2\pi/N)n} and sum over one period:

∑n=⟨N⟩x[n]e−jr(2π/N)n=∑n=⟨N⟩∑k=⟨N⟩akej(k−r)(2π/N)n=∑k=⟨N⟩ak∑n=⟨N⟩ej(k−r)(2π/N)n\begin{aligned} \sum_{n=\langle N\rangle}x[n]e^{-jr(2\pi/N)n} &= \sum_{n=\langle N\rangle}\sum_{k=\langle N\rangle}a_k e^{j(k-r)(2\pi/N)n} \\ &= \sum_{k=\langle N\rangle}a_k \sum_{n=\langle N\rangle}e^{j(k-r)(2\pi/N)n} \end{aligned}

The inner sum is NN when k=rk = r (within one period of kk) and 0 otherwise, so

∑n=⟨N⟩x[n]e−jr(2π/N)n=Nar\sum_{n=\langle N\rangle}x[n]e^{-jr(2\pi/N)n} = N a_r ak=1N∑n=⟨N⟩x[n] e−jk(2π/N)na_k = \frac{1}{N}\sum_{n=\langle N\rangle}x[n]\,e^{-jk(2\pi/N)n}

The DTFS pair

EquationFormula
Synthesisx[n]=∑k=⟨N⟩akejk(2π/N)nx[n] = \sum_{k=\langle N\rangle}a_k e^{jk(2\pi/N)n}
Analysisak=1N∑n=⟨N⟩x[n]e−jk(2π/N)na_k = \frac{1}{N}\sum_{n=\langle N\rangle}x[n]e^{-jk(2\pi/N)n}

Key points:

  • Both sums are finite (NN terms), so there is no convergence problem, unlike the CT Fourier series.
  • The coefficients are periodic: ak+N=aka_{k+N} = a_k.
  • a0=1N∑x[n]a_0 = \frac{1}{N}\sum x[n] is the average (dc) value.

Example: x[n]=sin⁡(2πn/5)x[n] = \sin(2\pi n/5), N=5N = 5: x[n]=12jej(2π/5)n−12je−j(2π/5)nx[n] = \frac{1}{2j}e^{j(2\pi/5)n} - \frac{1}{2j}e^{-j(2\pi/5)n}, so a1=12ja_1 = \frac{1}{2j}, a−1=−12ja_{-1} = -\frac{1}{2j}, and these repeat every 5 values of kk.

  • Asked 4 times
  • 2082 Chaitra · 4 marks
  • 2078 Chaitra · 2+3 marks
  • 2076 Baisakh · 2+4 marks
  • 2071 Bhadra · 6 marks

What is Parseval's relation? State and prove Parseval's relation for continuous-time periodic signals.

Answer

Parseval's relation states that the average power of a periodic signal, computed in the time domain over one period, equals the sum of the powers of all its harmonic components (the squared magnitudes of its Fourier series coefficients). Power is the same whether measured in time or frequency.

Statement (continuous-time periodic signals)

If x(t)x(t) has period TT and x(t)⟷FSakx(t) \overset{FS}{\longleftrightarrow} a_k, then

1T∫T∣x(t)∣2 dt=∑k=−∞∞∣ak∣2\frac{1}{T}\int_{T}|x(t)|^2\,dt = \sum_{k=-\infty}^{\infty}|a_k|^2

Proof

With ω0=2π/T\omega_0 = 2\pi/T, x(t)=∑kakejkω0tx(t) = \sum_k a_k e^{jk\omega_0 t} and ak=1T∫Tx(t)e−jkω0tdta_k = \frac{1}{T}\int_T x(t)e^{-jk\omega_0 t}dt.

1T∫T∣x(t)∣2dt=1T∫Tx(t) x∗(t) dt=1T∫Tx(t)[∑k=−∞∞ak∗e−jkω0t]dt=∑k=−∞∞ak∗[1T∫Tx(t)e−jkω0tdt]=∑k=−∞∞ak∗ ak=∑k=−∞∞∣ak∣2\begin{aligned} \frac{1}{T}\int_T |x(t)|^2dt &= \frac{1}{T}\int_T x(t)\,x^*(t)\,dt \\ &= \frac{1}{T}\int_T x(t)\left[\sum_{k=-\infty}^{\infty}a_k^* e^{-jk\omega_0 t}\right]dt \\ &= \sum_{k=-\infty}^{\infty}a_k^*\left[\frac{1}{T}\int_T x(t)e^{-jk\omega_0 t}dt\right] \\ &= \sum_{k=-\infty}^{\infty}a_k^*\,a_k = \sum_{k=-\infty}^{\infty}|a_k|^2 \end{aligned}

(the order of sum and integral is swapped, and the bracket is the analysis equation).

Interpretation

  • ∣ak∣2|a_k|^2 is the average power of the kkth harmonic, since 1T∫T∣akejkω0t∣2dt=∣ak∣2\frac{1}{T}\int_T |a_k e^{jk\omega_0 t}|^2dt = |a_k|^2.
  • A plot of ∣ak∣2|a_k|^2 against kω0k\omega_0 is the power spectrum of the signal.

Example: x(t)=Acos⁡ω0tx(t) = A\cos\omega_0 t has a±1=A/2a_{\pm1} = A/2. Frequency domain: (A/2)2+(A/2)2=A2/2(A/2)^2 + (A/2)^2 = A^2/2; time domain: 1T∫TA2cos⁡2ω0t dt=A2/2\frac{1}{T}\int_T A^2\cos^2\omega_0 t\,dt = A^2/2. They match.

  • Asked 3 times
  • 2081 Chaitra · 5 marks
  • 2079 Asoj · 5 marks
  • 2078 Chaitra · 5 marks

Derive the trigonometric Fourier series expansion of a continuous-time periodic signal x(t).

Answer

A periodic signal x(t)x(t) with period TT (ω0=2π/T\omega_0 = 2\pi/T) that meets the Dirichlet conditions can be written as a dc term plus sums of cosines and sines at multiples of ω0\omega_0. This is the trigonometric Fourier series:

x(t)=a0+∑n=1∞[ancos⁡nω0t+bnsin⁡nω0t]x(t) = a_0 + \sum_{n=1}^{\infty}\left[a_n\cos n\omega_0 t + b_n \sin n\omega_0 t\right]

Orthogonality relations (over one period TT)

For integers m,n≥1m, n \ge 1:

∫Tcos⁡mω0t dt=∫Tsin⁡mω0t dt=0∫Tcos⁡mω0t cos⁡nω0t dt={T/2,m=n0,m≠n∫Tsin⁡mω0t sin⁡nω0t dt={T/2,m=n0,m≠n∫Tsin⁡mω0t cos⁡nω0t dt=0for all m,n\begin{aligned} &\int_T \cos m\omega_0 t\,dt = \int_T \sin m\omega_0 t\,dt = 0 \\ &\int_T \cos m\omega_0 t\,\cos n\omega_0 t\,dt = \begin{cases} T/2, & m = n \\ 0, & m \ne n \end{cases} \\ &\int_T \sin m\omega_0 t\,\sin n\omega_0 t\,dt = \begin{cases} T/2, & m = n \\ 0, & m \ne n \end{cases} \\ &\int_T \sin m\omega_0 t\,\cos n\omega_0 t\,dt = 0 \quad \text{for all } m, n \end{aligned}

Finding a0a_0

Integrate both sides over one period. Every sine and cosine term integrates to zero:

∫Tx(t) dt=a0T⇒a0=1T∫Tx(t) dt\int_T x(t)\,dt = a_0 T \quad\Rightarrow\quad a_0 = \frac{1}{T}\int_T x(t)\,dt

Finding ana_n

Multiply both sides by cos⁡mω0t\cos m\omega_0 t and integrate over TT. Only the term with n=mn = m in the cosine sum survives:

∫Tx(t)cos⁡mω0t dt=amT2⇒an=2T∫Tx(t)cos⁡nω0t dt\int_T x(t)\cos m\omega_0 t\,dt = a_m\frac{T}{2} \quad\Rightarrow\quad a_n = \frac{2}{T}\int_T x(t)\cos n\omega_0 t\,dt

Finding bnb_n

Multiply by sin⁡mω0t\sin m\omega_0 t and integrate; only the n=mn = m sine term survives:

∫Tx(t)sin⁡mω0t dt=bmT2⇒bn=2T∫Tx(t)sin⁡nω0t dt\int_T x(t)\sin m\omega_0 t\,dt = b_m\frac{T}{2} \quad\Rightarrow\quad b_n = \frac{2}{T}\int_T x(t)\sin n\omega_0 t\,dt

Summary and notes

CoefficientFormula
a0a_0 (dc)1T∫Tx(t) dt\frac{1}{T}\int_T x(t)\,dt
ana_n2T∫Tx(t)cos⁡nω0t dt\frac{2}{T}\int_T x(t)\cos n\omega_0 t\,dt
bnb_n2T∫Tx(t)sin⁡nω0t dt\frac{2}{T}\int_T x(t)\sin n\omega_0 t\,dt
  • Even x(t)x(t): bn=0b_n = 0 (cosine terms only). Odd x(t)x(t): a0=an=0a_0 = a_n = 0 (sine terms only).
  • Compact (polar) form: x(t)=C0+∑n=1∞Cncos⁡(nω0t+θn)x(t) = C_0 + \sum_{n=1}^{\infty}C_n\cos(n\omega_0 t + \theta_n) with C0=a0C_0 = a_0, Cn=an2+bn2C_n = \sqrt{a_n^2 + b_n^2}, θn=−tan⁡−1(bn/an)\theta_n = -\tan^{-1}(b_n/a_n).
  • Link to exponential form: c0=a0c_0 = a_0, cn=an−jbn2c_n = \frac{a_n - jb_n}{2}, c−n=an+jbn2c_{-n} = \frac{a_n + jb_n}{2}.
  • Asked 3 times
  • 2082 Kartik · 5 marks
  • 2083 Bhadra (new course) · 5 marks
  • 2083 Baisakh (new course) · 4 marks

Find and plot the Fourier series coefficients (magnitude and associated phase) for the following signal with period N. x[n] = 1 + sin((2π/N)n) + 3cos((2π/N)n) + cos((4π/N)n + π/2)

Answer

Fundamental frequency ω0=2π/N\omega_0 = 2\pi/N. Expand each term with Euler's formula and read off the coefficients of ejkω0ne^{jk\omega_0 n}:

x[n]=1+12j[ejω0n−e−jω0n]+32[ejω0n+e−jω0n]+12[ej(2ω0n+π/2)+e−j(2ω0n+π/2)]=1+(32+12j)ejω0n+(32−12j)e−jω0n+12ejπ/2ej2ω0n+12e−jπ/2e−j2ω0n\begin{aligned} x[n] &= 1 + \frac{1}{2j}\left[e^{j\omega_0 n} - e^{-j\omega_0 n}\right] + \frac{3}{2}\left[e^{j\omega_0 n} + e^{-j\omega_0 n}\right] \\ &\quad + \frac{1}{2}\left[e^{j(2\omega_0 n + \pi/2)} + e^{-j(2\omega_0 n + \pi/2)}\right] \\ &= 1 + \left(\frac{3}{2} + \frac{1}{2j}\right)e^{j\omega_0 n} + \left(\frac{3}{2} - \frac{1}{2j}\right)e^{-j\omega_0 n} \\ &\quad + \frac{1}{2}e^{j\pi/2}e^{j2\omega_0 n} + \frac{1}{2}e^{-j\pi/2}e^{-j2\omega_0 n} \end{aligned}

Fourier series coefficients (one period, k=−2k = -2 to 22)

Using 12j=−j2\frac{1}{2j} = -\frac{j}{2}:

kkaka_k∣ak∣\lvert a_k\rvert∠ak\angle a_k
011100
132−j2\frac{3}{2} - \frac{j}{2}102=1.581\frac{\sqrt{10}}{2} = 1.581−18.43∘-18.43^\circ
−1-132+j2\frac{3}{2} + \frac{j}{2}1.581+18.43∘+18.43^\circ
2j2\frac{j}{2}0.5+90∘+90^\circ
−2-2−j2-\frac{j}{2}0.5−90∘-90^\circ
others in the period00-

Working: ∣a1∣=1.52+0.52=2.5=1.581|a_1| = \sqrt{1.5^2 + 0.5^2} = \sqrt{2.5} = 1.581, ∠a1=−tan⁡−1(0.5/1.5)=−18.43∘\angle a_1 = -\tan^{-1}(0.5/1.5) = -18.43^\circ (−0.322-0.322 rad).

The coefficients repeat with period NN: ak+N=aka_{k+N} = a_k (e.g. aN−1=a−1a_{N-1} = a_{-1}, aN−2=a−2a_{N-2} = a_{-2}). Here N≥5N \ge 5 is assumed so that the harmonics do not overlap.

Plots (one period, k=−2…2k = -2 \dots 2, repeating every NN)

 |a_k|  (repeats every N)
 1.581     |     |
     1     |  |  |
   0.5  |  |  |  |  |
       -+--+--+--+--+--> k
       -2 -1  0  1  2

 angle a_k (deg)
    90              |
 18.43     |        |
       -+--+--o--+--+--> k
-18.43  |        |
   -90  |
       -2 -1  0  1  2

Magnitude is even in kk and phase is odd in kk, as expected for a real signal (a−k=ak∗a_{-k} = a_k^*).

  • Asked 2 times
  • 2081 Chaitra · 5 marks
  • 2076 Baisakh · 4 marks

State and prove the time-scaling property for continuous-time (exponential) Fourier series pair.

Answer

Statement

If x(t)x(t) is periodic with period TT (fundamental frequency ω0=2π/T\omega_0 = 2\pi/T) and

x(t)=∑k=−∞∞akejkω0t,x(t) = \sum_{k=-\infty}^{\infty}a_k e^{jk\omega_0 t},

then for α>0\alpha > 0, the time-scaled signal x(αt)x(\alpha t) is periodic with period T/αT/\alpha and fundamental frequency αω0\alpha\omega_0, and its Fourier series coefficients are the same aka_k:

x(αt)=∑k=−∞∞ak ejk(αω0)tx(\alpha t) = \sum_{k=-\infty}^{\infty}a_k\, e^{jk(\alpha\omega_0)t}

Time scaling does not change the coefficients; it only changes the spacing of the spectral lines from ω0\omega_0 to αω0\alpha\omega_0.

Proof

Synthesis route: replace tt by αt\alpha t in the synthesis equation:

x(αt)=∑kakejkω0(αt)=∑kakejk(αω0)tx(\alpha t) = \sum_{k}a_k e^{jk\omega_0(\alpha t)} = \sum_k a_k e^{jk(\alpha\omega_0)t}

This is a Fourier series with fundamental frequency αω0\alpha\omega_0 and coefficients aka_k.

Analysis route (check): let y(t)=x(αt)y(t) = x(\alpha t), period T′=T/αT' = T/\alpha, ω0′=αω0\omega_0' = \alpha\omega_0:

bk=1T′∫0T′x(αt)e−jkαω0tdt(put τ=αt, dt=dτ/α)=αT∫0Tx(τ)e−jkω0τdτα=1T∫0Tx(τ)e−jkω0τdτ=ak\begin{aligned} b_k &= \frac{1}{T'}\int_0^{T'} x(\alpha t)e^{-jk\alpha\omega_0 t}dt \qquad (\text{put } \tau = \alpha t,\ dt = d\tau/\alpha) \\ &= \frac{\alpha}{T}\int_0^{T}x(\tau)e^{-jk\omega_0\tau}\frac{d\tau}{\alpha} \\ &= \frac{1}{T}\int_0^T x(\tau)e^{-jk\omega_0\tau}d\tau = a_k \end{aligned}

Hence bk=akb_k = a_k.

Interpretation

  • Compressing in time (α>1\alpha > 1) spreads the spectral lines further apart in frequency; expanding (α<1\alpha < 1) brings them closer.
  • The amplitude and phase of each harmonic are unchanged, only its frequency kω0k\omega_0 becomes kαω0k\alpha\omega_0.

Example: x(t)=cos⁡(2t)x(t) = \cos(2t) has a±1=1/2a_{\pm1} = 1/2 at ω0=2\omega_0 = 2. Then x(3t)=cos⁡(6t)x(3t) = \cos(6t) still has a±1=1/2a_{\pm1} = 1/2, now at ω0′=6\omega_0' = 6 rad/s.

  • Asked 2 times
  • 2081 Asoj · 6 marks
  • 2082 Bhadra (new course) · 4 marks

Consider a signal x[n] = sin(ω₀n) with period N = 6. Find the Fourier series coefficient of the given signal x[n]. Plot spectra. Also find the average power.

Answer

Since x[n]=sin⁡(ω0n)x[n] = \sin(\omega_0 n) has period N=6N = 6, the fundamental frequency is

ω0=2πN=2π6=π3\omega_0 = \frac{2\pi}{N} = \frac{2\pi}{6} = \frac{\pi}{3}

Fourier series coefficients

Using Euler's formula:

x[n]=sin⁡πn3=12jej(2π/6)n−12je−j(2π/6)nx[n] = \sin\frac{\pi n}{3} = \frac{1}{2j}e^{j(2\pi/6)n} - \frac{1}{2j}e^{-j(2\pi/6)n}

Comparing with x[n]=∑k=⟨6⟩akejk(2π/6)nx[n] = \sum_{k=\langle 6\rangle}a_k e^{jk(2\pi/6)n}:

a1=12j=−j2=12e−jπ/2,a−1=−12j=j2=12ejπ/2a_1 = \frac{1}{2j} = -\frac{j}{2} = \frac{1}{2}e^{-j\pi/2}, \qquad a_{-1} = -\frac{1}{2j} = \frac{j}{2} = \frac{1}{2}e^{j\pi/2}

All other aka_k in one period are zero. The coefficients repeat every 6: a−1=a5=a11a_{-1} = a_5 = a_{11}, a1=a7=a−5a_1 = a_7 = a_{-5}, etc.

kk012345
aka_k0−j/2-j/2000j/2j/2
∣ak∣\lvert a_k\rvert00.50000.5
∠ak\angle a_k-−90∘-90^\circ---+90∘+90^\circ

Spectra

 |a_k|  (repeats every 6)
   0.5  |           |     |           |
       -+--o--o--o--+--o--+--o--o--o--+--> k
       -5 -4 -3 -2 -1  0  1  2  3  4  5

 angle a_k (deg)
    90              |                 |
       -+--o--o--o--+--o--+--o--o--o--+--> k
   -90  |                 |
       -5 -4 -3 -2 -1  0  1  2  3  4  5

Lines of height 1/21/2 at k=±1,±5,±7,…k = \pm1, \pm5, \pm7, \dots; phase −90∘-90^\circ at k=1,7,−5,…k = 1, 7, -5, \dots and +90∘+90^\circ at k=−1,5,…k = -1, 5, \dots

Average power

By Parseval's relation:

P=∑k=⟨6⟩∣ak∣2=∣a1∣2+∣a5∣2=14+14=12P = \sum_{k=\langle 6\rangle}|a_k|^2 = |a_1|^2 + |a_5|^2 = \frac{1}{4} + \frac{1}{4} = \frac{1}{2}

Check in time: samples over one period are 0,32,32,0,−32,−320, \frac{\sqrt3}{2}, \frac{\sqrt3}{2}, 0, -\frac{\sqrt3}{2}, -\frac{\sqrt3}{2}, so P=16(4×34)=0.5P = \frac{1}{6}\left(4\times\frac{3}{4}\right) = 0.5.

Answer: a1=−j/2a_1 = -j/2, a−1=a5=j/2a_{-1} = a_5 = j/2, others 0; average power P=0.5P = 0.5 W.

  • Asked 2 times
  • 2079 Asoj · 4 marks
  • 2079 Chaitra · 4 marks

State and prove the time-scaling property for discrete-time Fourier series pair.

Answer

In discrete time, x[an]x[an] is not defined for non-integer anan, so time scaling is defined through the time-expanded signal (inserting zeros).

Statement

Let x[n]x[n] be periodic with period NN and DTFS coefficients aka_k. For a positive integer mm define

x(m)[n]={x[n/m],n a multiple of m0,otherwisex_{(m)}[n] = \begin{cases} x[n/m], & n \text{ a multiple of } m \\ 0, & \text{otherwise} \end{cases}

(m−1m - 1 zeros are inserted between samples). Then x(m)[n]x_{(m)}[n] is periodic with period mNmN and

x(m)[n]⟷FS1makx_{(m)}[n] \overset{FS}{\longleftrightarrow} \frac{1}{m}a_k

(with aka_k viewed as periodic with period NN, so the mNmN coefficients are ak/ma_k/m for k=0,1,…,mN−1k = 0, 1, \dots, mN-1).

Proof

The period of x(m)[n]x_{(m)}[n] is mNmN, so its coefficients are

bk=1mN∑n=0mN−1x(m)[n]e−jk2πmNnb_k = \frac{1}{mN}\sum_{n=0}^{mN-1}x_{(m)}[n]e^{-jk\frac{2\pi}{mN}n}

Only n=mrn = mr (r=0,1,…,N−1r = 0, 1, \dots, N-1) give non-zero terms, where x(m)[mr]=x[r]x_{(m)}[mr] = x[r]:

bk=1mN∑r=0N−1x[r] e−jk2πmNmr=1m[1N∑r=0N−1x[r]e−jk2πNr]=1mak\begin{aligned} b_k &= \frac{1}{mN}\sum_{r=0}^{N-1}x[r]\,e^{-jk\frac{2\pi}{mN}mr} \\ &= \frac{1}{m}\left[\frac{1}{N}\sum_{r=0}^{N-1}x[r]e^{-jk\frac{2\pi}{N}r}\right] \\ &= \frac{1}{m}a_k \end{aligned}

Interpretation

  • Expanding in time by mm multiplies the period by mm, so the fundamental frequency becomes ω0/m=2π/(mN)\omega_0/m = 2\pi/(mN) and the spectrum is packed mm times more densely.
  • The coefficient values are divided by mm, because the same energy per period is now averaged over mm times as many samples.
  • The coefficient pattern of x[n]x[n] is repeated mm times within one period of bkb_k.

Example: x[n]={1,1}x[n] = \{1, 1\} with N=2N = 2 has a0=1a_0 = 1, a1=0a_1 = 0. With m=2m = 2, x(2)[n]={1,0,1,0}x_{(2)}[n] = \{1, 0, 1, 0\} (N=4N = 4) gives bk=14(1+e−jπk)={0.5,0,0.5,0}b_k = \frac{1}{4}(1 + e^{-j\pi k}) = \{0.5, 0, 0.5, 0\}, which is 12{a0,a1,a0,a1}\frac{1}{2}\{a_0, a_1, a_0, a_1\}.

  • Asked 2 times
  • 2079 Asoj · 6 marks
  • 2079 Chaitra · 6 marks

Find the frequency domain representation (both magnitude and phase plot) of a periodic signal defined as x[n] = {1, −0.5, 0, 0, 0.5} (x[0] = 1), with period N = 5.

Answer

The signal has period N=5N = 5, so ω0=2π/5\omega_0 = 2\pi/5 and

ak=15∑n=04x[n] e−jk(2π/5)na_k = \frac{1}{5}\sum_{n=0}^{4}x[n]\,e^{-jk(2\pi/5)n}

One period: x[0]=1x[0] = 1, x[1]=−0.5x[1] = -0.5, x[2]=0x[2] = 0, x[3]=0x[3] = 0, x[4]=0.5x[4] = 0.5. Because of periodicity x[4]=x[−1]x[4] = x[-1], so it is easier to sum over n=−1,0,1n = -1, 0, 1 (note x[−1]=0.5x[-1] = 0.5):

ak=15[x[−1]ejk2π/5+x[0]+x[1]e−jk2π/5]=15[1+0.5ejk2π/5−0.5e−jk2π/5]=15[1+0.5 (2j)sin⁡2πk5]=15[1+jsin⁡2πk5]\begin{aligned} a_k &= \frac{1}{5}\left[x[-1]e^{jk2\pi/5} + x[0] + x[1]e^{-jk2\pi/5}\right] \\ &= \frac{1}{5}\left[1 + 0.5e^{jk2\pi/5} - 0.5e^{-jk2\pi/5}\right] \\ &= \frac{1}{5}\left[1 + 0.5\,(2j)\sin\frac{2\pi k}{5}\right] \\ &= \frac{1}{5}\left[1 + j\sin\frac{2\pi k}{5}\right] \end{aligned}

Magnitude and phase:

∣ak∣=151+sin⁡22πk5,∠ak=tan⁡−1(sin⁡2πk5)|a_k| = \frac{1}{5}\sqrt{1 + \sin^2\frac{2\pi k}{5}}, \qquad \angle a_k = \tan^{-1}\left(\sin\frac{2\pi k}{5}\right)

Values over one period

kksin⁡(2πk/5)\sin(2\pi k/5)aka_k∣ak∣\lvert a_k\rvert∠ak\angle a_k
000.20.20.2000∘0^\circ
10.9510.2+j0.1900.2 + j0.1900.27643.56∘43.56^\circ
20.5880.2+j0.1180.2 + j0.1180.23230.45∘30.45^\circ
3−0.588-0.5880.2−j0.1180.2 - j0.1180.232−30.45∘-30.45^\circ
4−0.951-0.9510.2−j0.1900.2 - j0.1900.276−43.56∘-43.56^\circ

(Since ak+5=aka_{k+5} = a_k: a−1=a4a_{-1} = a_4, a−2=a3a_{-2} = a_3.)

Magnitude and phase plots (repeat every 5)

 |a_k|  (repeats every 5)
 0.276     |     |
 0.232  |  |     |  |
   0.2  |  |  |  |  |
       -+--+--+--+--+--> k
       -2 -1  0  1  2

 angle a_k (deg)
 43.56           |
 30.45           |  |
       -+--+--o--+--+--> k
-30.45  |  |
-43.56     |
       -2 -1  0  1  2

The magnitude spectrum is even and the phase spectrum is odd in kk, as expected for a real x[n]x[n].

Check (synthesis at n=0n = 0): ∑k=04ak=5×0.2+j(0)=1=x[0]\sum_{k=0}^{4}a_k = 5\times 0.2 + j(0) = 1 = x[0]. Correct.

  • Asked 2 times
  • 2080 Asoj · 6 marks
  • 2070 Magh · 6 marks

Find the Fourier series coefficients of the signal x[n] = 1 + sin[(2π/N)n] + 3cos[(2π/N)n] + cos[(4π/N)n]

Answer

Fundamental frequency ω0=2π/N\omega_0 = 2\pi/N. Rewrite each term with Euler's formula:

x[n]=1+12j[ej2πNn−e−j2πNn]+32[ej2πNn+e−j2πNn]+12[ej4πNn+e−j4πNn]=1+(32+12j)ej(2π/N)n+(32−12j)e−j(2π/N)n+12ej2(2π/N)n+12e−j2(2π/N)n\begin{aligned} x[n] &= 1 + \frac{1}{2j}\left[e^{j\frac{2\pi}{N}n} - e^{-j\frac{2\pi}{N}n}\right] + \frac{3}{2}\left[e^{j\frac{2\pi}{N}n} + e^{-j\frac{2\pi}{N}n}\right] + \frac{1}{2}\left[e^{j\frac{4\pi}{N}n} + e^{-j\frac{4\pi}{N}n}\right] \\ &= 1 + \left(\frac{3}{2} + \frac{1}{2j}\right)e^{j(2\pi/N)n} + \left(\frac{3}{2} - \frac{1}{2j}\right)e^{-j(2\pi/N)n} + \frac{1}{2}e^{j2(2\pi/N)n} + \frac{1}{2}e^{-j2(2\pi/N)n} \end{aligned}

Comparing with x[n]=∑k=⟨N⟩akejk(2π/N)nx[n] = \sum_{k=\langle N\rangle}a_k e^{jk(2\pi/N)n} (and using 12j=−j2\frac{1}{2j} = -\frac{j}{2}):

kkaka_k∣ak∣\lvert a_k\rvert∠ak\angle a_k
011100
132−j2\frac{3}{2} - \frac{j}{2}102=1.581\frac{\sqrt{10}}{2} = 1.581−18.43∘-18.43^\circ
−1-132+j2\frac{3}{2} + \frac{j}{2}1.581+18.43∘+18.43^\circ
212\frac{1}{2}0.500
−2-212\frac{1}{2}0.500
other kk in one period00-

Working: ∣a1∣=1.52+0.52=2.5=1.581|a_1| = \sqrt{1.5^2 + 0.5^2} = \sqrt{2.5} = 1.581; ∠a1=−tan⁡−1(0.5/1.5)=−18.43∘\angle a_1 = -\tan^{-1}(0.5/1.5) = -18.43^\circ.

The coefficients are periodic with period NN: ak+N=aka_{k+N} = a_k, so for example aN−1=a−1a_{N-1} = a_{-1} and aN−2=a−2a_{N-2} = a_{-2} (assuming N≥5N \ge 5 so that the harmonics do not overlap).

 |a_k|  (repeats every N)
 1.581     |     |
     1     |  |  |
   0.5  |  |  |  |  |
       -+--+--+--+--+--> k
       -2 -1  0  1  2

Since x[n]x[n] is real, a−k=ak∗a_{-k} = a_k^*: the magnitude is even and the phase is odd.

  • Asked 2 times
  • 2078 Chaitra · 4 marks
  • 2075 Baisakh · 4 marks

State and prove the time-shifting property for discrete-time Fourier series pair.

Answer

Statement

If x[n]x[n] is periodic with period NN and x[n]⟷FSakx[n] \overset{FS}{\longleftrightarrow} a_k, then

x[n−n0]⟷FSak e−jk(2π/N)n0x[n - n_0] \overset{FS}{\longleftrightarrow} a_k\, e^{-jk(2\pi/N)n_0}

A shift in time multiplies each coefficient by a linear phase term; the magnitudes ∣ak∣|a_k| do not change.

Proof

x[n−n0]x[n - n_0] is also periodic with period NN. Its coefficients are

bk=1N∑n=⟨N⟩x[n−n0] e−jk(2π/N)nb_k = \frac{1}{N}\sum_{n=\langle N\rangle}x[n - n_0]\,e^{-jk(2\pi/N)n}

Put m=n−n0m = n - n_0, so n=m+n0n = m + n_0. As nn runs over one period, mm also runs over one period:

bk=1N∑m=⟨N⟩x[m] e−jk(2π/N)(m+n0)=e−jk(2π/N)n0[1N∑m=⟨N⟩x[m]e−jk(2π/N)m]=e−jk(2π/N)n0 ak\begin{aligned} b_k &= \frac{1}{N}\sum_{m=\langle N\rangle}x[m]\,e^{-jk(2\pi/N)(m + n_0)} \\ &= e^{-jk(2\pi/N)n_0}\left[\frac{1}{N}\sum_{m=\langle N\rangle}x[m]e^{-jk(2\pi/N)m}\right] \\ &= e^{-jk(2\pi/N)n_0}\,a_k \end{aligned}

Interpretation

  • ∣bk∣=∣ak∣|b_k| = |a_k|: the magnitude spectrum is unchanged by a shift.
  • ∠bk=∠ak−k2πNn0\angle b_k = \angle a_k - k\frac{2\pi}{N}n_0: the phase changes linearly with kk.

Example: x[n]=cos⁡(2πn/N)x[n] = \cos(2\pi n/N) has a±1=1/2a_{\pm1} = 1/2. Then x[n−1]x[n-1] has b1=12e−j2π/Nb_{1} = \frac{1}{2}e^{-j2\pi/N} and b−1=12ej2π/Nb_{-1} = \frac{1}{2}e^{j2\pi/N}, still of magnitude 1/21/2.

  • Asked 2 times
  • 2074 Bhadra · 4 marks
  • 2083 Bhadra (new course) · 4 marks

State and prove time shifting and time scaling properties of continuous time Fourier series.

Answer

Let x(t)x(t) be periodic with period TT, ω0=2π/T\omega_0 = 2\pi/T, and x(t)⟷FSakx(t) \overset{FS}{\longleftrightarrow} a_k, where ak=1T∫Tx(t)e−jkω0tdta_k = \frac{1}{T}\int_T x(t)e^{-jk\omega_0 t}dt.

Time shifting property

Statement:

x(t−t0)⟷FSak e−jkω0t0x(t - t_0) \overset{FS}{\longleftrightarrow} a_k\, e^{-jk\omega_0 t_0}

Proof: y(t)=x(t−t0)y(t) = x(t - t_0) has the same period TT. Its coefficients are

bk=1T∫Tx(t−t0)e−jkω0tdt(τ=t−t0)=1T∫Tx(τ)e−jkω0(τ+t0)dτ=e−jkω0t0 1T∫Tx(τ)e−jkω0τdτ=e−jkω0t0 ak\begin{aligned} b_k &= \frac{1}{T}\int_T x(t - t_0)e^{-jk\omega_0 t}dt \qquad (\tau = t - t_0) \\ &= \frac{1}{T}\int_T x(\tau)e^{-jk\omega_0(\tau + t_0)}d\tau \\ &= e^{-jk\omega_0 t_0}\,\frac{1}{T}\int_T x(\tau)e^{-jk\omega_0\tau}d\tau = e^{-jk\omega_0 t_0}\,a_k \end{aligned}

So ∣bk∣=∣ak∣|b_k| = |a_k| (magnitude unchanged) and the phase changes by −kω0t0-k\omega_0 t_0.

Time scaling property

Statement: for α>0\alpha > 0, x(αt)x(\alpha t) is periodic with period T/αT/\alpha and fundamental frequency αω0\alpha\omega_0, and

x(αt)=∑k=−∞∞akejk(αω0)tx(\alpha t) = \sum_{k=-\infty}^{\infty}a_k e^{jk(\alpha\omega_0)t}

The coefficients stay aka_k; only the fundamental frequency changes.

Proof: let y(t)=x(αt)y(t) = x(\alpha t), T′=T/αT' = T/\alpha:

bk=1T′∫0T′x(αt)e−jk(αω0)tdt(τ=αt, dt=dτ/α)=αT∫0Tx(τ)e−jkω0τdτα=ak\begin{aligned} b_k &= \frac{1}{T'}\int_0^{T'}x(\alpha t)e^{-jk(\alpha\omega_0)t}dt \qquad (\tau = \alpha t,\ dt = d\tau/\alpha) \\ &= \frac{\alpha}{T}\int_0^{T}x(\tau)e^{-jk\omega_0\tau}\frac{d\tau}{\alpha} = a_k \end{aligned}

Summary

PropertySignalCoefficientsFundamental
Time shiftx(t−t0)x(t - t_0)ake−jkω0t0a_k e^{-jk\omega_0 t_0}ω0\omega_0
Time scalex(αt)x(\alpha t)aka_kαω0\alpha\omega_0

Example: x(t)=cos⁡ω0tx(t) = \cos\omega_0 t (a±1=1/2a_{\pm1} = 1/2); x(t−t0)x(t - t_0) has a±1=12e∓jω0t0a_{\pm1} = \frac{1}{2}e^{\mp j\omega_0 t_0}; x(2t)=cos⁡2ω0tx(2t) = \cos 2\omega_0 t still has a±1=1/2a_{\pm1} = 1/2 but at fundamental 2ω02\omega_0.

  • 2082 Chaitra · 6 marks

Describe the time shifting and frequency shifting properties of continuous time Fourier series. Find Fourier series coefficient of signal x(t) = A cos(ω₀t).

Answer

Let x(t)x(t) have period TT, ω0=2π/T\omega_0 = 2\pi/T, and x(t)⟷FSakx(t) \overset{FS}{\longleftrightarrow} a_k.

Time shifting property

x(t−t0)⟷FSake−jkω0t0x(t - t_0) \overset{FS}{\longleftrightarrow} a_k e^{-jk\omega_0 t_0}

Proof sketch: with τ=t−t0\tau = t - t_0,

bk=1T∫Tx(t−t0)e−jkω0tdt=e−jkω0t01T∫Tx(τ)e−jkω0τdτ=e−jkω0t0akb_k = \frac{1}{T}\int_T x(t - t_0)e^{-jk\omega_0 t}dt = e^{-jk\omega_0 t_0}\frac{1}{T}\int_T x(\tau)e^{-jk\omega_0\tau}d\tau = e^{-jk\omega_0 t_0}a_k

A delay leaves ∣ak∣|a_k| unchanged and adds a phase −kω0t0-k\omega_0 t_0 that grows linearly with kk.

Frequency shifting property

ejMω0t x(t)⟷FSak−Me^{jM\omega_0 t}\,x(t) \overset{FS}{\longleftrightarrow} a_{k-M}

Proof:

bk=1T∫TejMω0tx(t)e−jkω0tdt=1T∫Tx(t)e−j(k−M)ω0tdt=ak−Mb_k = \frac{1}{T}\int_T e^{jM\omega_0 t}x(t)e^{-jk\omega_0 t}dt = \frac{1}{T}\int_T x(t)e^{-j(k-M)\omega_0 t}dt = a_{k-M}

Multiplying by a complex exponential at the MMth harmonic shifts the whole set of coefficients by MM places (this is the basis of modulation).

Fourier series coefficients of x(t)=Acos⁡(ω0t)x(t) = A\cos(\omega_0 t)

Using Euler's formula:

x(t)=Acos⁡ω0t=A2ejω0t+A2e−jω0tx(t) = A\cos\omega_0 t = \frac{A}{2}e^{j\omega_0 t} + \frac{A}{2}e^{-j\omega_0 t}

Comparing with x(t)=∑kakejkω0tx(t) = \sum_k a_k e^{jk\omega_0 t}:

a1=a−1=A2,ak=0 for k≠±1a_1 = a_{-1} = \frac{A}{2}, \qquad a_k = 0 \text{ for } k \ne \pm1

Check with the analysis equation for k=1k = 1:

a1=1T∫0TAcos⁡ω0t e−jω0tdt=A2T∫0T(1+e−j2ω0t)dt=A2a_1 = \frac{1}{T}\int_0^T A\cos\omega_0 t\,e^{-j\omega_0 t}dt = \frac{A}{2T}\int_0^T\left(1 + e^{-j2\omega_0 t}\right)dt = \frac{A}{2}
 a_k
 A/2      |       |
      ----+---+---+----> k
         -1   0   1

Answer: a±1=A/2a_{\pm1} = A/2 (real, zero phase), all other coefficients are zero.

  • 2081 Chaitra · 3 marks

State and prove frequency shifting property of discrete time Fourier series.

Answer

Statement

If x[n]x[n] is periodic with period NN and x[n]⟷FSakx[n] \overset{FS}{\longleftrightarrow} a_k, then for an integer MM:

ejM(2π/N)n x[n]⟷FSak−Me^{jM(2\pi/N)n}\,x[n] \overset{FS}{\longleftrightarrow} a_{k-M}

Multiplying by a complex exponential at the MMth harmonic shifts the coefficient sequence by MM.

Proof

y[n]=ejM(2π/N)nx[n]y[n] = e^{jM(2\pi/N)n}x[n] is also periodic with period NN (both factors are). Its coefficients are

bk=1N∑n=⟨N⟩ejM(2π/N)nx[n] e−jk(2π/N)n=1N∑n=⟨N⟩x[n] e−j(k−M)(2π/N)n=ak−M\begin{aligned} b_k &= \frac{1}{N}\sum_{n=\langle N\rangle}e^{jM(2\pi/N)n}x[n]\,e^{-jk(2\pi/N)n} \\ &= \frac{1}{N}\sum_{n=\langle N\rangle}x[n]\,e^{-j(k-M)(2\pi/N)n} \\ &= a_{k-M} \end{aligned}

Example: for N=4N = 4, multiplying by ejπn/2e^{j\pi n/2} (M=1M = 1) moves a0a_0 to position k=1k = 1, a1a_1 to k=2k = 2, and so on (cyclically, since aka_k has period 4).

  • 2081 Asoj · 6 marks

What is Fourier series representation of a signal? Derive an expression for analysis and synthesis equation of the Fourier Series for a CT periodic signal x(t).

Answer

The Fourier series representation of a periodic signal expresses it as a sum (linear combination) of harmonically related sinusoids or complex exponentials, whose frequencies are integer multiples of the fundamental frequency ω0=2π/T\omega_0 = 2\pi/T. It shows which frequencies the signal contains and how strongly.

Synthesis equation

Each ejkω0te^{jk\omega_0 t}, k=0,±1,±2,…k = 0, \pm1, \pm2, \dots, is periodic with period TT, so their combination is also periodic with TT:

x(t)=∑k=−∞∞ak ejkω0tx(t) = \sum_{k=-\infty}^{\infty}a_k\,e^{jk\omega_0 t}

This builds (synthesises) x(t)x(t) from its harmonics; aka_k are the Fourier series coefficients.

Analysis equation

Multiply both sides by e−jnω0te^{-jn\omega_0 t} and integrate over one period:

∫0Tx(t)e−jnω0tdt=∑k=−∞∞ak∫0Tej(k−n)ω0tdt\int_0^T x(t)e^{-jn\omega_0 t}dt = \sum_{k=-\infty}^{\infty}a_k\int_0^T e^{j(k-n)\omega_0 t}dt

For k≠nk \ne n, ej(k−n)ω0te^{j(k-n)\omega_0 t} completes (k−n)(k - n) whole cycles in TT, so

∫0Tej(k−n)ω0tdt=∫0Tcos⁡((k−n)ω0t)dt+j∫0Tsin⁡((k−n)ω0t)dt=0\int_0^T e^{j(k-n)\omega_0 t}dt = \int_0^T\cos((k-n)\omega_0 t)dt + j\int_0^T\sin((k-n)\omega_0 t)dt = 0

For k=nk = n the integrand is 1 and the integral is TT. So only k=nk = n survives:

∫0Tx(t)e−jnω0tdt=Tan\int_0^T x(t)e^{-jn\omega_0 t}dt = T a_n ak=1T∫Tx(t) e−jkω0t dta_k = \frac{1}{T}\int_T x(t)\,e^{-jk\omega_0 t}\,dt

This extracts (analyses) the amount of the kkth harmonic in x(t)x(t); the integral may be taken over any interval of length TT.

Summary

NameEquation
Synthesisx(t)=∑k=−∞∞akejkω0tx(t) = \sum_{k=-\infty}^{\infty}a_k e^{jk\omega_0 t}
Analysisak=1T∫Tx(t)e−jkω0tdta_k = \frac{1}{T}\int_T x(t)e^{-jk\omega_0 t}dt
DC terma0=1T∫Tx(t)dta_0 = \frac{1}{T}\int_T x(t)dt

Example: x(t)=sin⁡ω0t=12jejω0t−12je−jω0tx(t) = \sin\omega_0 t = \frac{1}{2j}e^{j\omega_0 t} - \frac{1}{2j}e^{-j\omega_0 t}, so a1=12ja_1 = \frac{1}{2j}, a−1=−12ja_{-1} = -\frac{1}{2j}, all others 0.

  • 2080 Asoj · 2+6 marks

What are the significance of Fourier series coefficients of a signal? State and prove time scaling and conjugation properties of continuous time Fourier series representation.

Answer

Significance of Fourier series coefficients

The coefficients ak=1T∫Tx(t)e−jkω0tdta_k = \frac{1}{T}\int_T x(t)e^{-jk\omega_0 t}dt describe the signal in the frequency domain:

  • ∣ak∣|a_k| gives the amplitude and ∠ak\angle a_k the phase of the kkth harmonic (frequency kω0k\omega_0); plotting them gives the magnitude and phase spectra.
  • a0a_0 is the dc (average) value of the signal.
  • ∣ak∣2|a_k|^2 is the power in the kkth harmonic, and by Parseval ∑∣ak∣2\sum|a_k|^2 is the total average power.
  • They show the bandwidth of the signal: how fast ∣ak∣|a_k| falls tells how many harmonics are needed (smooth signals need fewer).
  • An LTI system acts on each harmonic separately: the output coefficients are akH(jkω0)a_k H(jk\omega_0), which makes system analysis simple.

Time scaling property

Statement: if x(t)x(t) has period TT, ω0=2π/T\omega_0 = 2\pi/T, and x(t)⟷FSakx(t) \overset{FS}{\longleftrightarrow} a_k, then for α>0\alpha > 0, x(αt)x(\alpha t) has period T/αT/\alpha and

x(αt)=∑k=−∞∞akejk(αω0)tx(\alpha t) = \sum_{k=-\infty}^{\infty}a_k e^{jk(\alpha\omega_0)t}

i.e. the coefficients are unchanged but the fundamental frequency becomes αω0\alpha\omega_0.

Proof: let y(t)=x(αt)y(t) = x(\alpha t), period T′=T/αT' = T/\alpha; put τ=αt\tau = \alpha t:

bk=1T′∫0T′x(αt)e−jkαω0tdt=αT∫0Tx(τ)e−jkω0τdτα=ak\begin{aligned} b_k &= \frac{1}{T'}\int_0^{T'}x(\alpha t)e^{-jk\alpha\omega_0 t}dt \\ &= \frac{\alpha}{T}\int_0^{T}x(\tau)e^{-jk\omega_0\tau}\frac{d\tau}{\alpha} = a_k \end{aligned}

Conjugation property

Statement:

x∗(t)⟷FSa−k∗x^*(t) \overset{FS}{\longleftrightarrow} a_{-k}^*

Proof: take the conjugate of the synthesis equation:

x∗(t)=[∑k=−∞∞akejkω0t]∗=∑k=−∞∞ak∗e−jkω0t=∑m=−∞∞a−m∗ ejmω0t(m=−k)\begin{aligned} x^*(t) &= \left[\sum_{k=-\infty}^{\infty}a_k e^{jk\omega_0 t}\right]^* = \sum_{k=-\infty}^{\infty}a_k^* e^{-jk\omega_0 t} \\ &= \sum_{m=-\infty}^{\infty}a_{-m}^*\,e^{jm\omega_0 t} \qquad (m = -k) \end{aligned}

So the coefficient of ejmω0te^{jm\omega_0 t} in x∗(t)x^*(t) is a−m∗a_{-m}^*.

Consequence (conjugate symmetry): if x(t)x(t) is real, x(t)=x∗(t)x(t) = x^*(t), so ak=a−k∗a_k = a_{-k}^*. Hence ∣ak∣=∣a−k∣|a_k| = |a_{-k}| (even magnitude spectrum) and ∠a−k=−∠ak\angle a_{-k} = -\angle a_k (odd phase spectrum). If x(t)x(t) is also even, aka_k is real and even.

  • 2079 Jestha · 8 marks

Calculate and plot magnitude and phase spectrum of the following periodic signal. [Figure: sawtooth wave of period 1: x(t) rises linearly from 0 at t = 0 to 1 at t = 1 and then drops back to 0, i.e. x(t) = t for 0 ≤ t < 1, repeating every 1 s; shown from t = −3 to t = 4]

Answer

One period: x(t)=tx(t) = t for 0≤t<10 \le t < 1, with period T=1T = 1 s, so ω0=2π/T=2π\omega_0 = 2\pi/T = 2\pi rad/s.

ak=1T∫0Tx(t)e−jkω0tdt=∫01t e−j2πktdta_k = \frac{1}{T}\int_0^T x(t)e^{-jk\omega_0 t}dt = \int_0^1 t\,e^{-j2\pi kt}dt

DC term (k=0k = 0)

a0=∫01t dt=12a_0 = \int_0^1 t\,dt = \frac{1}{2}

Coefficients for k≠0k \ne 0

Integrate by parts with u=tu = t, dv=e−j2πktdtdv = e^{-j2\pi kt}dt, v=e−j2πkt−j2πkv = \frac{e^{-j2\pi kt}}{-j2\pi k}:

ak=[t e−j2πkt−j2πk]01−∫01e−j2πkt−j2πkdt=e−j2πk−j2πk+1j2πk[e−j2πkt−j2πk]01=1−j2πk+1j2πk⋅e−j2πk−1−j2πk\begin{aligned} a_k &= \left[\frac{t\,e^{-j2\pi kt}}{-j2\pi k}\right]_0^1 - \int_0^1\frac{e^{-j2\pi kt}}{-j2\pi k}dt \\ &= \frac{e^{-j2\pi k}}{-j2\pi k} + \frac{1}{j2\pi k}\left[\frac{e^{-j2\pi kt}}{-j2\pi k}\right]_0^1 \\ &= \frac{1}{-j2\pi k} + \frac{1}{j2\pi k}\cdot\frac{e^{-j2\pi k} - 1}{-j2\pi k} \end{aligned}

Since e−j2πk=1e^{-j2\pi k} = 1, the second term is zero:

ak=−1j2πk=j2πk,k≠0a_k = \frac{-1}{j2\pi k} = \frac{j}{2\pi k}, \qquad k \ne 0

Magnitude and phase

∣ak∣=12π∣k∣,∠ak={+90∘ (π/2),k>0−90∘ (−π/2),k<00,k=0|a_k| = \frac{1}{2\pi|k|}, \qquad \angle a_k = \begin{cases} +90^\circ\ (\pi/2), & k > 0 \\ -90^\circ\ (-\pi/2), & k < 0 \\ 0, & k = 0 \end{cases}
kk0±1\pm1±2\pm2±3\pm3±4\pm4
∣ak∣\lvert a_k\rvert0.50.1590.07960.05310.0398
∠ak\angle a_k (k>0k>0)0∘0^\circ90∘90^\circ90∘90^\circ90∘90^\circ90∘90^\circ
∠ak\angle a_k (k<0k<0)-−90∘-90^\circ−90∘-90^\circ−90∘-90^\circ−90∘-90^\circ

So the Fourier series is

x(t)=12+∑k≠0j2πkej2πkt=12−∑k=1∞1πksin⁡(2πkt)x(t) = \frac{1}{2} + \sum_{k\ne0}\frac{j}{2\pi k}e^{j2\pi kt} = \frac{1}{2} - \sum_{k=1}^{\infty}\frac{1}{\pi k}\sin(2\pi kt)

Spectra (lines at ω=2πk\omega = 2\pi k)

 |a_k|
   0.5           |
 0.159        |  |  |
0.0796     |  |  |  |  |
0.0531  |  |  |  |  |  |  |
       -+--+--+--+--+--+--+--> k
       -3 -2 -1  0  1  2  3

 angle a_k (deg)
    90              |  |  |
       -+--+--+--o--+--+--+--> k
   -90  |  |  |
       -3 -2 -1  0  1  2  3

The magnitude falls as 1/k1/k (because of the jump discontinuity each period), is even in kk, and the phase is odd: +90∘+90^\circ for positive kk and −90∘-90^\circ for negative kk.

Check: at t=0.25t = 0.25, the series 12−1π(1−13+15−… )=12−1π⋅π4=0.25\frac{1}{2} - \frac{1}{\pi}\left(1 - \frac{1}{3} + \frac{1}{5} - \dots\right) = \frac{1}{2} - \frac{1}{\pi}\cdot\frac{\pi}{4} = 0.25, which equals x(0.25)x(0.25).

  • 2078 Baisakh · 4 marks

"Even if a signal is periodic, its Fourier series might not exist." Justify the statement using three different examples.

Answer

A periodic signal has a Fourier series that converges to it only if it meets the Dirichlet conditions over one period: (1) absolutely integrable, (2) finite number of maxima and minima, (3) finite number of finite discontinuities. Periodicity alone is not enough. Three periodic signals (period T=1T = 1) that each break one condition:

Example 1: not absolutely integrable

x(t)=1t,0<t≤1,repeated with period 1x(t) = \frac{1}{t}, \quad 0 < t \le 1, \quad \text{repeated with period } 1 ∫01∣1t∣dt=[ln⁡t]01=∞\int_0^1\left|\frac{1}{t}\right|dt = \left[\ln t\right]_0^1 = \infty

The coefficients ak=∫011te−j2πktdta_k = \int_0^1 \frac{1}{t}e^{-j2\pi kt}dt are not finite, so the series does not exist.

Example 2: infinite number of maxima and minima

x(t)=sin⁡(2πt),0<t≤1,period 1x(t) = \sin\left(\frac{2\pi}{t}\right), \quad 0 < t \le 1, \quad \text{period } 1

It is bounded and absolutely integrable (∫01∣x∣ dt<1\int_0^1|x|\,dt < 1), but as t→0t \to 0 it oscillates infinitely fast, giving infinitely many maxima and minima in one period. The series does not converge to x(t)x(t).

Example 3: infinite number of discontinuities

A staircase signal over 0≤t<10 \le t < 1 that equals 1 on [0,1/2)[0, 1/2), 1/21/2 on [1/2,3/4)[1/2, 3/4), 1/41/4 on [3/4,7/8)[3/4, 7/8), and so on (each step half the height and half the length of the previous one), repeated every 1 s.

It is bounded and ∫01∣x∣ dt<1\int_0^1|x|\,dt < 1, but it has an infinite number of discontinuities within one period, so it violates the third condition.

ExampleCondition violated
1/t1/t on (0,1](0, 1]absolute integrability
sin⁡(2π/t)\sin(2\pi/t)finite maxima/minima
halving staircasefinite discontinuities

So even though all three are periodic, their Fourier series do not exist (or do not converge to them). Such signals are mathematical curiosities; practical signals usually meet the Dirichlet conditions.

  • 2078 Poush · 5+3 marks

Derive the formula for synthesis equation and analysis equation of Fourier series. Mention the three conditions that guarantee the existence of Fourier series.

Answer

A periodic signal x(t)x(t) with period TT and fundamental frequency ω0=2π/T\omega_0 = 2\pi/T can be written as a sum of harmonically related complex exponentials (the exponential Fourier series).

Synthesis equation

The functions ejkω0te^{jk\omega_0 t}, k=0,±1,±2,…k = 0, \pm1, \pm2, \dots, all repeat every TT, so

x(t)=∑k=−∞∞ak ejkω0tx(t) = \sum_{k=-\infty}^{\infty}a_k\,e^{jk\omega_0 t}

is periodic with period TT. This is the synthesis equation; aka_k are the Fourier series coefficients.

Analysis equation

Multiply both sides by e−jnω0te^{-jn\omega_0 t} and integrate over one period:

∫0Tx(t)e−jnω0tdt=∑k=−∞∞ak∫0Tej(k−n)ω0tdt\int_0^T x(t)e^{-jn\omega_0 t}dt = \sum_{k=-\infty}^{\infty}a_k\int_0^T e^{j(k-n)\omega_0 t}dt

Orthogonality of the exponentials:

∫0Tej(k−n)ω0tdt={T,k=n0,k≠n\int_0^T e^{j(k-n)\omega_0 t}dt = \begin{cases} T, & k = n \\ 0, & k \ne n \end{cases}

(for k≠nk \ne n, the cosine and sine parts complete a whole number of cycles and integrate to zero). Hence

∫0Tx(t)e−jnω0tdt=anT\int_0^T x(t)e^{-jn\omega_0 t}dt = a_n T ak=1T∫Tx(t) e−jkω0t dta_k = \frac{1}{T}\int_T x(t)\,e^{-jk\omega_0 t}\,dt

This is the analysis equation; a0=1T∫Tx(t)dta_0 = \frac{1}{T}\int_T x(t)dt is the dc value.

EquationFormula
Synthesisx(t)=∑kakejkω0tx(t) = \sum_k a_k e^{jk\omega_0 t}
Analysisak=1T∫Tx(t)e−jkω0tdta_k = \frac{1}{T}\int_T x(t)e^{-jk\omega_0 t}dt

Conditions for existence (Dirichlet conditions)

Over any one period:

  1. Absolutely integrable: ∫T∣x(t)∣ dt<∞\int_T |x(t)|\,dt < \infty, which guarantees every ∣ak∣|a_k| is finite.
  2. Finite number of maxima and minima (bounded variation) in one period.
  3. Finite number of discontinuities in one period, and each discontinuity is finite.

When these hold, the series converges to x(t)x(t) at every point where x(t)x(t) is continuous, and to the average 12[x(t−)+x(t+)]\frac{1}{2}[x(t^-) + x(t^+)] at a discontinuity. (A weaker condition, finite energy over one period, ∫T∣x(t)∣2dt<∞\int_T|x(t)|^2dt < \infty, guarantees convergence in the mean-square sense.)

  • 2078 Poush · 6 marks

State and explain four properties of discrete time Fourier series.

Answer

Let x[n]x[n] and y[n]y[n] be periodic with the same period NN, ω0=2π/N\omega_0 = 2\pi/N, with x[n]⟷FSakx[n] \overset{FS}{\longleftrightarrow} a_k and y[n]⟷FSbky[n] \overset{FS}{\longleftrightarrow} b_k, where ak=1N∑n=⟨N⟩x[n]e−jkω0na_k = \frac{1}{N}\sum_{n=\langle N\rangle}x[n]e^{-jk\omega_0 n}.

1. Linearity

A x[n]+B y[n]⟷FSAak+BbkA\,x[n] + B\,y[n] \overset{FS}{\longleftrightarrow} A a_k + B b_k

The analysis sum is linear, so the coefficients of a weighted sum are the same weighted sum of coefficients. Example: x[n]=2cos⁡(ω0n)+3x[n] = 2\cos(\omega_0 n) + 3 has a0=3a_0 = 3, a±1=1a_{\pm1} = 1.

2. Time shifting

x[n−n0]⟷FSake−jk(2π/N)n0x[n - n_0] \overset{FS}{\longleftrightarrow} a_k e^{-jk(2\pi/N)n_0}

Proof: with m=n−n0m = n - n_0, 1N∑nx[n−n0]e−jkω0n=e−jkω0n01N∑mx[m]e−jkω0m\frac{1}{N}\sum_n x[n-n_0]e^{-jk\omega_0 n} = e^{-jk\omega_0 n_0}\frac{1}{N}\sum_m x[m]e^{-jk\omega_0 m}. A shift leaves ∣ak∣|a_k| unchanged and adds a linear phase.

3. Frequency shifting

ejM(2π/N)nx[n]⟷FSak−Me^{jM(2\pi/N)n}x[n] \overset{FS}{\longleftrightarrow} a_{k-M}

Multiplying by the MMth harmonic exponential shifts the coefficients by MM (used in modulation).

4. Parseval's relation

1N∑n=⟨N⟩∣x[n]∣2=∑k=⟨N⟩∣ak∣2\frac{1}{N}\sum_{n=\langle N\rangle}|x[n]|^2 = \sum_{k=\langle N\rangle}|a_k|^2

The average power in one period equals the sum of the powers of the NN harmonics.

Other properties (for reference)

PropertySignalCoefficients
Time reversalx[−n]x[-n]a−ka_{-k}
Conjugationx∗[n]x^*[n]a−k∗a_{-k}^*
Periodic convolution∑r=⟨N⟩x[r]y[n−r]\sum_{r=\langle N\rangle}x[r]y[n-r]NakbkN a_k b_k
Multiplicationx[n]y[n]x[n]y[n]∑l=⟨N⟩albk−l\sum_{l=\langle N\rangle}a_l b_{k-l}
First differencex[n]−x[n−1]x[n] - x[n-1](1−e−jkω0)ak(1 - e^{-jk\omega_0})a_k
Real signalx[n]x[n] realak=a−k∗a_k = a_{-k}^*

A property special to the DTFS is that the coefficients are themselves periodic: ak+N=aka_{k+N} = a_k.

  • 2080 Chaitra · 3+5 marks

Write down the condition of convergence of continuous time Fourier series. Find the Fourier series coefficient a_k for the continuous time periodic signal x(t) = 1.5 for 0 ≤ t < 1; −1.5 for 1 ≤ t < 2.

Answer

Condition of convergence (Dirichlet conditions)

The Fourier series of a periodic x(t)x(t) converges if, over one period:

  1. x(t)x(t) is absolutely integrable: ∫T∣x(t)∣dt<∞\int_T|x(t)|dt < \infty.
  2. x(t)x(t) has a finite number of maxima and minima.
  3. x(t)x(t) has a finite number of finite discontinuities.

Then the series equals x(t)x(t) wherever x(t)x(t) is continuous and equals 12[x(t−)+x(t+)]\frac{1}{2}[x(t^-) + x(t^+)] at a jump. (Alternatively, finite energy in one period gives convergence in the mean-square sense.)

Fourier series coefficients

Period T=2T = 2 s, so ω0=2π/T=π\omega_0 = 2\pi/T = \pi rad/s.

x(t)={1.5,0≤t<1−1.5,1≤t<2x(t) = \begin{cases} 1.5, & 0 \le t < 1 \\ -1.5, & 1 \le t < 2 \end{cases}

DC term:

a0=12[∫011.5 dt+∫12(−1.5) dt]=12(1.5−1.5)=0a_0 = \frac{1}{2}\left[\int_0^1 1.5\,dt + \int_1^2(-1.5)\,dt\right] = \frac{1}{2}(1.5 - 1.5) = 0

For k≠0k \ne 0:

ak=12[∫011.5 e−jkπtdt−∫121.5 e−jkπtdt]=1.52[1−e−jkπjkπ−e−jkπ−e−j2kπjkπ]\begin{aligned} a_k &= \frac{1}{2}\left[\int_0^1 1.5\,e^{-jk\pi t}dt - \int_1^2 1.5\,e^{-jk\pi t}dt\right] \\ &= \frac{1.5}{2}\left[\frac{1 - e^{-jk\pi}}{jk\pi} - \frac{e^{-jk\pi} - e^{-j2k\pi}}{jk\pi}\right] \end{aligned}

Using e−j2kπ=1e^{-j2k\pi} = 1 and e−jkπ=(−1)ke^{-jk\pi} = (-1)^k:

ak=0.75jkπ[1−(−1)k−(−1)k+1]=0.75×2[1−(−1)k]jkπ=1.5[1−(−1)k]jkπ\begin{aligned} a_k &= \frac{0.75}{jk\pi}\left[1 - (-1)^k - (-1)^k + 1\right] = \frac{0.75 \times 2\left[1 - (-1)^k\right]}{jk\pi} \\ &= \frac{1.5\left[1 - (-1)^k\right]}{jk\pi} \end{aligned}

So

ak={3jkπ=−j3kπ,k odd0,k even (including k=0)a_k = \begin{cases} \dfrac{3}{jk\pi} = -j\dfrac{3}{k\pi}, & k \text{ odd} \\ 0, & k \text{ even (including } k = 0) \end{cases}
kk±1\pm1±2\pm2±3\pm3±5\pm5
aka_k∓j0.955\mp j0.9550∓j0.318\mp j0.318∓j0.191\mp j0.191
∣ak∣\lvert a_k\rvert0.95500.3180.191

Phase: −90∘-90^\circ for positive odd kk and +90∘+90^\circ for negative odd kk (the signal is odd, so the aka_k are purely imaginary).

Equivalent trigonometric form:

x(t)=6π[sin⁡πt+13sin⁡3πt+15sin⁡5πt+… ]x(t) = \frac{6}{\pi}\left[\sin\pi t + \frac{1}{3}\sin 3\pi t + \frac{1}{5}\sin 5\pi t + \dots\right]

Answer: a0=0a_0 = 0; ak=3jkπa_k = \dfrac{3}{jk\pi} for odd kk; ak=0a_k = 0 for even kk.

  • 2079 Chaitra · 5+3 marks

Derive the expression of Fourier Series Coefficient in exponential form for Continuous time periodic signal. Describe the condition for existence of Fourier series with suitable examples.

Answer

Exponential Fourier series coefficients

Let x(t)x(t) be periodic with period TT, ω0=2π/T\omega_0 = 2\pi/T. It is expressed as a sum of harmonically related complex exponentials (synthesis equation):

x(t)=∑k=−∞∞akejkω0tx(t) = \sum_{k=-\infty}^{\infty}a_k e^{jk\omega_0 t}

Multiply both sides by e−jnω0te^{-jn\omega_0 t} and integrate over one period:

∫0Tx(t)e−jnω0tdt=∑k=−∞∞ak∫0Tej(k−n)ω0tdt\int_0^T x(t)e^{-jn\omega_0 t}dt = \sum_{k=-\infty}^{\infty}a_k\int_0^T e^{j(k-n)\omega_0 t}dt

The right-hand integral:

∫0Tej(k−n)ω0tdt=[ej(k−n)ω0tj(k−n)ω0]0T=ej(k−n)2π−1j(k−n)ω0=0(k≠n)\int_0^T e^{j(k-n)\omega_0 t}dt = \left[\frac{e^{j(k-n)\omega_0 t}}{j(k-n)\omega_0}\right]_0^T = \frac{e^{j(k-n)2\pi} - 1}{j(k-n)\omega_0} = 0 \quad (k \ne n)

and equals TT when k=nk = n. Only the k=nk = n term remains:

∫0Tx(t)e−jnω0tdt=anT\int_0^T x(t)e^{-jn\omega_0 t}dt = a_n T ak=1T∫Tx(t)e−jkω0tdta_k = \frac{1}{T}\int_T x(t)e^{-jk\omega_0 t}dt

This is the analysis equation. a0=1T∫Tx(t)dta_0 = \frac{1}{T}\int_T x(t)dt is the average value, and ∣ak∣|a_k|, ∠ak\angle a_k give the magnitude and phase spectra.

Conditions for existence (Dirichlet conditions)

Over one period:

  1. Absolutely integrable: ∫T∣x(t)∣dt<∞\int_T|x(t)|dt < \infty.
    • Violated by x(t)=1/tx(t) = 1/t on 0<t≤10 < t \le 1 (period 1): ∫01dt/t=∞\int_0^1 dt/t = \infty.
  2. Finite number of maxima and minima in a period.
    • Violated by x(t)=sin⁡(2π/t)x(t) = \sin(2\pi/t) on 0<t≤10 < t \le 1: it oscillates infinitely often near t=0t = 0.
  3. Finite number of discontinuities, each of finite size, in a period.
    • Violated by a staircase that steps down by half on [0,12),[12,34),[34,78),…[0, \frac12), [\frac12, \frac34), [\frac34, \frac78), \dots: infinitely many jumps in one period.

A signal that satisfies them: a square wave has two finite jumps per period, is bounded and has no oscillations, so its Fourier series exists; it converges to the midpoint value at each jump.

  • 2078 Chaitra · 2+3 marks

State and prove convolution property of discrete time periodic signals.

Answer

Statement

Let x[n]x[n] and y[n]y[n] be periodic with the same period NN, with x[n]⟷FSakx[n] \overset{FS}{\longleftrightarrow} a_k and y[n]⟷FSbky[n] \overset{FS}{\longleftrightarrow} b_k. Their periodic convolution

z[n]=∑r=⟨N⟩x[r] y[n−r]z[n] = \sum_{r=\langle N\rangle}x[r]\,y[n - r]

is periodic with period NN, and

z[n]⟷FSck=Nakbkz[n] \overset{FS}{\longleftrightarrow} c_k = N a_k b_k

Convolution in time becomes multiplication of coefficients.

Proof

With ω0=2π/N\omega_0 = 2\pi/N:

ck=1N∑n=⟨N⟩z[n]e−jkω0n=1N∑n=⟨N⟩∑r=⟨N⟩x[r] y[n−r] e−jkω0n=1N∑r=⟨N⟩x[r] e−jkω0r∑n=⟨N⟩y[n−r] e−jkω0(n−r)\begin{aligned} c_k &= \frac{1}{N}\sum_{n=\langle N\rangle}z[n]e^{-jk\omega_0 n} = \frac{1}{N}\sum_{n=\langle N\rangle}\sum_{r=\langle N\rangle}x[r]\,y[n-r]\,e^{-jk\omega_0 n} \\ &= \frac{1}{N}\sum_{r=\langle N\rangle}x[r]\,e^{-jk\omega_0 r}\sum_{n=\langle N\rangle}y[n-r]\,e^{-jk\omega_0(n-r)} \end{aligned}

Put m=n−rm = n - r; since yy is periodic, mm still runs over one full period, so the inner sum is ∑m=⟨N⟩y[m]e−jkω0m=Nbk\sum_{m=\langle N\rangle}y[m]e^{-jk\omega_0 m} = N b_k. Then

ck=1N∑r=⟨N⟩x[r]e−jkω0r⋅Nbk=N ak bkc_k = \frac{1}{N}\sum_{r=\langle N\rangle}x[r]e^{-jk\omega_0 r}\cdot N b_k = N\,a_k\,b_k

Hence ∑r=⟨N⟩x[r]y[n−r]⟷FSNakbk\sum_{r=\langle N\rangle}x[r]y[n-r] \overset{FS}{\longleftrightarrow} N a_k b_k. (The dual property: multiplication x[n]y[n]x[n]y[n] gives the periodic convolution of coefficients ∑l=⟨N⟩albk−l\sum_{l=\langle N\rangle}a_l b_{k-l}.)

  • 2077 Chaitra · 2 marks

A continuous time Fourier Series is evaluated over infinite duration while a discrete time Fourier series is evaluated over one time period only. Why?

Answer

The difference comes from how many distinct harmonics each signal can have:

  • Continuous time: the exponentials ejkω0te^{jk\omega_0 t} are all different for every integer kk; a higher kk is always a faster oscillation. A CT periodic signal (e.g. a square wave with sharp jumps) may need infinitely many harmonics, so the synthesis sum runs from k=−∞k = -\infty to ∞\infty, and convergence must be examined (Dirichlet conditions).
  • Discrete time: ej(k+N)(2π/N)n=ejk(2π/N)nej2πn=ejk(2π/N)ne^{j(k+N)(2\pi/N)n} = e^{jk(2\pi/N)n}e^{j2\pi n} = e^{jk(2\pi/N)n}, because nn is an integer. Only NN distinct exponentials exist, so both the synthesis and analysis sums need only one period (NN terms) and the coefficients repeat, ak+N=aka_{k+N} = a_k. The DTFS is a finite sum, so it always exists and has no convergence problem.

In both cases the analysis is done over one period of the signal; the "infinite" range in CT refers to the infinite number of harmonic terms needed.

  • 2077 Chaitra · 2+2+2 marks

Obtain fourier series coefficients of the signal x[n] = 1 + 3cos(2πn/N + π/3), where N is the fundamental time period of x[n]. Also, plot magnitude and phase spectrum of the coefficients.

Answer

With fundamental period NN, ω0=2π/N\omega_0 = 2\pi/N. Expand the cosine using Euler's formula:

x[n]=1+3cos⁡(2πnN+π3)=1+32ejπ/3ej(2π/N)n+32e−jπ/3e−j(2π/N)n\begin{aligned} x[n] &= 1 + 3\cos\left(\frac{2\pi n}{N} + \frac{\pi}{3}\right) \\ &= 1 + \frac{3}{2}e^{j\pi/3}e^{j(2\pi/N)n} + \frac{3}{2}e^{-j\pi/3}e^{-j(2\pi/N)n} \end{aligned}

Fourier series coefficients

Comparing with x[n]=∑k=⟨N⟩akejk(2π/N)nx[n] = \sum_{k=\langle N\rangle}a_k e^{jk(2\pi/N)n}:

a0=1,a1=32ejπ/3=0.75+j1.299,a−1=32e−jπ/3=0.75−j1.299a_0 = 1, \qquad a_1 = \frac{3}{2}e^{j\pi/3} = 0.75 + j1.299, \qquad a_{-1} = \frac{3}{2}e^{-j\pi/3} = 0.75 - j1.299

All other coefficients in one period are zero, and ak+N=aka_{k+N} = a_k (so aN−1=a−1a_{N-1} = a_{-1}, aN+1=a1a_{N+1} = a_1).

kkaka_k∣ak∣\lvert a_k\rvert∠ak\angle a_k
−1-11.5e−jπ/31.5e^{-j\pi/3}1.5−60∘-60^\circ
0110∘0^\circ
11.5ejπ/31.5e^{j\pi/3}1.5+60∘+60^\circ

Magnitude spectrum

 |a_k|  (one period shown; repeats every N)
   1.5     |     |
     1     |  |  |
       -o--+--+--+--o--> k
       -2 -1  0  1  2

Phase spectrum

 angle a_k (deg)
    60           |
       -o--+--o--+--o--> k
   -60     |
       -2 -1  0  1  2

The magnitude is even and the phase odd in kk (real signal), and both repeat with period NN.

  • 2076 Bhadra · 8 marks

Derive the expression to compute Fourier series coefficients of exponential Fourier series. Explain Gibbs phenomenon.

Answer

Exponential Fourier series coefficients

A periodic signal x(t)x(t) with period TT and ω0=2π/T\omega_0 = 2\pi/T is written as

x(t)=∑k=−∞∞akejkω0t(synthesis)x(t) = \sum_{k=-\infty}^{\infty}a_k e^{jk\omega_0 t} \qquad \text{(synthesis)}

Multiply both sides by e−jnω0te^{-jn\omega_0 t} and integrate over one period TT:

∫0Tx(t)e−jnω0tdt=∑k=−∞∞ak∫0Tej(k−n)ω0tdt\int_0^T x(t)e^{-jn\omega_0 t}dt = \sum_{k=-\infty}^{\infty}a_k\int_0^T e^{j(k-n)\omega_0 t}dt

For k≠nk \ne n:

∫0Tej(k−n)ω0tdt=ej(k−n)ω0T−1j(k−n)ω0=ej2π(k−n)−1j(k−n)ω0=0\int_0^T e^{j(k-n)\omega_0 t}dt = \frac{e^{j(k-n)\omega_0 T} - 1}{j(k-n)\omega_0} = \frac{e^{j2\pi(k-n)} - 1}{j(k-n)\omega_0} = 0

For k=nk = n the integral is ∫0T1 dt=T\int_0^T 1\,dt = T. Therefore

∫0Tx(t)e−jnω0tdt=anT⇒ak=1T∫Tx(t)e−jkω0tdt(analysis)\int_0^T x(t)e^{-jn\omega_0 t}dt = a_n T \quad\Rightarrow\quad a_k = \frac{1}{T}\int_T x(t)e^{-jk\omega_0 t}dt \qquad \text{(analysis)}
  • a0=1T∫Tx(t)dta_0 = \frac{1}{T}\int_T x(t)dt is the dc value.
  • aka_k is complex: ∣ak∣|a_k| is the magnitude spectrum and ∠ak\angle a_k the phase spectrum.
  • For real x(t)x(t), a−k=ak∗a_{-k} = a_k^*.

Example: a square wave of height 1 for ∣t∣<T1|t| < T_1 in period TT gives a0=2T1/Ta_0 = 2T_1/T and ak=sin⁡(kω0T1)kπa_k = \frac{\sin(k\omega_0 T_1)}{k\pi}.

Gibbs phenomenon

When a signal with a jump discontinuity (e.g. a square wave) is approximated by a finite number of Fourier terms,

xN(t)=∑k=−NNakejkω0t,x_N(t) = \sum_{k=-N}^{N}a_k e^{jk\omega_0 t},

the partial sum shows ripples near the discontinuity, with an overshoot of about 9% of the jump height on each side of it.

Key features:

  • As NN increases, the ripples become narrower and crowd closer to the discontinuity, but the peak overshoot does not decrease; it stays at about 9% of the jump (about 1.09 times for a unit jump).
  • At the discontinuity itself, the series converges to the midpoint 12[x(t−)+x(t+)]\frac{1}{2}[x(t^-) + x(t^+)].
  • The energy of the error ∫T∣x(t)−xN(t)∣2dt→0\int_T|x(t) - x_N(t)|^2dt \to 0 as N→∞N \to \infty (mean-square convergence), because the area under the ripples shrinks.
  • It was explained by J. W. Gibbs (1899) after Michelson noticed the overshoot in a mechanical harmonic analyser.
 x_N(t) near a jump (square wave)
       overshoot ~9%
        _/\_/\/\/\________
       |
       |   finite N ripples
 ______|
  \/\/\
       ^ jump at t = T1

Practical meaning: truncating a Fourier series (or ideal filtering, which cuts off high harmonics) produces ringing near sharp edges; window functions are used to reduce it.

  • 2076 Bhadra · 5 marks

A periodic square wave as shown in figure below is defined over one period as: x(t) = A for |t| < T₁; 0 for T₁ < |t| < T/2. Determine the Fourier series coefficients for x(t). [Figure: periodic rectangular pulse train of height A and period T; each pulse is centred on t = 0, ±T, … and extends from −T₁ to T₁]

Answer

The signal has period TT, so ω0=2π/T\omega_0 = 2\pi/T. Over one period (−T/2-T/2 to T/2T/2) it equals AA only for ∣t∣<T1|t| < T_1.

ak=1T∫−T/2T/2x(t)e−jkω0tdt=1T∫−T1T1A e−jkω0tdta_k = \frac{1}{T}\int_{-T/2}^{T/2}x(t)e^{-jk\omega_0 t}dt = \frac{1}{T}\int_{-T_1}^{T_1}A\,e^{-jk\omega_0 t}dt

DC term (k=0k = 0)

a0=1T∫−T1T1A dt=2AT1Ta_0 = \frac{1}{T}\int_{-T_1}^{T_1}A\,dt = \frac{2AT_1}{T}

This is the average value (height times duty cycle).

For k≠0k \ne 0

ak=AT[e−jkω0t−jkω0]−T1T1=AT⋅ejkω0T1−e−jkω0T1jkω0=2Akω0Tsin⁡(kω0T1)=Asin⁡(kω0T1)kπ(ω0T=2π)\begin{aligned} a_k &= \frac{A}{T}\left[\frac{e^{-jk\omega_0 t}}{-jk\omega_0}\right]_{-T_1}^{T_1} = \frac{A}{T}\cdot\frac{e^{jk\omega_0 T_1} - e^{-jk\omega_0 T_1}}{jk\omega_0} \\ &= \frac{2A}{k\omega_0 T}\sin(k\omega_0 T_1) \\ &= \frac{A\sin(k\omega_0 T_1)}{k\pi} \qquad (\omega_0 T = 2\pi) \end{aligned}

Result

ak={2AT1T,k=0Asin⁡(kω0T1)kπ=2AT1T sinc(2kT1T),k≠0a_k = \begin{cases} \dfrac{2AT_1}{T}, & k = 0 \\[2mm] \dfrac{A\sin(k\omega_0 T_1)}{k\pi} = \dfrac{2AT_1}{T}\,\text{sinc}\left(\dfrac{2kT_1}{T}\right), & k \ne 0 \end{cases}

(with sinc(x)=sin⁡(πx)/(πx)\text{sinc}(x) = \sin(\pi x)/(\pi x)). The coefficients are real and even because x(t)x(t) is real and even.

Special case, 50% duty cycle (T=4T1T = 4T_1): ω0T1=π/2\omega_0 T_1 = \pi/2, so a0=A/2a_0 = A/2 and ak=Asin⁡(kπ/2)kπa_k = \frac{A\sin(k\pi/2)}{k\pi}: a±1=A/πa_{\pm1} = A/\pi, a±2=0a_{\pm2} = 0, a±3=−A/3πa_{\pm3} = -A/3\pi, a±5=A/5πa_{\pm5} = A/5\pi, i.e. even harmonics are zero.

 a_k / A  (T = 4T1)
   0.5                 |
 0.318              |  |  |
0.0637  |           |  |  |           |
       -+--o--+--o--+--+--+--o--+--o--+--> k
-0.106        |                 |
       -5 -4 -3 -2 -1  0  1  2  3  4  5

The envelope of aka_k follows a sinc shape; making the pulse narrower (smaller T1/TT_1/T) spreads the spectrum wider.

  • 2075 Bhadra · 2+5 marks

State Parseval's relation for discrete time Fourier series. Find the average power of discrete time periodic signal x[n] = sin(πn/3) using Fourier series coefficient.

Answer

Parseval's relation for the DTFS

For a periodic x[n]x[n] with period NN and coefficients aka_k:

1N∑n=⟨N⟩∣x[n]∣2=∑k=⟨N⟩∣ak∣2\frac{1}{N}\sum_{n=\langle N\rangle}|x[n]|^2 = \sum_{k=\langle N\rangle}|a_k|^2

The average power over one period equals the sum of the powers ∣ak∣2|a_k|^2 of the NN harmonic components.

Average power of x[n]=sin⁡(πn/3)x[n] = \sin(\pi n/3)

Period: ω0=π/3\omega_0 = \pi/3, ω02π=16\frac{\omega_0}{2\pi} = \frac{1}{6}, so N=6N = 6 and the fundamental is 2π/6=π/32\pi/6 = \pi/3.

Coefficients:

x[n]=12jej(2π/6)n−12je−j(2π/6)nx[n] = \frac{1}{2j}e^{j(2\pi/6)n} - \frac{1}{2j}e^{-j(2\pi/6)n} a1=12j=−j2,a−1=a5=−12j=j2a_1 = \frac{1}{2j} = -\frac{j}{2}, \qquad a_{-1} = a_5 = -\frac{1}{2j} = \frac{j}{2}

and a0=a2=a3=a4=0a_0 = a_2 = a_3 = a_4 = 0 in one period.

Power by Parseval:

P=∑k=05∣ak∣2=∣a1∣2+∣a5∣2=(12)2+(12)2=12P = \sum_{k=0}^{5}|a_k|^2 = |a_1|^2 + |a_5|^2 = \left(\frac{1}{2}\right)^2 + \left(\frac{1}{2}\right)^2 = \frac{1}{2}

Check in the time domain: x[n]x[n] for n=0…5n = 0 \dots 5 is 0,32,32,0,−32,−320, \frac{\sqrt3}{2}, \frac{\sqrt3}{2}, 0, -\frac{\sqrt3}{2}, -\frac{\sqrt3}{2}:

P=16(4×34)=36=0.5P = \frac{1}{6}\left(4\times\frac{3}{4}\right) = \frac{3}{6} = 0.5

Answer: average power P=0.5P = 0.5 W.

  • 2075 Bhadra · 6 marks

Show that fourier series coefficient of discrete time signal is periodic in nature.

Answer

For a discrete-time periodic signal x[n]x[n] with period NN, the DTFS coefficients are

ak=1N∑n=⟨N⟩x[n] e−jk(2π/N)na_k = \frac{1}{N}\sum_{n=\langle N\rangle}x[n]\,e^{-jk(2\pi/N)n}

To show: ak+N=aka_{k+N} = a_k

Replace kk by k+Nk + N:

ak+N=1N∑n=⟨N⟩x[n] e−j(k+N)(2π/N)n=1N∑n=⟨N⟩x[n] e−jk(2π/N)n e−j2πn\begin{aligned} a_{k+N} &= \frac{1}{N}\sum_{n=\langle N\rangle}x[n]\,e^{-j(k+N)(2\pi/N)n} \\ &= \frac{1}{N}\sum_{n=\langle N\rangle}x[n]\,e^{-jk(2\pi/N)n}\,e^{-j2\pi n} \end{aligned}

Since nn is an integer, e−j2πn=cos⁡(2πn)−jsin⁡(2πn)=1e^{-j2\pi n} = \cos(2\pi n) - j\sin(2\pi n) = 1. Hence

ak+N=1N∑n=⟨N⟩x[n] e−jk(2π/N)n=aka_{k+N} = \frac{1}{N}\sum_{n=\langle N\rangle}x[n]\,e^{-jk(2\pi/N)n} = a_k

Repeating the argument, ak+mN=aka_{k+mN} = a_k for any integer mm. So the sequence of DTFS coefficients is periodic in kk with period NN.

Reason

The root cause is that discrete-time complex exponentials whose frequencies differ by 2π2\pi are identical:

ϕk+N[n]=ej(k+N)(2π/N)n=ejk(2π/N)n=ϕk[n]\phi_{k+N}[n] = e^{j(k+N)(2\pi/N)n} = e^{jk(2\pi/N)n} = \phi_k[n]

So there are only NN distinct harmonics, and the "(k+N)(k+N)th harmonic" is the same signal as the kkth; its coefficient must be the same. This is why the synthesis equation sums over only NN consecutive kk:

x[n]=∑k=⟨N⟩akejk(2π/N)nx[n] = \sum_{k=\langle N\rangle}a_k e^{jk(2\pi/N)n}

This is different from the CT Fourier series, where every ejkω0te^{jk\omega_0 t} is distinct and the coefficients are not periodic.

Example

x[n]=cos⁡(πn/2)x[n] = \cos(\pi n/2), N=4N = 4: x[n]=12ej(2π/4)n+12e−j(2π/4)nx[n] = \frac{1}{2}e^{j(2\pi/4)n} + \frac{1}{2}e^{-j(2\pi/4)n}, so a1=a−1=12a_1 = a_{-1} = \frac12. Computing from the analysis equation with samples {1,0,−1,0}\{1, 0, -1, 0\}:

kk−3-3−2-2−1-1012345
aka_k0.500.500.500.500.5

The pattern {0,0.5,0,0.5}\{0, 0.5, 0, 0.5\} repeats every 4 values of kk, e.g. a5=a1a_5 = a_1 and a3=a−1a_3 = a_{-1}.

 a_k
   0.5  |     |     |     |     |
       -+--o--+--o--+--o--+--o--+--> k
       -3 -2 -1  0  1  2  3  4  5
  • 2075 Baisakh · 5 marks

Find Fourier series coefficients and associated phase of the signal with fundamental frequency ω₀: x(t) = 1 + sin(ω₀t) + 2cos(ω₀t) + cos(2ω₀t + π/4)

Answer

Write each term using Euler's formula, with fundamental frequency ω0\omega_0:

x(t)=1+12j[ejω0t−e−jω0t]+[ejω0t+e−jω0t]+12[ej(2ω0t+π/4)+e−j(2ω0t+π/4)]=1+(1+12j)ejω0t+(1−12j)e−jω0t+12ejπ/4ej2ω0t+12e−jπ/4e−j2ω0t\begin{aligned} x(t) &= 1 + \frac{1}{2j}\left[e^{j\omega_0 t} - e^{-j\omega_0 t}\right] + \left[e^{j\omega_0 t} + e^{-j\omega_0 t}\right] + \frac{1}{2}\left[e^{j(2\omega_0 t + \pi/4)} + e^{-j(2\omega_0 t + \pi/4)}\right] \\ &= 1 + \left(1 + \frac{1}{2j}\right)e^{j\omega_0 t} + \left(1 - \frac{1}{2j}\right)e^{-j\omega_0 t} + \frac{1}{2}e^{j\pi/4}e^{j2\omega_0 t} + \frac{1}{2}e^{-j\pi/4}e^{-j2\omega_0 t} \end{aligned}

Comparing with x(t)=∑kakejkω0tx(t) = \sum_k a_k e^{jk\omega_0 t} and using 12j=−j2\frac{1}{2j} = -\frac{j}{2}:

a0=1,a1=1−j2,a−1=1+j2,a2=12ejπ/4=24(1+j),a−2=12e−jπ/4a_0 = 1, \quad a_1 = 1 - \frac{j}{2}, \quad a_{-1} = 1 + \frac{j}{2}, \quad a_2 = \frac{1}{2}e^{j\pi/4} = \frac{\sqrt2}{4}(1 + j), \quad a_{-2} = \frac{1}{2}e^{-j\pi/4}

ak=0a_k = 0 for ∣k∣>2|k| > 2.

Magnitude and phase

∣a±1∣=12+0.52=52=1.118,∠a±1=∓tan⁡−1(0.5)=∓26.57∘|a_{\pm1}| = \sqrt{1^2 + 0.5^2} = \frac{\sqrt5}{2} = 1.118, \qquad \angle a_{\pm1} = \mp\tan^{-1}(0.5) = \mp 26.57^\circ
kkaka_k∣ak∣\lvert a_k\rvert∠ak\angle a_k
−2-20.354−j0.3540.354 - j0.3540.5−45∘-45^\circ (−π/4-\pi/4)
−1-11+j0.51 + j0.51.118+26.57∘+26.57^\circ (0.464 rad)
0110∘0^\circ
11−j0.51 - j0.51.118−26.57∘-26.57^\circ (−0.464-0.464 rad)
20.354+j0.3540.354 + j0.3540.5+45∘+45^\circ (π/4\pi/4)
 |a_k|
 1.118     |     |
     1     |  |  |
   0.5  |  |  |  |  |
       -+--+--+--+--+--> k
       -2 -1  0  1  2

 angle a_k (deg)
    45              |
 26.57     |        |
       -+--+--o--+--+--> k
-26.57  |        |
   -45  |
       -2 -1  0  1  2

As expected for a real signal, a−k=ak∗a_{-k} = a_k^*: the magnitude is even and the phase is odd.

  • 2074 Bhadra · 3 marks

Explain the Dirichlet conditions for convergence of Fourier series.

Answer

The Dirichlet conditions are a set of sufficient conditions on a periodic signal x(t)x(t) (period TT) which guarantee that its Fourier series converges to x(t)x(t) at every point of continuity.

  1. Absolute integrability over one period: the signal must be absolutely integrable over a period,
∫T∣x(t)∣ dt<∞\int_{T} |x(t)|\, dt < \infty

This makes every coefficient finite, since ∣ak∣≤1T∫T∣x(t)∣ dt|a_k| \le \frac{1}{T}\int_T |x(t)|\,dt.

  • Violating example: x(t)=1/tx(t) = 1/t for 0<t≤10 < t \le 1, repeated with period 1.
  1. Finite number of maxima and minima: in any one period, x(t)x(t) has only a finite number of maxima and minima (bounded variation).
    • Violating example: x(t)=sin⁡(2π/t)x(t) = \sin(2\pi/t) for 0<t≤10 < t \le 1.
  2. Finite number of discontinuities: in any finite interval there are only a finite number of discontinuities, and each one is finite (a finite jump).
    • Violating example: a staircase signal whose steps get halved again and again within one period, giving infinitely many jumps.

Result: if these conditions hold, the Fourier series equals x(t)x(t) wherever x(t)x(t) is continuous. At a jump it converges to the average of the left and right limits, 12[x(t−)+x(t+)]\frac{1}{2}[x(t^-) + x(t^+)]. Near a jump there is still an overshoot of about 9% (the Gibbs phenomenon).

The conditions are sufficient but not necessary. Signals that break them are mostly artificial, so almost every practical signal has a convergent Fourier series.

  • 2074 Bhadra · 4 marks

Find Fourier series coefficient of discrete time periodic signal x[n] = sin ω₀n where ω₀ = 2π/N. Plot the signal when N = 5.

Answer

Given x[n]=sin⁡ω0nx[n] = \sin \omega_0 n with ω0=2π/N\omega_0 = 2\pi/N, the signal is periodic with fundamental period NN. Its DTFS coefficients can be found by inspection, using Euler's formula.

Fourier series coefficients

x[n]=sin⁡(2πNn)=12jej(2π/N)n−12je−j(2π/N)nx[n] = \sin\left(\frac{2\pi}{N}n\right) = \frac{1}{2j}e^{j(2\pi/N)n} - \frac{1}{2j}e^{-j(2\pi/N)n}

Compare this with the synthesis equation x[n]=∑k=⟨N⟩akejk(2π/N)nx[n] = \sum_{k=\langle N\rangle} a_k e^{jk(2\pi/N)n}:

a1=12j=−j2,a−1=−12j=j2a_1 = \frac{1}{2j} = -\frac{j}{2}, \qquad a_{-1} = -\frac{1}{2j} = \frac{j}{2}

All other coefficients in one period are zero. DTFS coefficients repeat with period NN, so ak+N=aka_{k+N} = a_k:

  • a1=a1+N=a1−N=⋯=−j/2a_1 = a_{1+N} = a_{1-N} = \dots = -j/2
  • a−1=aN−1=a2N−1=⋯=j/2a_{-1} = a_{N-1} = a_{2N-1} = \dots = j/2

Over k=0,1,…,N−1k = 0, 1, \dots, N-1: a1=−j/2a_1 = -j/2, aN−1=j/2a_{N-1} = j/2, and every other ak=0a_k = 0.

For N = 5: a1=−j/2a_1 = -j/2, a4=j/2a_4 = j/2, and a0=a2=a3=0a_0 = a_2 = a_3 = 0, repeating every 5 values of kk. In magnitude form ∣a1∣=∣a4∣=1/2|a_1| = |a_4| = 1/2, with ∠a1=−π/2\angle a_1 = -\pi/2 and ∠a4=+π/2\angle a_4 = +\pi/2.

Plot of the signal for N = 5

x[n]=sin⁡(2πn/5)x[n] = \sin(2\pi n/5):

n01234
x[n]00.9510.588−0.588−0.951
 x[n]
  0.95 |     *              *
  0.59 |     |  *           |  *
       |     |  |           |  |
  0  --*-----+--+--+--+--*--+--+--+--+-- n
       0     1  2  3  4  5  6  7  8  9
 -0.59 |           |  |           |  |
 -0.95 |           *  *           *  *

The pattern repeats every 5 samples.

  • 2073 Magh · 6 marks

Derive the expression of continuous time exponential Fourier series x[n] = A sin ω₀n

Answer

The signal is read as the continuous-time sinusoid x(t)=Asin⁡ω0tx(t) = A\sin\omega_0 t (the question writes x[n]x[n], but asks for the continuous-time series). The derivation of the general coefficient formula comes first, then the formula is applied.

Exponential Fourier series and its coefficients

A periodic signal with period TT and ω0=2π/T\omega_0 = 2\pi/T can be written as a sum of harmonically related complex exponentials:

x(t)=∑k=−∞∞akejkω0tx(t) = \sum_{k=-\infty}^{\infty} a_k e^{jk\omega_0 t}

To find ana_n, multiply both sides by e−jnω0te^{-jn\omega_0 t} and integrate over one period:

∫0Tx(t)e−jnω0tdt=∑k=−∞∞ak∫0Tej(k−n)ω0tdt\int_0^T x(t)e^{-jn\omega_0 t}dt = \sum_{k=-\infty}^{\infty} a_k \int_0^T e^{j(k-n)\omega_0 t}dt

Orthogonality of the exponentials gives

∫0Tej(k−n)ω0tdt={T,k=n0,k≠n\int_0^T e^{j(k-n)\omega_0 t}dt = \begin{cases} T, & k = n \\ 0, & k \ne n \end{cases}

so only the k=nk = n term remains, and

ak=1T∫Tx(t)e−jkω0tdta_k = \frac{1}{T}\int_T x(t)e^{-jk\omega_0 t}dt

Applying it to x(t) = A sin ω₀t

Using Euler's identity:

x(t)=A ejω0t−e−jω0t2j=A2jejω0t−A2je−jω0tx(t) = A\,\frac{e^{j\omega_0 t} - e^{-j\omega_0 t}}{2j} = \frac{A}{2j}e^{j\omega_0 t} - \frac{A}{2j}e^{-j\omega_0 t}

Comparing with ∑akejkω0t\sum a_k e^{jk\omega_0 t}:

a1=A2j=−jA2,a−1=−A2j=jA2,ak=0 (k≠±1)a_1 = \frac{A}{2j} = -j\frac{A}{2}, \qquad a_{-1} = -\frac{A}{2j} = j\frac{A}{2}, \qquad a_k = 0 \ (k \ne \pm 1)

The integral formula gives the same result. For example,

a1=1T∫0TA ejω0t−e−jω0t2je−jω0tdt=A2j⋅1T[T−0]=A2ja_1 = \frac{1}{T}\int_0^T A\,\frac{e^{j\omega_0 t}-e^{-j\omega_0 t}}{2j}e^{-j\omega_0 t}dt = \frac{A}{2j}\cdot\frac{1}{T}\left[T - 0\right] = \frac{A}{2j}

Spectrum

  • Magnitude: ∣a1∣=∣a−1∣=A/2|a_1| = |a_{-1}| = A/2
  • Phase: ∠a1=−π/2\angle a_1 = -\pi/2, ∠a−1=+π/2\angle a_{-1} = +\pi/2
 |a_k|               angle(a_k)
  A/2  |     |        +pi/2 |
       |     |              |
 ------+--+--+--  k   --+---+---+--  k
      -1  0  1         -1   0   1
                             |
                       -pi/2 |

The spectrum has only two lines, at ±ω0\pm\omega_0. The DT signal Asin⁡ω0nA\sin\omega_0 n with ω0=2π/N\omega_0 = 2\pi/N gives the same coefficients a±1=±A/(2j)a_{\pm1} = \pm A/(2j), repeated every NN.

  • 2073 Magh · 5 marks

Find Fourier series coefficients of signal and plot the coefficients. [The signal is not printed on the paper.]

Answer

The paper does not print the signal. This answer takes the standard textbook case: a periodic square wave of period TT with x(t)=1x(t) = 1 for ∣t∣<T1|t| < T_1 and x(t)=0x(t) = 0 for T1<∣t∣<T/2T_1 < |t| < T/2.

 x(t)
  1     ___       ___       ___
       |   |     |   |     |   |
 ______|   |_____|   |_____|   |___ t
     -T1   T1  T-T1  T+T1
       (pulse centred at t = 0)

Coefficients

With ω0=2π/T\omega_0 = 2\pi/T, the DC term is

a0=1T∫−T1T11 dt=2T1Ta_0 = \frac{1}{T}\int_{-T_1}^{T_1} 1\,dt = \frac{2T_1}{T}

For k≠0k \ne 0:

ak=1T∫−T1T1e−jkω0tdt=1T[e−jkω0t−jkω0]−T1T1=2kω0T⋅ejkω0T1−e−jkω0T12j=sin⁡(kω0T1)kπ\begin{aligned} a_k &= \frac{1}{T}\int_{-T_1}^{T_1} e^{-jk\omega_0 t}dt = \frac{1}{T}\left[\frac{e^{-jk\omega_0 t}}{-jk\omega_0}\right]_{-T_1}^{T_1} \\ &= \frac{2}{k\omega_0 T}\cdot\frac{e^{jk\omega_0 T_1} - e^{-jk\omega_0 T_1}}{2j} = \frac{\sin(k\omega_0 T_1)}{k\pi} \end{aligned}

The coefficients are real, because x(t)x(t) is real and even.

Example values for T = 4T₁ (50% duty cycle)

Here ω0T1=π/2\omega_0 T_1 = \pi/2, so ak=sin⁡(kπ/2)kπa_k = \frac{\sin(k\pi/2)}{k\pi}:

k0±1±2±3±4±5
aka_k1/21/π1/\pi ≈ 0.3180−1/(3π)-1/(3\pi) ≈ −0.10601/(5π)1/(5\pi) ≈ 0.064

Plot of the coefficients

 a_k
  0.50                   |
  0.32                |  |  |
  0.06    |           |  |  |           |
  0  -----+--+--+--+--+--+--+--+--+--+--+--- k
         -5 -4 -3 -2 -1  0  1  2  3  4  5
 -0.11          |                 |

The coefficients are samples, at ω=kω0\omega = k\omega_0, of the sinc-shaped envelope ω0π⋅sin⁡(ωT1)ω\frac{\omega_0}{\pi}\cdot\frac{\sin(\omega T_1)}{\omega}. Even harmonics vanish for a 50% duty cycle, and the odd ones decrease as 1/k1/k with alternating sign.

  • 2073 Magh · 3+3 marks

State and prove the frequency shifting and convolution properties for discrete time Fourier series pair.

Answer

Let x[n]x[n] and y[n]y[n] be periodic with the same period NN, and ω0=2π/N\omega_0 = 2\pi/N, with x[n]↔FSakx[n] \xleftrightarrow{FS} a_k and y[n]↔FSbky[n] \xleftrightarrow{FS} b_k, where

ak=1N∑n=⟨N⟩x[n]e−jkω0na_k = \frac{1}{N}\sum_{n=\langle N\rangle} x[n]e^{-jk\omega_0 n}

Frequency shifting property

Statement: multiplying a signal by a harmonic exponential shifts its coefficients:

ejMω0nx[n]↔FSak−Me^{jM\omega_0 n}x[n] \xleftrightarrow{FS} a_{k-M}

Proof: let z[n]=ejM(2π/N)nx[n]z[n] = e^{jM(2\pi/N)n}x[n]. Its coefficients are

ck=1N∑n=⟨N⟩ejMω0nx[n]e−jkω0n=1N∑n=⟨N⟩x[n]e−j(k−M)ω0n=ak−M\begin{aligned} c_k &= \frac{1}{N}\sum_{n=\langle N\rangle} e^{jM\omega_0 n}x[n]e^{-jk\omega_0 n} \\ &= \frac{1}{N}\sum_{n=\langle N\rangle} x[n]e^{-j(k-M)\omega_0 n} = a_{k-M} \end{aligned}

So multiplying by ejMω0ne^{jM\omega_0 n} moves the whole spectrum MM places to the right. This is the DT form of modulation.

Convolution property (periodic convolution)

Statement: periodic convolution in time corresponds to multiplication of the coefficients, scaled by NN:

∑r=⟨N⟩x[r] y[n−r]↔FSNakbk\sum_{r=\langle N\rangle} x[r]\,y[n-r] \xleftrightarrow{FS} N a_k b_k

Proof: let z[n]=∑r=⟨N⟩x[r]y[n−r]z[n] = \sum_{r=\langle N\rangle} x[r]y[n-r], which is also periodic with period NN. Then

ck=1N∑n=⟨N⟩∑r=⟨N⟩x[r] y[n−r] e−jkω0n=1N∑r=⟨N⟩x[r]e−jkω0r∑n=⟨N⟩y[n−r]e−jkω0(n−r)\begin{aligned} c_k &= \frac{1}{N}\sum_{n=\langle N\rangle}\sum_{r=\langle N\rangle} x[r]\,y[n-r]\,e^{-jk\omega_0 n} \\ &= \frac{1}{N}\sum_{r=\langle N\rangle} x[r]e^{-jk\omega_0 r}\sum_{n=\langle N\rangle} y[n-r]e^{-jk\omega_0 (n-r)} \end{aligned}

Put m=n−rm = n - r. Because y[m]e−jkω0my[m]e^{-jk\omega_0 m} is periodic with period NN, the inner sum over any NN consecutive mm is

∑m=⟨N⟩y[m]e−jkω0m=Nbk\sum_{m=\langle N\rangle} y[m]e^{-jk\omega_0 m} = N b_k

Therefore

ck=1N⋅Nbk∑r=⟨N⟩x[r]e−jkω0r=Nbk⋅ak=Nakbkc_k = \frac{1}{N}\cdot N b_k\sum_{r=\langle N\rangle} x[r]e^{-jk\omega_0 r} = N b_k\cdot a_k = N a_k b_k

Dual (multiplication) property: x[n]y[n]↔FS∑l=⟨N⟩albk−lx[n]y[n] \xleftrightarrow{FS} \sum_{l=\langle N\rangle} a_l b_{k-l}, which is periodic convolution of the coefficients.

  • 2072 Magh · 5 marks

Find the Fourier Series coefficient of the following continuous time periodic signal x(t) = 1.5 for 0 ≤ t < 1; −1.5 for 1 ≤ t < 2, with fundamental frequency π.

Answer

The fundamental frequency is ω0=π\omega_0 = \pi, so the period is T=2π/ω0=2T = 2\pi/\omega_0 = 2. Over one period, x(t)=1.5x(t) = 1.5 on 0≤t<10 \le t < 1 and x(t)=−1.5x(t) = -1.5 on 1≤t<21 \le t < 2.

 x(t)
  1.5  ____      ____      ____
      |    |    |    |    |
 -----+----+----+----+----+----- t
      0    1    2    3    4
           |____|    |____|
 -1.5

DC coefficient

a0=1T∫02x(t) dt=12[1.5(1)−1.5(1)]=0a_0 = \frac{1}{T}\int_0^2 x(t)\,dt = \frac{1}{2}\left[1.5(1) - 1.5(1)\right] = 0

Coefficients for k ≠ 0

ak=1T∫0Tx(t)e−jkω0tdt=12[∫011.5 e−jkπtdt−∫121.5 e−jkπtdt]a_k = \frac{1}{T}\int_0^T x(t)e^{-jk\omega_0 t}dt = \frac{1}{2}\left[\int_0^1 1.5\,e^{-jk\pi t}dt - \int_1^2 1.5\,e^{-jk\pi t}dt\right]

Evaluate each integral, using e−jkπ=(−1)ke^{-jk\pi} = (-1)^k and e−j2kπ=1e^{-j2k\pi} = 1:

∫01e−jkπtdt=1−e−jkπjkπ=1−(−1)kjkπ∫12e−jkπtdt=e−jkπ−e−j2kπjkπ=(−1)k−1jkπ\begin{aligned} \int_0^1 e^{-jk\pi t}dt &= \frac{1 - e^{-jk\pi}}{jk\pi} = \frac{1-(-1)^k}{jk\pi} \\ \int_1^2 e^{-jk\pi t}dt &= \frac{e^{-jk\pi} - e^{-j2k\pi}}{jk\pi} = \frac{(-1)^k - 1}{jk\pi} \end{aligned}

Substituting:

ak=1.52⋅[1−(−1)k]−[(−1)k−1]jkπ=1.52⋅2[1−(−1)k]jkπ=1.5 [1−(−1)k]jkπ\begin{aligned} a_k &= \frac{1.5}{2}\cdot\frac{[1-(-1)^k] - [(-1)^k - 1]}{jk\pi} \\ &= \frac{1.5}{2}\cdot\frac{2[1-(-1)^k]}{jk\pi} = \frac{1.5\,[1-(-1)^k]}{jk\pi} \end{aligned}

So

ak={3jkπ=−j3kπ,k odd0,k evena_k = \begin{cases} \dfrac{3}{jk\pi} = -j\dfrac{3}{k\pi}, & k \text{ odd} \\ 0, & k \text{ even} \end{cases}
k±1±2±3±5
∣ak∣\lvert a_k\rvert3/π3/\pi ≈ 0.95501/π1/\pi ≈ 0.3183/(5π)3/(5\pi) ≈ 0.191
∠ak\angle a_k∓π/2—∓π/2∓π/2

The coefficients are purely imaginary and odd (a−k=−aka_{-k} = -a_k). This is expected, because x(t)x(t) is real and has odd symmetry about t=0t = 0 (half-wave symmetry also removes the even harmonics). In trigonometric form:

x(t)=6π[sin⁡πt+13sin⁡3πt+15sin⁡5πt+… ]x(t) = \frac{6}{\pi}\left[\sin\pi t + \frac{1}{3}\sin 3\pi t + \frac{1}{5}\sin 5\pi t + \dots\right]

Answer: a0=0a_0 = 0, ak=−j 3/(kπ)a_k = -j\,3/(k\pi) for odd kk, and ak=0a_k = 0 for even kk.

  • 2072 Magh · 2+4 marks

State and prove the convolution property of continuous-time Fourier series.

Answer

Statement

Let x(t)x(t) and y(t)y(t) be periodic with the same period TT (ω0=2π/T\omega_0 = 2\pi/T), with x(t)↔FSakx(t) \xleftrightarrow{FS} a_k and y(t)↔FSbky(t) \xleftrightarrow{FS} b_k. The periodic convolution of the two signals has Fourier series coefficients TakbkT a_k b_k:

z(t)=∫Tx(τ) y(t−τ) dτ  ↔FS  ck=Takbkz(t) = \int_T x(\tau)\,y(t-\tau)\,d\tau \;\xleftrightarrow{FS}\; c_k = T a_k b_k

So convolution in time becomes multiplication of the coefficients.

Proof

z(t)z(t) is periodic with period TT, because y(t−τ)y(t - \tau) is periodic in tt. Its coefficients are

ck=1T∫Tz(t)e−jkω0tdt=1T∫T[∫Tx(τ)y(t−τ)dτ]e−jkω0tdtc_k = \frac{1}{T}\int_T z(t)e^{-jk\omega_0 t}dt = \frac{1}{T}\int_T\left[\int_T x(\tau)y(t-\tau)d\tau\right]e^{-jk\omega_0 t}dt

Swap the order of integration, and write e−jkω0t=e−jkω0τe−jkω0(t−τ)e^{-jk\omega_0 t} = e^{-jk\omega_0 \tau}e^{-jk\omega_0 (t-\tau)}:

ck=1T∫Tx(τ)e−jkω0τ[∫Ty(t−τ)e−jkω0(t−τ)dt]dτc_k = \frac{1}{T}\int_T x(\tau)e^{-jk\omega_0 \tau}\left[\int_T y(t-\tau)e^{-jk\omega_0 (t-\tau)}dt\right]d\tau

Put σ=t−τ\sigma = t - \tau. The integrand is periodic with period TT, so integrating over any interval of length TT gives the same result:

∫Ty(σ)e−jkω0σdσ=Tbk\int_T y(\sigma)e^{-jk\omega_0 \sigma}d\sigma = T b_k

Therefore

ck=1T⋅Tbk∫Tx(τ)e−jkω0τdτ=Tbk⋅akc_k = \frac{1}{T}\cdot T b_k\int_T x(\tau)e^{-jk\omega_0 \tau}d\tau = T b_k\cdot a_k ck=Takbk\boxed{c_k = T a_k b_k}

Use: the output of an LTI system for a periodic input can be found by multiplying spectra instead of convolving. The dual result is the multiplication property: x(t)y(t)↔FS∑lalbk−lx(t)y(t) \xleftrightarrow{FS} \sum_{l} a_l b_{k-l}.

  • 2072 Asoj · 6 marks

State and prove the linearity and time shift properties of discrete time Fourier series representation.

Answer

Let x[n]x[n] and y[n]y[n] be periodic with period NN (ω0=2π/N\omega_0 = 2\pi/N), with x[n]↔FSakx[n] \xleftrightarrow{FS} a_k and y[n]↔FSbky[n] \xleftrightarrow{FS} b_k, where

ak=1N∑n=⟨N⟩x[n]e−jkω0na_k = \frac{1}{N}\sum_{n=\langle N\rangle} x[n]e^{-jk\omega_0 n}

Linearity

Statement: for any constants AA and BB,

z[n]=Ax[n]+By[n]↔FSck=Aak+Bbkz[n] = A x[n] + B y[n] \xleftrightarrow{FS} c_k = A a_k + B b_k

Proof: z[n]z[n] is also periodic with period NN. Then

ck=1N∑n=⟨N⟩(Ax[n]+By[n])e−jkω0n=A⋅1N∑n=⟨N⟩x[n]e−jkω0n+B⋅1N∑n=⟨N⟩y[n]e−jkω0n=Aak+Bbk\begin{aligned} c_k &= \frac{1}{N}\sum_{n=\langle N\rangle}\left(A x[n] + B y[n]\right)e^{-jk\omega_0 n} \\ &= A\cdot\frac{1}{N}\sum_{n=\langle N\rangle} x[n]e^{-jk\omega_0 n} + B\cdot\frac{1}{N}\sum_{n=\langle N\rangle} y[n]e^{-jk\omega_0 n} \\ &= A a_k + B b_k \end{aligned}

Example: if x[n]=cos⁡(2πn/N)x[n] = \cos(2\pi n/N) (a±1=1/2a_{\pm1} = 1/2) and y[n]=1y[n] = 1 (b0=1b_0 = 1), then 3x[n]+2y[n]3x[n] + 2y[n] has c0=2c_0 = 2 and c±1=3/2c_{\pm1} = 3/2.

Time shifting

Statement: a delay of n0n_0 samples multiplies each coefficient by a linear-phase term:

x[n−n0]↔FSake−jkω0n0=ake−jk(2π/N)n0x[n - n_0] \xleftrightarrow{FS} a_k e^{-jk\omega_0 n_0} = a_k e^{-jk(2\pi/N)n_0}

Proof: let z[n]=x[n−n0]z[n] = x[n - n_0]. Its coefficients are

ck=1N∑n=⟨N⟩x[n−n0]e−jkω0nc_k = \frac{1}{N}\sum_{n=\langle N\rangle} x[n-n_0]e^{-jk\omega_0 n}

Put m=n−n0m = n - n_0, so n=m+n0n = m + n_0. As nn runs over one period, so does mm:

ck=1N∑m=⟨N⟩x[m]e−jkω0(m+n0)=e−jkω0n0⋅1N∑m=⟨N⟩x[m]e−jkω0m=e−jkω0n0ak\begin{aligned} c_k &= \frac{1}{N}\sum_{m=\langle N\rangle} x[m]e^{-jk\omega_0 (m+n_0)} \\ &= e^{-jk\omega_0 n_0}\cdot\frac{1}{N}\sum_{m=\langle N\rangle} x[m]e^{-jk\omega_0 m} = e^{-jk\omega_0 n_0}a_k \end{aligned}

Meaning: ∣ck∣=∣ak∣|c_k| = |a_k|, so a time shift does not change the magnitude spectrum. It only adds a phase of −kω0n0-k\omega_0 n_0, which is linear in kk.

  • 2071 Magh · 5+1 marks

Find Fourier series coefficients for periodic rectangular pulses with unity amplitude. Draw the magnitude spectrum.

Answer

Take a periodic train of rectangular pulses with unity amplitude, period TT and pulse width 2T12T_1, centred at t=0t = 0:

x(t)={1,∣t∣<T10,T1<∣t∣<T/2x(t+T)=x(t),ω0=2πTx(t) = \begin{cases} 1, & |t| < T_1 \\ 0, & T_1 < |t| < T/2 \end{cases} \qquad x(t+T) = x(t), \quad \omega_0 = \frac{2\pi}{T}
 x(t)
  1     ___       ___       ___
       |   |     |   |     |   |
 ______|   |_____|   |_____|   |___ t
     -T1   T1  T-T1  T+T1
       (width 2T1, period T)

Fourier series coefficients

DC term (the average value):

a0=1T∫−T1T1dt=2T1Ta_0 = \frac{1}{T}\int_{-T_1}^{T_1} dt = \frac{2T_1}{T}

For k≠0k \ne 0:

ak=1T∫−T1T1e−jkω0tdt=1T⋅e−jkω0T1−ejkω0T1−jkω0=2kω0Tsin⁡(kω0T1)=sin⁡(kω0T1)kπ\begin{aligned} a_k &= \frac{1}{T}\int_{-T_1}^{T_1} e^{-jk\omega_0 t}dt = \frac{1}{T}\cdot\frac{e^{-jk\omega_0 T_1} - e^{jk\omega_0 T_1}}{-jk\omega_0} \\ &= \frac{2}{k\omega_0 T}\sin(k\omega_0 T_1) = \frac{\sin(k\omega_0 T_1)}{k\pi} \end{aligned}

This can also be written ak=2T1T sinc ⁣(2kT1T)a_k = \frac{2T_1}{T}\,\mathrm{sinc}\!\left(\frac{2kT_1}{T}\right), with sinc(x)=sin⁡πxπx\mathrm{sinc}(x) = \frac{\sin \pi x}{\pi x}.

Example: T = 4T₁ (duty cycle 50%)

ak=sin⁡(kπ/2)kπa_k = \frac{\sin(k\pi/2)}{k\pi}: a0=0.5a_0 = 0.5, a±1=1/π≈0.318a_{\pm1} = 1/\pi \approx 0.318, a±2=0a_{\pm2} = 0, a±3=−1/(3π)≈−0.106a_{\pm3} = -1/(3\pi) \approx -0.106, a±5≈0.064a_{\pm5} \approx 0.064.

Magnitude spectrum

 |a_k|
 0.5              |
 0.32          |  |  |
 0.11       |  |  |  |  |
 0.06    |  |  |  |  |  |  |
 --------+--+--+--+--+--+--+--- k
        -5 -3 -1  0  1  3  5
     (|a_k| = 0 for even k != 0)

The lines sit at kω0k\omega_0 under a sinc-shaped envelope. A wider pulse (T1T_1 larger) gives a narrower envelope. A longer period (TT larger) gives more closely spaced lines, and as T→∞T \to \infty they merge into the continuous Fourier transform 2sin⁡(ωT1)/ω2\sin(\omega T_1)/\omega.

  • 2071 Bhadra · 6 marks

Discuss the following properties of continuous time Fourier series. (a) Time shifting (b) Time scaling (c) Conjugation

Answer

Let x(t)x(t) be periodic with period TT (ω0=2π/T\omega_0 = 2\pi/T) and x(t)↔FSakx(t) \xleftrightarrow{FS} a_k, where ak=1T∫Tx(t)e−jkω0tdta_k = \frac{1}{T}\int_T x(t)e^{-jk\omega_0 t}dt.

(a) Time shifting

x(t−t0)↔FSake−jkω0t0x(t - t_0) \xleftrightarrow{FS} a_k e^{-jk\omega_0 t_0}

Proof: put τ=t−t0\tau = t - t_0:

bk=1T∫Tx(t−t0)e−jkω0tdt=e−jkω0t01T∫Tx(τ)e−jkω0τdτ=e−jkω0t0akb_k = \frac{1}{T}\int_T x(t-t_0)e^{-jk\omega_0 t}dt = e^{-jk\omega_0 t_0}\frac{1}{T}\int_T x(\tau)e^{-jk\omega_0 \tau}d\tau = e^{-jk\omega_0 t_0}a_k

The magnitudes ∣bk∣=∣ak∣|b_k| = |a_k| stay the same. Only a linear phase −kω0t0-k\omega_0 t_0 is added.

(b) Time scaling

If α>0\alpha > 0, then x(αt)x(\alpha t) is periodic with period T/αT/\alpha and fundamental frequency αω0\alpha\omega_0. Its Fourier series is

x(αt)=∑k=−∞∞akejk(αω0)tx(\alpha t) = \sum_{k=-\infty}^{\infty} a_k e^{jk(\alpha\omega_0)t}

Proof: with period T′=T/αT' = T/\alpha, put τ=αt\tau = \alpha t:

bk=αT∫T/αx(αt)e−jk(αω0)tdt=1T∫Tx(τ)e−jkω0τdτ=akb_k = \frac{\alpha}{T}\int_{T/\alpha} x(\alpha t)e^{-jk(\alpha\omega_0)t}dt = \frac{1}{T}\int_T x(\tau)e^{-jk\omega_0\tau}d\tau = a_k

The coefficients do not change. What changes is the spacing of the harmonics, which becomes αω0\alpha\omega_0 instead of ω0\omega_0. Compressing a signal in time spreads its spectral lines further apart.

(c) Conjugation

x∗(t)↔FSa−k∗x^*(t) \xleftrightarrow{FS} a_{-k}^*

Proof: take the conjugate of the synthesis equation:

x∗(t)=∑kak∗e−jkω0t=∑ka−k∗ejkω0t(k→−k)x^*(t) = \sum_{k} a_k^* e^{-jk\omega_0 t} = \sum_{k} a_{-k}^* e^{jk\omega_0 t} \quad (k \to -k)

So the coefficients of x∗(t)x^*(t) are a−k∗a_{-k}^*.

Consequence (conjugate symmetry): if x(t)x(t) is real, then x(t)=x∗(t)x(t) = x^*(t), so ak=a−k∗a_k = a_{-k}^*. This gives ∣ak∣=∣a−k∣|a_k| = |a_{-k}| (even magnitude) and ∠a−k=−∠ak\angle a_{-k} = -\angle a_k (odd phase). If x(t)x(t) is also even, the aka_k are real and even.

  • 2070 Magh · 4+4 marks

How could you represent a signal x(t) with harmonically related exponentials? State and prove conjugation and conjugate symmetry property of CTFS.

Answer

Representation with harmonically related exponentials

Take a periodic signal x(t)x(t) with fundamental period TT and ω0=2π/T\omega_0 = 2\pi/T. The set of harmonically related complex exponentials

ϕk(t)=ejkω0t,k=0,±1,±2,…\phi_k(t) = e^{jk\omega_0 t}, \quad k = 0, \pm1, \pm2, \dots

are all periodic with period TT, since each frequency is an integer multiple of ω0\omega_0. The term for k=0k = 0 is the DC part, k=±1k = \pm1 are the fundamental components, and k=±Nk = \pm N are the NNth harmonics. A linear combination of them is periodic with period TT:

x(t)=∑k=−∞∞akejkω0t(synthesis equation)x(t) = \sum_{k=-\infty}^{\infty} a_k e^{jk\omega_0 t} \quad \text{(synthesis equation)}

The functions are orthogonal over one period: ∫Tejkω0te−jnω0tdt\int_T e^{jk\omega_0 t}e^{-jn\omega_0 t}dt equals TT when k=nk = n and 0 otherwise. Multiplying the synthesis equation by e−jnω0te^{-jn\omega_0 t} and integrating over TT therefore leaves only one term. This gives

ak=1T∫Tx(t)e−jkω0tdt(analysis equation)a_k = \frac{1}{T}\int_T x(t)e^{-jk\omega_0 t}dt \quad \text{(analysis equation)}

The series converges to x(t)x(t) when x(t)x(t) meets the Dirichlet conditions. Example: cos⁡ω0t=12ejω0t+12e−jω0t\cos\omega_0 t = \frac{1}{2}e^{j\omega_0 t} + \frac{1}{2}e^{-j\omega_0 t}, so a±1=1/2a_{\pm1} = 1/2.

Conjugation property

Statement: if x(t)↔FSakx(t) \xleftrightarrow{FS} a_k, then

x∗(t)↔FSa−k∗x^*(t) \xleftrightarrow{FS} a_{-k}^*

Proof: let bkb_k be the coefficients of x∗(t)x^*(t):

bk=1T∫Tx∗(t)e−jkω0tdt=[1T∫Tx(t)ejkω0tdt]∗=[1T∫Tx(t)e−j(−k)ω0tdt]∗=a−k∗b_k = \frac{1}{T}\int_T x^*(t)e^{-jk\omega_0 t}dt = \left[\frac{1}{T}\int_T x(t)e^{jk\omega_0 t}dt\right]^* = \left[\frac{1}{T}\int_T x(t)e^{-j(-k)\omega_0 t}dt\right]^* = a_{-k}^*

Conjugate symmetry property

Statement: if x(t)x(t) is real, its coefficients are conjugate symmetric:

a−k=ak∗a_{-k} = a_k^*

Proof: for a real signal x(t)=x∗(t)x(t) = x^*(t), so the two signals have the same coefficients: ak=bk=a−k∗a_k = b_k = a_{-k}^*. Taking the conjugate of both sides gives ak∗=a−ka_k^* = a_{-k}.

Consequences for real x(t)x(t):

  • ∣a−k∣=∣ak∣|a_{-k}| = |a_k|, so the magnitude spectrum is even.
  • ∠a−k=−∠ak\angle a_{-k} = -\angle a_k, so the phase spectrum is odd.
  • Re{ak}\mathrm{Re}\{a_k\} is even and Im{ak}\mathrm{Im}\{a_k\} is odd.
  • Real and even x(t)x(t) gives real, even aka_k. Real and odd x(t)x(t) gives purely imaginary, odd aka_k.

For example, sin⁡ω0t\sin\omega_0 t has a1=1/(2j)a_1 = 1/(2j) and a−1=−1/(2j)=a1∗a_{-1} = -1/(2j) = a_1^*.

  • 2070 Bhadra · 1+6 marks

What is the information provided by Fourier series coefficients of a signal? State and prove time shifting and conjugation properties of continuous time Fourier series representation.

Answer

Information provided by Fourier series coefficients

The coefficient aka_k gives the amount (magnitude ∣ak∣|a_k|) and the phase (∠ak\angle a_k) of the kkth harmonic ejkω0te^{jk\omega_0 t} in the signal. This makes the coefficients the frequency-domain picture of the signal. a0a_0 is the DC (average) value, ∣ak∣2|a_k|^2 is the power at frequency kω0k\omega_0 (by Parseval, P=∑∣ak∣2P = \sum|a_k|^2), and the spread of the aka_k shows the bandwidth of the signal.

Time shifting property

Statement: if x(t)↔FSakx(t) \xleftrightarrow{FS} a_k with period TT and ω0=2π/T\omega_0 = 2\pi/T, then

x(t−t0)↔FSe−jkω0t0akx(t - t_0) \xleftrightarrow{FS} e^{-jk\omega_0 t_0}a_k

Proof: the coefficients of y(t)=x(t−t0)y(t) = x(t - t_0) are

bk=1T∫Tx(t−t0)e−jkω0tdtb_k = \frac{1}{T}\int_T x(t-t_0)e^{-jk\omega_0 t}dt

Put τ=t−t0\tau = t - t_0 (so dt=dτdt = d\tau). The range of τ\tau is still one full period:

bk=1T∫Tx(τ)e−jkω0(τ+t0)dτ=e−jkω0t0⋅1T∫Tx(τ)e−jkω0τdτ=e−jkω0t0ak\begin{aligned} b_k &= \frac{1}{T}\int_T x(\tau)e^{-jk\omega_0(\tau + t_0)}d\tau \\ &= e^{-jk\omega_0 t_0}\cdot\frac{1}{T}\int_T x(\tau)e^{-jk\omega_0 \tau}d\tau = e^{-jk\omega_0 t_0}a_k \end{aligned}

Interpretation: ∣bk∣=∣ak∣|b_k| = |a_k|, so the magnitude spectrum does not change. The phase changes by −kω0t0-k\omega_0 t_0, which is proportional to the harmonic number.

Conjugation property

Statement:

x∗(t)↔FSa−k∗x^*(t) \xleftrightarrow{FS} a_{-k}^*

Proof: the coefficients of x∗(t)x^*(t) are

bk=1T∫Tx∗(t)e−jkω0tdt=[1T∫Tx(t)e+jkω0tdt]∗=[1T∫Tx(t)e−j(−k)ω0tdt]∗=a−k∗\begin{aligned} b_k &= \frac{1}{T}\int_T x^*(t)e^{-jk\omega_0 t}dt = \left[\frac{1}{T}\int_T x(t)e^{+jk\omega_0 t}dt\right]^* \\ &= \left[\frac{1}{T}\int_T x(t)e^{-j(-k)\omega_0 t}dt\right]^* = a_{-k}^* \end{aligned}

Consequence: if x(t)x(t) is real, then x(t)=x∗(t)x(t) = x^*(t), which gives a−k=ak∗a_{-k} = a_k^* (conjugate symmetry). The magnitude spectrum is then even and the phase spectrum is odd.

  • 2070 Bhadra · 5+2 marks

Find out Fourier series coefficients of a periodic discrete time signal described over a period as x[n] = 2 for |n| ≤ 1; 0 for 1 < |n| ≤ 3. Using the Fourier series coefficients calculated above, find the Fourier series coefficients of the signal e^(j4πn/7) x[n].

Answer

The signal is x[n]=2x[n] = 2 for n=−1,0,1n = -1, 0, 1 and x[n]=0x[n] = 0 for n=±2,±3n = \pm2, \pm3. One period runs over n=−3n = -3 to 33, so N=7N = 7 and ω0=2π/7\omega_0 = 2\pi/7.

 x[n]
  2        *  *  *              *  *  *
           |  |  |              |  |  |
--*--*--*--+--+--+--*--*--*--*--+--+--+--- n
 -4 -3 -2 -1  0  1  2  3  4  5  6  7  8

Fourier series coefficients of x[n]

ak=1N∑n=−112 e−jk(2π/7)n=27[1+e−j2πk/7+ej2πk/7]=27[1+2cos⁡2πk7]a_k = \frac{1}{N}\sum_{n=-1}^{1} 2\,e^{-jk(2\pi/7)n} = \frac{2}{7}\left[1 + e^{-j2\pi k/7} + e^{j2\pi k/7}\right] = \frac{2}{7}\left[1 + 2\cos\frac{2\pi k}{7}\right]

In closed form (the standard result for a DT rectangular pulse with N1=1N_1 = 1):

ak=27⋅sin⁡(3πk/7)sin⁡(πk/7),k≠0,±7,… ;a0=2(2N1+1)N=67a_k = \frac{2}{7}\cdot\frac{\sin(3\pi k/7)}{\sin(\pi k/7)}, \quad k \ne 0, \pm7, \dots; \qquad a_0 = \frac{2(2N_1+1)}{N} = \frac{6}{7}
k0123456
aka_k0.8570.6420.159−0.229−0.2290.1590.642

The coefficients are real and even (ak=a−k=a7−ka_k = a_{-k} = a_{7-k}), because x[n]x[n] is real and even. They repeat with period 7.

Coefficients of y[n] = e^{j4πn/7} x[n]

ej4πn/7=ej2(2π/7)n=ejMω0nwith M=2e^{j4\pi n/7} = e^{j2(2\pi/7)n} = e^{jM\omega_0 n} \quad \text{with } M = 2

By the frequency shifting property, ejMω0nx[n]↔FSak−Me^{jM\omega_0 n}x[n] \xleftrightarrow{FS} a_{k-M}, so

bk=ak−2=27⋅sin⁡(3π(k−2)/7)sin⁡(π(k−2)/7),b2=67b_k = a_{k-2} = \frac{2}{7}\cdot\frac{\sin\left(3\pi(k-2)/7\right)}{\sin\left(\pi(k-2)/7\right)}, \qquad b_2 = \frac{6}{7}
k0123456
bk=ak−2b_k = a_{k-2}0.1590.6420.8570.6420.159−0.229−0.229

Answer: ak=27sin⁡(3πk/7)sin⁡(πk/7)a_k = \frac{2}{7}\frac{\sin(3\pi k/7)}{\sin(\pi k/7)} with a0=6/7a_0 = 6/7, and bk=ak−2b_k = a_{k-2}. The spectrum is the same shape, moved 2 places to the right (peak at k=2k = 2).

  • 2069 Bhadra · 6 marks

Find the Fourier series representation of the signal x[n] = Σ_{ℓ=−∞}^{∞} δ[n − ℓN].

Answer

The signal is a periodic impulse train: a unit impulse at every multiple of NN.

x[n]=∑ℓ=−∞∞δ[n−ℓN]x[n] = \sum_{\ell=-\infty}^{\infty}\delta[n - \ell N]
 x[n]
  1  *              *              *
     |              |              |
 ----+--*--*--...---+--*--*--...---+---- n
     0              N              2N

It is periodic with period NN, so ω0=2π/N\omega_0 = 2\pi/N.

Coefficients

ak=1N∑n=⟨N⟩x[n] e−jk(2π/N)na_k = \frac{1}{N}\sum_{n=\langle N\rangle} x[n]\,e^{-jk(2\pi/N)n}

Choose the period n=0,1,…,N−1n = 0, 1, \dots, N-1. Within it x[n]=δ[n]x[n] = \delta[n], which is 1 at n=0n = 0 and 0 elsewhere:

ak=1N∑n=0N−1δ[n] e−jk(2π/N)n=1Ne0=1Nfor all ka_k = \frac{1}{N}\sum_{n=0}^{N-1}\delta[n]\,e^{-jk(2\pi/N)n} = \frac{1}{N}e^{0} = \frac{1}{N} \quad \text{for all } k

Fourier series representation

x[n]=∑k=0N−11Nejk(2π/N)n=1N∑k=⟨N⟩ejk(2π/N)nx[n] = \sum_{k=0}^{N-1}\frac{1}{N}e^{jk(2\pi/N)n} = \frac{1}{N}\sum_{k=\langle N\rangle} e^{jk(2\pi/N)n}

Check: for n=0,±N,…n = 0, \pm N, \dots every exponential equals 1, so the sum is N⋅1N=1N\cdot\frac{1}{N} = 1. For other nn the sum is a full geometric series over the NNth roots of unity, which is 1N⋅1−ej2πn1−ej2πn/N=0\frac{1}{N}\cdot\frac{1 - e^{j2\pi n}}{1 - e^{j2\pi n/N}} = 0. This matches x[n]x[n].

Spectrum: all NN harmonics have the same amplitude 1/N1/N and zero phase. The spectrum is flat. A train of narrow impulses contains all frequencies equally, just as the CT impulse train ∑δ(t−kT)\sum\delta(t - kT) has ak=1/Ta_k = 1/T. This result is the basis of the sampling theorem.

  • 2083 Baisakh (new course) · 4 marks

Describe the continuous-time Fourier series' time reversal, time scaling, and frequency shifting properties.

Answer

Let x(t)x(t) be periodic with period TT, ω0=2π/T\omega_0 = 2\pi/T, and x(t)↔FSakx(t) \xleftrightarrow{FS} a_k.

Time reversal

x(−t)↔FSa−kx(-t) \xleftrightarrow{FS} a_{-k}

Put k=−mk = -m in the synthesis equation:

x(−t)=∑kake−jkω0t=∑ma−mejmω0tx(-t) = \sum_k a_k e^{-jk\omega_0 t} = \sum_m a_{-m}e^{jm\omega_0 t}

Reversing the signal in time reverses its coefficient sequence. If x(t)x(t) is even, a−k=aka_{-k} = a_k. If it is odd, a−k=−aka_{-k} = -a_k.

Time scaling

x(αt), α>0:x(αt)=∑kakejk(αω0)tx(\alpha t),\ \alpha > 0: \quad x(\alpha t) = \sum_k a_k e^{jk(\alpha\omega_0)t}

x(αt)x(\alpha t) has period T/αT/\alpha and fundamental frequency αω0\alpha\omega_0. The coefficients aka_k stay the same, and only the harmonic spacing changes from ω0\omega_0 to αω0\alpha\omega_0. Compressing the signal (α>1\alpha > 1) spreads the spectral lines apart.

Frequency shifting

ejMω0tx(t)↔FSak−Me^{jM\omega_0 t}x(t) \xleftrightarrow{FS} a_{k-M}

The coefficients of the product are

1T∫TejMω0tx(t)e−jkω0tdt=1T∫Tx(t)e−j(k−M)ω0tdt=ak−M\frac{1}{T}\int_T e^{jM\omega_0 t}x(t)e^{-jk\omega_0 t}dt = \frac{1}{T}\int_T x(t)e^{-j(k-M)\omega_0 t}dt = a_{k-M}

Multiplying by a harmonic exponential moves the whole line spectrum MM places to the right. This is the dual of the time-shift property.

Questions from Old Question Collection (Signal Analysis EX 651) (IOE Signal Analysis (EX 651) exam papers from 2069 to 2082) and 2080 course paper (ENEX 255) (IOE Signals and Systems (ENEX 255, new course) papers: 2082 Bhadra, 2083 Baisakh, 2083 Bhadra). Answers are written for this site; check them against your class notes.

Chapter titles and hours from the IOE syllabus ↗