Esta página aún no se ha traducido para esta versión. Puede ver la versión más reciente de esta página en inglés.

czt

Sintaxis

y = czt(x,m,w,a)
y = czt(x)

Descripción

y = czt(x,m,w,a) devuelve el chirp Z-transform de la señal.x La transformación Z chirp es la transformación Z de a lo largo de un contorno espiral definido por y .  es un escalar que especifica la longitud de la transformación, es la relación entre los puntos a lo largo del contorno de espiral -plano de interés, y escalar es el punto de partida complejo en ese contorno.xwamwza El contorno, una espiral o "chirp" en el plano, es dado porz

z = a*(w.^-(0:m-1)) 

y = czt(x) utiliza los siguientes valores predeterminados:

  • =mlength(x)

  • =wexp(-j*2*pi/m)

  • =a1

Con estos valores predeterminados, devuelve la transformación Z de en puntos igualmente espaciados alrededor del círculo de la unidad.cztxm Esto es equivalente a la transformación discreta de Fourier de , o .xfft(x) La matriz vacía especifica el valor predeterminado para un parámetro.[]

Si es una matriz, transforma las columnas de .xczt(x,m,w,a)x

Ejemplos

contraer todo

Cree un vector aleatorio, , de longitud 1013.x Calcular su DFT utilizando .czt

rng default x = randn(1013,1); y = czt(x);

Se utiliza para acercar una sección de banda estrecha de la respuesta de frecuencia de un filtro.czt

Diseñe un filtro FIR de paso bajo de 30o orden utilizando el método window. Especifique una frecuencia de muestreo de 1 kHz y una frecuencia de corte de 125 Hz. Utilice una ventana rectangular. Busque la función de transferencia del filtro.

fs = 1000; d = designfilt('lowpassfir','FilterOrder',30,'CutoffFrequency',125, ...     'DesignMethod','window','Window',@rectwin,'SampleRate',fs); h = tf(d);

Calcular el DFT y el CZT del filtro. Restringir el rango de frecuencia del CZT a la banda entre 100 y 150 Hz. Genere 1024 muestras en cada caso.

m = 1024; y = fft(h,m);  f1 = 100; f2 = 150; w = exp(-j*2*pi*(f2-f1)/(m*fs)); a = exp(j*2*pi*f1/fs); z = czt(h,m,w,a);

Trazar las transformaciones. Amplíe el área de interés.

fn = (0:m-1)'/m; fy = fs*fn; fz = (f2-f1)*fn + f1;  subplot(2,1,1) plot(fy,abs(y)) axis([f1 f2 0 1.2]) title('FFT')  subplot(2,1,2) plot(fz,abs(z)) axis([f1 f2 0 1.2]) title('CZT') xlabel('Frequency (Hz)')

Diagnóstico

Si , , o no es un escalar, aparece el siguiente mensaje de error:mwaczt

Las entradas M, W y A deben ser escalares. 

Algoritmos

utiliza la siguiente longitud de potencia de 2 FFT para realizar una convolución rápida al calcular la transformación -en un contorno chirp especificado.cztz[1]

Referencias

[1] Rabiner, Lawrence R., and Bernard Gold. Theory and Application of Digital Signal Processing. Englewood Cliffs, NJ: Prentice-Hall, 1975, pp. 393–399.

Capacidades ampliadas

Consulte también

|

Introducido antes de R2006a