function [X,f] = fftej5(N,D)

% fftej5(N,D)
% Ejemplo de una FFT de una se¤al exponencial modulada en 
% frecuecia con una frecuencia portadora de 200Hz.
% N es el n£mero de puntos de muestreo durante D s de señal. 
% Elegir preferiblemente valores que sean potencias de 2
% Hacer notar que al ser la frecuencia portadora de 200Hz,
% necisitamos una frecuencia de muestreo de por lo menos
% 400Hz, por lo que N/D > 400

ts=D/N;
d=ts/2;
t=0:ts:D-d;
x=sin(2*pi*200*t-5*exp(-2*t));
X=fft(x);

%Reordenar X
M=N/2;
Xaux=X;X(M+1:N)=Xaux(1:M);
X(1:M)=Xaux(M+1:N);

Xm=abs(X)*ts;
Xf=unwrap(angle(X))*180/pi;

%Reordenar los ¡ndices k
faux(M+1:N)=0:M-1;faux(1:M)=-M:-1;
f=faux/D;

figure;plot(t,x,'-g');zoom;
title('x(t)=sin(2·pi·200·t-5·exp(-2t))');
xlabel('Tiempo (t)');ylabel('x(t)');
figure;lines(f,Xm,'oc3','-c3');zoom;
title('Módulo de los coeficientes espectrales de x(t)');
xlabel('Frecuencia (Hz)');ylabel('|X[k]|');
figure;lines(f,Xf,'oc3','-c3');zoom;
title('Fase de los coeficientes espectrales X[k]');
xlabel('Frecuencia (Hz)');ylabel('Fase X[k]');