MATLAB Answers

Area under multiple peaks (Exponentially Modified Gaussians)

9 views (last 30 days)
FW
FW on 3 Apr 2017
Edited: KSSV on 4 Apr 2017
Hello, I have numerical data from an instrument (time vs. absorbance) consisting of three overlapping peaks (attached Excel). Is there a way to numerically integrate this data from time t1(0.45 s) to t2 (2 s) to find the total area under the three peaks. Once the total area is estimated, I want to make a single exponentially modified Gaussian of the same area as the real data. Could anyone assist in the functions to be used for such a situation in MATLAB? Amplitude=A; mu= mean; lambda= variable (0.1 to 20); standard deviation =s; EMG=A*s*lambda*sqrt(pi/2)*exp(0.5*(s*lambda)^2-lambda.*(t-mu)).*erfc((1/sqrt(2))*(s*lambda-((t-mu)/s)))
Thanks for your time.

  0 Comments

Sign in to comment.

Answers (1)

KSSV
KSSV on 3 Apr 2017
Edited: KSSV on 3 Apr 2017
[num,txt,raw] = xlsread('Fast Chromatogram.xls') ;
t = num(:,1) ;
a = num(:,2) ;
idx = t>=0.45 & t <= 2 ;
ti = t(idx) ;
ai = a(idx) ;
Int = trapz(ti,ai) ;
area(ti,ai)

  2 Comments

FW
FW on 3 Apr 2017
Thanks you. Is it also possible to integrate the exponential modified Gaussian when the functional form is known?
Amplitude=A; mu= mean; lambda= variable (0.1 to 20); standard deviation =s; EMG=A*s*lambda*sqrt(pi/2)*exp(0.5*(s*lambda)^2-lambda.*(t-mu)).*erfc((1/sqrt(2))*(s*lambda-((t-mu)/s)))
Thanks once again.

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!