Main Content

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.

Frecuencia instantánea de quirp complejo

En este ejemplo se muestra cómo calcular la frecuencia instantánea de una señal mediante la transformación sincronera de Fourier.

Generar un chirrido con contenido de frecuencia sinusoidalmente variable. La señal está incrustada en ruido gaussiano blanco y se muestrea a 3 kHz durante 1 segundo.

fs = 3000; t = 0:1/fs:1-1/fs;  x = exp(2j*pi*100*cos(2*pi*2*t)) + randn(size(t))/100;

Calcular y trazar la transformación sincronaexprimida de Fourier de la señal. Visualice la hora en el eje y la frecuencia en el eje - .xy

fsst(x,fs,'yaxis')

Encuentra la frecuencia instantánea de la señal extrayendo la cresta de frecuencia de tiempo de máxima energía de la transformación Sincroprimida de Fourier.

[sst,f,tfs] = fsst(x,fs);  fridge = tfridge(sst,f);

Superponga la cresta en la gráfica de transformación. Convierta el tiempo a milisegundos y la frecuencia a kHz.

hold on plot(t*1000,fridge/1000,'r') hold off

Para una señal real, puede encontrar la frecuencia instantánea más fácilmente utilizando la función.instfreq Por ejemplo, muestre la frecuencia instantánea de la parte real del chirp complejo calculando la señal analítica y diferenciando su fase.

ax = real(x);  instfreq(ax,fs,'Method','hilbert')

Consulte también

| | | | |

Temas relacionados