Pulses and Spectra

From LNTwww
Revision as of 17:48, 25 March 2023 by Noah (talk | contribs)

Open Applet in new Tab   Deutsche Version Öffnen


Applet Description


Time-limited symmetric signals   ⇒   "pulses"  $x(t)$  and the corresponding spectral functions  $X(f)$  are considered, namely

  • Gaussian pulse, 
  • rectangular pulse,  
  • triangular pulse, 
  • trapezoidal pulse, 
  • raised cosine pulse,
  • cosine square pulse.


Further it is to be noted:

  • The functions  $x(t)$  resp.  $X(f)$  are shown for up to two parameter sets in one diagram each.
  • The red curves and numbers apply to the left parameter set, the blue ones to the right parameter set.
  • The abscissas  $t$  (time) and  $f$  (frequency) as well as the ordinates  $x(t)$  (signal values) and  $X(f)$  (spectral values) are normalized.


Theoretical Background


Relationship $x(t)\Leftrightarrow X(f)$

  • The relationship between the time function  $x(t)$  and the spectrum  $X(f)$  is given by the  "first Fourier integral" :
$$X(f)={\rm FT} [x(t)] = \int_{-\infty}^{+\infty}x(t)\cdot {\rm e}^{-{\rm j}2\pi f t}\hspace{0.15cm} {\rm d}t\hspace{1cm} \rm FT\hspace{-0.1cm}: \ Fourier \ transform.$$
  • In order to calculate the time function  $x(f)$  from the spectral function  $x(t)$  one needs the  "second Fourier integral":
$$x(t)={\rm IFT} [X(f)] = \int_{-\infty}^{+\infty}X(f)\cdot {\rm e}^{+{\rm j}2\pi f t}\hspace{0.15cm} {\rm d}f\hspace{1cm} {\rm IFT}\hspace{-0.1cm}: \rm Inverse \ Fourier \ transform.$$
  • In all examples we use real and even functions.  Thus:
$$x(t)=\int_{-\infty}^{+\infty}X(f)\cdot \cos(2\pi ft) \hspace{0.15cm} {\rm d}f \ \ \circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\ \ \ X(f)=\int_{-\infty}^{+\infty}x(t)\cdot \cos(2\pi ft) \hspace{0.15cm} {\rm d}t .$$
  • $x(t)$  and  $X(f)$  have different units, for example  $x(t)$  in  $\rm V$,  $X(f)$  in  $\rm V/Hz$.
  • The relationship between this module and the similarly constructed applet  "Frequency response & Impulse response"  is based on the  "Duality Theorem".
  • All times are normalized to a time  $T$  and all frequencies are normalized to  $1/T$   ⇒   the spectral values  $X(f)$  still have to be multiplied by the normalization time  $T$ .


$\text{Example:}$   If one sets a rectangular pulse with amplitude  $A_1 = 1$  and equivalent pulse duration  $\Delta t_1 = 1$  then  $x_1(t)$  in the range  $-0.5 < t < +0. 5$  equal to one and outside this range equal to zero.  The spectral function  $X_1(f)$  proceeds  $\rm si$–shaped with  $X_1(f= 0) = 1$  and the first zero at  $f=1$.

  • If a rectangular pulse with  $A = K = 3 \ \rm V$  and  $\delta t = T = 2 \ \rm ms$  is to be simulated with this setting, then all signal values with  $K = 3 \ \rm V$  and all spectral values with  $K \cdot T = 0. 006 \ \rm V/Hz$  to be multiplied by.
  • The maximum spectral value is then  $X(f= 0) = 0.006 \ \rm V/Hz$  and the first zero is at  $f=1/T = 0.5 \ \rm kHz$.


Gaussian Pulse

  • Die Zeitfunktion des Gaußimpulses mit der Höhe  $K$  und der (äquivalenten) Dauer  $\Delta t$  lautet:
