¿Existe algún código fuente de transformada de Fourier escrito en matlab?

Encuentre 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;%Wmax se divide uniformemente en 500 puntos, W es el factor de rotación de discretización

Xa=xa*exp(-j*t'*W)*Dt ;

Xa=real(Xa); %Xa=real(Xa) es en realidad el módulo (amplitud) de cada elemento de Xa

Xa=real(Xa). modo (amplitud)

% Transformada de Fourier de 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 unidad: milisegundos');

ylabel('xa(t)');

gtext('señal analógica');

subplot(2 ,1 ,2);

plot(W/(2*pi*1000),xa*1000,'.') )

xlabel('Frecuencia en KHz'); p>

ylabel('Xa(jw)*1000');

gtext('Transformada de Fourier en tiempo continuo');