Peaks search error in a gaussian Comb

1 view (last 30 days)
MartinM
MartinM on 11 Oct 2022
Hello eveyone
I try to annalyse a comb of gaussian.
But the peaksearch function cant.
c=299792458;
LambdaC=1030; % [nm]
FreqC=LambdaC.*1e-9./c; %[Hz]
Freq1=c./(1014.5.*1e-9).*1e-12; %THz
Freq2=c./(1065.58.*1e-9).*1e-12;%THz
FreqTHZ=Freq2:0.0005:Freq1;
FreqTHZ=round(FreqTHZ*1000)/1000;
Lambda=(fliplr(c./(FreqTHZ.*1e12))).*1e9;% [nm]
Largeur= 0.1;% Largeur de la porte
DL = 0.3; % Espacement des portes
%% Gaussienne
Lc(1)=1014;
Peigne = zeros(length(Lambda),1);
for ii = 1 : length(Lambda)
Tempo=(exp(-log(2).*((Lambda-Lc(ii))/((Largeur)/(2))).^2));
Peigne=Peigne+Tempo';
Lc(ii+1)= Lc(ii)+DL;
end
[P_REF_Lambda,Idx_REF_Lambda] = findpeaks(Peigne);
I am not able to find the peak index
Do you have an idea why?
Best regards

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!