Encuentra la transformada de Fourier de xa=exp(-1000*abs(t)) en t=[-0.005,0.005].
Dt=0.00005;
t=-0.005:Dt:0.005;
xa=exp(-1000*abs(t));
Wmax=2*pi*2000; %Dt=0.00005 por lo que el período es 2*pi*2000
K=500;k=0:1:K
;W=k*Wmax/K % Divida Wmax en 500 puntos igualmente espaciados, W es el factor de rotación discretizado
Xa=xa*exp(-j*t'*W )*Dt;
Xa=real(Xa); %Xa=real(Xa) en realidad toma el módulo (amplitud) de cada elemento de Xa
%Transformada de Fourier en tiempo continuo
W=[-fliplr(W),W(2:501)];
% frecuencia de -Wmax a Wmax
Xa=[fliplr(Xa ), Xa (2:501)];
% Xa rango-Wmax a Wmax
figura(1)
subplot(2,1,1) ; p>
plot(t*1000,xa,'.');
xlabel('t en msec');
ylabel('xa(t )') ;
gtext('señal analógica');
subplot(2,1,2);
plot(W/(2*pi *1000) ,Xa*1000,'.');
xlabel('Frecuencia en KHz');
ylabel('Xa(jw)*1000');
p >gtext('Transformada de Fourier en tiempo continuo');