Experimento MATLAB para extracción de frecuencia

La parte que escribí puede dibujar el diagrama de señal original y el diagrama de espectro después de la transformación fft:

fs = 1000 % Establecer la frecuencia de muestreo

n = 1024;

N = 0:N-1;

t = n/(fs * 50);

SF = 100 % establece la frecuencia de la señal sinusoidal

%Señal generada

x = sin(2 * pi * SF * t)+sin(6 * SF * pi * t);

Figura;

Gráfico de rama (221);

plot(t, x);% se utiliza como forma de onda en el dominio del tiempo de la señal.

xlabel(' t ');

y label(' y ');

Título ("Forma de onda en el dominio del tiempo"); p>Cuadrícula;

Gráfico de rama (222);

t = 1;

t=linspace(0, T, N);

x = sin(2 * pi * SF * t)+cos(6 * SF * pi * t);

dt = t(2)-t(1);

f = 1/dt;

x = FFT(x);

f = X(1:N/2+1);

f = f *(0:N/2)/N;

Gráfico (F, abs(F)*2/N)

Xlabel('Frecuencia (Hz)' ) ;

Ylabel('amplitud');

Título('Espectrograma N = 128 ');