Skip to main content

Fourier Transform of some common analog signals


  • Fourier Transform of a delta function

If x(t) = \( \delta(t) \), then Fourier transform,

X(\( \omega \)) = \( \int_{-\infty}^{\infty} \delta(t)e^{-j\omega t} dt \)

= \( \int_{-\epsilon}^{\epsilon} \delta(t)e^{-j\omega t} dt \)

= \( \int_{-\epsilon}^{\epsilon} \delta(t)e^{-j\omega 0} dt \)

= 1

Thus, the Fourier transform of a delta/impulse is a constant equal to 1, independent of frequency. Remember that derivation used the shifting property of the impulse to eliminate the integral.

 

 

 

  • Fourier transform of a unit step function

 

\( \gamma(t) \) = 0 for t < 0

                                                            = 1 for t \( \ge \) 1

We know that a unit-step function is an integration of a delta function. So for a unit step function,

\( \gamma(t) \) = \( \int \delta(t) dt \)

So, X(\( \omega \)) = \( \frac{1}{j\omega} F(\delta(t)) + \pi F(\delta(t))\delta(\omega) \)

[See the property of integration above]

           = \( \frac{1}{j\omega} + \pi\delta(\omega) \)       [as F(\( \delta(t) \)) = 1]

When a function (i.e., x(t)) is not an energy function and hence the Fourier transform of [\( \int x(t) dt \)] includes an impulse function.

 

 

 

  • Fourier Transform of a rectangular / unit pulse function

A pulse function can be represented as,

x(t)=Π(t) = \( \gamma(t + \frac{1}{2}) - \gamma(t - \frac{1}{2}) \)

For a function rect(t) = Π(t) = 1 for |t| \( \le \frac{1}{2} \)

                                             = 0 otherwise

Given that

x(t) = Π(t)

Hence from the definition of the Fourier transform we have

F (Π(t)) = X(\( \omega \)) = \( \int_{-\infty}^{\infty} \Pi(t)e^{-j\omega t} dt \)

                           = \( \int_{-1/2}^{1/2} e^{-j\omega t} dt \)      [as Π(t) = 1 for |t| \( \le \frac{1}{2} \)]

                            = [(e –jωt)/-jω]-1/21/2

                                           = [e –jω/2 - e jω/2] / -jω

                            = [e jω/2 - e -jω/2] / jω

                            = 2/ω . {[e jω/2 - e -jω/2] / 2j}

                            = 2/ω . sin(ω/2)

                            = {sin(ω/2) / (ω/2)}

                           = {sin(\( \pi \)(ω/2\( \pi \))) / (\( \pi \)(ω/2\( \pi \)))}

                            = sinc(ω/2\( \pi \))

 

For the above case, the rectangular function has a pulse width value of 1 over the interval of [-½, ½]; 0 otherwise.

Now we’ll discuss a rectangular pulse that has a width of T

Then,    rect(t/T) = Π(t/T) = 1 for |t| \( \le \) T/2

                                             = 0 otherwise

Given that

x(t/T) = Π(t/T)

Hence from the definition of the Fourier transform we have

F (Π(t/T)) = \( \int_{-\infty}^{\infty} \Pi(t/T)e^{-j\omega t} dt \)

                           = \( \int_{-T/2}^{T/2} e^{-j\omega t} dt \)      [as Π(t/T) = 1 for |t| \( \le \) T/2]

                            = [(e –jωt)/-jω]-T/2T/2

                                           = [e –jωT/2 - e jωT/2] / -jω

                            = [e jωT/2 - e -jωT/2] / jω

                            = 2/ω . {[e jωT/2 - e -jωT/2] / 2j}

                            = 2/ω . sin(ω(T/2))

                            = {sin(ω(T/2)) / (ω/2)}

                            = {sin(\( \pi \)ωT/2\( \pi \))) / (\( \pi \)ω/2\( \pi \))}

                            = {sin(\( \pi \)ωT/2\( \pi \))) / (\( \pi \)ωT/2\( \pi \))}.T

                            = T. sinc(ωT/2\( \pi \))

 

 

 

 

  • Fourier Transform of a unit triangle pulse

A unit triangle pulse is simply the convolution of a unit pulse function with itself.

Here, Λ(t) = Π(t) * Π(t)               

[Π(t) is a unit pulse function & ‘*’ denotes convolution]

So, Λ(\( \omega \)) = sinc(ω/2\( \pi \)) . sinc(ω/2\( \pi \)) = sinc2(ω/2\( \pi \))

 

 

 

  • Fourier Transform of a Sawtooth function

1738659393_fourier-transform-40.png

s(t) = 0, for t < 0 and t > 1

                                                  = t, for 0 \( \le \) t \( \le \) 1

We can represent sawtooth as the integral of shifted unit pulse function (to give the ramp) and a negative impulse (delayed by one second) to give the discontinuity at the end of the ramp

