Repositorio de tesis abierto de la PUCP

La Pontificia Universidad Católica de Peru (PUCP) es una de las 20 mejores universidades de Latinoamérica. Recientemente me entero que tienen un enorme repositorio de tesis de acceso público que van desde nivel bachillerato hasta doctorado en un amplio número de áreas como Ciencias Sociales, Filosofía e Ingeniería y Ciencias Exactas. Me ha ayudado bastante a informarme sobre el estado del arte y para citar trabajos en mi tesis de licenciatura. Definitivamente un gran aporte al conocimiento para el mundo hispanohablante.

http://tesis.pucp.edu.pe/repositorio/

Small Ku Band Radio Telescope (1m)

Early this year I started to build a small radio telescope for my undergraduate thesis.The idea came to me as a smaller and open (without Radio Eyes) version of the Robert Arrowsmith's SDR radio telescope. In order to make a contribution to other students and educators all over the world, the main goals of the project are the following:

• Create an open hardware and affordable altazimuth mount  
• Create Python modules for control and coordinate conversions
• Start a complete open hardware/software astronomical platform for educational proposes in long term.  

This week I've already finished the mount and I began to test the movement of the antenna. I'm using a Raspberry Pi as main control unit. The digital controllers for the motors are implemented in a microcontroller board called Miuva with a Pololu Qik 2s12v10 dual motor driver:

 And this is how the system looks like when is assembled:

Once the mount is calibrated after add a magnetometer and an IMU, I will take some data using a FunCube and a RSA306B USB Spectrum Analyzer (which is not mine, sadly, but my adviser is letting me use it).

This post is only a brief teaser, but I promise I will release everything as soon as I get my degree!

Morelia Vol. 2

Yo aparentando saber de lo que hablo.

A principios de agosto regresé a Morelia para presentar lo que en esos días era mi propuesta de proyecto de tesis (un radiotelescopio de un plato de un metro del que escribiré pronto). Lo que más me tenía preocupado era tener algún episodio de pánico escénico de los cuales había creído erradicados hasta un pequeño insistente a mediados de junio, una semana después de la escuela de verano. Me tomé muy enserio el no repetirlo y estuve ensayando bastante mi corta presentación de 10 minutos. Me da mucho gusto contar que logré hacerlo bien esta vez. Toda la semana del pequeño congreso en Morelia recibí muy buenas criticas, incluida la de un investigador de del IRyA, Stan, que admiro bastante y a quien había invitado a mi charla. Disfrute mucho este viaje. Me reencontré con varios amigos e hice algunos nuevos más. Fue muy curioso que uno de mis compañeros de cuarto fue el amigo de N. de la UNAM al que había conocido un mes antes en CDMX. Conocí a una chica norteña que era asistente de verano de Stan por el programa Delfín y esta trabajando con SRT del Haystak. Me dijo que le gustaba mucho la instrumentación astronómica, me pidió algunos consejos pero después de darme cuenta lo mucho que sabía no fue tanto el extra que pude recomendarle estudiar. Hubo un día que perdí el autobús que nos llevaba del instituto al hotel por ayudarle con un problema en la comunicación serial con el analizador de espectros que estaban utilizando pero fue muy amable en ayudarme a regresar a Morelia. Ese día en el camino de regreso le contaba que a mi jamás me habían aceptado en el programa Delfín por no tener el promedio suficiente, o al menos que mi universidad no había respondido a mis solicitud de excepción. Le pareció gracioso e injusto porque yo le había dado la impresión de conocer bastante. El último día tuve que cambiar los planes de pasar unos días en Guanajuato para visitar a un amigo por no encontrar corridas de autobuses para allá después de las 4 pm y tuve que ir a la CDMX con un amigo en coche. Llegamos muy tarde a la ciudad por perdernos en Santa Fe. El carro en el íbamos de una amiga de mi amigo y nos dejaron en Taxqueña (nos habían dejado en una parada del trolebús pero nunca pasó y tuvo que regresar a rescatarnos). Era tan tarde que no alcanzamos boletos para el último autobús para Puebla, así que nos arriesgamos a cruzar media ciudad en metro para llegar a la TAPO. Al final lo logramos y llegamos en la madrugada. A mi amigo lo había conocido en la escuela de verano de junio. Me la pasé hablandole de lo mucho que sigo pensando en N.

