%% 需要文件夹格式一致
clc;
close all;
clear all;
f = 92.08*1e6; % 激光器的重复频率
delta_f = 580; % 激光器的重频差
M = f/delta_f; % 异步光采样放大倍数
load('t.mat');
load('wave_all.mat');
N = length(t);
wave_all = wave_all(end:-1:1);
wave_all = wave_all-min(wave_all);
wave_all = wave_all./max(wave_all);
wave_all = wave_all-0.24;
for j = 1:N
if wave_all(j)<0
wave_all(j)=0;
end
end
global Ts ;
Ts = t(2)-t(1);
v = [(0:N/2-1),(-N/2:-1)]'/(N*Ts);
% filter out the biggest envelop
f3dB = 2e6; % 单边带带宽
wave_f = fft(abs(wave_all));
% filter
filter_f = 2.^(-0.5.*(v./(f3dB)).^10);
wave_ff = wave_f.*filter_f';
wave_af = ifft(wave_ff);
wave_af = wave_af./max(wave_af);
figure(1)
plot(t./M*1e12,wave_af,'LineWidth',1);
xlabel('time(ps)','Fontname', 'Times New Roman','fontsize',17);
ylabel('Intensity(a.u)','Fontname', 'Times New Roman','fontsize',17);
set(gca, 'Fontname', 'Times New Roman','FontSize',18);
title('waveform');
axis([0 1050 0.2 1])
% print myfigure1.jpg -djpeg -r2400 直接保存的,最后那个-r2400代表分辨率,可以调
mkdir('graphene');
path = ['graphene' '\figure1'];
saveas(1,path,'jpg');
input_mattrix=imread('C:\Users\麦兜兜是天才\Documents\Tencent Files\526320241\FileRecv\Tec_01_22_22_52_09\graphene\figure1.jpg');
fid = fopen('myfigure.dat','wt'); %写入文件途径
[m,n]=size(input_mattrix);
for i=1:1:m
for j=1;1:n
if j==n
fprintf(fid,'%g\n',input_mattrix(i,j));
else
fprintf(fid,'%g\t',input_mattrix(i,j));
end
end
end
fclose(fid);
%两种方法 你自己选