TEMA: Introducción a la simulación de señales y sistemas

OBJETIVOS:

• Simular los conocimientos teóricos adquiridos sobre señales y sistemas lineales y el análisis de señales
determinísticas.
• Calcular transformadas de Fourier.
• Obtener e interpretar espectros de señales.
• Analizar y diferenciar señales de energía y señales de potencia
• Obtener e interpretar espectros de densidad de potencia.

MARCO TEÓRICO:

Transformada de Fourier
La transformada de Fourier nos es una herramienta que nos ayuda a simplificar el estudio de la solución
de determinado tipo de ecuaciones diferenciales, de esta forma se convierte el problema de la solución
de una ED en un problema que se puede solucionar a través de ecuaciones algebraicas. Además, posee
distintas propiedades algebraicas a la hora de aplicar a las derivadas de una señal, o trasladar la señal.
(Almira, 2009)
Esta operación es aplicable a cualquier tipo de función, ya sea periódica o no, a través de las siguientes
funciones.
+∞
𝑇𝑟𝑎𝑛𝑠𝑓𝑜𝑟𝑚𝑎𝑑𝑎 𝑑𝑒 𝐹𝑜𝑢𝑟𝑖𝑒𝑟 𝐹(𝜔) = 𝔍{𝑓(𝑡)} = ∫ 𝑓(𝑡)𝑒 −𝑖𝜔𝑡 𝑑𝑡
−∞

1 +∞
𝑇𝑟𝑎𝑛𝑠𝑓𝑜𝑟𝑚𝑎𝑑𝑎 𝑖𝑛𝑣𝑒𝑟𝑠𝑎 𝑑𝑒 𝐹𝑜𝑢𝑟𝑖𝑒𝑟 𝑓(𝑡) = 𝔍−1 {𝑓(𝑡)} = ∫ 𝐹(𝜔)𝑒 −𝑖𝜔𝑡 𝑑𝜔
2𝜋 −∞
Donde 𝑡 representa el tiempo en 𝑠 y 𝜔 la frecuencia en 𝑟𝑎𝑑/𝑠.
Con la fórmula de la transformada se obtiene una función continua en el dominio de la frecuencia 𝜔 en
rad/s. De esta forma 𝐹(𝜔) se llama espectro continuo de 𝑓(𝑡), o Transformada de Fourier de 𝑓(𝑡).
(Carrillo, 2003)
Señales determinísticas
Este tipo de señales son un tipo especial de señales estacionarias y poseen un contenido de frecuencia y
de nivel relativamente constante por un largo periodo de tiempo. Son conocidas debido a que su
descripción física es completamente conocida por su forma matemática o gráfica. Pueden ser
especificadas para cualquier instante de tiempo a través de funciones o sistemas de ecuaciones. Estas son
generadas por maquinaria rotativa, instrumentos musicales, y generadores de funciones eléctricas. Se
puede dividir en señales periódicas, estas tienen formas de ondas con un patrón que se repite a igual
distancia en el tiempo; y casi periódicas, que son producidas por maquinaria rotativa, especialmente por
equipo activado por banda. (Señales deterministicas, s.f.)
Señal de energía
Una señal de energía posee energía finita, por lo general toda señal tiene duración finita, la energía se
define como se muestra a continuación.
∞
𝐸 = ∫ 𝑥 2 (𝑡)𝑑𝑡
−∞

No se considera únicamente como energía el área bajo la curva, o integral de la señal, esto se debe a que
estas áreas pueden ser de signo negativo que pueden cancelar la media. Como ejemplo de señal de
energía tenemos un audio limitado en el tiempo, su amplitud tiende a 0 cuando |𝑡| → ∞. (Rosario, 2017)

Fig.1. Señal de energía

Si en 𝑥(𝑡) la amplitud no tiende → 0 al mismo tiempo que |𝑡| → ∞, por lo que la energía de la señal será
infinita.
Señal de potencia
Una mejor medida es el promedio de energía en el intervalo de tiempo 𝑇, si la señal es periódica, existe
𝑇.
Esta se define por
1 𝑇/2
𝑃𝑥 = lim ∫ |𝑥(𝑡)|2 𝑑𝑡
𝑇→∞ 𝑇 −𝑇/2

Como ejemplo tenemos una señal periódica, integrando en un determinado periodo, por lo que se puede
calcular la energía y la potencia. (Rosario, 2017)
 Fig.2. Señal de potencia.
Usando la fórmula de la potencia se obtiene:

Densidad espectral de potencia

Si se cambia el dominio de la función de tiempo a frecuencia, se conocería como periodo grama, llegando
así a determinar la densidad espectral de potencia de 𝑥(𝑡). (Rosario, Procesos estocásticos, 2017)
1
𝑆𝑥 (𝑓) = lim 𝐸[|𝑥̃(𝑓)|2 ]
𝑇→∞ 𝑇

