Source: http://docplayer.es/12635952-Modulo-de-calculo-de-la-transformada-rapida-de-fourier-para-analizador-de-espectros-en-tiempo-real-en-fpga.html
Timestamp: 2018-10-22 17:53:22+00:00

Document:
Módulo de cálculo de la Transformada Rápida de Fourier para analizador de espectros en tiempo real en FPGA - PDF
Download "Módulo de cálculo de la Transformada Rápida de Fourier para analizador de espectros en tiempo real en FPGA"
Rocío Cano Valverde
1 Módulo de cálculo de la Transformada Rápida de Fourier para analizador de espectros en tiempo real en FPGA Víctor Yelpo; Diego Costa; Carlos Sosa Páez (*) Laboratorio de Electrónica, Investigación y Servicios Facultad de Ciencias Físico, Matemáticas y Naturales / Universidad Nacional de San Luis San Luis, Argentina {vayelpo, dec, Resumen Un analizador de espectros digital en tiempo real requiere un módulo de cálculo de la Transformada Discreta de Fourier y otros bloques de procesamiento de señales y control para manejo de los conversores A/D y D/A, enventanado, cálculo del módulo, permutación, generador de disparo y circuitería para manejo de reloj, reinicio y habilitación. En el presente trabajo se hizo una implementación mediante el uso de nuevas herramientas de diseño y simulación desde el modelado, que permiten diseñar sistemas digitales con gran versatilidad describiéndolos mediante un diagrama de bloques esquemáticos prediseñados con parámetros ajustables. La depuración se realiza con rapidez porque se puede simular desde el modelo de bloques (donde se genera una tabla de señales de prueba para las entradas y salidas) hasta la etapa de posicionamiento y ruteo (usando dicha tabla y comparando los resultados). La descripción se realizó en VHDL y se implementó en FPGA obteniendo resultados satisfactorios. Palabras clave: VHDL, FPGA, modelado, DFT, analizador de espectros, tiempo real, instrumentación virtual. I. INTRODUCCIÓN. Un analizador de espectros digital en tiempo real es un sistema que realiza el cálculo de la Transformada Discreta de Fourier (DFT) de una señal digitalizada y presenta el resultado en una pantalla de cristal líquido, un osciloscopio o un monitor. La señal de entrada ingresa a un conversor A/D (ADC) para el muestreo y la cuantización. Después debe ser enventanada ya que el cálculo se realiza para un número definido de muestras de la señal, conformando tramas cuya duración, en general, no coincide con un período de la misma. Luego se transfiere al bloque central que instrumenta un algoritmo eficiente de cálculo como la Transformada Rápida de Fourier (FFT) [1] y finalmente se calcula el módulo de dicho resultado. Para obtener una correcta presentación de la salida en pantalla, es preciso realizar una permutación del resultado anterior que debe egresar del sistema a través del conversor D/A (DAC). En el presente trabajo, se optó por presentar la salida en la pantalla de un osciloscopio, para lo cual se generó una señal de disparo a fin de lograr la sincronización en el barrido de frecuencias También resulta imprescindible contar con circuitería para enlazar los bloques y hacer un correcto manejo de reloj, reinicio y habilitación de cada sección. Figura 1. Ambientes de software de diseño desde el modelado. II. HERRAMIENTAS DE DISEÑO DESDE EL MODELADO. En este trabajo se hizo una implementación mediante una herramienta de diseño desde el modelado (**) que permite trabajar en forma integrada de principio al fin del diseño, con gran versatilidad y rapidez, realizando la verificación del funcionamiento en cada etapa de la realización y en cada uno de los ambientes del software como muestra la Figura 1. Las etapas de un diseño completo se resumen a continuación. A. Modelado. Como primer paso del procedimiento, se describe el diseño con bloques esquemáticos pre-diseñados de parámetros ajustables mediante una herramienta de modelado como Matlab/Simulink [2]. Esto es mostrado en la Figura 2. En este paso, para testear el algoritmo, se puede hacer una simulación en punto flotante de doble precisión. Luego se puede simular y optimizar el diseño trabajando en punto fijo, para evaluar los efectos de longitud de palabra finita. Este procedimiento de simulación-corrección en el nivel de abstracción mayor, configura una ventaja clave respecto a otras metodologías de diseño porque permite una depuración temprana. (*) Pertenecientes al Proyecto Instrumentación Virtual Reconfigurable financiado por Ciencia y Técnica de la UNSL. (**) Las herramientas de software fueron recibidas como premio Designer Forum Mention del III Southern Conference on Programable Logic. 23
2 Figura 2. Diagrama de bloques esquemáticos en el ambiente de modelado. B. Conversión de modelo a código HDL. Luego, la herramienta de traducción del modelo de alto nivel a RTL (Register Transfer Level) genera un código en lenguaje de descripción de circuitos (HDL) a partir del diagrama en bloques. Para esta tarea se usó Synplify DSP [3]. Esta herramienta de traducción provee una librería de bloques de diseño DSP que se adiciona al ambiente de la herramienta de modelado. También genera automáticamente un código de banco de pruebas de simulación o test bench en HDL y una tabla de valores de señales de entrada y las salidas correspondientes, para ser utilizado en el siguiente paso. Una ventaja adicional, es la opción de optimización de la descripción por área, para lo cual, la herramienta verifica la aparición de estructuras repetidas y hace un re-uso de las mismas mediante un mecanismo de plegado en la etapa temprana del diseño. Esto redunda en un ahorro de celdas, agregando sólo circuitería para arbitraje. También hay una opción de optimización por velocidad que hace el trabajo inverso. C. Síntesis, posicionamiento-ruteo y programación. El código HDL generado en la traducción, finalmente se importa al ambiente de la herramienta usual de depuración, síntesis, posicionamiento-ruteo y programación. Luego se realizan los pasos habituales con sus correspondientes simulaciones (pre-síntesis, post-síntesis, y postposicionamiento-ruteo). Usando la tabla de señales del paso anterior, la herramienta de simulación compara automáticamente los resultados para detectar errores. En el presente trabajo se utilizó el Libero 8.6 para depuración, Synplify DSP AE para la síntesis, Designer para posicionamiento-ruteo, Flash Pro para programación y ModelSim para simulación en todos lo niveles de la etapa [4]. Todas estas herramientas demostraron ser aptas para trabajar integradamente brindando total portabilidad de los archivos generados en el proyecto entre una y otra etapa. Como lenguaje, se eligió el VHDL [5] [6] por su amplia difusión en el ambiente académico. Además, existe gran disponibilidad de códigos abiertos en este lenguaje, que describen IP cores de sistemas variados, los cuales pueden sintetizarse en el mismo FPGA para implementar funciones más complejas haciendo re-uso de diseños probados. III. DISPOSITIVOS UTILIZADOS. Se utilizó la placa de desarrollo RVI Prototype Board junto con la placa de expansión LP Data Conversion Daughter Board del ICTP (***). La placa principal mencionada contiene un FPGA ProASIC3E (A3PE1500) de Actel [7], conectores para incorporar placas de expansión, circuito de configuración, extensión para memoria SDRAM externa, indicadores y pulsadores para control y depuración, como así también varios puertos de comunicación (paralelo, RS232, USB y ethernet 10/100). La placa de expansión posee un conversor A/D (AD9201) y un conversor D/A (LTC1654), cada uno con conectores y selector de rango, además de sensores de temperatura. Esta plataforma fue diseñada para aplicaciones abiertas de instrumentación reconfigurable [8]. Los bloques constitutivos y el conexionado se muestran en la Figura 3. IV. DESCRIPCIÓN DEL DISEÑO. El diseño cuenta con los bloques descriptos a continuación, los cuales son mostrados en la Figura 4. (***) Esta placa ha sido provista por el ICTP (International Center For Theoretical Physics) y fue desarrollada por A. Cicuttin, M. L. Crespo y A. Shapiro. 24
3 Figura 3. Diagrama de bloques de la placa utilizada. A. Bloque generador de reloj. Toma como entrada el oscilador de 25 MHz provisto en la placa. Mediante un PLL (lazo de enganche de fase) y un divisor de frecuencia, se obtiene la señal de reloj de 20 MHz, adecuada para alimentar al bloque DAC. B. Bloque DAC. Este bloque constituye la salida del sistema. Inicializa el chip del conversor D/A y transforma el dato de entrada paralelo al formato serie del chip. También genera un pulso que es usado como reloj lento (clock out) para el resto del sistema. Cabe mencionar que todas las entradas y salidas tienen formato paralelo con excepción de las del chip del conversor D/A. Por eso, este chip es quien impone la primera limitación en la frecuencia de trabajo. Este conversor D/A disponible en la placa posee una frecuencia de reloj máxima admitida de 22,5 MHz. Con esta frecuencia, la tasa máxima de conversión es de 0,9 MHz porque el protocolo de comunicación serie del chip requiere 25 ciclos para poner un dato de 14 bits de carga útil en la salida. Tomando un margen de seguridad, se fijó la frecuencia en 20 MHz generando la señal clock out de 800 khz (20 MHz/25). C. Bloque ADC. El bloque ADC constituye la entrada del sistema. El chip del conversor A/D contenido en la placa permite una frecuencia máxima de muestreo de 20 MHz, pero debido a las limitaciones impuestas por conversor D/A mencionadas, el conversor A/D se alimenta con un reloj de 800 khz. La entrada tiene un ancho de palabra de 10 bits en paralelo. Dado que el bloque que calcula la DFT realiza internamente una serie de productos y sumas en cascada, cada bloque que lo compone es susceptible de sufrir desbordamiento que se arrastra hasta el final produciendo errores en el resultado. Por este motivo, debe restringirse el ancho de palabra de entrada a un valor apropiado para el ancho de palabra de la salida. Esto se verificó mediante simulación y prueba en la placa, fijando la salida del bloque analizador en 14 bits, se comparó una señal sinusoidal de 10 bits y otra de 8 bits de ancho de palabra. En el primer caso, la salida saturó al llegar a su valor máximo y no permitió buena discriminación entre el pico y el ruido de fondo. En el segundo caso, hubo una mejor discriminación y no saturó. Por este motivo, en el diseño se descartaron los 2 bits menos significativos, de los 10 bits de entrada. Esto representa un filtrado del contenido de baja amplitud de la señal. La conversión A/D introduce un offset a la señal de entrada para que todos los valores de tensión queden comprendidos en el rango de 0 a 1 V. Por esto, en el bloque ADC, se invierte el bit más significativo para representar los datos de la señal de salida del bloque en complemento a 2 y eliminar el offset mencionado. En el caso de 8 bits, el máximo valor positivo resultante es , que corresponde a 1 V, y el máximo negativo es , que corresponde a 0 V. D. Bloque analizador. Se compone de los siguientes sub-bloques: Puerto de entrada: Registra la salida del bloque ADC con un ancho de palabra de 8 bits en complemento a 2. Sub-bloque generador de ventana: Genera una secuencia de 1024 puntos de una función para enventanado. Esta secuencia es multiplicada con la señal de entrada y el resultado alimenta el sub-bloque de cálculo de la FFT. Se eligió la ventana de Blackman por ser la que produce menor derrame espectral. Sub-bloque de cálculo de la FFT: Implementa el algoritmo de la FFT. Presenta las partes real e imaginaria del resultado, en tramas de 1024 muestras, dadas en dos vectores de 13 bits con formato entero y con signo. Tiene una latencia de 2061 ciclos de reloj. Para indicar que hay un resultado válido a la salida, el sub-bloque cuenta con una salida ready y otra salida valid en las que se generan un pulso y un 1 lógico respectivamente. Sub-bloque de cálculo de Módulo: Calcula el valor absoluto de la salida del bloque de cálculo de la FFT, en un vector con formato entero sin signo de 14 bits. Para optimizar los recursos, se utiliza un algoritmo CORDIC (COordinate Rotation DIgital Computer) ya que requiere sólo sumas, restas, desplazamientos y búsqueda en tablas. Sub-bloque de Permutación: Deja inalteradas las primeras 512 muestras e invierte el orden de las últimas 512. De esta forma, no se descartan valores y se observa la representación del espectro sólo en el rango de las frecuencias positivas. Esto aprovecha la propiedad de que, para una secuencia x[n] real, el módulo de la DFT, X[k], tiene simetría par. Por definición, la DFT se calcula para los valores del índice k que corresponden al intervalo ω=[0, 2π-2π/N], donde los valores de frecuencia angular digital ω se corresponden con los valores de frecuencia analógica F, con F creciente para ω de 0 a π, y F decreciente para ω de π a 2π (donde ω=0 corresponde a la continua F=0, ω=π a la mitad de la frecuencia de muestreo F=F m /2 y ω=2π nuevamente a la continua F=0). La salida de este bloque tiene 14 bits en formato entero y sin signo. Posee una entrada de reset que se conecta a la salida ready del sub-bloque de cálculo de la FFT, para limpiar todos los registros antes de iniciar el procedimiento de permutación. Puerto de Salida: Registra la salida del bloque anterior. Está configurado para trabajar con datos sin signo. La salida alimenta la entrada de bloque DAC que tiene este formato. 25
4 Figura 4. Diagrama de bloques del diseño final. Sub-bloque de disparo: Genera una onda rectangular de 512 puntos de período. Esto permite que la representación del módulo de la DFT (dada en las tramas de 512 muestras que se suceden cronológicamente), se pueda visualizar en el osciloscopio como una imagen estable. Posee una salida de 1 bit que es la señal de disparo para conectar al osciloscopio, y una entrada enable de habilitación de 1 bit la cual se conecta a la salida valid del sub-bloque de cálculo de la FFT. E. Bloque generador de habilitación y reinicio local La herramienta de diseño que realiza la traducción del modelo de bloques esquemáticos a código HDL permite trabajar con dos tipos de señales de re-inicialización para llevar a un estado inicial a todos los registros que forman parte del hardware: Un reset global y un reset local. El primero es generado por la herramienta de traducción y no es visible en el ambiente de la herramienta de modelado. Puede configurarse como sincrónico o asincrónico. En este diseño se eligió trabajar en forma sincrónica por compatibilidad con el resto de los bloques cuyos códigos se generaron manualmente. El Segundo es opcional. Es visible en el ambiente de la herramienta de modelado y puede habilitarse en los subbloques de lógica secuencial del diseño de esa herramienta. Siempre es sincrónico. En el presente trabajo, se habilitó un reset local en los sub-bloques del generador de ventana, de cálculo de la FFT y el generador de disparo. La herramienta de traducción también genera una señal de entrada de habilitación que se activa por alto y no es visible en el ambiente de modelado. El bloque generador de habilitación y reinicio local se encarga de hacer una adecuada sincronización que requieren estas señales. Este bloque tiene una entrada de clock conectada al reloj de 800 khz y otra entrada de reset vinculada a un pulsador de la placa. Posee también una salida reset local que se conecta al reset local de los sub-bloques ya mencionados y la salida de habilitación enable que se conecta a la entrada enable del bloque analizador. El diseño de este bloque fue realizado directamente en VHDL describiéndolo mediante una máquina de estados. V. CONCLUSIONES. A. Prestaciones. El rango dinámico del diseño es ajustable según la configuración de puentes de la placa. Se seleccionó un rango dinámico de entrada de -0.5 a 0.5 V y de salida de 0 a 3.3 V. El ancho de banda está limitado principalmente por restricciones de los componentes que posee la placa utilizada, es decir que son factores ajenos al diseño. En primer término, está restringido por la tasa máxima de reconstrucción de 0,9 MHz del chip del conversor D/A serial que posee la placa. Tomando un margen de seguridad, se configuró una frecuencia de muestreo de 0,8 MHz dando un máximo de 0,4 MHz para la señal analizada, cumpliendo el teorema del muestreo. Para mejorar el ancho de banda, debería cambiarse el chip del conversor D/A por uno que trabaje en paralelo. Superada esta restricción, el ancho de banda estaría limitado a 10 MHz, dado que el chip del conversor A/D trabaja con una frecuencia de muestreo máxima de 20 MHz. Cambiando este conversor, la limitación estaría dada por la máxima frecuencia soportada por el bloque Analizador que según el reporte es de 188 Mhz, por lo que se tendría un ancho de banda de 94 MHz. La resolución en amplitud está dada por el ancho de palabra de 14 bits del chip del conversor D/A de la placa. La resolución interna se ajustó para obtener la mejor precisión. El análisis de una secuencia de longitud finita genera derrame espectral. La envolvente del espectro tiene forma de seno cardinal donde el ancho del lóbulo principal y la amplitud de los lóbulos laterales dependen de la longitud de la secuencia y del tipo de ventana de ponderación. En el diseño realizado, el rango de frecuencias de 400 khz dado por el ancho de banda logrado, se divide en 512 muestras. 26
5 TABLA I. REPORTE DE RECURSOS CONSUMIDOS Recurso Detalle Cantidad % Celdas lógicas Bloques E/S Entradas 11 Salidas 10 Memoria RAM TABLA II. REPORTE DE TEMPORIZACIÓN Nombre Período Frecuencia Máx Clock Min Clock (ns) (MHz) - Out (ns) Out (ns) Ext Clk UO/U1/Core: GLA U2/present_state[1]:Q U4/Clk_0:Q Figura 5. Simulación de una señal de entrada con forma de onda cuadrada, envantando de la misma y salida de las tramas con el módulo de la DFT. Para mejorar la resolución en frecuencia se usó la ventana de Blackman ya que, si bien el ancho del lóbulo principal de su DFT es mayor que el de la rectangular, la amplitud de sus lóbulos laterales es menor logrando mejor resolución entre frecuencias adyacentes. El enventanado se muestra en la Figura 5. Con estas premisas, se logró que la discriminación entre frecuencias consecutivas, sin excesiva distorsión, sea menor al 1% del ancho de banda. En este trabajo, con un ancho de banda de 0,4 MHz se logra una resolución menor a 3,5 khz. En las pruebas realizadas inicialmente, se observó la linealidad de la amplitud en el 30% del rango. Para mejorar esta prestación, se modificaron las salidas del bloque FFT para trabajar en Full Precision. Esto produjo internamente un ensanchamiento de palabra a 41 bits. Además, se introdujo un bloque de conversión que hace extrapolación inversa para llevar la salida del bloque de cálculo de módulo a formato de 14 bits, sin signo, apropiado para el conversor D/A. Con esto, se amplió el rango de linealidad al 70% de la excursión. El reporte generado por la herramienta SmartPower de Libero muestra que el consumo de potencia estático es de 18.12mW y el consumo dinámico es de mw, de los cuales el consumo debido a las líneas de conexionado es de 0,007mW y el de bloques E/S es de 0.491mW. Esto suma 18,618mW de consumo total. B. Recursos utilizados. En la Tabla I se detalla el reporte de recursos. En la Tabla II se consigna el informe de tiempos. El bloque analizador, etiquetado como U0, muestra que la limitación de ese bloque es de 188MHz. C. Cumplimiento de los objetivos. Se ha logrado implementar un analizador de espectros compacto, totalmente implementado en hardware que puede ser utilizado en cualquier laboratorio que cuente con un osciloscopio. En la Figura 6 se puede visualizar la salida mostrada en la pantalla del osciloscopio y en la Figura 7 se muestra el equipo completo. El diseño con herramienta de conversión de modelo a código HDL resulta más versátil que el procedimiento realizado con la metodología usual de descripción manual del código [9]. La integración de los pasos de diseño (partiendo del modelado, la descripción en código HDL, la síntesis, y demás pasos hasta la programación del dispositivo) permite realizar un diseño óptimo logrando el total cumplimiento de las premisas especificadas en el más alto nivel de abstracción que es el modelo inicial. El proceso de validación se acelera al simular en alto nivel durante el modelado y a bajo nivel usando la tabla de señales patrón de entrada y sus correspondientes salidas que permiten el contraste. 27
6 Figura 6. Pantalla del osciloscopio mostrando el espectro de una señal de amplitud modulada en un canal y la señal de disparo en el otro. El software utilizado para la traducción del modelo a código HDL, permite optimizar la descripción por área o por velocidad al momento de generar el código, haciendo re-uso de estructuras por un mecanismo de plegado. Es posible realizar diseños para plataformas de diferentes empresas o diseñar cores para FPGAs o CPLDs genéricos. D. Posibles implementaciones futuras. La instrumentación virtual reconfigurable es una tendencia creciente [10] [11]. Por eso se pretende realizar una conexión a computadora por USB y generar una interfaz de software para ver el espectro. Esto desafectaría la limitación en frecuencia del chip del conversor D/A. También se requiere un amplificador de entrada con bajo ruido y control calibrado de atenuación. Igualmente, sería conveniente contar con filtros anti-solapamiento y reconstructor realizado con un banco de capacitores para ajustar la frecuencia de corte por programación. Asimismo, se podrían agregar otras opciones de enventanado, ajustar la frecuencia de muestreo para regular ancho de banda y resolución o incluir una sección que implemente el desplazamiento espectral de la entrada, para analizar señales con alto contenido armónico. Por último, se podría integrar este trabajo con otros instrumentos virtuales como osciloscopios, generadores de señales, etc., para brindar un set completo de instrumentación para laboratorio. Figura 7. Placa de desarrollo funcionando con el instrumental de prueba. [8] A. Cicuttin, M. L. Crespo, A. Shapiro y N. Abdallah, A block-bassed open source approach for a reconfigurable virtual implementation platform ussing FPGA technology, IEEE International Conference on Reconfigurable Computing and FPGA s ReConFig 2006, 2006, pp 1 a 8. [9] J. F.Wakeerly, Diseño digital, principios y prácticas, Prentice-Hall, [10] C. Quintáns, M. Moure, R. García-Valladares, M. Valdés, E. Mandado, A virtual instrumentation lab based on a reconfigurable coprocessor, IEEE Trans. Instrumentation and Measurement, Apr. 2006, pp [11] H. Goldberg, What is virtual instrumentation?, IEEE Instrumentation & Measurement Magazine, vol. 3, no. 4, 2000, pp 10 a 13. VI. REFERENCIAS BIBLIOGRÁFICAS. [1] J. G. Pproakis y D. G. Manolakis, Tratamiento digital de señales, 3ra ed., Prentice-Hall, 2003, pp 457 a 507. [2] S. Nakamura, Análisis numérico y visualización gráfica con MATLAB, Prentice-Hall Hispanoamericana S.A. [3] Synplify DSP User Guide [4] Manuales de uso del Libero, Designer, Flash Pro y ModelSim. [5] K. Skahill, VHDL for programmable logic, Addison-Wesley, [6] J. Bhasker, VHDL primer, 3th edition, Prentice Hall, [7] Documentos técnicos de Actel Corporation, 28
Laboratorio 2: Mediciones Digitales
Objetivos: Laboratorio 2: Mediciones Digitales Conocer y utilizar con propiedad un osciloscopio de señal mixta. Manejar los conceptos de sincronización, disparo, nivel de disparo, y base de tiempo de un
Análisis y Síntesis. Introducción a los Sistemas Lógicos y Digitales 2008
Introducción a los Sistemas Lógicos y Digitales 2008 Métodos de análisis: Tabla de verdad. Heurístico. Diagramas de estado. Simulación. Test del hardware. etc... Ejemplo de análisis heurístico (planteo
DISEÑO E IMPLANTACIÓN DE UN FRECUENCÍMETRO UTILIZANDO EL CPLD EPM7128SLC84-15.
DISEÑO E IMPLANTACIÓN DE UN FRECUENCÍMETRO UTILIZANDO EL CPLD EPM7128SLC84-15. J. Valeriano, C. Ojeda, J. Castillo, S. Quintana Centro de Instrumentos UNAM. Laboratorio de Electrónica. Circuito Exterior

References: resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución