Transformada Z de un muestreador

12

En este artículo o filtrado multirrata, el autor establece la siguiente relación matemática. Deje ser la salida de un muestreador de modo tal queyD

yD[n]=x[Mn]

donde es el factor de disminución de resolución. En otras palabras, conservamos cada -ésima muestra de la señal original. El autor continúa diciendo lo siguiente:MMM

... la transformación z de viene dada poryD[n]

YD[z]=1Mk=0M1X[z1/MWk]

donde es el núcleo de la Transformada de Fourier discreta de punto , es decir, . M e ( - j 2 π k ) / MWkMe(j2πk)/M

¿Cómo podemos pasar de la primera expresión a la segunda? ¿Cuál es la relación entre DFT y la transformación Z que permite tal transición?

Phonon
fuente

Respuestas:

9

Esta derivación es complicada. El enfoque sugerido anteriormente tiene un defecto. Déjame demostrarte esto primero; entonces daré la solución correcta.

Queremos relacionar la de la señal , , con la de la señal original .Y D ( z ) = Z { x [ M n ] } Z X ( z ) = Z { x [ n ] }ZYD(z)=Z{x[Mn]}ZX(z)=Z{x[n]}

La forma incorrecta

Se podría pensar en simplemente conectar la expresión para la señal muestreada a la expresión de la transformación :Z

YD(z)=n=+x[Mn]zn

Un cambio de la variable parece obvio:n=Mn

YD(z)=nMZx[n]zn/M

Sin embargo, es importante darse cuenta de que aunque el nuevo índice de suma todavía se ejecuta desde a , la suma ahora supera 1 de M números enteros . En otras palabras, - n

nMZ={...,2M,M,0,M,2M,...} ,

mientras que la definición de la requiereZ

n{...,2,1,0,1,2,...} .

Como esto ya no es una , no podemos escribir:Z

YD(z)=X(z1/M)

La direccion correcta

Primero definamos una señal de tren de impulso 'auxiliar' como:tM[n]

