022 Matlab Fourier transform of experiments

Fourier transform of experiments

clear all
figure
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; % 시간의 범위의 간격의 개수
 n=1.52; % 렌즈의 굴절률 1.5~1.7 정도범위
 d=1e-3; % 렌즈의 두께 1mm

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

 for p=1:dscale
     ft = (1/(2*pi)).*w_intensity(p).*cos(w(p).*t+((w(p)/c)*n*d)).*(w(2)-w(1));
     ftt=ftt+ft;
 end

 %subplot(3,1,3);
 plot(t,ftt)