021 MTALB Data Fourier transform

Data Fourier transform

clear all

ti_laser=xlsread('fouriertransform.xlsx');
ti_laser=ti_laser';
lamda=ti_laser(1,:).*1e-9;
intensity=ti_laser(2,:)*1e-3;
subplot(3,1,1)
plot(lamda,intensity)
c=3e8;
w=[];
for m=1:length(lamda);
    w(length(lamda)-m+1)=(2*pi*c)/lamda(m);
    w_intensity(length(lamda)-m+1)=intensity(m);
end
subplot(3,1,2)
plot(w,w_intensity)

 dscale = length(lamda); % 범위안의 간격의 개수
 dtimescale = 2000; % 시간의 범위의 간격의 개수

 t = linspace(-80e-15,80e-15,dtimescale);
 ftt=linspace(0,0,dtimescale);

 for p=1:dscale
     ft = (1/(2*pi)).*w_intensity(p).*cos(w(p).*t).*(w(2)-w(1));
     ftt=ftt+ft;
 end
 subplot(3,1,3);
 plot(t,ftt)