(-1.0, 4.0)
(-1.0, 4.0)
Ingeniería Biomédica
2026-03-26
Euler’s formula shows the deep connection between trigonometry and complex exponential functions:
\[e^{i\alpha} = \cos(\alpha) + i \sin(\alpha)\]
To prove this, we use the Taylor Series method at \(a = 0\). This allows us to write complex functions as infinite sums of simple terms.
First, we define the standard power series for the three functions involved:
Exponential Function: \[e^x = \sum_{n=0}^{\infty} \frac{x^n}{n!} = 1 + x + \frac{x^2}{2!} + \frac{x^3}{3!} + \frac{x^4}{4!} + \dots\]
Trigonometric Functions: \[\cos(\alpha) = 1 - \frac{\alpha^2}{2!} + \frac{\alpha^4}{4!} - \frac{\alpha^6}{6!} + \dots\] \[\sin(\alpha) = \alpha - \frac{\alpha^3}{3!} + \frac{\alpha^5}{5!} - \frac{\alpha^7}{7!} + \dots\]
Next, we replace the variable \(x\) in the exponential series with the imaginary term \(i\alpha\):
\[e^{i\alpha} = 1 + (i\alpha) + \frac{(i\alpha)^2}{2!} + \frac{(i\alpha)^3}{3!} + \frac{(i\alpha)^4}{4!} + \dots\]
We must apply the powers of the imaginary unit \(i\): * \(i^2 = -1\) * \(i^3 = -i\) * \(i^4 = 1\) * \(i^5 = i\)
After applying the powers of \(i\), the series looks like this:
\[e^{i\alpha} = 1 + i\alpha - \frac{\alpha^2}{2!} - i\frac{\alpha^3}{3!} + \frac{\alpha^4}{4!} + i\frac{\alpha^5}{5!} - \dots\]
Now, we group the Real parts and the Imaginary parts:
\[e^{i\alpha} = \underbrace{\left( 1 - \frac{\alpha^2}{2!} + \frac{\alpha^4}{4!} - \dots \right)}_{\text{Real Part}} + i \underbrace{\left( \alpha - \frac{\alpha^3}{3!} + \frac{\alpha^5}{5!} - \dots \right)}_{\text{Imaginary Part}}\]
We can now see that the two groups match the original Taylor series for Sine and Cosine:
Final Identity: \[e^{i\alpha} = \cos(\alpha) + i \sin(\alpha)\]
A complex number \(z\) is typically written in Cartesian coordinates as:
\[z = a + bi\]
Where:
We can treat \((a, b)\) as coordinates on a 2D plane (the Argand plane).
To convert to polar form, we define:
From the geometry of a right triangle, we know: * \(a = r \cos(\alpha)\) * \(b = r \sin(\alpha)\)
Substituting these into \(z = a + bi\): \[z = r \cos(\alpha) + i(r \sin(\alpha))\] \[z = r (\cos(\alpha) + i \sin(\alpha))\]
Using Euler’s Formula (\(e^{i\alpha} = \cos\alpha + i\sin\alpha\)), we simplify this to: \[z = r e^{i\alpha}\]
In biomedical engineering, the form \(z = r e^{i\alpha}\) is superior for calculations:
1. Multiplication: \[z_1 z_2 = (r_1 e^{i\alpha_1})(r_2 e^{i\alpha_2}) = r_1 r_2 e^{i(\alpha_1 + \alpha_2)}\]
2. Division: \[\frac{z_1}{z_2} = \frac{r_1}{r_2} e^{i(\alpha_1 - \alpha_2)}\]
3. Signal Analysis: It allows us to represent a rotating phasor, which is the basis for the Fourier Transform.
| Feature | Rectangular | Exponential |
|---|---|---|
| Notation | \(a + bi\) | \(r e^{i\alpha}\) |
| Best for | Addition / Subtraction | Mult. / Div. / Powers |
| Biomedical | Signal Amplitude | Phase Shift Analysis |
\[e^{i\alpha} = \cos(\alpha) + i \sin(\alpha)\]
A continuous-time LTI system is defined by:
\[ y(t) = x(t) * h(t) = \int_{-\infty}^{\infty} h(\tau)\,x(t - \tau)\,d\tau \]
where
Let the input be a sinusoid:
\[ x(t) = A \cos(\omega_0 t + \phi) \]
Using Euler’s identity:
\[ x(t) = \frac{A}{2} \left[e^{j(\omega_0 t + \phi)} + e^{-j(\omega_0 t + \phi)}\right] \]
Because the system is linear, we can treat each exponential term separately:
\[ y(t) = \frac{A}{2}\left[ y_1(t) + y_2(t) \right] \] where \[ y_1(t) = \frac{A}{2} e^{j\phi} \int_{-\infty}^{\infty} h(\tau)\, e^{j\omega_0(t-\tau)}\, d\tau \] \[ y_2(t) = \frac{A}{2} e^{-j\phi} \int_{-\infty}^{\infty} h(\tau)\, e^{-j\omega_0(t-\tau)}\, d\tau \]
We can factor out the term \(e^{j\omega_0 t}\):
\[ y_1(t) = e^{j\phi} e^{j\omega_0 t} \int_{-\infty}^{\infty} h(\tau)\, e^{-j\omega_0 \tau}\, d\tau \]
Let \[ H_1 = \int_{-\infty}^{\infty} h(\tau)\, e^{-j\omega_0 \tau}\, d\tau \]
Then, \[ y_1(t) = e^{j\phi} e^{j\omega_0 t} H_1 \]
Similarly, \[ y_2(t) = e^{-j\phi} e^{-j\omega_0 t} H_1 \]
The total output is:
\[ y(t) = \frac{A}{2}\left[e^{j\phi} H_1 e^{j\omega_0 t} + e^{-j\phi} H_1^* e^{-j\omega_0 t}\right] \]
If we express \(H_1\) in polar form: \[ H_1 = |H_1| e^{j\theta} \]
Then, \[ y(t) = A |H_1| \cos(\omega_0 t + \phi + \theta) \]
\[ \boxed{y(t) = A |H_1| \cos(\omega_0 t + \phi + \theta)} \]
A discrete-time LTI system is defined by the convolution sum:
\[ y[n] = \sum_{k=-\infty}^{\infty} h[k]\,x[n - k] \]
where
Let the input be a discrete sinusoid:
\[ x[n] = A \cos(\omega_0 n + \phi) \]
Using Euler’s identity:
\[ x[n] = \frac{A}{2}\left[e^{j(\omega_0 n + \phi)} + e^{-j(\omega_0 n + \phi)}\right] \]
Because the system is linear, each exponential term can be treated independently:
\[ y[n] = \frac{A}{2}\left[y_1[n] + y_2[n]\right] \]
where \[ y_1[n] = e^{j\phi} \sum_{k=-\infty}^{\infty} h[k]\, e^{j\omega_0 (n - k)} \] \[ y_2[n] = e^{-j\phi} \sum_{k=-\infty}^{\infty} h[k]\, e^{-j\omega_0 (n - k)} \]
We can factor out \(e^{j\omega_0 n}\):
\[ y_1[n] = e^{j\phi} e^{j\omega_0 n} \sum_{k=-\infty}^{\infty} h[k]\, e^{-j\omega_0 k} \]
Let \[ H_1 = \sum_{k=-\infty}^{\infty} h[k]\, e^{-j\omega_0 k} \]
Then, \[ y_1[n] = e^{j\phi} H_1 e^{j\omega_0 n} \]
Similarly, \[ y_2[n] = e^{-j\phi} H_1^* e^{-j\omega_0 n} \]
The total output becomes:
\[ y[n] = \frac{A}{2}\left[e^{j\phi} H_1 e^{j\omega_0 n} + e^{-j\phi} H_1^* e^{-j\omega_0 n}\right] \]
If we write \(H_1\) in polar form: \[ H_1 = |H_1| e^{j\theta} \]
Then, \[ y[n] = A |H_1| \cos(\omega_0 n + \phi + \theta) \]
\[ \boxed{y[n] = A |H_1| \cos(\omega_0 n + \phi + \theta)} \]
\[y(t) = x(t) * h(t) = \int_{-\infty}^{\infty} x(\tau) h(t - \tau) d\tau\]
\[y[n] = \sum_{k=-\infty}^{\infty} x[k] h[n-k]\]
\[X(f) H(f) = Y(f)\]
(-1.0, 4.0)
(-1.0, 4.0)
\[x(t) = \sum_{n=-\infty}^{\infty} C_n e^{jn\omega_0 t}\]
where \(C_n\) are the Fourier coefficients.
\[C_n = \frac{1}{T} \int_{0}^{T} x(t) e^{-jn\omega_0 t} dt\]
For a periodic signal \(x(t)\) with period \(T\):
\[ x(t) = \sum_{k=-\infty}^{\infty} c_k e^{j k \Omega_0 t}, \]
where:
\[ c_k = \frac{1}{T} \int_{T_0}^{T_0 + T} x(t)e^{-jk\Omega_0 t}\,dt, \quad \Omega_0 = \frac{2\pi}{T}. \]
We define a periodic signal \(x(t)\) with period \(2\pi\):
\[ x(t) = \begin{cases} 0, & -\pi \le t < 0,\\ 1, & 0 \le t < \pi. \end{cases} \]
This corresponds to a 50% duty cycle pulse.
Expected result:
\[ c_k = \begin{cases} \dfrac{1}{2}, & k=0,\\[4pt] \dfrac{j\left((-1)^k - 1\right)}{2\pi k}, & k \ne 0. \end{cases} \]
import numpy as np
import pandas as pd
def ck_numeric(k):
k = np.asarray(k, dtype=float)
return np.where(
k != 0,
1j*(np.exp(1j*np.pi*k)-1)/(2*np.pi*k),
0.5
)
K = 10
ks = np.arange(-K, K+1)
cks = ck_numeric(ks)
pd.DataFrame({
"k": ks,
"Re{c_k}": np.real(cks),
"Im{c_k}": np.imag(cks),
"|c_k|": np.abs(cks)
}) k Re{c_k} Im{c_k} |c_k|
0 -10 1.949086e-17 -0.000000 1.949086e-17
1 -9 -1.949086e-17 0.035368 3.536777e-02
2 -8 1.949086e-17 -0.000000 1.949086e-17
3 -7 -1.949086e-17 0.045473 4.547284e-02
4 -6 1.949086e-17 -0.000000 1.949086e-17
5 -5 -1.949086e-17 0.063662 6.366198e-02
6 -4 1.949086e-17 -0.000000 1.949086e-17
7 -3 -1.949086e-17 0.106103 1.061033e-01
8 -2 1.949086e-17 -0.000000 1.949086e-17
9 -1 -1.949086e-17 0.318310 3.183099e-01
10 0 5.000000e-01 0.000000 5.000000e-01
11 1 -1.949086e-17 -0.318310 3.183099e-01
12 2 1.949086e-17 0.000000 1.949086e-17
13 3 -1.949086e-17 -0.106103 1.061033e-01
14 4 1.949086e-17 0.000000 1.949086e-17
15 5 -1.949086e-17 -0.063662 6.366198e-02
16 6 1.949086e-17 0.000000 1.949086e-17
17 7 -1.949086e-17 -0.045473 4.547284e-02
18 8 1.949086e-17 0.000000 1.949086e-17
19 9 -1.949086e-17 -0.035368 3.536777e-02
20 10 1.949086e-17 0.000000 1.949086e-17
import numpy as np
import matplotlib.pyplot as plt
def x_numeric(t):
# Original pulse over [-π, π): 0 on [-π,0), 1 on [0,π).
tw = ((t + np.pi) % (2*np.pi)) - np.pi
return np.where((tw >= 0) & (tw < np.pi), 1.0, 0.0)
def partial_sum(t, N):
ks = np.arange(-N, N+1)
cks = ck_numeric(ks)
return np.sum(cks[:, None] * np.exp(1j*np.outer(ks, t)), axis=0)
tgrid = np.linspace(-2*np.pi, 2*np.pi, 4000, endpoint=False)
xg = x_numeric(tgrid)(-6.283185307179586, 6.283185307179586)
(-6.283185307179586, 6.283185307179586)
(-6.283185307179586, 6.283185307179586)
| Concept | Expression |
|---|---|
| Fundamental frequency | \(\Omega_0 = \frac{2\pi}{T}\) |
| Coefficient \(c_k\) | \(\frac{1}{T}\int x(t)e^{-jk\Omega_0 t}dt\) |
| Reconstructed signal | \(S_N(t) = \sum_{k=-N}^{N} c_k e^{jk\Omega_0 t}\) |
| Example \(x(t)\) | Half-wave pulse in \([-\pi, \pi)\) |
(array([-5., 0., 5., 10., 15., 20., 25., 30., 35.]), [Text(-5.0, 0, '−5'), Text(0.0, 0, '0'), Text(5.0, 0, '5'), Text(10.0, 0, '10'), Text(15.0, 0, '15'), Text(20.0, 0, '20'), Text(25.0, 0, '25'), Text(30.0, 0, '30'), Text(35.0, 0, '35')])
(array([4.5, 5. , 5.5, 6. , 6.5, 7. , 7.5, 8. , 8.5]), [Text(0, 4.5, '4.5'), Text(0, 5.0, '5.0'), Text(0, 5.5, '5.5'), Text(0, 6.0, '6.0'), Text(0, 6.5, '6.5'), Text(0, 7.0, '7.0'), Text(0, 7.5, '7.5'), Text(0, 8.0, '8.0'), Text(0, 8.5, '8.5')])
(array([-40., -30., -20., -10., 0., 10., 20., 30., 40.]), [Text(-40.0, 0, '−40'), Text(-30.0, 0, '−30'), Text(-20.0, 0, '−20'), Text(-10.0, 0, '−10'), Text(0.0, 0, '0'), Text(10.0, 0, '10'), Text(20.0, 0, '20'), Text(30.0, 0, '30'), Text(40.0, 0, '40')])
(array([4.5, 5. , 5.5, 6. , 6.5, 7. , 7.5, 8. , 8.5]), [Text(0, 4.5, '4.5'), Text(0, 5.0, '5.0'), Text(0, 5.5, '5.5'), Text(0, 6.0, '6.0'), Text(0, 6.5, '6.5'), Text(0, 7.0, '7.0'), Text(0, 7.5, '7.5'), Text(0, 8.0, '8.0'), Text(0, 8.5, '8.5')])
DTFS Coefficients:
C[0] = 2.0000+0.0000j
C[1] = -0.6036-0.6036j
C[2] = 0.0000+0.0000j
C[3] = 0.1036-0.1036j
C[4] = 0.0000+0.0000j
C[5] = 0.1036+0.1036j
C[6] = 0.0000+0.0000j
C[7] = -0.6036+0.6036j
For a continuous-time signal \(x(t)\), its Fourier transform is
\[ X(\omega)=\int_{-\infty}^{\infty} x(t)e^{-j\omega t}\,dt \]
For a discrete-time signal \(x[n]\), its Fourier transform is
\[ X(e^{j\omega})=\sum_{n=-\infty}^{\infty} x[n]e^{-j\omega n} \]
We will prove:
Assume \(x(t)\in\mathbb{R}\). Then
\[ X(\omega)=\int_{-\infty}^{\infty} x(t)e^{-j\omega t}\,dt \]
Take complex conjugate:
\[ X^*(\omega) = \left(\int_{-\infty}^{\infty} x(t)e^{-j\omega t}\,dt\right)^* \]
Because \(x(t)\) is real, \(x^*(t)=x(t)\). Therefore
\[ X^*(\omega) = \int_{-\infty}^{\infty} x(t)e^{j\omega t}\,dt \]
But
\[ X(-\omega) = \int_{-\infty}^{\infty} x(t)e^{-j(-\omega)t}\,dt = \int_{-\infty}^{\infty} x(t)e^{j\omega t}\,dt \]
Hence
\[ \boxed{X(-\omega)=X^*(\omega)} \]
This is called conjugate symmetry.
From
\[ X(-\omega)=X^*(\omega) \]
we get two useful results:
\[ |X(-\omega)|=|X(\omega)| \]
If \(X(\omega)=|X(\omega)|e^{j\phi(\omega)}\), then
\[ \phi(-\omega)=-\phi(\omega) \]
So, for a real signal:
The same proof works for the DTFT if \(x[n]\) is real.
Start from the DTFT:
\[ X(e^{j\omega})=\sum_{n=-\infty}^{\infty} x[n]e^{-j\omega n} \]
Now evaluate it at \(\omega+2\pi\):
\[ X(e^{j(\omega+2\pi)}) = \sum_{n=-\infty}^{\infty} x[n]e^{-j(\omega+2\pi)n} \]
Separate the exponential:
\[ X(e^{j(\omega+2\pi)}) = \sum_{n=-\infty}^{\infty} x[n]e^{-j\omega n}e^{-j2\pi n} \]
Because \(n\) is an integer,
\[ e^{-j2\pi n}=1 \]
for every integer \(n\). Then
\[ X(e^{j(\omega+2\pi)}) = \sum_{n=-\infty}^{\infty} x[n]e^{-j\omega n} = X(e^{j\omega}) \]
So,
\[ \boxed{X(e^{j(\omega+2\pi)})=X(e^{j\omega})} \]
More generally,
\[ \boxed{X(e^{j(\omega+2\pi k)})=X(e^{j\omega}),\quad k\in\mathbb{Z}} \]
For the continuous-time Fourier transform,
\[ X(\omega)=\int_{-\infty}^{\infty} x(t)e^{-j\omega t}\,dt \]
Now test \(\omega+2\pi\):
\[ X(\omega+2\pi) = \int_{-\infty}^{\infty} x(t)e^{-j(\omega+2\pi)t}\,dt = \int_{-\infty}^{\infty} x(t)e^{-j\omega t}e^{-j2\pi t}\,dt \]
Here is the key point:
Therefore, in general,
\[ \boxed{X(\omega+2\pi)\neq X(\omega)} \]
So the CTFT is not periodic in general.
Let
\[ x[n]=\delta[n]+\delta[n-1] \]
Its DTFT is
\[ X(e^{j\omega})=1+e^{-j\omega} \]
Factor it:
\[ X(e^{j\omega}) = e^{-j\omega/2}\left(e^{j\omega/2}+e^{-j\omega/2}\right) = 2\cos\left(\frac{\omega}{2}\right)e^{-j\omega/2} \]
Magnitude:
\[ |X(e^{j\omega})|=2\left|\cos\left(\frac{\omega}{2}\right)\right| \]
This magnitude is:
So this example shows both properties at the same time.
If \(x(t)\) is real:
\[ \boxed{X(-\omega)=X^*(\omega)} \]
So the spectrum is symmetric in the conjugate sense, but not periodic in general.
For any \(x[n]\):
\[ \boxed{X(e^{j(\omega+2\pi)})=X(e^{j\omega})} \]
So the spectrum is periodic.
If \(x[n]\) is also real, then:
\[ \boxed{X(e^{-j\omega})=X^*(e^{j\omega})} \]
So the DTFT spectrum is both periodic and symmetric.