1738659393_fourier-transform-41.png

 

s(t) = \( \int_{-\infty}^{t} \Pi(\tau)d\tau - \int_{-\infty}^{t} \delta(\tau-1)d\tau = \int_{-\infty}^{t} y(\tau)d\tau \)

y(t) = \( \Pi(t) - \delta(t-1) \)

Now, we’ve to find the Fourier transform of y(t),

Y(\( \omega \)) = sinc(ω/2\( \pi \))e-jω/2 - e-jω

We can now apply integral property with Y(0) = 0, to find S(\( \omega \))

S(\( \omega \)) = F(\( \int_{-\infty}^{t} y(\tau)d\tau \)) = Y(\( \omega \))/jω - \( \pi \)Y(0)\( \delta \)(0) = Y(\( \omega \))/jω

= {(sinc(ω/2\( \pi \))e-jω/2 - e-jω) / jω}

= {((sin(π . ω/2\( \pi \)) / (π . ω/2\( \pi \)))e-jω/2 - e-jω) / jω}

=  (((sin(Ï€ . ω/2\( \pi \)) / (Ï€ . ω/2\( \pi \)))e-jω/2) / jω) – (e-jω / jω)

= (2(sin(ω/2)e-jω/2) / jω2) – (je-jω / j2ω)

= (2((ejω/2 - e-jω/2) / 2j)e-jω/2) / jω2) + (je-jω / ω)     [as j2 = -1]

= (((ejω/2 - e-jω/2)e-jω/2) / j2ω2) + (je-jω / ω)

= (((e-jω/2 - ejω/2)e-jω/2) / ω2) + (je-jω / ω)

= (((e-jω/2 - ejω/2)e-jω/2) / ω2) + (je-jω / ω)

= ((((e-jω/2 - ejω/2)e-jω/2) + jωe-jω) / ω2)

= ((e-jω - 1 + jωe-jω) / ω2)

