creo que le vale esta progrma
%Estudio de una señal no estacionaria con contenido de 10,25,50 y 100 Hz
%y su respectivo espectro de frecuencia
clear all
close all
npuntos=1024;
w=2*pi;
t=0:1/npuntos:1;
t1=t(1:npuntos*0.3);
t2=t(npuntos*0.3:npuntos*0.6);
t3=t(npuntos*0.6:npuntos*0.8);
t4=t(npuntos*0.8:npuntos);
tiempo=1000*t;
x1=sin(w*100*t1);
x2=sin(w*50*t2);
x3=sin(w*25*t3);
x4=sin(w*10*t4);
x=[x1 x2 x3 x4];
figure(5)
plot(tiempo,x)
xlabel('Tiempo [ms]');
title('Señal estacionaria no estacionaria con contenido de 10, 25, 50 y 100 Hz');
grid on
ft=fft(x);
figure(6)
subplot(2,1,1)
plot(abs(ft))
xlabel('Frecuencia [Hz]');
title('Espectro de frecuencias de la señal de 10,25,50 y 100 Hz no estacionaria','fontweight','bold');
axis([0,npuntos,0,200]);
grid on
subplot(2,1,2)
plot(abs(ft))
xlabel('Frecuencia [Hz]');
axis([0,npuntos/3,0,200]);
grid on