你好
试试这个
当然我不能画出光谱，因为这里没有提供
我只是让你们看看如何区分谐波和非谐波数据。
顺便说一下，数据的排序
findpeaks
不会给这里带来任何附加价值——事实上不需要;
这里我从数据开始:
轮([pkt lct '])
%轮([pkt lct '])
data = (
78 240
18 60
180年14
960年14
1050年13
12 480];
PKT = data(:，1);
LCT = data(:，2);
%谐波/非谐波群偏析
F_base = min(lct);
F_ratio = (lct./f_base);
f_residual = abs(f_ratio-round(f_ratio));
Tol = 1e-2;% 1%公差
Ind_harmo = find(f_residual <=tol);
Ind_non_harmo = find(f_residual >tol);
%谐波群
Pkt_harmo = data(ind_harmo,1);
Lct_harmo = data(ind_harmo,2);
%非谐波群
Pkt_non_harmo = data(ind_non_harmo,1);
Lct_non_harmo = data(ind_non_harmo,2);
图(1)
情节(lct_harmo pkt_harmo,“博士”lct_non_harmo pkt_non_harmo,‘* k”）;
传奇(“harmo峰”，“非谐波峰”）;