La potencia promedio de 𝑥(𝑡)

 ∞
𝐸[𝑥 2 (𝑡)] = 𝑅𝑥 (0) = ∫ 𝑆𝑥 (𝑓)𝑑𝑓
−∞

La densidad espectral de potencia también se relaciona con la autocorrelación y autocovarianza a través

de la transformada de Fourier
𝑆𝑥 (𝑓) = 𝐹𝑜𝑢𝑟𝑖𝑒𝑟{𝐶𝑥 (𝜏)} + 𝑚𝑥2 𝛿(𝑓)
Y potencia cruzada 𝑆𝑥,𝑦 (𝑓) como
𝑆𝑥,𝑦 (𝑓) = 𝐹𝑜𝑢𝑟𝑖𝑒𝑟{𝑅𝑥,𝑦 (𝜏)}

DESARROLLO:

SEÑAL: 𝒙𝟐 (𝒕) = (𝒔𝒊𝒏(𝟐𝝅𝒕) + 𝒋𝒄𝒐𝒔(𝟒𝟎𝝅𝒕)

𝑇
𝑥2 (𝑡) = ∫ (𝑠𝑖𝑛(2𝜋𝑡) + 𝑗𝑐𝑜𝑠(40𝜋𝑡)) 𝑑𝑡
−𝑇
𝑇 𝑇
𝑥2 (𝑡) = ∫ (𝑠𝑖𝑛(2𝜋𝑡)) 𝑑𝑡 + 𝑗 ∫ (𝑐𝑜𝑠(40𝜋𝑡)) 𝑑𝑡
−𝑇 −𝑇

1 1 𝑇
𝑥2 (𝑡) = [ ∗ −𝑐𝑜𝑠(20𝜋𝑡) + 𝑗 ∗ −𝑠𝑒𝑛(40𝜋𝑡)]
20𝜋 40𝜋 −𝑇
1 1 1 1
𝑥2 (𝑡) = [− 𝑐𝑜𝑠(20𝜋𝑇) − (− 𝑐𝑜𝑠(−20𝜋𝑇)) + 𝑗 ( 𝑠𝑒𝑛(40𝜋𝑇) − 𝑠𝑒𝑛(−40𝜋𝑇))]
20𝜋 20 40𝜋 40𝜋

1 1 1 1
− 𝑐𝑜𝑠(20𝜋𝑇) − (− 𝑐𝑜𝑠(−20𝜋𝑇)) = − 𝑐𝑜𝑠(20𝜋) + ( 𝑐𝑜𝑠(−20𝜋))
20 20 20 20
1 1 1 1
𝑗( 𝑠𝑒𝑛(40𝜋𝑇) − 𝑠𝑒𝑛(40𝜋𝑇)) = 𝑗 ( 𝑠𝑒𝑛(40𝜋𝑇) + 𝑠𝑒𝑛(40𝜋𝑇))
40𝜋 40𝜋 40𝜋 40𝜋
1 1
= 𝑗 (2 ( 𝑠𝑒𝑛(40𝜋𝑇))) = 𝑗 ∗ 𝑠𝑒𝑛(40𝜋𝑇)
40𝜋 20𝜋

1
𝑥2 (𝑡) = 𝑠𝑒𝑛(40𝜋𝑡) → 𝑵𝒐 𝒆𝒔 𝒖𝒏𝒂 𝒔𝒆ñ𝒂𝒍 𝒅𝒆 𝒆𝒏𝒆𝒓𝒈í𝒂
20𝜋
POTENCIA

𝟏 𝑻/𝟐
𝑷𝒙 = 𝐥𝐢𝐦 ∫ [𝒙(𝒕)]𝟐 𝒅𝒕
𝑻→∞ 𝑻 −𝑻/𝟐

1 𝑇/2
𝑃𝑥 = ∫ [𝑠𝑒𝑛2 (20𝜋𝑡) + 𝑐𝑜𝑠 2 (40𝜋𝑡)]𝑑𝑡
𝑇 −𝑇/2

𝑃𝑥 = 1
La señal es periódica de periodo T=1/10

1 𝑗20𝜋𝑡 1 −𝑗20𝜋𝑡 1 𝑗40𝜋𝑡 1 −𝑗40𝜋𝑡

𝑥2 (𝑡) = 𝑒 − 𝑒 + 𝑒 − 𝑒
2𝑗 2𝑗 2𝑗 2𝑗

𝑪𝟎 = 𝟎
1
𝑪𝟏 = −𝑪−𝟏 =
2𝑗
𝑗
𝑪𝟐 = 𝑪−𝟐 =
2
FUNCIÓN DE AUTOCORRELACIÓN
∞ 𝟏 𝒋𝟐𝟎𝝅𝝉 𝟏 −𝒋𝟐𝟎𝝅𝝉 𝟏 𝒋𝟒𝟎𝝅𝝉 𝟏 −𝒋𝟒𝟎𝝅𝝉
𝟐𝝅𝒌𝝉
𝒓𝒙𝒙(𝝉) = ∑ |𝑪𝒌 |𝟐 𝒆𝒋 𝑻 𝒆 =+ 𝒆 + 𝒆 + 𝒆
𝒏=−∞ 𝟒 𝟒 𝟒 𝟒
𝒄𝒐𝒔(𝟐𝟎𝝅𝝉) + 𝒄𝒐𝒔(𝟒𝟎𝝅𝝉)
=
𝟐

CÓDIGO MATLAB
clear all
clc
fs = 1000; %Frecuencia de muestreo
Ts = 1/fs; %Periodo de muestreo
N = 500; %Número de muestras
t = 0:Ts:(N-1)*Ts; %Vector de tiempo
fc = 200; %Frecuencia de la señal
x2 = sin(20*pi*t)+(1j*cos(40*pi*t));
x=cos(40*pi*t);
figure(1)
plot(t,x2,'b');
hold on;
plot(t,x);

title('Señal x_2(t) = sin(20*pi*t)+(1j*cos(40*pi*t))');

ylabel('x_2(t)');
xlabel('Tiempo (s)');
%axis([0 0.5 -1 1.15]);
% Transformada de Fourier

F = fs*(0:(N/2))/N;
Y = fft(x2);
P = abs(Y/N);
P2 = P(1:N/2+1);
figure(2)
stem(F,P2,'r','LineWidth',1);
title ('Transformada de Fourier')
xlabel('Frecuencia(Hz)');
ylabel('Amplitud');

% Analisis espectral
% Espectro de amplitud
figure(3)
stem(abs(Y/N),'k','LineWidth',1);
title ('Análisis espectral');

% Densidad espectral de potencia

r = xcorr(x2);
w = 0:fs/N:(fs-fs/N);
DPS = (Y/N).*conj(Y/N);
figure(4)
stem(w,DPS,'g','LineWidth',1 );
%axis ([0 1000 0 1]);
title ('Densidad espectral de potencia');
xlabel('Frecuencia(Hz)');
ylabel('Potencia');

Figura 1. Grafica señal original

Figura 2. Transformada de Fourier

Figura 3. Análisis espectral

 Figura 4. Densidad espectral de potencia.

DISCUSION:

En comparativa de los resultados obtenidos de manera analítica y mediante el software de simulación

Matlab encontramos similitud en los valores obtenidos, pero cuesta diferenciar la parte negativa en la
interfaz debido a que graficamos en valores positivos.

CONCLUSIONES:

Hemos descrito de una forma general y práctica la Transformada de Fourier para su uso correcto mediante
la implementación con ordenador, mencionando el efecto que nos produce los distintos parámetros que
intervienen en el cálculo de los espectros variantes en el tiempo. La potencia media (o la energía en un
período) de la señal puede computarse en el dominio frecuencial.

Para realizar la práctica utilizamos los conceptos aprendidos, Si una señal x(t) tiene Energía Total (E) finita
y mayor que cero, se clasifica como una Señal de Energía. Estas señales tienen, además, una Potencia
Promedio igual a cero. Si la señal x(t) tiene Potencia Promedio (P) finita y mayor que cero, se clasifica
como una Señal de Potencia.

RECOMENDACIONES:

Longitud de transformada, especificada un escalar entero no negativo. Especificar un escalar entero

positivo para la longitud de transformada puede mejorar el rendimiento de fft. La longitud se especifica
normalmente como una potencia de 2 o un valor que puede factorizarse en un producto de números
primos pequeños. Si n es menor que la longitud de la señal, fft ignora el resto de valores de señal más allá
de la n-ésima entrada y devuelve el resultado truncado. Si n es 0, fft devuelve una matriz vacía.
BIBLIOGRAFÍA:

Almira, J. M. (2009). Transformada de Fourier. Jaén: Universidad de Jaén.

Carrillo, C. J. (2003). Fundamentos del análisis de Fourier. Vigo: Universidad de Vigo.
Rosario, E. d. (7 de agosto de 2017). Procesos estocásticos. Obtenido de Densidad Espectral de
Potencia Concepto: http://blog.espol.edu.ec/estg1003/densidad-espectral-de-potencia-
concepto/
Rosario, E. d. (27 de Febrero de 2017). Señales y sistemas. Obtenido de 1.9 Señales de Energía y
Potencia: http://blog.espol.edu.ec/telg1001/senales-de-energia-y-potencia/
Señales deterministicas. (s.f.). Obtenido de https://azimadli.com/vibman-
spanish/sealesdeterministas.htm