= ((e-jω(1+jω) - 1) / ω2


Further Reading

  1.  

People are good at skipping over material they already know!

View Related Topics to







Contact Us

Name

Email *

Message *

Popular Posts

Online Simulator for ASK, FSK, and PSK

Try our new Digital Signal Processing Simulator!   Start Simulator for binary ASK Modulation Message Bits (e.g. 1,0,1,0) Carrier Frequency (Hz) Sampling Frequency (Hz) Run Simulation Simulator for binary FSK Modulation Input Bits (e.g. 1,0,1,0) Freq for '1' (Hz) Freq for '0' (Hz) Sampling Rate (Hz) Visualize FSK Signal Simulator for BPSK Modulation ...

Constellation Diagrams of ASK, PSK, and FSK

📘 Overview of Energy per Bit (Eb / N0) 🧮 Online Simulator for constellation diagrams of ASK, FSK, and PSK 🧮 Theory behind Constellation Diagrams of ASK, FSK, and PSK 🧮 MATLAB Codes for Constellation Diagrams of ASK, FSK, and PSK 📚 Further Reading 📂 Other Topics on Constellation Diagrams of ASK, PSK, and FSK ... 🧮 Simulator for constellation diagrams of m-ary PSK 🧮 Simulator for constellation diagrams of m-ary QAM BASK (Binary ASK) Modulation: Transmits one of two signals: 0 or -√Eb, where Eb​ is the energy per bit. These signals represent binary 0 and 1.    BFSK (Binary FSK) Modulation: Transmits one of two signals: +√Eb​ ( On the y-axis, the phase shift of 90 degrees with respect to the x-axis, which is also termed phase offset ) or √Eb (on x-axis), where Eb​ is the energy per bit. These signals represent binary 0 and 1.  BPSK (Binary PSK) Modulation: Transmits one of two signals...

Power Spectral Density Calculation Using FFT in MATLAB

📘 Overview 🧮 Steps to calculate the PSD of a signal 🧮 MATLAB Codes 📚 Further Reading Power spectral density (PSD) tells us how the power of a signal is distributed across different frequency components, whereas Fourier Magnitude gives you the amplitude (or strength) of each frequency component in the signal. Steps to calculate the PSD of a signal Firstly, calculate the first Fourier transform (FFT) of a signal Then, calculate the Fourier magnitude of the signal The power spectrum is the square of the Fourier magnitude To calculate power spectrum density (PSD), divide the power spectrum by the total number of samples and the frequency resolution. {Frequency resolution = (sampling frequency / total number of samples)} Sampling frequency (fs): The rate at which the continuous-time signal is sampled (in Hz). ...

FFT Magnitude and Phase Spectrum using MATLAB

📘 Overview & Theory 🧮 MATLAB Code 1 🧮 MATLAB Code 2 📚 Further Reading   MATLAB Code  % Developed by SalimWireless.Com clc; clear; close all; % Configuration parameters fs = 10000; % Sampling rate (Hz) t = 0:1/fs:1-1/fs; % Time vector creation % Signal definition x = sin(2 * pi * 100 * t) + cos(2 * pi * 1000 * t); % Calculate the Fourier Transform y = fft(x); z = fftshift(y); % Create frequency vector ly = length(y); f = (-ly/2:ly/2-1) / ly * fs; % Calculate phase while avoiding numerical precision issues tol = 1e-6; % Tolerance threshold for zeroing small values z(abs(z) < tol) = 0; phase = angle(z); % Plot the original Signal figure; subplot(3, 1, 1); plot(t, x, 'b'); xlabel('Time (s)'); ylabel('|y|'); title('Original Messge Signal'); grid on; % Plot the magnitude of the Fourier Transform subplot(3, 1, 2); stem(f, abs(z), 'b'); xlabel('Frequency (Hz)'); ylabel('|y|'); title('Magnitude o...

MATLAB Code for Rms Delay Spread

RMS delay spread is crucial when you need to know how much the signal is dispersed in time due to multipath propagation, the spread (variance) around the average. In high-data-rate systems like LTE, 5G, or Wi-Fi, even small time dispersions can cause ISI. RMS delay spread is directly related to the amount of ISI in such systems. RMS Delay Spread [↗] Delay Spread Calculator Enter delays (ns) separated by commas: Enter powers (dB) separated by commas: Calculate   The above calculator Converts Power to Linear Scale: It correctly converts the power values from decibels (dB) to a linear scale. Calculates Mean Delay: It accurately computes the mean excess delay, which is the first moment of the power delay profile. Calculates RMS Delay Spread: It correctly calculates the RMS delay spread, defined as the square root of the second central moment of the power delay profile.   MATLAB Code  clc...

Coherence Bandwidth and Coherence Time

🧮 Coherence Bandwidth 🧮 Coherence Time 🧮 MATLAB Code s 📚 Further Reading For Doppler Delay or Multi-path Delay Coherence time T coh ∝ 1 / v max (For slow fading, coherence time T coh is greater than the signaling interval.) Coherence bandwidth W coh ∝ 1 / Ï„ max (For frequency-flat fading, coherence bandwidth W coh is greater than the signaling bandwidth.) Where: T coh = coherence time W coh = coherence bandwidth v max = maximum Doppler frequency (or maximum Doppler shift) Ï„ max = maximum excess delay (maximum time delay spread) Notes: The notation v max −1 and Ï„ max −1 indicate inverse proportionality. Doppler spread refers to the range of frequency shifts caused by relative motion, determining T coh . Delay spread (or multipath delay spread) determines W coh . Frequency-flat fading occurs when W coh is greater than the signaling bandwidth. Coherence Bandwidth Coherence bandwidth is...

MATLAB Code for ASK, FSK, and PSK

📘 Overview & Theory 🧮 MATLAB Code for ASK 🧮 MATLAB Code for FSK 🧮 MATLAB Code for PSK 🧮 Simulator for binary ASK, FSK, and PSK Modulations 📚 Further Reading ASK, FSK & PSK HomePage MATLAB Code MATLAB Code for ASK Modulation and Demodulation % The code is written by SalimWireless.Com % Clear previous data and plots clc; clear all; close all; % Parameters Tb = 1; % Bit duration (s) fc = 10; % Carrier frequency (Hz) N_bits = 10; % Number of bits Fs = 100 * fc; % Sampling frequency (ensure at least 2*fc, more for better representation) Ts = 1/Fs; % Sampling interval samples_per_bit = Fs * Tb; % Number of samples per bit duration % Generate random binary data rng(10); % Set random seed for reproducibility binary_data = randi([0, 1], 1, N_bits); % Generate random binary data (0 or 1) % Initialize arrays for continuous signals t_overall = 0:Ts:(N_bits...

BER vs SNR for M-ary QAM, M-ary PSK, QPSK, BPSK, ...

📘 Overview of BER and SNR 🧮 Online Simulator for BER calculation of m-ary QAM and m-ary PSK 🧮 MATLAB Code for BER calculation of M-ary QAM, M-ary PSK, QPSK, BPSK, ... 📚 Further Reading 📂 View Other Topics on M-ary QAM, M-ary PSK, QPSK ... 🧮 Online Simulator for Constellation Diagram of m-ary QAM 🧮 Online Simulator for Constellation Diagram of m-ary PSK 🧮 MATLAB Code for BER calculation of ASK, FSK, and PSK 🧮 MATLAB Code for BER calculation of Alamouti Scheme 🧮 Different approaches to calculate BER vs SNR What is Bit Error Rate (BER)? The abbreviation BER stands for Bit Error Rate, which indicates how many corrupted bits are received (after the demodulation process) compared to the total number of bits sent in a communication process. BER = (number of bits received in error) / (total number of tran...