Estimación de la transferencia del canal

Estimación de la transferencia del canal

de Julian Vera Piano -
Número de respuestas: 10

Hola, estamos con lo del taller, en la parte s, que tenemos que ver que tan buena es la estimación de la transferencia del canal.

Para eso, como tenemos los taps en 1, 0.2 y 0.1, consideramos que nuestra respuesta a impulso es de la forma h[n]=δ[n]-0.2δ[n-1]-0.1δ[n-2] (con δ[n] la delta).

A partir de eso, obtuvimos una transferencia que es de la forma H(e^{jw})=1-0.2e^{-jw}-0.1e^{-2jw}.

Al compararlo lo nuestro con la transferencia estimada en Gnu, vemos que son relativamente parecidas, excepto en la zona del medio de la transferencia, donde nosotros tenemos un pico pronunciado hacia abajo que en Gnu no ocurre. Pero hacia los lados son relativamente parecidas.

No estamos seguros si vamos por buen camino por como consideramos la transferencia, en particular la manera que llegamos a la H.
Tampoco sabemos a que podría deberse ese pico que nos aparece en el medio.

Saludos y gracias,
Julián

En respuesta a Julian Vera Piano

Re: Estimación de la transferencia del canal

de Federico "Larroca" La Rocca -

Hola,

Mandá una captura de pantalla o algo así. No entiendo bien a qué te referís con "la zona del medio".

saludos

En respuesta a Federico "Larroca" La Rocca

Re: Estimación de la transferencia del canal

de Julian Vera Piano -

Te mando una imagen.
La de la izquierda es lo que hicimos nosotros en Octave, y la de la derecha es la estimada por Gnu.

A los costados son relativamente parecidas, pero el bajón que ocurre en el medio no sabemos a que puede deberse,
si es en la forma que escribimos la transferencia o donde podria estar el problema.

saludos

Adjunto Transferencia.png
En respuesta a Julian Vera Piano

Re: Estimación de la transferencia del canal

de Federico "Larroca" La Rocca -

Hola,

La verdad no los veo muy parecidos. Cómo hicieron para hacer la gráfica con octave? Veo que el eje de las x va hasta 1, y debería estar entre -pi y pi (o 0 y 2pi).

saludos

En respuesta a Federico "Larroca" La Rocca

Re: Estimación de la transferencia del canal

de Julian Vera Piano -
Hola, nosotros por la manera que tomamos la transferencia, lo hacemos en función de la frecuencia (y no del tiempo como esta en Gnu)

Primero definimos un vector "f=[0:0.00001:1]" y luego graficamos la función "H=1-0.2*e.^(-j*2*pi*f)-0.1*e.^(-j*4*pi*f)" en función del vector f con "plot(f,H)"
El vector f va entre 0 y 1, por que luego lo que aparece es la misma formada que se ve ahí periodizada.

El eje va entre 0 y 1 por que lo graficamos directamente en función de f. Si hubiéramos puesto w=2*pi*f y colocado w en los exponentes ahí hubieran quedado los ejes como mencionas.

No sabemos si esta bien esa manera de pensarlo, pero fue la idea que surgio.

Saludos
En respuesta a Julian Vera Piano

Re: Estimación de la transferencia del canal

de Federico "Larroca" La Rocca -

Hola,

Sobre el eje x en GNU Radio, como el bloque saca el vector del canal estimado, la manera que encontré de graficarla fue usando el QT Time Sink, que en el eje muestra tiempo, pero le fijo un total de puntos igual al largo del vector. En definitiva, el eje está mal, no le den corte, y responde a la manera que encontré de graficar la estimación.

Sobre lo que hicieron, no entiendo porqué usan ese cambio de variable con f que les complica el razonamiento (o al menos a mí me complica). De todas formas, por lo que entendí del cambio de variable, estarían graficando entre 0 y 2\pi, por lo que sería menos parecido todavía a lo que está mostrando.

De todas formas, cuando plotean H, deberían usar el módulo. No sé qué hace cuando le pasas complejos al Octave (me parece que te muestra solo la parte real).

saludos

En respuesta a Federico "Larroca" La Rocca

Re: Estimación de la transferencia del canal

de Julian Vera Piano -

Hola,
primero que nada, hicimos el cambio en el eje x para tener w, y que la grafica vaya desde 0 hasta 2pi.

Despues, probamos comparar graficando el modulo de la transferencia (en azul, con "plot(w,abs(H))" ), con lo que teniamos graficado originalmente (en rojo, "plot(w, H)" )
El cambio entre una y otra es mínimo.

Aun graficando el módulo tenemos el comportamiento ese en el medio que no sabemos porque es. Tenemos algo mal en la manera que consideramos la transferencia? No se me ocurre donde podría estar el problema.

Saludos

Adjunto Transferencia.png
En respuesta a Julian Vera Piano

Re: Estimación de la transferencia del canal

de Federico "Larroca" La Rocca -

Hola,

Me parece raro que cambie tan poco al tomar el módulo.

De todas formas, el diagrama en GNU Radio saca la estimación del canal entre -pi y pi, por lo que se parecen cada vez menos las gráficas. porqué no usan la fft a la respuesta al impulso y se fijan qué da?

saludos

En respuesta a Federico "Larroca" La Rocca

Re: Estimación de la transferencia del canal

de Julian Vera Piano -

Hola
Intentamos escribir la respuesta al impulso, pero no encontramos (no sabemos ni si existe) forma de escribir las deltas en octave para escribir h[n]=δ[n]-0.2δ[n-1]-0.1δ[n-2] como habíamos pensado. Alguna idea como se podría expresarla en Octave?
Supongo que luego de hacer eso, hacerle la fft a h[n] no seria muy complicado.

Saludos

En respuesta a Julian Vera Piano

Re: Estimación de la transferencia del canal

de Federico "Larroca" La Rocca -
Hola,
Se le pasa la h[n] como vector. O sea, fft([1, -0.2, -0.1, 0, ... ,0]) con varios ceros al final para que la resolución en frecuencia sea razonable.
saludos
En respuesta a Federico "Larroca" La Rocca

Re: Estimación de la transferencia del canal

de Julian Vera Piano -

Hola,
Probe hacer lo que me sugeriste, y el resultado es exactamente el mismo que obtenía por el método que usaba al principio.
No le estaria encontrando la diferencia entre los dos casos.
Saludos

Adjunto FFT.png