clear all
close all
k=2;
data=[];
data2=[];%取0.75THz处的点的值
for i1=1:31
stringA=strcat('1x',num2str(i1),'y.txt');
a1=load(stringA);
x=a1(:,1);
y=a1(:,2)-mean(a1(:,2));
fy=fft(y,256);
A=abs(fy);
time=[0:255]*0.02e-3*2/3e8;
freq=eps:1/time(length(time)):(length(time)-1)/time(length(time));
% plot(freq(1:90),A(1:90))
data(:,1)=freq(1:90);
data(:,k)=A(1:90);
data2(i1,1)=5+(i1-1)*0.5;
data2(i1,2)=A(10); %freq(10)代表0.26THz
i1
k=k+1;
end
save data.txt -ascii data
save data2.txt -ascii data2
plot(data2(:,1),data2(:,2)*1e9,'linewidth',2.5)
title('0.26THz')
set(gca,'fontsize',20)
% save('imageset.mat','imageset')
% hh=sum(sum(imageset));
% hh=reshape(hh,126,1);%改点数
% plot(hh)
