Utilice matlab para verificar las propiedades de la transformada de Fourier, cómo escribir un programa

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) ;

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');