$$x(t)=K\cdot {\rm e}^{-\pi\hspace{0.05cm}\cdot \hspace{0.05cm}(t/\Delta t)^2}.$$
  • Die äquivalente Zeitdauer  $\Delta t$  ergibt sich aus dem flächengleichen Rechteck.
  • Der Wert bei  $t = \Delta t/2$  ist um den Faktor  $0.456$  kleiner als der Wert bei  $t=0$.
  • Für die Spektralfunktion erhält man gemäß der Fouriertransformation:
$$X(f)=K\cdot \Delta t \cdot {\rm e}^{-\pi(f\hspace{0.05cm}\cdot \hspace{0.05cm} \Delta t)^2} .$$
  • Je kleiner die äquivalente Zeitdauer  $\Delta t$  ist, um so breiter und niedriger ist das Spektrum   ⇒   Reziprozitätsgesetz von Bandbreite und Impulsdauer.
  • Sowohl  $x(t)$  als auch  $X(f)$  sind zu keinem  $f$–  bzw.  $t$–Wert exakt gleich Null.
  • Für praktische Anwendungen kann der Gaußimpuls jedoch in Zeit und Frequenz als begrenzt angenommen werden.  Zum Beispiel ist  $x(t)$  bereits bei  $t=1.5 \Delta t$  auf weniger als  $0.1\% $  des Maximums abgefallen.


Rectangular Pulse

  • Die Zeitfunktion des Rechteckimpulses mit der Höhe  $K$  und der (äquivalenten) Dauer  $\Delta t$  lautet:
