clc
clear
filename='20221012203847+0.005cc+10ps+2553ppmv-1.gcidx.mass.csv';
Data = readmatrix(filename);%从文件中读取矩阵
A=Data(:,2);%质荷比
B=Data(:,3);%强度
data=Data(:,2:3);%质荷比+强度
x1=find(A>=27.8&A<=28.2);%元素是N2
[pks1,locs1,widths1,proms1] = findpeaks(B(x1),A(x1),'MinPeakDistance',0.2);
subplot(2,2,1)
findpeaks(B(x1),A(x1),'MinPeakProminence',4,'Annotate','extents','MinPeakDistance',0.2)
title('N2');
xlabel('质荷比');
ylabel('强度');
x2=find(A>=17.8&A<=18.2);%元素是H20
[pks2,locs2,widths2,proms2] = findpeaks(B(x2),A(x2),'MinPeakDistance',0.2);
subplot(2,2,2)
findpeaks(B(x2),A(x2),'MinPeakProminence',4,'Annotate','extents','MinPeakDistance',0.2)
title('H20');
xlabel('质荷比'