% Supongamos que tienes una señal temporal llamada 'senal' y una frecuencia de muestreo 'fs'
% Cambia esto con tus datos reales
% Parámetros para la STFT
ventana = hamming(256); % Puedes cambiar el tipo de ventana según tus necesidades
noverlap = 128; % Solapamiento entre ventanas
% Calcular la STFT
[frecuencias, tiempos, espectro] = spectrogram(senal, ventana, noverlap, [], fs);
% Calcular la frecuencia media en cada ventana
frecuencia_media = sum(bsxfun(@times, espectro, frecuencias')) ./ sum(espectro);
% Graficar la frecuencia media a lo largo del tiempo
figure;
plot(tiempos, frecuencia_media);
xlabel('Tiempo (s)');
ylabel('Frecuencia Media (Hz)');
title('Frecuencia Media a lo Largo del Tiempo');