tM[n]=k=+δ[nkM]={1:nMZ0:nMZ

Esta función es en una de cada muestras, y cero en cualquier otro lugar.m1M

De manera equivalente, la función del tren de pulsos se puede escribir como:

tM[n]=1Mk=0M1ej2πkn/M

Prueba: debemos considerar por separado los casos y : n M ZnMZnMZ

nMZ

tM[n]=1Mk=0M1ej2πkn/M={1Mk=0M11:nMZ1M1ej2πkn1ej2πkn/M:nMZ={1MM:nMZ1M111ej2πkn/M:nMZ={1:nMZ0:nMZ
En el caso ,nMZ

Ahora volvamos a nuestro problema original de encontrar la de un muestreador inferior:Z

YD(z)=n=+x[Mn]zn

Aplicamos la sustitución , teniendo en cuenta que esto hace que la suma se ejecute solo sobre múltiplos enteros de M:n=Mn

YD(z)=nMZx[n]zn/M

Ahora podemos usar la función de tren de impulsos anterior para reescribir esto de manera segura como una suma de todos :nZ

YD(z)=n=+tM[n]x[n]zn/M

Usando la formulación anterior para la función de tren de impulsos como una suma finita de exponenciales, obtenemos:

YD(z)=n=+(1Mk=0M1ej2πkn/M)x[n]zn/M=1Mk=0M1n=+ej2πkn/Mx[n]zn/M=1Mk=0M1n=+x[n](ej2πk/Mz1/M)n

La suma de la derecha es una suma de todos los enteros y, por lo tanto, es una válida de en términos de . Por lo tanto, podemos escribir:Zz=ej2πk/Mz1/M

YD(z)=1Mk=0M1X(ej2πk/Mz1/M)

Esta es la fórmula para la de un reductor de muestras.Z

usuario9037
fuente
1
Muy agradable. Mientras leía mi respuesta anterior, también noté la misma falla que tú.
Jason R
5

No he visto esta notación antes. Sin embargo, parece tener sentido. El -downsampler se define por la ecuación:M

yD[n]=x[Mn]

Su transformación está definida por la ecuación:z

YD(z)=n=yD[n]zn=n=x[Mn]zn

Aplicar un cambio de variable, dejando . Los rangos de la suma no se ven afectados por el cambio de variable, ya que se extienden hasta el infinito.n=Mn

YD(z)=n=x[n]zn/M

Esto se parece a la transformación de sí. Recuerde que se define como:zx[n]

X(z)=n=x[n]zn

Por inspección, podemos concluir la siguiente relación entre las transformaciones de y :zx[n]yD[n]

YD(z)=X(z1/M)

Por lo tanto, la transformación de la salida del muestreador inferior está estrechamente relacionada con la transformación de la señal de entrada, que es de esperar. En el dominio de la frecuencia, esto da como resultado un estiramiento pliegue del contenido de frecuencia de la señal.zzM

Pero, ¿cómo se pasa de la ecuación anterior a la referenciada en el documento? Da una definición de en términos de solamente, mientras que la expresión que derivamos es una función de . Entonces, para un valor particular de que le gustaría evaluar , primero debe calcular (es decir, tomar la raíz -ésima de ) y luego sustituirla por . Sin embargo, todos los distintos de cero tienen raíces -th distintas :YD(z)zz1/MzYD(z)z1/MMzX(z)zCMM

{rp, rpej2πM, rpej2π2M,  , rpej2π(M1)M}

={rp, rpW, rpW2,  , rpWM1}

donde es el valor del núcleo DFT referencia en su pregunta, y es lo que defino como la raíz -ésima principal del valor complejo :Wkej2πk/MrpMz

rp=|z|MejzM

Esto es, 's principal de -ésima raíz se obtiene mediante la conversión de a la forma polar, teniendo el raíz-ésimo de ' s magnitud (que es un número real), y dividiendo 's ángulo por . Los valores resultantes expresan en forma polar.zMrpzMzzMrp

¿Por qué ir a todos estos problemas? Porque, como señalé antes, la asignación del de al dominio de no es uno a uno. Ahora comenzaré a agitar las manos. Para cualquier valor particular de que le gustaría evaluar , hay puntos correspondientes en que podría asignar. Por lo tanto, cada uno de esos puntos en contribuye al valor correspondiente de . Luego terminas con una suma como la que se muestra en el documento:YD(z)X(z1/M)zYD(z)MX(z1/M)MX(z1/M)YD(z)

YD(z)=1Mk=0M1X(rp(z)Wk)

donde refiere al cálculo de la raíz -ésima principal que mostré anteriormente. En realidad, podría elegir cualquiera de las raíces de como la principal; Elegí esta definición porque es la más sencilla. Si tuviera que derivar esta relación de manera adecuada y rigurosa, creo que el factor entra en juego debido a una derivada de .rp(z)MzM1Mz1/M

En el lenguaje matemático, creo que esto se denominaría una composición de funciones; , donde y . Para desenrollar la composición de la función y escribir como una función de , debe cortar el dominio de en fragmentos que son uno a uno, invertir la función en esos intervalos y luego sumar Los resultados con factores de escala apropiados. He usado esta técnica antes para calcular la función de distribución de probabilidad de una función de una variable aleatoria dado el pdf original de la variable aleatoria (por ejemplo, para derivar el pdf de dadof ( z ) = X ( z ) g ( z ) = z 1 / M Y D ( z ) z Y D ( z ) YD(z)=f(g(z))f(z)=X(z)g(z)=z1/MYD(z)zYD(z) XXX's pdf), pero el nombre de la técnica se me escapa.

Jason R
fuente
Muy buena respuesta.
Spacey
Gracias. Cualquier matemático con licencia se estremecería ante mi intento de una descripción (obviamente soy un ingeniero). No creo que esté muy claro, pero tal vez alguien más pueda sugerir una explicación más limpia, o tal vez piense en una mejor manera de decirlo.
Jason R
Entiendo la primera mitad, pero las cosas se ponen borrosas hacia el final para mí.
Spacey
Debería reescribir la segunda mitad cuando tenga la oportunidad. Realmente es solo una técnica estándar para derivar una expresión para la composición de dos funciones. Necesito recordar los detalles de cómo hacerlo.
Jason R