fs = 16000; % frecuencia de muestreo (Hz)
fcarr = 280; % f_c (Hz)
fmod = 200; % f_m (Hz)
Ac = 1; % amplitud
Am = 2; % índice de modulación
n = 0:1/fs:5; % vector de tiempo en segundos

% Síntesis FM
y = Ac * sin( 2*pi*n.*fcarr + Am*sin(2*pi*fmod*n)).*exp(-5*n);

% Grafico el espectro
Nfft = 1048576;
df = fs / Nfft;
f = (0:(Nfft-1))*df;
Y = abs( fft(y,Nfft))/length(y);
plot(f(1:round(length(Y)/2)), Y(1:round(length(Y)/2)), 'k','Linewidth',2);
title('Magnitude spectrum');
xlabel('frequency');
axis([0 2500 0 0.015])
soundsc(y, fs); % Para escuchar el sonido