$$x(t) = \left\{ \begin{array}{l} \hspace{0.25cm}K \\ K /2 \\ \hspace{0.25cm} 0 \\ \end{array} \right.\quad \quad \begin{array}{*{20}c} {\rm{f\ddot{u}r}} \\ {\rm{f\ddot{u}r}} \\ {\rm{f\ddot{u}r}} \\ \end{array}\begin{array}{*{20}c} {\left| \hspace{0.05cm} t\hspace{0.05cm} \right| < T/2,} \\ {\left| \hspace{0.05cm}t\hspace{0.05cm} \right| = T/2,} \\ {\left|\hspace{0.05cm} t \hspace{0.05cm} \right| > T/2.} \\ \end{array}$$
  • Der  $\pm \Delta t/2$–Wert liegt mittig zwischen links- und rechtsseitigem Grenzwert.
  • Für die Spektralfunktion erhält man entsprechend den Gesetzmäßigkeiten der Fouriertransformation (1. Fourierintegral):
$$X(f)=K\cdot \Delta t \cdot {\rm si}(\pi\cdot \Delta t \cdot f) \quad \text{mit} \ {\rm si}(x)=\frac{\sin(x)}{x}.$$
  • Der Spektralwert bei  $f=0$  ist gleich der Rechteckfläche der Zeitfunktion.
  • Die Spektralfunktion besitzt Nullstellen in äquidistanten Abständen  $1/\Delta t$.
  • Das Integral über der Spektralfunktion  $X(f)$  ist gleich dem Signalwert zum Zeitpunkt  $t=0$, also der Impulshöhe  $K$.


Triangular Pulse

  • Die Zeitfunktion des Dreieckimpulses mit der Höhe  $K$  und der (äquivalenten) Dauer  $\Delta t$  lautet:
$$x(t) = \left\{ \begin{array}{l} \hspace{0.25cm}K\cdot (1-|t|/{\Delta t}) \\ \hspace{0.25cm} 0 \\ \end{array} \right.\quad \quad \begin{array}{*{20}c} {\rm{f\ddot{u}r}} \\ {\rm{f\ddot{u}r}} \\ \end{array}\begin{array}{*{20}c} {\left| \hspace{0.05cm} t\hspace{0.05cm} \right| < \Delta t,} \\ {\left| \hspace{0.05cm}t\hspace{0.05cm} \right| \ge \Delta t.} \\ \end{array}$$
  • Die absolute Zeitdauer ist  $2 \cdot \Delta t$;  diese ist doppelt so groß als die des Rechtecks.
  • Für die Spektralfunktion erhält man gemäß der Fouriertransformation:
$$X(f)=K\cdot \Delta f \cdot {\rm si}^2(\pi\cdot \Delta t \cdot f) \quad \text{mit} \quad {\rm si}(x)=\frac{\sin(x)}{x}.$$
  • Obige Zeitfunktion ist gleich der Faltung zweier Rechteckimpulse, jeweils mit Breite  $\Delta t$.
  • Daraus folgt:  $X(f)$  beinhaltet anstelle der  ${\rm si}$-Funktion die  ${\rm si}^2$-Funktion.
  • $X(f)$  weist somit ebenfalls Nullstellen im äquidistanten Abständen  $1/\Delta f$  auf.
  • Der asymptotische Abfall von  $X(f)$  erfolgt hier mit  $1/f^2$, während zum Vergleich der Rechteckimpuls mit  $1/f$  abfällt.


Trapezoidal Pulse

Die Zeitfunktion des Trapezimpulses mit der Höhe  $K$  und den Zeitparametern  $t_1$  und  $t_2$  lautet:

$$x(t) = \left\{ \begin{array}{l} \hspace{0.25cm}K \\ K\cdot \frac{t_2-|t|}{t_2-t_1} \\ \hspace{0.25cm} 0 \\ \end{array} \right.\quad \quad \begin{array}{*{20}c} {\rm{f\ddot{u}r}}\quad \\ {\rm{f\ddot{u}r}}\quad \\ {\rm{f\ddot{u}r}} \quad \\ \end{array}\begin{array}{*{20}c} {\left| \hspace{0.05cm} t\hspace{0.05cm} \right| \le t_1,} \\ {t_1\le \left| \hspace{0.05cm}t\hspace{0.05cm} \right| \le t_2,} \\ {\left|\hspace{0.05cm} t \hspace{0.05cm} \right| \ge t_2.} \\ \end{array}$$
  • Für die äquivalente Impulsdauer (flächengleiches Rechteck) gilt:   $\Delta t = t_1+t_2$.
  • Der Rolloff-Faktor (im Zeitbereich) kennzeichnet die Flankensteilheit:
$$r=\frac{t_2-t_1}{t_2+t_1}.$$
  • Der Sonderfall  $r=0$  entspricht dem Rechteckimpuls und der Sonderfall  $r=1$  dem Dreieckimpuls.
  • Für die Spektralfunktion erhält man gemäß der Fouriertransformation:
$$X(f)=K\cdot \Delta t \cdot {\rm si}(\pi\cdot \Delta t \cdot f)\cdot {\rm si}(\pi \cdot r \cdot \Delta t \cdot f) \quad \text{mit} \quad {\rm si}(x)=\frac{\sin(x)}{x}.$$
  • Der asymptotische Abfall von  $X(f)$  liegt zwischen  $1/f$  $($für Rechteck,  $r=0)$  und  $1/f^2$  $($für Dreieck,  $r=1)$.


Cosine roll-off Pulse

Die Zeitfunktion des Cosinus-Rolloff-Impulses mit der Höhe  $K$  und den Zeitparametern  $t_1$  und  $t_2$  lautet:

$$x(t) = \left\{ \begin{array}{l} \hspace{0.25cm}K \\ K\cdot \cos^2\Big(\frac{|t|-t_1}{t_2-t_1}\cdot {\pi}/{2}\Big) \\ \hspace{0.25cm} 0 \\ \end{array} \right.\quad \quad \begin{array}{*{20}c} {\rm{f\ddot{u}r}}\quad \\ {\rm{f\ddot{u}r}}\quad \\ {\rm{f\ddot{u}r}}\quad \\ \end{array}\begin{array}{*{20}c} {\left| \hspace{0.05cm} t\hspace{0.05cm} \right| \le t_1,} \\ {t_1\le \left| \hspace{0.05cm}t\hspace{0.05cm} \right| \le t_2,} \\ {\left|\hspace{0.05cm} t \hspace{0.05cm} \right| \ge t_2.} \\ \end{array}$$
  • Für die äquivalente Impulsdauer (flächengleiches Rechteck) gilt:   $\Delta t = t_1+t_2$.
  • Der Rolloff-Faktor (im Zeitbereich) kennzeichnet die Flankensteilheit:
$$r=\frac{t_2-t_1}{t_2+t_1}.$$
  • Der Sonderfall  $r=0$  entspricht dem Rechteckimpuls und der Sonderfall  $r=1$  dem Cosinus-Quadrat-Impuls.
  • Für die Spektralfunktion erhält man gemäß der Fouriertransformation:
$$X(f)=K\cdot \Delta t \cdot \frac{\cos(\pi \cdot r\cdot \Delta t \cdot f)}{1-(2\cdot r\cdot \Delta t \cdot f)^2} \cdot {\rm si}(\pi \cdot \Delta t \cdot f).$$
  • Je größer der Rolloff-Faktor  $r$  ist, desto schneller nimmt  $X(f)$  asymptotisch mit  $f$  ab.


Cosine square Pulse

  • Dies ist ein Sonderfall des Cosinus-Rolloff-Impulses und ergibt sich für  $r=1 \hspace{0.3cm} \Rightarrow \hspace{0.3cm}t_1=0, \ t_2= \Delta t$:
$$x(t) = \left\{ \begin{array}{l} \hspace{0.25cm}K\cdot \cos^2\Big(\frac{|t|\hspace{0.05cm}\cdot \hspace{0.05cm} \pi}{2\hspace{0.05cm}\cdot \hspace{0.05cm} \Delta t}\Big) \\ \hspace{0.25cm} 0 \\ \end{array} \right.\quad \quad \begin{array}{*{20}c} {\rm{f\ddot{u}r}} \\ {\rm{f\ddot{u}r}} \\ \end{array}\begin{array}{*{20}c} {\left| \hspace{0.05cm} t\hspace{0.05cm} \right| < \Delta t,} \\ {\left| \hspace{0.05cm}t\hspace{0.05cm} \right| \ge \Delta t.} \\ \end{array}$$
  • Für die Spektralfunktion erhält man gemäß der Fouriertransformation:
$$X(f)=K\cdot \Delta f \cdot \frac{\pi}{4}\cdot \big [{\rm si}(\pi(\Delta t\cdot f +0.5))+{\rm si}(\pi(\Delta t\cdot f -0.5))\big ]\cdot {\rm si}(\pi \cdot \Delta t \cdot f).$$
  • Wegen der letzten  ${\rm si}$-Funktion ist  $X(f)=0$  für alle Vielfachen von  $F=1/\Delta t$.  Die äquidistanten Nulldurchgänge des Cos-Rolloff-Impulses bleiben erhalten.
  • Aufgrund des Klammerausdrucks weist  $X(f)$  nun weitere Nulldurchgänge bei  $f=\pm1.5 F$,  $\pm2.5 F$,  $\pm3.5 F$, ... auf.
  • Für die Frequenz  $f=\pm F/2$  erhält man die Spektralwerte  $K\cdot \Delta t/2$.
  • Der asymptotische Abfall von  $X(f)$  verläuft in diesem Sonderfall mit  $1/f^3$.

Exercises


  • First select the number  $(1,\text{...}, 7)$  of the exercise.  The number  $0$  corresponds to a "Reset":  Same setting as at program start.
  • A task description is displayed.  The parameter values ​​are adjusted.  Solution after pressing "Show solution".
  • "Red" refers to the first parameter set ⇒ $x_1(t) \circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\ X_1(f)$,  "Blue" refers to the second parameter set ⇒ $x_2(t) \circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\ X_2(f)$.
  • Values with magnitude less than  $0.0005$  are output in the program as "zero".


(1)   Compare the  red Gaussian pulse  $(A_1 = 1, \Delta t_1 = 1)$  with the  blue rectangular pulse  $(A_2 = 1, \Delta t_2 = 1)$ ⇒ default setting.
          What are the differences in the time and frequency domain?

  • The Gaussian pulse theoretically reaches infinity in the time as well as in the frequency domain.
  • Practically  $x_1(t)$  for  $|t| > 1.5$  and  $X_1(f)$  for  $|f| > 1.5$  are almost zero.
  • The rectangle is strictly limited in time:  $x_2(|t| > 0.5) \equiv 0$.  $X_2(f)$  has shares in a much larger range than  $X_1(f)$.
  • It holds  $X_1(f = 0) = X_2(f = 0)$  since the integral over the Gaussian pulse  $x_1(t)$  is equal to the integral over the rectangular pulse  $x_2(t)$.


(2)   Compare the  red Gaussian pulse  $(A_1 = 1,  \Delta t_1 = 1)$ with the  blue rectangular pulse  $(A_2 = 1,  \Delta t_2)$.
          Vary the equivalent pulse duration  $\Delta t_2$  between  $0.5$  and  $2$.  Interpret the displayed graphs.

  • One can recognize the reciprocity law of bandwidth and pulse duration.  The greater  $\Delta t_2$, the higher and narrower the spectral function  $X_2(f)$.
  • For each setting of  $\Delta t_2$,  $x_1(t=0)$  and  $x_2(t=0)$  are equal   ⇒   Also, the integrals over  $X_1(f)$  and  $X_2(f)$  are identical.


(3)   Compare the  red Gaussian pulse  $(A_1 = 1,  \Delta t_1 = 1)$ with the  blue rectangular pulse  $(A_2 = 1,  \Delta t_2 = 0.5)$.
          Vary  $\Delta t_2$  between  $0.05$  and  $2$.  Interpret the displayed graphs and extrapolate the result.

  • The blue spectrum is now twice as wide as the red one, but only half as high.  First zero of  $X_1(f)$  at  $f = 1$, of  $X_2(f)$  at  $f = 2$.
  • Reduction of  $\Delta t_2$:  $X_2(f)$  lower and wider.  Very flat course at  $\Delta t_2 = 0.05$:  $X_2(f = 0)= 0.05$,  $X_2(f = \pm 3)= 0.048$.
  • If one choose  $\Delta t_2 = \varepsilon \to 0$  (not possible in the program),  the result would be the almost constant, very small spectrum  $X_2(f)=A \cdot \varepsilon \to 0$.
  • Increasing the amplitude to  $A=1/\varepsilon$  results in the constant spectral function  $X_2(f) = 1$  of the Dirac function  $\delta(t)$.  That means:
  • $\delta(t)$  is approximated by a rectangle  $($width  $\Delta t = \varepsilon \to 0$,  height  $A = 1/\varepsilon \to \infty)$.  The weight of the Dirac function is one:  $x(t) = 1 \cdot \delta (t)$.


(4)   Compare the  rectangular pulse  $(A_1 = 1,   \Delta t_1 = 1)$  with the  triangular pulse  $(A_2 = 1,   \Delta t_2 = 1)$.  Interpret the spectral functions.

  • The (normalized) spectrum of the rectangle  $x_1(t)$  with the (normalized) parameters  $A_1 = 1, \ \ \Delta t_1 = 1$  is:  $X_1(f)= {\rm si}(\pi\cdot f)= {\rm sinc}(f)$.
  • The convolution of the rectangle  $x_1(t)$  with itself gives the triangle  $x_2(t) = x_1(t) \star x_1(t)$.  By the convolution theorem:   $X_2(f) = X_1(f)^2 $.
  • By squaring the  ${\rm sinc}(f)$–shaped spectral function  $X_1(f)$  the zeros of  $X_2(f)$  remain unchanged.  But now it holds that: $X_2(f) \ge 0$.


(5)   Compare the  trapezoidal pulse  $(A_1 = 1, \Delta t_1 = 1, r_1 = 0.5)$  with the  triangular pulse $(A_2 = 1,   \Delta t_2 = 1)$.
         Vary  $r_1$  between  $0$  and  $1$.  Interpret the spectral function  $X_1(f)$.

  • The trapezoidal pulse with roll–off factor  $r_1= 0$  is identical to the rectangular pulse.  The "normalized spectrum" is  $X_1(f)= {\rm sinc}(f)$.
  • The trapezoidal pulse with roll–off factor  $r_1= 1$  is identical to the triangular pulse.  The "normalized spectrum" is  $X_1(f)= {\rm sinc}^2(f)$.
  • In both cases  $X_1(f)$  has equidistant zeros at  $\pm 1$,  $\pm 2$, ...  (none else);   $0 < r_1 < 1$:  depending on  $r_1$  further zeros.


(6)   Compare this  trapezoidal pulse  with the cosine roll-off pulse  $(A_2 = 1,\ \Delta t_2 = 1.0,\ r_2 = 0.5)$.
         Vary  $r_2$  between  $0$  and  $1$.  Interpret the spectral function  $X_2(f)$  for  $r_2 = 0.7$.

  • With the same  $r= 0.5$  the cosine roll-off pulse  $X_2(f)$ is for  $f > 1$  greater in magnitude than the trapezoidal pulse.
  • With the same roll-off factor  $(r_1 = r_2= 0.5)$  the drop of  $X_2(f)$  around the frequency  $f = 0.5$  is steeper than the drop of  $X_1(f)$.
  • With  $r_1 = 0.5$  and  $r_2 = 0.7$  $x_1(t) \approx x_2(t)$  is valid and therefore also  $X_1(f) \approx X_2(f)$.  Comparable edge steepness.


(7)   Compare the  red trapezoidal pulse  $(A_1 = 1, \Delta t_1 = 1, \ r_1 = 1)$  with the  blue cosine roll-off pulse  $(A_2 = 1,\ \Delta t_2 = 1.0, \ r_2 = 1)$.
          Interpret the time function  $x_2(t)$  and the spectral function  $X_2(f)$  system theoretically.

  • $x_2(t) = \cos^2(|t|\cdot \pi/2) \ \ \ \text{for} \ |t| \le 1$  is the  cosine square pulse.  Zeros at  $f = \pm 1$,  $\pm 2$, ...
  • For the frequency  $f=\pm 0.5$  one obtains the spectral values  $X_2(f)=0.5$.  The asymptotic decline is shown here with  $1/f^3$.


Applet Manual


Screenshot

    (A)     Theme (changeable graphical user interface design)

  • Dark:   dark background  (recommended by the authors)
  • Bright:   white background  (recommended for beamers and printouts)
  • Deuteranopia:   for users with pronounced green visual impairment
  • Protanopia:   for users with pronounced red visual impairment

    (B)     Preselection for pulse shape  $x_1(t)$  (red curve)

    (C)     Parameter definition for  $x_1(t)$ 

    (D)     Numeric output for  $x_1(t_*)$  and  $X_1(f_*)$

    (E)     Preselection for pulse shape  $x_2(t)$  (blue curve)

    (F)     Parameter definition for  $x_2(t)$ 

    (G)     Numeric output for  $x_2(t_*)$  and  $X_2(f_*)$

    (H)     Setting the time  $t_*$  for the numeric output

    (I)     Setting the frequency  $f_*$  for the numeric output

    (J)     Graphic field for the time domain

    (K)     Graphic field for the frequency domain

    (L)     Selection of the exercise according to the numbers

    (M)     Task description and questions

    (N)     Show and hide sample solution


About the Authors


This interactive calculation tool was designed and implemented at the  Institute for Communications Engineering  at the  Technical University of Munich.

  • The first version was created in 2005 by  Ji Li  as part of her diploma thesis with “FlashMX – Actionscript”  (Supervisor:  Günter Söder).
  • In 2017 the program was redesigned by  David Jobst  (Ingenieurspraxis_Math, Supervisor: Tasnád Kernetzky ) via "HTML5".
  • Last revision and English version 2020 by  Carolin Mirschina  in the context of a working student activity. 


Once again: Open Applet in new Tab

Open Applet in new Tab   Deutsche Version Öffnen