Libro Prácticas de Laboratorio Comunicaciones Analógicas
Libro Prácticas de Laboratorio Comunicaciones Analógicas
Libro Prácticas de Laboratorio Comunicaciones Analógicas
Cada capítulo al final de la práctica identifica los enlaces o referencias que fueron utilizados
para acompañar el desarrollo de la guía, se respeta los derechos de autor, los cuales son citados
a través de sus referencias donde se utilizan y se agradece al profesor Gil Hernández, Helen y
John C. Hartmann del Departamento de Ingeniería Eléctrica y de computación, del Instituto
New Jersey de Tecnología, University Heights, Newark, New Jersey, por el acceso a recursos
didácticos particularmente utilizando Radio Definida por Software (SDR) para acompañar los
ejercicios de simulación, con prácticas que utilizan dispositivos reales para la transmisión y
recepción de señales en tiempo real.
Este material es para uso exclusivo y como apoyo a la docencia en el curso de Comunicaciones
Análogas impartido en la Universidad Pedagógica y Tecnológica de Colombia – Facultad Sede
Seccional Sogamoso, su uso está permitido siempre y cuando se haga el respectivo
reconocimiento a los autores.
INTRODUCCIÓN .................................................................................................................................. VI
GLOSARIO .......................................................................................................................................... VII
PRÓLOGO ............................................................................................................................................ X
LABORATORIO 1: ANÁLISIS DE SEÑALES, POTENCIA Y ESPECTRO ........................................................ 11
INTRODUCCIÓN ........................................................................................................................................... 11
OBJETIVOS ................................................................................................................................................. 11
METODOLOGÍA ........................................................................................................................................... 12
REFERENCIAS .............................................................................................................................................. 28
LABORATORIO 2: RUIDO Y RELACIÓN SEÑAL A RUIDO........................................................................ 29
INTRODUCCIÓN ........................................................................................................................................... 29
OBJETIVOS ................................................................................................................................................. 29
METODOLOGÍA ........................................................................................................................................... 29
REFERENCIAS .............................................................................................................................................. 37
LABORATORIO 3: FILTRADO DE SEÑALES ............................................................................................ 38
INTRODUCCIÓN ........................................................................................................................................... 38
OBJETIVO................................................................................................................................................... 38
METODOLOGÍA ........................................................................................................................................... 38
REFERENCIAS .............................................................................................................................................. 56
LABORATORIO 4: MODULACIÓN Y DEMODULACIÓN DE AMPLITUD - PARTE A .................................. 57
INTRODUCCIÓN ........................................................................................................................................... 57
OBJETIVOS ................................................................................................................................................. 57
METODOLOGÍA ........................................................................................................................................... 57
REFERENCIAS .............................................................................................................................................. 74
LABORATORIO 4: MODULACIÓN Y DEMODULACIÓN DE AMPLITUD – PARTE B ................................... 75
INTRODUCCIÓN ........................................................................................................................................... 75
OBJETIVOS ................................................................................................................................................. 75
METODOLOGÍA ........................................................................................................................................... 75
REFERENCIAS .............................................................................................................................................. 88
LABORATORIO 5: DISEÑO Y ANÁLISIS DE MODULADOR Y DEMODULADOR DE FRECUENCIA ............... 89
INTRODUCCIÓN ........................................................................................................................................... 89
OBJETIVOS ................................................................................................................................................. 89
REFERENCIAS ............................................................................................................................................ 104
BIBLIOGRAFÍA .................................................................................................................................. 106
UPTC Ingeniería Electrónica Prácticas De Laboratorio Para Comunicaciones ii
Análogas
LISTA DE FIGURAS
Teniendo en cuenta que esta guía de prácticas sirve como apoyo a la docencia es indispensable
contar con una herramienta que permita determinar, analizar y evaluar los parámetros de
diseño de los sistemas de comunicaciones análogos en amplitud y frecuencia, dado que son
esenciales en la formación de ingenieros electrónicos y de telecomunicaciones.
El desarrollo de habilidades que se pretende adquirir con estas prácticas de laboratorio, les
permitirá a los participantes mejorar sus destrezas en programación y análisis de señales en
cursos posteriores como los cursos de comunicaciones digitales, radio-comunicaciones,
comunicaciones ópticas, etc. Para ello se ha incorporado ejercicios donde se pide generar
señales, aplicarlas para encontrar respuestas a sistemas, operar con estas señales aplicando
principios y teoremas matemáticos fundamentales, y la teoría de comunicaciones análogas para
la modulación y demodulación de amplitud y frecuencia en sistemas analógicos.
El texto inicia en el capítulo 1 con una explicación básica de sistemas y señales, energía y
potencia, convolución, y autocorrelación entre otros, con ejercicios que permiten entender
estos conceptos desde la implementación y generación gráfica de resultados para su análisis. En
el capítulo 2, se enfoca la práctica a la generación del ruido, la relación señal a ruido y el efecto
particular en las señales cuando se contaminan con ruido, que es propio a los sistemas de
comunicaciones. El capítulo 3, aborda los conceptos de filtros analógicos y hace uso de las
herramientas de Matlab para su diseño e implementación. Por practicidad, se propone diseños
específicos de filtros tanto analógicos como digitales y algunos ejercicios de aplicación para
observar el resultado de esta etapa particularmente en la selección de un canal en sistema
multiportadora. Los últimos tres capítulos, introducen los conceptos de la modulación en
amplitud y fase, y luego mediante radios digitales usando el concepto de la Radio Definida por
Software (SDR), se implementan bloques de transmisión y recepción de señales usando técnicas
de modulación analógicas, modulación en amplitud y frecuencia. El uso de dispositivos
periféricos de radio software universal (USRP) ha tomado fuerza, dado que la digitalización de
señales habilitó la migración de los servicios analógicos a digitales en la mayoría de sistemas de
telecomunicaciones. Estos radios incorporan procesadores digitales de señal, con interfaces
transmisor/receptor que permiten traslado de señales de banda base a pasabanda, y viceversa,
filtrado, modulación en fase y cuadratura, etc., con lo cual se puede diseñar, evaluar y emular
sistemas de comunicaciones analógicos y digitales.
Amplificador: Circuito electrónico que multiplica por un factor la señal de entrada para
entregar a la salida una versión escalada y proporcional al factor de amplificación.
Ancho de Banda: Rango de frecuencias con límite inferior y superior que ocupa algún sistema
de comunicaciones en el espectro electromagnético.
Atenuación: Efecto producido sobre la amplitud de una señal por la pérdida de energía durante
la propagación a través del canal de comunicaciones.
Banda Base: Ancho de banda de la señal de información pura, sin modular, tal y como se
genera desde la fuente.
Canal: Referente al medio de comunicación, transmisión o al ancho de banda del canal de voz.
Distorsión: Cambio que sufre una señal respecto a otra originalmente transmitida debido al
proceso de generación, transmisión, propagación, detección y recuperación de la señal de
información.
Frecuencia: Valor numérico que expresa el número de veces que se repite un ciclo de onda en
un periodo por segundo.
Índice de modulación: Es una medida de la variación de amplitud que rodea una portadora no
modulada.
Interferencia: Fenómeno físico en el que dos o más ondas se superponen para formar una onda
resultante de mayor o menor amplitud.
Modulación: Técnicas que se usan para trasportar información mediante el uso de una señal
portadora a una frecuencia más alta que la señal de información, generalmente es de tipo
sinusoidal.
Oscilador Local: Generador de la señal de referencia o de la señal portadora que se mezcla para
entregar otra señal de mayor o menor frecuencia.
Radio receptor: dispositivo del lado del receptor de un sistema de comunicación por radio que
demodula las señales transmitidas, particularmente en amplitud y frecuencia.
Sistema analógico: Sistema de comunicaciones que usa señales analógicas en los procesos de
modulación y transmisión.
Sistema no lineal: sistema que tiene un comportamiento que no puede expresarse en términos
de sus descriptores. Es decir, que estos comportamientos no pueden definirse mediante
dinámicas de modelos lineales.
Tono de prueba: Señal de frecuencia pura usada para calibrar los niveles de voltaje o potencia
en dispositivos de los sistemas de comunicaciones.
Introducción
Objetivos
Crear y analizar señales continuas en tiempo, señales discretas, señales aleatorias, y determinar
parámetros de interés como potencia, ancho de banda, periodo, etc.
Las señales se definen como magnitudes físicas que varían con el tiempo, el espacio, o cualquier
otra variable independiente. Es decir, que son patrones de variación que representan
información codificada. En comunicaciones son muy utilizadas desde la generación de la
información, su procesamiento para adatarlo a las características del canal, y luego para la
reconstrucción y recuperación en el lado del transmisor.
1. Generación de señales
1 1
0.5 0.5
0 0
-0.5 -0.5
-1 -1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 -10 -8 -6 -4 -2 0 2 4 6 8 10
Tiempo [Segundos] n [Muestras]
Fig. 1.1. Tren de pulso rectangulares Fig. 1.2. Tren de pulso rectangulares
continuo en tiempo. discreto en tiempo.
Fig. 1.4. Función exponencial decreciente. Fig. 1.5. Función exponencial creciente.
Fig. 1.6. Onda sinusoidal continua en Fig. 1.7. Onda sinusoidal discreta en
tiempo. tiempo.
Los comandos usados entonces para generar una función sinusoidal amortiguada
exponencialmente son como sigue:
Funciones escalón unitario, impulso o delta de Dirac y rampa: La función escalón unitario es
generada a partir de los comandos:
% Escalón unitario
u=[zeros(1,50),ones(1,50)];
plot(u)
Fig. 1.9. Función Escalón unitario. Fig. 1.10. Función delta de Dirac discreto.
UPTC Ingeniería Electrónica Ruido y Relación Señal a Ruido 15
% Secuencia discreta
n=-10:10;
x=5*exp(-0.2*n);
stem(n,x)
z=x.*y; % y es un coseno
Amplitud
stem(n,z)
Matlab tiene implementadas dos funciones para determinar la Transformada directa de Fourier
(DTF) y la Inversa (IDTF) [1]:
X fft ( x) y x ifft ( X ) (2)
x y X son vectores de N elementos cada una, el primero con las muestras temporales y el
segundo con los resultados de la DTF, la DTF aproxima el espectro de una forma de onda w(t).
Donde los cruces por cero de la función sinc(x) se dan para To f n , de donde f n To .
Utilice Matlab para determinar una aproximación al espectro del pulso rectangular, a
continuación, se presenta un ejemplo posible
En el código de ejemplo anterior se tomó M = 7, lo que da lugar a una FFT (Transformada rápida
de Fourier) de 128 puntos. Igualmente, se calculó el espectro para un pulso con duración de un
segundo, muestreando sobre una duración total de 10 segundos. En la figura 13, se observa que
los nulos o cruces por cero se encuentran en frecuencias de 1 Hz, 2 Hz, … n Hz, dado el tiempo
T0 que se encuentra en un nivel de amplitud unitario la señal. Por otra parte, el espectro de
magnitud de un pulso rectangular es teóricamente infinito, aunque a partir del 5º nulo se
puede considerar que sus valores son prácticamente despreciables.
Teniendo en cuenta esto, y dado que el 5º nulo se da para f = 5/T0, habría que muestrear la
señal a una tasa lo suficientemente rápido de tal forma que el criterio de Nyquist [1] se cumpla:
f s 2B (4)
donde B, es el ancho de banda de la señal que se desea muestrear. Se pretende demostrar que
la selección apropiada de la frecuencia de muestreo es un parámetro muy importante para
evitar que se tomen muy pocas muestras y no sea posible tener una versión aproximada de la
señal que se quiere digitalizar, o también no cambiar la forma de onda debido a cambios suaves
que puedan inducirnos en un error. Con lo anterior se hace referencia a señales cosenoidales
que pueden recuperarse, pero se altera el periodo de la señal.
Continuando con el análisis del ejercicio, si se toma B ≈ 5/T0, la frecuencia de muestreo debería
ser fs ≥ 10/T0, es decir, requeriríamos aproximadamente 10 muestras del pulso.
Por favor, comente de forma concreta respecto al efecto en el cambio del tamaño de la
transformada de Fourier o del periodo total de duración de la señal, esto dado que al cambiar
de dominio existe un compromiso importante que debe ser interiorizado.
3. Estimación de la Potencia
El espectro de una senoide se puede obtener a través del teorema de Parseval [2], por
ejemplo, si se desea determinar el espectro de la señal dada por:
w(t ) 3sen( wot 20) (5)
Se puede representar la anterior expresión en términos de la exponencial compleja como:
3 3
w(t ) e j 20 e jwot e j ( 20) e jwot (6)
2j 2j
Donde los coeficientes de Fourier vienen dados por
3
c1 e j 20 1.5 70
2j
(7)
3
c2 e j ( 20) 1.5 70
2j
Y los demás coeficientes son cero.
Fig. 1.14. Análisis de señal senosoidal para obtener espectro de magnitud y fase.
4. Autocorrelación
La autocorrelación se define [2] como la correlación cruzada de la señal consigo misma. Esta
función, ofrece gran utilidad para encontrar patrones repetitivos dentro de una señal, por
ejemplo, para determinar inicio de una secuencia discreta en un vector de datos binarios
que es típico en comunicaciones digitales para encontrar sincronismo.
Cargue las señales usando el comando load SensorData, observe que estas señales están
muestreadas a una tasa fs = 11025 Hz. Cree el vector de tiempo y grafique las señales, debe
obtener unas señales como se muestra en la figura 16. ¿Comente respecto a cada señal
(amplitud, variabilidad, etc.), temporalmente, cual señal ocurrirá primero?
La correlación cruzada de las dos señales de los sensores es máxima en la muestra igual al
retardo. Para ello se debe expresar el retardo como el número de muestras y en segundos.
Explique y comente el código siguiente:
[acor,Retardo] = xcorr(s2,s1);
[~,I] = max(abs(acor));
DifRetardo = lag(I)
¿Cómo se puede obtener el tiempo de diferencia entre las dos señales? ¿Cuál es el tiempo de diferencia
entre las dos señales?
DifTiempo = ________/Fs;
figure
plot(Retardo,acor)
a3 = gca;
a3.XTick = sort([-3000:1000:3000 DifRetardo]);
Ahora debe alinear las dos secuencias y dibújelas usando subplot, deberá obtener la señal que se
muestra en la figura 17. Adicione 1 a la diferencia de retardo para tener en cuenta el sistema indexación
usado por Matlab, y obtendrás las dos secuencias alineadas en tiempo, eliminando el retardo existente
entre las mediciones de los dos sensores.
5. Convolución
La convolución [2] es una operación matemática lineal que combina dos señales y las
transforma para producir una tercera señal de salida. Tiene gran importancia en los
sistemas de comunicaciones donde usualmente se multiplican (proceso de convolución)
señales en el dominio del tiempo, pero por dualidad en el dominio de la frecuencia la
convolución corresponde al producto de las transformadas de Fourier de cada señal, lo que
simplifica considerablemente el procedimiento operacional.
Convolución de pulsos rectangulares: Suponga que las secuencias x[n] y v[n] definen un
pulso rectangular dado por:
1, 0n9
p n (8)
0, otro valor de n
i) Encuentre el producto de convolución para las secuencias x[n] y v[n] utilizando la función
de Matlab conv. Los comandos en Matlab para calcular la convolución para este caso son:
% Convolución Continua
clc
clear
pause on
x=[0 1 2 3 4 3 2 1 0]; %Respuesta al impulso unitario
h=[1 1 1 1 1 1 1 1 ]; %Señal de entrada
m=length(x);
n=length(h);
%invierte el vector h
hi=fliplr(h);
k=20;
X=[x,zeros(1,40-m)];
X= X([ end-k+1:end 1:end-k ]); % Refleja la señal y la atrasa
H=[h,zeros(1,40-n)];
H= H([ end-k+1:end 1:end-k ]); % Respuesta al impulso
xn=-20:20-1;
Y=zeros(1,40);
p=zeros(1,40);
h1=subplot(3,1,1);
plot(xn,X,'-b')
ylabel('h(t)')
title('Respuesta al impulso unitario')
% colocar(h1,'YLim',[a b])
h2=subplot(3,1,2);
plot(xn,H,'-k')
% colocar(h2,'YLim',[a b])
pause(3)
h2=subplot(3,1,3);
plot(xn,Y,'-m')
%ciclo
Hi=[hi,zeros(1,40-n)];
for i=1:40-n
p=X.*Hi;
Y(i+n-1)=sum(p);
subplot(3,1,2);
plot(xn,Hi,'-r')
ylabel('x(t)')
title('Señal de entrada')
subplot(3,1,3)
plot(xn,Y,'-m')
Mientras que la convolución discreta en tiempo se obtiene a partir del siguiente script.
xn=-20:20-1;
Y=zeros(1,40);
p=zeros(1,40);
h1=subplot(3,1,1);
stem(xn,X,'MarkerFaceColor','red')
ylabel('h[n]')
title('Respuesta al impulso unitario')
% set(h1,'YLim',[a b])
h2=subplot(3,1,2);
stem(xn,H,'MarkerFaceColor','blue')
% set(h2,'YLim',[a b])
pause(3)
h2=subplot(3,1,3);
stem(xn,Y,'MarkerFaceColor','green')
%
Hi=[hi,zeros(1,40-n)];
UPTC Ingeniería Electrónica Ruido y Relación Señal a Ruido 26
for i=1:40-n
p=X.*Hi;
Y(i+n-1)=sum(p);
subplot(3,1,2);
stem(xn,Hi,'MarkerFaceColor','blue')
ylabel('x[n]')
title('Señal de entrada')
subplot(3,1,3)
stem(xn,Y,'r','MarkerFaceColor','green')
xlabel('Tiempo [n]')
ylabel('y[n]')
title('Señal de salida')
Hi= Hi([ end 1:end-1 ]);
pause(0.5)
end
h[n]
x[n]
y[n]
[1] Haykin, Simon., and Barry Van Veen. “Señales y Sistemas”. 1a. ed., 1a. reimp. Mexico: Limusa Wiley,
2006.
[2] Oppenheim, A. V., Willsky, A. S., & Young, I. T.. Signals and systems. Englewood Cliffs, N.J. Prentice-
Hall, 1983.
Introducción
Objetivos
Utilizar las diferentes funciones de Matlab relacionadas con el ruido aplicadas en el análisis del
nivel de potencia de señal y de ruido.
Metodología
a) Revise la ayuda de Matlab usando help randn, rand, randi, rng, RandStream y aplique en un
ejemplo donde el efecto sea evidente en señales típicas. Debe graficar rotulando
apropiadamente los ejes y comentando el código.
b) Vamos a calcular la SNR de un pulso rectangular de duración 20ms, muestreado a 10kHz en
la presencia de ruido Gaussiano. Ajuste el generador aleatorio a la configuración por
defecto para datos reproducibles.
rng default
Tpulse = 20e-3;
Fs = 10e3;
t = -1:1/Fs:1;
x = rectpuls(t,Tpulse);
y =
0.00001*randn(size(x));
s = x + y;
pulseSNR = snr(x,s-x)
pulseSNR =
80.0818
Fig. 2.1. Código Matlab y pulso rectangular ideal y contaminado con ruido.
d) Calcule la SNR de una sinusoide con frecuencia de operación 2.5 kHz muestreada a 48 kHz.
Adicione ruido con una varianza de 0.001². En ese caso el argumento de la función
sinusoidal será de la forma wo*Fo/Fs*(1:N). ¿Cuál es el nivel de relación SNR? Dibuje el
espectro de un solo lado (DEP positiva – use la función pwelch - onesided) y adicione
f) SNR con y sin armónicos aliased (señales continuas distintas que se tornan indistinguibles
cuando se muestrean digitalmente). Genere una señal que parezca la salida de un
amplificador no lineal débilmente con un tono a 2.1 kHz como entrada. La señal es
muestreada por 1 segundo a 10 kHz. Calcule y dibuje el espectro de la señal. Use el
periodograma con la ventana Kaiser del punto e. Identifique claramente los picos de la
UPTC Ingeniería Electrónica Ruido y Relación Señal a Ruido 32
frecuencia de interés y los armónicos que deberían aparecer con marcadores y adicione
rótulos para marcarlos. Calcule la SNR de la señal con snr(x,Fs,7); y con snr(x,Fs,7,'aliased').
¿Qué observa? Explique y comente al respecto.
g) La potencia de ruido corresponde a la varianza también tal como se planteó anteriormente.
Cree una sinusoide muestreada a 18kHz, configure una amplitud unitaria y frecuencia
fundamental de 1kHz. Además, contiene un armónico con media amplitud y ruido aditivo
con varianza 0.12. Calcule la potencia de ruido de la señal y verifique que concuerda con la
definición teórica.
h) Reconocidas las diferentes funciones y aplicada a señales sinusoidales vamos a aplicar ruido
aditivo blanco Gaussiano (Additive White Gaussian Noise - AWGN) para evaluar el impacto
en la generación de una señal rampa. Tome un segmento de tiempo de 10 s en incrementos
de 0.1.
Potencia (dB)
Mediante la función sawtooth cree una rampa y adicione ruido con la función awgn con
niveles de snr = [0 10 20 30] dBW usando ‘measured’ como signalpower. Grafique y rotule la
señal para cada uno de los cuatro casos utilizando subplot. ¿Qué puede decir de los cuatro
niveles de SNR?
A continuación, se muestra el efecto de la adición de ruido para una señal rampa con una
SNR=10dB.
3. Antes de entrar a evaluar algunos efectos de los sistemas, es conveniente definir lo que es
un sistema lineal y uno no lineal [3]:
y t k1 x t k2 x 2 t k3 x 3 t (5)
a) Distorsión de un tono: Sea x(t)=x1(t) siendo A1=1.5 y f1=20Hz. Obtenga una figura con
subplot para i) la señal en el dominio del tempo, ii) su espectro potencia, iii) la salida del
sistema y el espectro de potencia de la salida del sistema. Se recomienda hacer la
gráfica en el dominio del tiempo y luego comparar usando stem que permitiría una
mejor visualización para contar armónicos. ¿Qué diferencias se observan entre la señal
de salida y la señal original? ¿Qué tipo de componentes aparecen?
b) Distorsión de dos tonos: Sea x(t)=x1(t) + x2(t), siendo A1=1.3, A2=0.5, f1=5Hz y f2=5f1.
Obtenga una figura usando subplot de forma similar al punto a. ¿Qué diferencias se
observan entre la señal de salida y la señal original? ¿Qué tipo de componentes
aparecen en el espectro? Diferencie con colores cada una de las señales y rotule
apropiadamente los ejes.
c) Para los dos casos anteriores adicione AWGN para SNR = [0 10 20 30] dBW y grafique
usando subplot las figuras resultantes en cada caso. ¿Cómo afecta el ruido la forma de
onda y la potencia de la señal deseada?
Para ilustrar el fenómeno de la distorsión no lineal [1,3], se tiene una señal cosenoidal con
frecuencia 60 Hz y amplitud 2, esta señal se pasa a través de un sistema no lineal con salida
y(t)=5x(t)-2x2(t)+4x3(t). Se observa en la figura 2.6 que la señal presenta una compresión en
el semiciclo positivo lo que cambia la forma de onda y que corresponde al efecto de la suma
de componentes lineales, cuadráticas e incluso cúbicas.
En la figura 2.7, se superpone los dos tonos (color rojo) que se usan como entrada al sistema no
lineal, con frecuencias de portadora a 1.5 kHz y 25 kHz. Y los tonos en color azul, corresponden
a la salida del sistema no lineal que incluye armónicos de tercer orden. Replique el espectro de
salida para los valores presentados anteriormente y comente qué más señales aparecen, en
que frecuencias ocurren y desde el análisis matemático a que se deben estos picos adicionales.
Potencia/Frecuencia (dB/Hz)
Potencia (dB)
Fig. 2.8. Distorsión armónica total para dos tonos y sus productos de intermodulación.
Referencias
[1] Carlson, Bruce, and Paul B. Crilly. “Communication Systems”, Ed. Mc Graw Hill 5th Ed. New
York 2009.
[2] The Mathworks Inc® Signal to Noise Ratio. (Enlace recuperado el 28 de abril de
https://la.mathworks.com/help/signal/ref/snr.html)
[3] The Mathworks Inc® Total Harmonic Distortion (Enlace recuperado el 28 de abril de
https://la.mathworks.com/help/signal/ref/thd.html)
Introducción
Los filtros existen en todos los sistemas de comunicaciones, sirven para diseñar y permitir el
paso de una banda de frecuencias específicas, mientras atenúan las demás señales fuera de
esta banda. Los filtros son usualmente clasificados de acuerdo al rango de filtrado, respuesta en
frecuencia en la banda de paso, y componentes de circuitos. En la clasificación por el rango de
filtrado, existen cuatro tipos de filtros: filtro paso bajo, paso alto, pasabanda y filtros rechaza
banda.
Por otra parte, de acuerdo a la respuesta en frecuencia en la banda de paso, existen dos tipos
de filtros: activos y pasivos. Los filtros pasivos, están constituidos por circuitos que contienen
solamente componentes pasivos (resistores, inductores y capacitores) conectados en tal forma
que permitirán pasar ciertas frecuencias mientras que rechazan otras. Los filtros activos por su
parte emplean dispositivos activos (transistores o amplificadores operacionales) combinados
con resistencias capacitores e inductores. Los filtros activos son los más utilizados en los
sistemas de comunicaciones modernos, debido a que tienen las siguientes ventajas:
1. Debido a que la función de transferencia con características inductivas puede obtenerse con
un diseño particular, se puede usar resistencias en lugar de inductancias.
3. Dado que los componentes activos proporcionan amplificación, entonces los filtros activos
pueden incorporar ganancias.
Objetivo
Metodología
El diseño e implementación de filtros analógicos y digitales es una tarea importante que debe
aprenderse y aplicarse en la vida práctica pues sistemas de instrumentación, en bioingeniería,
electrónica de potencia, y telecomunicaciones entre otros, de forma constante los utilizan en
procesamiento analógico y digital. Los filtros son un arreglo de dispositivos que permiten
fs 2 fm (1)
Fig. 3.1. Respuesta en amplitud de Filtro ideal vs. Filtros Prácticos [1].
Los filtros ideales paso bajo deberían tener una banda de transición y una ganancia plana
perfecta en la pasabanda como se muestra por el pulso cuadrado en línea punteada en la figura
3.1. Además, la atenuación de la banda de parada debería reducir cualquier señal fuera de
banda a cero. En la figura 3.1 también se observa la respuesta de algunos filtros analógicos
prácticos comúnmente, identificando las diferentes zonas de transición, como la banda de
paso, la caída de -3 dB, la banda de transición y la banda de para. La atenuación de la banda de
parada no es aquella que limita el ruido fuera de banda, sino que puede tener una banda de
transición sin caída abrupta, cuyo ruido degrada la atenuación alrededor de la frecuencia de
corte. Adicionalmente, todos los filtros no ideales, introducen un retardo de fase o de grupo.
Fig. 3.2. Banda de transición de filtros Butterworth ideales para diferentes órdenes [1].
Pero exige un diseño complejo debido a la alta sensibilidad en la ganancia de acople para ser
prácticos solamente algunos ordenes de atenuación de magnitud. Además, cualquier
componente tal como una resistencia o amplificador operacional, adicionado en la cadena de
señal introducirá ruido en banda. De acuerdo a la figura 3.2, un filtro RC de primer orden podría
ser usado en la protección de canales, mientras que un filtro RC de segundo orden proporciona
una mejor planicidad en banda y un factor de caída más agresivo en la transición para la medida
de los canales.
Función de transferencia
Los coeficientes de la función de transferencia del filtro, retornados como vectores fila de
longitud n+1 para filtros paso bajo y paso alto, y longitud 2n +1 para filtros pasabanda y
rechaza-banda.
H ( z)
A( z ) a 1 a 2 z 1 a n 1 z n
H (s)
A( s ) a 1 s n a 2 s n1 a n 1
Ceros, polos y ganancia
Los ceros, polos y la ganancia de un filtro, retornados como dos vectores columna de longitud n
(2n para diseños pasabanda y rechaza-banda) y un escalar.
Representación en espacio de estados del filtro, retornado como matrices. Si m=n para diseños
paso bajo y paso alto, y m=2n, para filtros pasabanda y rechaza-banda, entonces, A es mm, B
es m1, C es 1, y D es 11.
Para filtros digitales, las matrices en espacio de estados relacionan el vector de estado x, la
entrada u, y la salida y, a través de:
x k 1 Ax( k ) Bu ( k )
y (k ) Cx(k ) Du
De acuerdo a lo anterior, diferentes comandos en Matlab serán usados para cada caso de
acuerdo a las necesidades. A continuación, se ilustra el diseño de filtros paso bajo y pasabanda
usando los modelos más comunes como lo son Chebyshev y Butterworth.
Respuesta en frecuencia y al impulso de filtro paso bajo FIR e IIR usando objetos y fvtool
Cree un filtro paso bajo FIR de orden mínimo para datos muestreados a 44.1 kHz. Especifique
una frecuencia pasabanda de 8 kHz, una frecuencia en la banda de parada de 12 kHz, un ripple
en la banda de paso de 0.1 dB, y atenuación en la banda de parada de 80 dB. Use para ello la
herramienta fvtool [2].
Diseñe un Filtro paso bajo IIR de orden mínimo con las mismas propiedades del filtro FIR.
Grafique la respuesta al impulse del filtro FIR paso bajo. El coeficiente de orden cero está
retrasado 19 muestras, que es igual al retardo de grupo del filtro. El filtro paso bajo FIR es un
filtro FIR causal.
fvtool(FIRLPF,'Analysis','impulse')
fvtool(FIRLPF,'Analysis','freq')
Y la respuesta en frecuencia para la magnitud y fase del filtro paso bajo IIR
UPTC Ingeniería Electrónica Filtrado de Señales 43
fvtool(IIRLPF,'Analysis','impulse')
fvtool(IIRLPF,'Analysis','freq')
Determine el costo computacional de implementación de los filtros FIR e IIR, para ello use el
comando cost, además, usando también grpdelay puede determinar el retardo de grupo o de
fase de cada filtro. Luego se podrá también observar usando la herramienta fdaTool o filter
Designer en el ícono i de información del filtro.
0 -0.022
-10 -3.544
-20 -7.066
-30 -10.588
-40 -14.11
-50 -17.632
-60 -21.154
-70 -24.676
-80 -28.198
-90 -31.72
0 5 10 15 20
Frecuencia (kHz)
Fig. 3.4. Respuesta en magnitud y fase para el filtro paso bajo diseñado.
Diseñe un filtro paso bajo Chebyshev tipo I de orden 6, con 10 dB de rizado (ripple) de banda de
paso y una frecuencia de transición de 300 Hz, que, para datos muestreados corresponde a
1000 Hz, o equivalentemente 0.6 rad/muestra (wc=2fc=2300), que normalizado en
rad/muestra es 1000/600=0.6. ¿Cómo obtiene esta frecuencia angular muestreada? Dibuje sus
respuestas en magnitud y fase. Use el filtro para una señal aleatoria de 1000 muestras. Revise
help cheby1 para mayor detalle en los parámetros de entrada y salida. ¿Para qué se utiliza el
comando freqz? Usando filter puede filtrar los datos en el vector de entrada X con los
parámetros descritos por el vector de coeficientes A y B, revise la ayuda del comando para
tener mayor detalle.
Usando hold y formateando la figura muestre la señal aleatoria generada y filtrada, ¿y comente
si observa alguna particularidad, que se aprecia? ¿Modifique ahora el rizado a 5dB (Mantenga
la gráfica superpuesta con la señal original y la filtrada de rizado 10dB) y analice nuevamente,
que cambia? Por tanto, concluya respecto al efecto del rizado en la señal filtrada.
Amplitud
Amplitud
Fig. 3.5. Señal con muestras de entrada aleatorias y filtradas paso bajo.
Diseñe un filtro pasabanda Chebyshev Tipo I de orden 20 con una frecuencia pasabanda inferior
de 500Hz y una frecuencia pasabanda mayor de 560Hz. Especifique un rizado pasabanda de 3dB
y una tasa de muestreo de 1500Hz. Use la representación en espacio de estados y diseñe un
filtro idéntico usando la función designfilt [2]. Revise la ayuda de Matlab para reconocer los
parámetros de entrada de la función designfilt.
sos = ss2sos(A,B,C,D);
fvt = fvtool(sos,d,'Fs',1500);
legend(fvt,'Filtro Pasabanda Chebyshev Tipo 1','designfilt')
Diseñe un filtro paso bajo Buttherworth con frecuencia de corte a 300 Hz, para datos
muestreados a 1000 Hz (0.6 rad/muestra). El diseño no debe tener más de 3dB de ripple en la
banda de 0 a 40 Hz, y al menos 60 dB de atenuación en la banda de parada. Encuentre el orden
del filtro y la frecuencia de corte. Para ello revise la ayuda del comando buttord.
Magnitud (dB)
Fase (Grados)
Fig. 3.6. Respuesta en magnitud y fase para el filtro paso bajo diseñado.
[n,Wn] = buttord(Wp,Ws,3,60)
n = 5
Wn = 0.0810
[z,p,k] = butter(n,Wn);
sos = zp2sos(z,p,k);
Diseñe un filtro pasabanda con una banda de paso de 100 a 200 Hz, con al menos 3 dB de ripple
en la banda de paso y al menos 40 dB de atenuación en la banda de parada. Especifique una
tasa de muestreo de 1 kHz. Ajuste el ancho de la banda de parada a 50 Hz a los dos lados de la
pasabanda. Encuentre el orden del filtro y las frecuencias de corte.
Wp = [100 200]/500;
Ws = [50 250]/500;
Rp = 3;
Rs = 40;
[n,Wn] = buttord(Wp,Ws,Rp,Rs)
[z,p,k] = butter(n,Wn);
sos = zp2sos(z,p,k);
Ajuste a los parámetros del diseño del filtro Chebyshev [3] pasabanda anterior y compare el
desempeño superponiendo las dos respuestas de amplitud frecuencia y de fase. Reflexione al
respecto y explique lo que analiza. La respuesta en magnitud y fase se muestra en la figura 3.7.
-100
Magnitud (dB)
-200
-300
-400
-500
0 50 100 150 200 250 300 350 400 450 500
Frecuencia (Hz)
1000
500
Fase (Grados)
-500
-1000
0 50 100 150 200 250 300 350 400 450 500
Frecuencia (Hz)
Fig. 3.7. Respuesta en magnitud y fase para el filtro paso bajo diseñado.
Ejercicio: Diseñe un filtro pasabanda Butterworth de orden 20 con una frecuencia de corte
inferior de 500Hz y de corte superior 560Hz. Especifique una tasa de muestreo de 1500 Hz. Use
la representación en espacio de estados y diseñe un filtro idéntico usando designfilt. Convierta
la representación en espacio de estados a secciones de segundo orden, y visualice las
respuestas en frecuencia usando fvtool.
La aplicación filterDesigner incorpora características del toolbox de Sistema DSP [2-3]. Esto
permite diseñar y analizar filtros cuantizados, así como convertir filtros cuantizados a diferentes
estructuras, transformaciones de filtros, diseños de filtros multi-tasa, e implementar modelos
de filtros.
Reconozca las diferentes estructuras que hay explorando la barra lateral y las opciones en la
paleta superior tal como se muestra en la figura.
Crear un filtro multi-tasa (multirate filter) — Dando clic a este botón se cambia el filtro
diseñador a modo filtro multi-tasa de forma que puede diseñarse filtros interpoladores,
decimadores, y filtros de cabio de tasa fraccional.
Realización del modelo (Realize Model) — inicia desde su filtro FIR, de forma directa y
cuantizada, dando clic en este botón se crea un modelo en Simulink de la estructura del filtro
en una nueva ventana.
Además, si selecciona Store Filter, podrá guardar el diseño en el administrador de filtros. Por
otra parte, podrá exportar al espacio de trabajo (Workspace) los coeficientes del filtro, objeto y
objeto del sistema para poder cargar y operar con las señales que desea acondicionar, debe
asignar un nombre en Variables Names – Numerator: sugerido por defecto num, pero se
recomienda siempre dar un nombre relacionado con el diseño, por ejemplo,
LPF_n7_fc500_fs_10k, de esa forma además tendrá presente los parámetros usados en el
diseño. De igual forma puede exportar los coeficientes como Objects y el nombre por defecto
será Hd.
1. Diseñe un filtro Butterworth paso bajo, estable y de orden mínimo, con una frecuencia
pasabanda de 1kHz y una frecuencia de banda de paso de 1.4kHz. Introduzca una
atenuación de 1dB en la frecuencia pasabanda y de 80dB en la frecuencia de banda de
parada.
2. Diseñe un filtro Chebyshev tipo I, paso bajo, estable y de orden mínimo, con las mismas
características del punto 1.
UPTC Ingeniería Electrónica Filtrado de Señales 50
3. Diseñe un filtro FIR paso bajo usando la ventana Blackman con una frecuencia de corte de
1kHz. Especifique el orden del filtro, tal que el primer mínimo en la banda de parada
(precediendo el primer lóbulo) es tan cercano a 1.4kHz como sea posible sin excederlo.
4. Diseñe un filtro pasabanda (BPF) que limite la banda de audio digitalizado a una tasa de
muestreo de 48kHz, y que se extiende de 1800-12000 Hz. Asuma una Fstop1=1500 Hz,
Fpass1=1800 Hz, Fpass2= 12000 Hz, Fstop2=13000 Hz, Astop1=20dB, Apass=3dB y
Astop2=20dB. Observe en particular que para este diseño si se considera un filtro elíptico,
¿cuál sería el orden del filtro? En general, si no hay banda de paso o ripple en la banda de
rechazo, se debe seleccionar el diseño Butterworth (maximally flat). Si se requiere un
diseño con retardo constante sobre la banda de paso (desplazamiento de fase lineal) se
requiere eliminar la distorsión de retardo de la forma de onda en banda, entonces, se debe
seleccionar el diseño Bessel. Ahora, en el menú análisis seleccione Group Delay Response, y
obtenga la gráfica con los diferentes retardos a diferentes frecuencias en la banda de paso
de 1800 a 12000 Hz. Esto significa que el filtro introduce alguna distorsión en la forma de
onda de las señales en banda multi-frecuencia, puesto que diferentes componentes en
frecuencia dentro de la banda de paso llegarán a salida del filtro en diferentes instantes de
tiempo. Esto no es problema cuando se filtra señales de audio, puesto que el oído humano
tiene una propiedad de ser esencialmente “insensitivo a la fase”. Finalmente, obtenga la
gráfica de polos y ceros, es decir, la función de transferencia del filtro H(z).
Preguntas:
1. ¿Cuál es el orden del filtro paso bajo Butterworth que diseño anteriormente?, y, ¿cuál el
orden del filtro Chebyshev?
2. Compare el uso de memoria de cada uno de los tres primeros filtros asumiendo se empleó
una realización de Forma Directa II. ¿Cuán ineficiente es la técnica de ventaneo? ¿Cuánto
más costoso es la técnica de ventaneo respecto al mejor filtro IIR?
3. ¿Cuáles son los retardos de grupo para cada uno de los tres primeros diseños de filtros
utilizados anteriormente?
Este ejercicio muestra cómo diseñar e implementar un filtro FIR usando dos funciones en líneas
de comando, fir1 y designfilt, y la aplicación interactiva Filter Designer [3].
Cree una señal para usar en los ejemplos. La señal es una onda coseno de 100Hz con ruido
aditivo blanco Gaussiano normal N(0,0.25). Configure el número generador aleatorio (randn) al
El diseño de filtro es un FIR paso bajo con orden 20 y una frecuencia de corte de 150Hz. Use una
ventana Kaiser con longitud de una muestra mayor que el orden del filtro y =3. Vea la ayuda
de la función kaiser para más detalles sobre los parámetros de entrada y salida.
Use fir1 para diseñar el filtro [3], esta función requiere frecuencias normalizadas en el intervalo
[0,1], donde 1 corresponde a rad/muestra. Use fir1, debe convertir todas las especificaciones
de frecuencia a frecuencias normalizadas.
fc = 150;
Wn = (2/Fs)*fc;
b = fir1(20,Wn,'low',kaiser(21,3));
fvtool(b,1,'Fs',Fs)
Ahora, aplique el filtro a la señal y obtenga el resultado para los primeros diez periodos de la
sinusoide de 100Hz.
y = filter(b,1,x);
plot(t,x,t,y)
xlim([0 0.1])
xlabel('Time (s)')
ylabel('Amplitude')
legend('Original Signal','Filtered Data')
Diseñe el mismo filtro usando designfilt. Observe la respuesta del filtro a ‘lowpassfir’ e ingrese
las especificaciones como los pares Name, Value. Con designfilt, usted debe especificar el
diseño del filtro en Hz.
Fs = 1000;
Hd = designfilt('lowpassfir','FilterOrder',20,'CutoffFrequency',150, ...
'DesignMethod','window','Window',{@kaiser,3},'SampleRate',Fs);
y1 = filter(Hd,x);
plot(t,x,t,y1)
xlim([0 0.1])
Se tiene un sistema de comunicaciones que transmite por dos portadoras, una centrada a 8
MHz y la otra en 12 MHz. Se pide seleccionar el canal de 8 MHz, por tanto, deberá diseñar y
validar un filtro pasabanda para la frecuencia de interés.
Fig. 3.12. Señal con portadoras en 8 y 12 Fig. 3.13. Respuesta en magnitud del filtro.
MHz.
Configure Design Method a FIR y seleccione el método Equirriple, para el sistema la frecuencia
de muestreo es Fs = 40 MS/s.
Bajo Filter order, seleccione Minimum order e ingrese los siguientes parámetros:
Fstop1: 6e6
Fpass1: 7e6
Fstop1: 9e6
Fpass1: 10e6
Astop1: 80
Apass: 1
Astop2: 80
y = filter(Hd,r);
pwelch(r,[],[],[],1/Ts,'onesided');
-70
-80
-90
-100
-110
-120
-130
-140
-150
Portadoras a 8 y 12MHz
Canal Filtrado fc =8 MHz
-160
0 2 4 6 8 10 12 14 16 18 20
Frecuencia (MHz)
Fig. 3.14. Canal a 8 MHz filtrado paso banda.
De clic en Design Filter, de esa forma se genera el diseño específico, se obtienen los coeficientes
y puede hacer una revisión y análisis de las respuestas en magnitud y fase, la respuesta al
impulso, retardo de grupo, etc.
Referencias
[1] Steven, Xie. “Practical Filter Design Challenges and Considerations for Precision ADCs” en
Analog Dialogue 2016. Documento recuperado de Internet el 29 de abril de 2018, en
http://www.analog.com/en/analog-dialogue/articles/practical-filter-design-precision-
adcs.html.
[2] Filter Designer with DSP System Toolbox Software Documento recuperado del internet el 29
de abril de 2018, en https://la.mathworks.com/help/dsp/ug/use-fdatool-with-dsp-system-
toolbox-software.html y https://la.mathworks.com/help/signal/ug/filtering-data-with-signal-
processing-toolbox.html
[3] Otros enlaces de interés para aplicaciones prácticas pueden ser encontrados en
https://www.electronicshub.org/butterworth-filter/ y https://www.electronics-
tutorials.ws/filter/filter_8.html.
Introducción
Objetivos
Metodología
Se propone un análisis matemático de los diferentes diagramas de bloque para desarrollar los
módulos que modelan la arquitectura del sistema de comunicaciones. Por tanto, se hace una
Las técnicas de modulación en amplitud [1] también se conocen como de envolvente variable y
pueden sintetizarse en:
AM con dos bandas laterales y vestigio o piloto de portadora (Vestigial Sideband - VSB)
Banda lateral única superior o inferior (Upper Sideband – USB o lower sideband - LSB)
Por otra parte, las técnicas de modulación analógica se clasifican de acuerdo a la figura 4.2.a.
1. El tamaño de la antena se reduce cuando una señal es modulada por una señal
portadora de frecuencia mayor
2. Se usa modulación para transmitir la señal a través del espacio a largas distancias. Por
tanto, las técnicas de comunicaciones inalámbricas han dado origen a nuevos
estándares que evolucionan constantemente por mayor ancho de banda, tasa de datos
y cobertura.
Como ejemplo, la figura 4.3.a muestra el esquema de modulación en amplitud con la señal
mensaje o de información m(t)=sen(2 t), amplitud de la portadora Ac=1, índice de modulación
=0.9, y frecuencia de la señal portadora fc=10Hz.
En color rojo se grafica la señal de información, luego en color morado la señal portadora con
una frecuencia lógicamente más alta que la señal mensaje, en color café la señal modulada y
finalmente en color azul, la señal demodulada y recibida.
0.5
Amplitud
0
-0.5
-1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Señal Modulada
5
Amplitud
-5
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Mensaje Recibido
0.5
Amplitud
-0.5
1
m(t ) cos 2 f c t M ( f f c ) M ( f f c )
2
-50
-100
-150
-200
-250
-300
-350
-500 -400 -300 -200 -100 0 100 200 300 400 500
Frecuencia
Como función portadora, puede usarse una función seno o coseno, no hay mayor diferencia,
excepto un corrimiento de fase de -90°.
Observación:
Cualquier señal se suma por un valor constante significa que esta señal se eleva por el mismo
valor constante con respecto al eje vertical en el dominio del tiempo. En el dominio de
frecuencia, el valor constante está representado por un impulso en f=0 Hz.
Demodulación
La señal de alta frecuencia es removida luego del filtro paso bajo (LPF).
Demodulador Síncrono
Entonces el filtro paso bajo remueve las componentes de alta frecuencia (en este caso las
componentes a 2fc), y se puede recuperar el mensaje m(t).
Parámetros:
Haga doble clic en el generador de señal, y luego configure la frecuencia como 1 kHz con
una forma de onda de seno.
Ajuste la frecuencia de la onda sinusoidal de la portadora en 20 kHz.
Configure el tiempo de simulación como 0.01 para observar las señales claramente.
Ejecuta tu simulación.
Para observar el analizador de espectro, aumente el tiempo de simulación a 1 o 2
segundos.
En el filtro paso bajo Butterworth escoja la frecuencia del borde de parada como 2X, y conecte
los modelos de Modulador y demodulador para probar el sistema de transmisión de doble
banda lateral (DSB-AM), que debe quedar como se muestra en la figura 4.10.a.
En la figura se modula una señal a una frecuencia de 1 KHz, señal entregada por un generador y
que se traslada a una frecuencia de 20 KHz, que en el espacio de trabajo puede observarse
usando el comando pwelch y seleccionando la propiedad para tener espectro de doble lado
usando “twosided”. En la figura 4.12.a, se demodula la señal y se trasladada a bandabase
recuperando el mensaje en la frecuencia de 1 kHz.
Transmitir y recibir un archivo multimedia usando DSB-AM a través de una tarjeta Periférico de
Software Radio Universal (USRP) [3].
Ahora, se avanza un paso más para transmitir el archivo de música y luego se recibirá a través
de una tarjeta USRP. En este caso, la transmisión es en tiempo real, por lo tanto, a diferencia de
las simulaciones, observará el ruido a través del aire. El modelo se expresa como se muestra a
continuación. En la figura 4.17.a, se muestra el transmisor que utiliza el modulador en Simulink
como entrada a la USRP.
En la siguiente etapa, la señal pasabanda recibida sc(t) luego del proceso de demodulación es
trasladada a banda base.
Fig. 4.23.a. Modulación (señal trasladada a RF), y demodulación (señal trasladada a banda base).
5. Tareas de Laboratorio
1. Detector síncrono
Construya el modelo que se muestra a continuación tomado de [3], y luego configure los
parámetros del bloque como sigue:
a. Observe los 3 analizadores de espectro, luego explique las formas de onda desde el punto
de vista de la frecuencia (Sugerencia: recuerda la propiedad de modulación). Comenta tu
resultado.
b. Cambia el tiempo de simulación a 500 / 100k (para ver claramente la onda sinusoidal).
¿Comparar las señales en el alcance del tiempo? ¿Recuperaste m (t)? ¿Hay alguna demora
entre dos señales? Si es así, explica, ¿por qué?
Conecte el canal AWGN. Establezca la varianza de la máscara como 0.01, 0.05, 0.1 y 0.5
respectivamente. ¿Qué observas en cada caso? Comenta tu resultado.
4. Establezca el índice de modulación μ como -10, -5, -0.9, -0.1, 0.5, 0.9, 5, y 10,
respectivamente.
a. ¿Qué sucede con la forma de onda de la señal modulada para cada caso?
6. Pasos: USRP
6. Autoestudio:
>> dspenvdet
Este código abrirá el modelo Simulink del modulador DSB-AM y las técnicas de demodulación
basadas en detector de envolvente [4] con potencia cuadrática y transformador de Hilbert.
Pasos a realizar:
2. Con la ayuda del generador de funciones obtenga una señal senoidal de 1 Vpp con una
frecuencia de 700 kHz, esta señal será utilizada como portadora. Con otro generador,
obtenga una señal senoidal de 1 Vpp y frecuencia 5 kHz, la cual representa la señal
mensaje dentro del circuito.
Nota: La señal portadora y la señal mensaje pueden ser obtenidas de circuitos externos
(osciladores y señales de audio), tenga en cuenta los acoples necesarios para poder
conectar dichos circuitos al modulador AM.
3. Una vez conectado y configurado el modulador AM, visualice la señal modulada con la
ayuda del osciloscopio, obtendrá una señal como la que se muestra en la figura 4.25.a.,
donde se aprecia la envolvente:
Referencias
[1] Carlson, Bruce, and Paul B. Crilly. “Communication Systems”, Ed. Mc Graw Hill 5th Ed. New
York 2009.
[2] Silage, Denis. “Digital Communication systems using Matlab and Simulink”, Ed. Bookstand
Publishing, Gilroy, USA - 2009.
[3] Simulink model of Perfect Modulation and Demodulation, Software Defined Radio using
MATLAB & Simulink and the RTL-SDR, Strathclyde Academic Media, 2015.
[5] Li Tan and Jean Jiang. “Digital Signal Processing Fundamentals and Applications” Ed.
Academic Press, 2nd Edition. London 2013.
Introducción
La forma práctica de implementar un demodulador es usando un circuito con diodo y filtro paso
bajo RC, el diodo rectifica los picos negativos, y el filtro mantiene el nivel de voltaje, que para el
caso corresponde a la envolvente. La constante de tiempo del filtro RC debe cumplir con un
compromiso, no debe ser demasiado baja, dado que no conseguiría filtrar las componentes de
alta frecuencia (se verían los picos de la portadora), ni demasiado alta, ya que no sería capaz de
seguir las variaciones rápidas de la señal demodulada.
Objetivos
Metodología
Un filtro en cuadratura es una red pasa todo que desplaza la fase de componentes de
frecuencia positiva por -90° y componentes de frecuencia negativa por +90° [1]. Dado que ese
desplazamiento de fase de 90° es equivalente a multiplicar por ej90°=j, la función de
transferencia puede ser escrita en términos de la función signo.
Por tanto, estando x(t) una seña arbitraria a la entrada del filtro en cuadratura, la señal de
salida de x(t) = x(t) hQ(t), será definida como transformador de Hilbert de x(t), denotado por
xˆ (t ) . Es importante observar que la transformación Hilbert es una convolución que no cambia
el domino, de forma que tanto x(t) como xˆ(t ) son funciones del tiempo.
1 x d
1
xˆ (t ) x(t ) * t (4)
t t
Por otra parte, el espectro de xˆ (t ) puede obtenerse como:
F xˆ (t ) j sign f X ( f ) (5)
Puesto que un desplazamiento de fase produce el espectro de salida HQ( f )X( f ).
1. Observe que al modular una señal real x(t) por una exponencial compleja, se requiere la
multiplicación por las funciones coseno y seno. Esto produce un canal real (cos) y uno
imaginario (sin). La demodulación implica multiplicar por una exponencial compleja que
requerirá cuatro múltiplos reales y dos sumas. La parte imaginaria del resultado debería ser
cero. Por consiguiente, necesita calcular la parte real y dejar que w(t) sea la parte real de la
multiplicación compleja. Sin embargo, calcule la parte imaginaria también y muestre que es
cero o un valor muy pequeño.
2. Use los bloques To Workspace para exportar las señales y(t) y w(t) al espacio de trabajo.
3. Genere las gráficas espectrales de las tres señales x(t), y(t), y w(t), y obtenga conclusiones al
respecto.
1. Adicione los bloques como se muestra en la figura 4.6.b. y ejecute la simulación para ver la
operación del detector de envolvente.
2. Los parámetros para el bloque de Función de Matlab son como se muestra en la figura
4.5.b.
1.5
0.5
9. La figura inferior es una versión magnificada alrededor de la componente DC. Puede ver
que la salida del detector de envolvente z(t)-A difiere de la señal moduladora x(t) por
algunas componentes de alta frecuencia y una componente en DC. Esas gráficas espectrales
dan idea de cómo diseñar el filtro paso bajo y remover las componentes de alta frecuencia
en z(t)-A. Para ilustración, la figura 8 muestra x(t) y z(t)-A. ¿Cuál señal es cuál?
10. Complete el demodulador asíncrono diseñando el filtro paso bajo y filtrando z(t)-A para
producir w(t). La siguiente figura muestra x(t) y w(t). ¿Identifique las señales? ¿Es w(t) lo
suficientemente cercano a x(t) para ser aceptable?
11. Imprima las gráficas de las señales x(t), y(t), z(t), y w(t). También imprima gráficas del
espectro de x(t) y z(t)-A. Dibuje la respuesta de magnitud del filtro paso bajo aplicado a z(t)-
A para producir w(t). Explique cómo seleccionó los parámetros del filtro (frecuencia de
corte, ganancia, etc.).
3. Transformada Hilbert
N = 15;
h = (2/pi)*sin(pi*[-N:N]/2).^2./[-N:-1,1,1:N];
stem(h)
plot(h,N)
En la figura 4.9.b. se muestra la respuesta al impulso del transformador Hilbert, como se observa
corresponde a una señal con simetría impar, es decir, simetría respecto al origen.
0.6
0.4
0.2
-0.2
-0.4
-0.6
-0.8
0 5 10 15 20 25 30
Tiempo [segundos]
Si intenta reproducir las gráficas de espectro usando la función FFT en Matlab, la gráfica de fase
resultante probablemente parecerá a una línea recta con pendiente hacia abajo. La pendiente
descendente es debida al término de fase lineal introducido en el dominio de la frecuencia debido al
retardo en tiempo de la respuesta al impulso en el dominio del tiempo. No se preocupe por reproducir
esas gráficas exactamente. Usted no ha aprendido adecuadamente cómo usar la función FFT aún.
Una manera más fácil de diseñar un transformador Hilbert es usar la función remez de Matlab y su
resultado se observa en la figura 4.11.b.
load x
2. Esto cargará una variable llamada x en el espacio de trabajo de Matlab. Filtre x con su
transformador Hilbert escribiendo el siguiente comando en el prompt de Matlab.
xp=filter(h,1,x);
L = 31;
delayfilter = zeros(L,1).';
delayfilter(16) = 1;
x = filter(delayfilter,1,x);
3. A continuación, en la figura 12 se grafica las partes real e imaginaria del espectro de xp, en lugar de
tan solo la magnitud. Puede usar la función pwelch o implementar un código que permita obtener
la magnitud y densidad espectral de potencia de la señal.
Como se observa para la componente real, el espectro de magnitud es simétrico respecto al eje y,
mientras que la componente imaginaria es simétrica respecto al origen. A continuación, se observa en la
figura 4.13.b la respuesta al impulso del transformador Hilbert utilizando la función remez.
4. Usando una función, grafique las partes real e imaginaria del espectro de x y xp. Escriba una
descripción de los dos espectros. ¿Se parece a lo que estaba esperando? ¿El espectro de xp parece
al espectro de x multiplicado por la respuesta de un transformador Hilbert? ¿Es la parte real una
función par y la parte imaginaria una función impar? En las figuras 4.12.b y 4.13.b se muestra lo que
debería obtener.
Ahora que se tiene un poco de experiencia en la transformada de Hilbert se está listo para entrar en la
modulación de banda lateral única.
Un modulador de banda lateral única (sin portadora) en diagramas de bloque se muestra en la figura
4.15.b. El bloque rotulado con H corresponde al transformador de Hilbert. El bloque rotulado con D
(Delay) es un bloque de retardo puro. La señal x(t) experimenta un retardo mientras pasa por el
transformador Hilbert. Un correspondiente retardo debe ser insertado en la rama superior si está
usando un transformador Hilbert con L=31 derivadores, por tanto, debe usar un retardo de 15
muestras.
1. Construya un modelo en Simulink del diagrama de la figura 4.16.b. Para la entrada, use la misma
señal que utilizó para probar el transformador Hilbert anterior. Escriba "load x" en Matlab para
refrescar la variable x. Para la frecuencia de portadora use 0.25 y un tiempo de muestreo de 1. En
la caja de diálogo de los parámetros de Simulación, use Start Time y Stop Time configurados en 0 y
65535. El modelo Simulink debe parecerse al mostrado en la figura 16. Observe que este diagrama
Fig. 4.16.b. Modelo Simulink del Modulador de banda lateral única usando Transformador Hilbert
3. ¿Cuál es la tasa de muestreo más baja posible que podría evitar la interferencia intersímbolo (ISI)
debida a aliasing en todos los puntos de entrada en el proceso de la modulación?
4. En el modelo Simulink, inserte un bloque de exportar señal al espacio de trabajo para capturar las
señales y0(t), y1(t) y z(t). Ejecute la simulación.
Nota: Cargue la variable x nuevamente antes de ejecutar la simulación.
5. Dibuje la magnitud (lineal no escala dB) y fase del espectro de las señales y0(t), y1(t) y z(t). Escriba
una descripción que explique lo que observa.
¿Es z(t) una señal de banda lateral única? ¿Qué banda lateral es obtenida? ¿Son el espectro de y0(t)
y y1(t) lo que esperaba?
Fig. 4.18.b. Espectro de señal original y señal remodulada - Banda lateral única
Referencias
[1] Carlson, Bruce, and Paul B. Crilly. “Communication Systems”, Ed. Mc Graw Hill 5th Ed. New York 2009.
[2] The Mathworks Inc. “Analytic Signal and Hilbert Transform” (Enlace recuperado el 28 de abril de
https://la.mathworks.com/help/signal/ug/analytic-signal-and-hilbert-transform.html)
Introducción
Objetivos
1. Fundamentos teóricos
La modulación de amplitud fue el primer tipo de modulación que se consideró en los sistemas
de comunicación analógicos. La modulación de amplitud tiene la ventaja obvia de ser simple y
relativamente eficiente en cuanto a ancho de banda. Las desventajas de la modulación de
amplitud son [1]:
Cuando el mensaje pasa por un período de silencio en los sistemas de doble banda
lateral (Double Sideband - DSB) o banda lateral única (Single Sideband - SSB), se
transmiten señales de portadora muy pequeñas. La ausencia de la señal tiende a
acentuar el ruido.
Modulación de ángulo
Modulación de frecuencia
i (t ) 2 fC t (t )
1 d 1 d (t )
fi (t ) 2 fC t (t ) fC (8)
2 dt 2 dt
d (t )
Dónde es llamada desviación de fase.
dt
La desviación de fase de la portadora i(t) está relacionado con el mensaje de banda base m(t),
entonces
d (t )
K f m(t )
dt
t
(t ) K f m( )d (9)
Finalmente, la señal modulada en frecuencia se expresa como dominio de tiempo:
t
s (t ) AC cos 2 fC t K f m( )d (10)
La regla de Carson se usa para determinar el ancho de banda de la onda FM. De acuerdo con la
regla de Carson, el ancho de banda viene dado por:
BT 2W (1 f )
Kf
f max m(t )
2 W
Demodulación de frecuencia
ds (t ) t
AC 2 f C k f m(t ) sen 2 f C t 2 K f m( ) d (11)
dt
La señal diferenciada se modula tanto en amplitud como en frecuencia, la envolvente
AC [2 fC k f m(t )] está relacionada linealmente con la señal mensaje (componente de
t
amplitud) y sen[2 fC t 2 K f m( )d ] es un componente de alta frecuencia. Por lo tanto,
b. Demodulación FM no coherente
La señal FM compleja recibida tiene componentes tanto reales como imaginarios. Esta señal
tiene la forma:
j2 fC t F M ( t )
s(t ) AC e (12)
La señal compleja recibida es la entrada a dos bloques. El bloque Desplazador de Fase, que se
usa para tomar el conjugado de la señal para cambiar la fase, y el bloque Retardo, z-1, que
agrega un retardo a la señal.
j2 fC t F M ( t )
sconj (n) AC e
(13)
sdelay (t ) AC e
j 2 fC ( t ) F M ( t )
S s pd (2 fC t FM (t )) (2 fC (t ) FM (t ))
(15)
2 fC t 2 fC (t ) FM (t ) FM (t )
Τ se supone que es un valor muy pequeño, se puede considerar un operador diferenciador
d d
sd s pd (2 fC t ) FM (t ) 2 fC t 2 K f m(t ) (16)
dt dt
Usted puede reconocer fácilmente que este resultado es muy similar a la señal mediante la
demodulación utilizando el método del diferenciador.
sd sip n sqp n sip n 1 sqp n 1 (17)
Un PLL es un circuito capaz de generar una oscilación cuya fase con respecto a una señal de
entrada se mantiene acotada, contando para ello con una realimentación que compara la fase
de las dos señales y actúa modificando la frecuencia de la oscilación generada. Por su parte, el
VCO es un dispositivo electrónico que usa amplificación, realimentación y circuitos resonantes
que entrega a su salida una señal eléctrica de frecuencia proporcional a la tensión de entrada.
Entonces el PLL incorporado en el circuito para demodular la señal FM hace uso de la fuerza de
retroalimentación de un Oscilador Controlado por Voltaje (VCO) para permanecer en fase con la
portadora de la señal entrante. El mensaje se recupera como la entrada de control del VCO [4].
Modelo 1
Como se observa se mantiene un desplazamiento de fase de 90°, que se calcula mediante una
función trigonométrica y el tiempo de muestra corresponde a 1.0204e-5 segundos por
condiciones propias de la señal generada. De igual forma, no se define ninguna condición
inicial, ni límites de saturación inferior o superior. Finalmente, el método de integración
utilizado corresponde a Euler directo (Forward Euler).
Ajustando el tiempo de simulación a 0.2 segundos se puede observar las formas de onda con la
señal original y la señal modulada FM. Durante el semiciclo positivo la frecuencia de la señal es
mayor, y decrece en el semiciclo negativo.
Debe modificar los parámetros del filtro Butterworth paso bajo y del VCO como se muestra en
las figuras de configuración de parámetros 5.8, en particular el filtro tiene un orden de 10 y la
frecuencia la banda de paso se define para una frecuencia de 10 Hz. Para el oscilador se define
una sensitividad de entrada de 50 Hz/V y una fase inicial de 0 Hz.
Configurando el tiempo de simulación en 0.3 segundos, se obtienen las gráficas para la señal
sinusoidal original (onda sinusoidal en color azul), la señal modulada en frecuencia, y el
UPTC Ingeniería Electrónica Bibliografía 97
estimado de la señal demodulada (onda sinusoidal en color rojo), donde se puede apreciar en la
figura 6.8 que la sinusoide es recuperada en el receptor. Como se observa, hay un retardo entre
la señal original transmitida y la onda demodulada, aunque con propósitos didácticos como se
pretende mostrar debe prestar atención pues generalmente el sincronismo es una tarea
complicada y ardua, especialmente en el dominio digital de la señal.
Modelo 2
En el modelo 1, se aprendió los fundamentos teóricos de FM. En este segundo modelo, en lugar
de usar una estructura compleja de modulador/demodulador, se implementará un sistema FM
usando directamente los bloques de modulador/demodulador definidos en Simulink. La
entrada, es un bloque de onda seno discreto (luego deberá cambiarlo por un generador para
probar tres diferentes formas: onda senoidal, tren de pulsos rectangular y ondas triangulares,
así será capaz de observar las variaciones de frecuencia usando una gran variedad de entradas).
El modelo 2 es expresado como:
Nota para tarea: El bloque generador de señal es simplemente una entrada análoga. Para usar
este bloque como una entrada del modulador FM, se necesita digitalizarla. El bloque de tasa de
transición (retenedor de orden cero - ZOH) muestreará la información análoga muestreando el
periodo, en este caso equivale a Ts=1/100e3.
Tome la frecuencia de la señal de entrada como 5 Hz (use una onda seno), y configure la
desviación de frecuencia del modulador y demodulador como 100 Hz, respectivamente.
También, ajuste el osciloscopio con tres capas verticales para analizar m(t) vs. la señal
modulada m(t).
Sin embargo, los encontrará útiles cuando comprenda la tasa de muestreo, tasa de conversión,
respuesta al impulso finita (FIR), decimación, interpolación, etc. La salida resultante del
osciloscopio será:
1
-1
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
Tiempo (s)
-1
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
Tiempo (s)
-1
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
Tiempo (s)
Fig. 5.11. Formas de onda de en sistema de FM usando bloques Simulink.
5. Instrucciones de Prelaboratorio
6. Tareas de Laboratorio
1. Realice las siguientes tareas: Construya el modelo 1 (figura 6.1). Explique la parte de la
modulación como en la figura 6.1, describiendo cada paso.
b. Nos gustaría tener un tren de pulsos periódicos en ese caso, por tanto, debe
modificar la fuente cambiando su magnitud y desplazando el eje vertical. Los
siguientes pasos deben ser realizados desde el inicio:
4. Pasos: USRP
Referencias
[2] Hwei P. Hsu, Schaum’s Outlines of Theory and Problems of Signals and Systems, pp.1-5,
1995, McGraw-Hill.
[3] Robert W. Stewart, Kenneth W. Barlee, Dale S.W. Atkinson, and Louise H. Crockett,
Software Defined Radio using MATLAB & Simulink and the RTL-SDR, pp. 355-358 and
figure 9.17, Strathclyde Academic Media, 2015.
[6] Robert W. Stewart, Kenneth W. Barlee, Dale S.W. Atkinson, and Louise H. Crockett,
Software Defined Radio using MATLAB & Simulink and the RTL-SDR, pp. 370-375,
Strathclyde Academic Media, 2015.
[9] Robert W. Stewart, Kenneth W. Barlee, Dale S.W. Atkinson, and Louise H. Crockett,
Software Defined Radio using MATLAB & Simulink and the RTL-SDR, pp. 371, Strathclyde
Academic Media, 2015
[1] Carlson, Bruce, and Paul B. Crilly. “Communication Systems”, Ed. Mc Graw Hill 5th Ed. New
York 2009.
[2] Filter Designer with DSP System Toolbox Software Documento recuperado del internet el 29
de abril de 2018, en https://la.mathworks.com/help/dsp/ug/use-fdatool-with-dsp-system-
toolbox-software.html y https://la.mathworks.com/help/signal/ug/filtering-data-with-signal-
processing-toolbox.html
[3] Haykin, Simon and Barry Van Veen. “Señales y Sistemas”. 1a. ed., 1a. Reimpreso en México.
Limusa Wiley, 2006.
[4] Hwei P. Hsu, Schaum’s Outlines of Theory and Problems of Signals and Systems, pp.1-5,
1995, McGraw-Hill.
[5] Li Tan and Jean Jiang. “Digital Signal Processing Fundamentals and Applications” Ed.
Academic Press, 2nd Edition. London 2013.
[11] Oppenheim, A. V., Willsky, A. S., & Young, I. T.. Signals and systems. Englewood Cliffs, N.J.
Prentice-Hall, 1983.
[13] Robert W. Stewart, Kenneth W. Barlee, Dale S.W. Atkinson, and Louise H. Crockett,
Software Defined Radio using MATLAB & Simulink and the RTL-SDR, pp. 355-358 and figure
9.17, Strathclyde Academic Media, 2015.
[14] Silage, Denis. “Digital Communication systems using Matlab and Simulink”, Ed. Bookstand
Publishing, Gilroy, USA - 2009.
[15] Simulink model of Perfect Modulation and Demodulation, Software Defined Radio using
MATLAB & Simulink and the RTL-SDR, Strathclyde Academic Media, 2015.
[16] Steven, Xie. “Practical Filter Design Challenges and Considerations for Precision ADCs” en
Analog Dialogue 2016. Documento recuperado de Internet el 29 de abril de 2018, en
http://www.analog.com/en/analog-dialogue/articles/practical-filter-design-precision-
adcs.html.
[17] The Mathworks Inc® Total Harmonic Distortion (Enlace recuperado el 28 de abril de
https://la.mathworks.com/help/signal/ref/thd.html)
[19] The Mathworks Inc® Signal to Noise Ratio. (Enlace recuperado el 28 de abril de
https://la.mathworks.com/help/signal/ref/snr.html)