Identificación de la función de transferencia de un motor de Lego Mindstorms NXT

Utilizando la App de System Identification Toolbox de Matlab es muy sencillo estimar experimentalmente la función de transferencia de un sistema dinámico. Básicamente, existen 3 formas de identificar la planta de un sistema dinámico: el modelo de caja negra, el modelo de caja blanca y el modelo de caja gris. En el modelo de caja blanca se obtiene la función de transferencia a partir de los principios físicos conocidos sobre el sistema (que es lo que suele hacerse en un primer curso de sistemas de control al obtener la planta de un sistema masa-resorte-amortiguador). Sin embargo, cuando se trata de sistemas reales esto o suele ser muy complicado. El método de modelado por caja negra se encarga de estimar la función de transferencia con herramientas estadísticas utilizando únicamente una función de entrada de prueba y la función de respuesta de la planta ante dicha función:

Para la que la identificación sea certera es un requisito que u(t) tengo un ancho de banda muy amplio. Esto le permite al algoritmo de estimación tener más información disponible para la estimación del modelo. Idealmente la función que cumple con ancho de banda infinito es la delta de Dirac, sin embargo, al ser imposible de generar físicamente se opta en la practica por usar secuencias pseudoaleatorias binarias (PRBS's). Estas secuencias pueden generarse con un registro de corrimiento de retroalimentación lineal (LFSR) y tienen características espectrales de ruido blanco lo que las hace ideales para ser usadas como señales de prueba. La descripción de los generadores de PRBS's es todo un tema aparte del que espero escribir a detalle más adelante, pero pueden darse muy buena idea de como funcinan con este documento.

Para esta práctica se requiere tener instalar el paquete de soporte de  Lego Mindstorms para Simulink. Pueden instalarlo ejecutando el comando supportPackageInstaller en la ventana de comandos de Matlab y buscando Lego Mindstorms NXT. El modelo implementado en el NXT es el siguiente:

 








Este diagrama nos permite adquirir los datos tanto de la función de entrada como la de salida y los envía al workspace. Para poder hacer la adquisición a la PC debe ejecutarse en modo external. [Aquí pueden ver un video tutorial para hacer esto]. Los datos obtenidos son los siguentes:

Estas dos señales son las que utilizaremos en al app para realizar la identificación. Esta aplicación tiene una interfaz gráfica muy sencilla de utilizar.[Aquí hay un videotutorial en español para aprender a usarla]. Después de experimentar con varias combinaciones de polos y ceros, el modelo preliminar que logró un acierto del 96.12 % fue el siguiente:

Este aún es el modelo crudo encontrado por el identificador y es necesario ajustarlo. Lo primero que salta a la vista el el termino 5,434e-11 que es despreciable y eliminaremos en el modelo final. Lo que resta es hacer un ajuste de ganancia. Para ello requerimos la respuesta experimental en lazo cerrado del motor ante un escalón unitario. Ejecutamos entonces el siguiente diagrama a bloques para poder obtenerlo:

Teniendo los datos de la respuesta experimental, debemos ahora compararlo con la respuesta teórica de nuestro modelo ajustando la ganancia hasta observar resultados favorables (debe haber algún procedimiento analítico para esto, pero lo he realizado de forma empírica). El diagrama del modelo teórico, junto con el modelo de caja blanca realizado por la Universidad  Politécnica de Madrid[1] como segunda referencia, es el siguiente:

 Las respuestas fueron las siguientes:

   Con esto queda verificado que la función de transferencia final de posición angular (en grados) contra potencia (%) nos queda:

 La ventaja de este proceso es que puede repetirse para caracterizar la función de transferencia de los motores con carga como puede ser un brazo robótico o un robot móvil, algo que representaría un proceso tedioso con el enfoque de caja blanca. Adelantantome a la petición de los archivos:

Descargar Modelos

Reptile - The Church

Concepto de frecuecia en tiempo discreto

Este es un tema que si no se entiende bien al principio de un curso de procesamiento digital de señales puede acarrear muchos problemas al ir avanzando en temas más complicados. En esta entrada trataré de explicar lo más claro posible lo que es la frecuencia de una señal discreta. Primero, comencemos describiendo a una función cosenoidal en tiempo continuo con la notación del Proakis [1]:

El subindice a se utiliza para denotar que x_a(t) es una señal analógica y no confundirla con su contraparte discreta x[n]. Ahora, tomemos esta misma señal discreta y multipliquémosla por tren de pulsos unitarios con un periodo de muestreo Ts (proceso conocido como muestreo).  Gráficamente, el resulto obtenido es el siguiente:

Tenemos ahora una nueva función que depende de una variable entera n multiplicada por el periodo de muestreo Ts. Es este escalamiento por Ts lo que permite que la señal discreta x[n] coincida en el tiempo con la señal x_a(t). x[n] es una nueva señal que está definida en función de una secuencia de números enteros que representan una posición en un arreglo de datos (estas posiciones pueden tomar valores negativos si se toman con respecto a un origen arbitrario en el arreglo):

Como consecuencia del teorema de muestreo de Nyquist-Shannon, sólo es posible obtener información sobre una señal o componente a una frecuencia dada sin riesgo de ambigüedades (conflictos con sus alias) mientras su frecuencia sea menor o igual a la mitad de la frecuencia de muestreo (que es simplemente el inverso de su periodo de muestreo Fs = 1/Ts). Por lo tanto, la frecuencia relativa (también llamada frecuencia normalizada) esta limitada al intervalo:

Si se grafica un ejemplo numérico, resulta bastante intuitivo visualizar que cuando un coseno discreto tiene una frecuencia relativa de pi radianes por muestra se alcanza el periodo mínimo distinguible para una función periódica discreta. Los sinusoides en tiempo discreto muestran periodicidad en frecuencia. Esto significa que los sinusoides discretos que guarden una diferencia en frecuencia de 2*pi son idénticos. Esto puede verse más claramente en esta animación que hice en Matlab:

Naturalmente, lo que nos interesan son los parámetros físicos a los que las señales en tiempo discreto hacen referencia. El periodo de muestreo Ts es el puente clave para darle significado en tiempo continuo a estas señales. Resumiendo las transformaciones de frecuencia:

 Referencias:
Sec. 1.3 "Concepto de frecuencia en señales continuas y discretas en el tiempo",  "Tratamiento Digital de Señales", John G. Proakis
Teorema del Muestreo, L. Javier Morales Mendoza, UGto

¿Qué tan lejos están las estrellas que podemos ver a simple vista?

Hace unos años se había abierto una discusión en reddit (también en Quora)a causa de la tira de xkcd de arriba. Hoy me encontré con una discusión similar ahora por un post en Facebook. Como ni las respuestas en reddit ni Quora me convencieron, decidí encontrar la respuesta yo mismo analizando algunos datos.

Primero, descargué la base de datos HYG (Hipparcos, Yale Bright Star, and Gliese) desde aquí, que nos da la información combinada de los 3 catálogos estelares mencionados, una muestra de cerca de 120,000 estrellas. Teniendo el archivo, escribí un programa en Python usando la librería Pandas [pueden bajarlo aquí] para analizar los datos. Y esto fue lo que encontré para la todas las estrellas resueltas [que pueden distinguirse como fuentes puntuales individuales] con una magnitud menor a 6:

De acuerdo al catálogo, el 98% de las estrellas que podemos ver a simple vista en la Tierra (muy lejos de las ciudades y sin luna) se encuentran a menos de 1000 años luz. Aún más interesante, las distancias de estas estrellas tienen una media de sólo 145 años luz con una desviación estándar de 136 años luz . Considerando los porcentajes en una distribución normal:

estamos hablando de que el 95% de esas estrellas se encuentran a menos de 417 años luz. Debe notarse también de que el 2% de las estrellas que se encuentran a más de 1000 años luz en este catalogo tienen una incertidumbre considerable y todas ellas se redondean a 10,000 años luz. En conclusión, cuando uno mira una estrella en la noche, lo más probable es que esté a menos de 400 años luz de distancia.

Convolución de dos señales en Matlab

La convolución continua de dos funciones esta dada por:

De esta forma, lo que esperaríamos ver al convolucionar dos pulsos cuadrados de amplitud y ancho unitario se muestra en la siguiente animación:

 Para implementar la convolución de dos señales en Matlab se debe considerar que la función conv() [ver documentación] está realizando la operación de convolución discreta que está expresada como:

Si uno utiliza esta función esperando los resultados de una convolución continua, se lleva la sorpresa de que los valores de la convolución resultante son muchos más altos que los valores correctos. Esto se ocurre porque hace falta considerar que una convolución continua se aproxima como el producto del periodo de muestreo por la convolución discreta de las mismas señales evaluadas en tiempo discreto. Esto es:

Así, el programa en Matlab nos queda:

Ts = 0.01;
t = -2:Ts:2;
f = (t>-0.5) - (t>0.5);
g = f;

%Convolución
cnv = Ts*conv(f,g,'same');

subplot(3,1,1),plot(t,f), ylim([0 1.5]),title('f(t)'),grid on
subplot(3,1,2),plot(t,g,'r'), ylim([0 1.5]),title('g(t)'),grid on
subplot(3,1,3),plot(t,cnv,'g'), ylim([0 1.5]),title('conv(f(t),g(t))'),grid on

Referencias:
Using MATLAB with the Convolution Method, California State University Northridge

Language Barriers

I hadn't written in English here before. Now, regarding that I need to take the GRE this year, it's something that I enormously regret. Mostly what I post here are technical notes of problems I solved and that there is no documentation online about it in Spanish. Someone commented once that I should write my technical notes in English so people all over the world can use them. What I always say is that there's plenty of resources written in English while good ones written in Spanish are scarce. That's why I prefer to write in my mother language. I feel I'm contributing to the develop of my region, not only Mexico but all the Hispanic world. But this blog is not only about engineering and science stuff (things I really passionate about), although it seems to be the case. It's about my life too, the whole of it. I think writing about ourself and our lives is not narcissism. We need to know other peoples thoughts and passions in order to make a better view of the world. So I know that writing always in Spanish have trapped me in bubble too. I haven't done blogger friends outside Hispanic America. But the blogosphere is kind of dead these day anyway, isn't it? Or may be it's just in the Spanish Internet! I don't even know it for the same reason. Well, I'll give it shot.

Duloxetina

A Quay at Caen - Stanislas Lépine

Algo muy extraño que ocurrió cuando tomaba Cymbalta fue perder sólo durante esos días el gusto por las pinturas oscuras de Stanislas Lépine, Christian Dahl y Atkinson Grimshaw. No es un suceso que me haya caído en completa sorpresa, pero me parece muy curioso que el gusto en arte esté tan directamente condicionado por la química cerebral. Otro efecto interesante fue perder temporalmente la capacidad de acceder a memorias asociadas con ansiedad o tristeza. ¿Es este mismo mecanismo el mismo que me alejaba del gusto por las pinturas oscuras? Tan pronto dejé el antidepresivo, hace ya algunas semanas, mi gusto volvió. No lo dejé por orden médica ni porque no me sintiera cómodo con los resultados, sino porque me incomodaba el hecho de de tener que vivir con una tranquilidad artificial.

Atkinson Grimshaw

Debo aclarar que no padezco de depresión sino que he sufrido de un desorden de ansiedad desde que tengo memoria. Un problema que mis psicólogos no supieron identificar hasta que no pude soportarlo más y tuve que visitar a un psiquiatra. Ahora lo único que lamento es no haberme tratado este problema antes. Pienso en todo el tiempo que la ansiedad me arrebató de mi vida sin poder impedirlo a pesar de mis intentos. En el fondo sabía que en mí había problema más grande que la timidez. Había sido como una batalla en un cuarto competentemente oscuro. Tarde, si, pero finalmente me siento verdaderamente bien conmigo mismo, como no lo sentía desde niño. Sé ahora que esa presión en el pecho y esa angustia constante que me envolvía no eran parte de mí. He aprendido tanto en estas últimas semanas.