Scribd
Upload
5 views

Chapter6_Prob28

The document presents a method to model bacteria growth using the equation ln Nt - ln N0 = μ (t - t0), where μ is the growth rate constant. It provides a solution using both a user-defined function LinReg and MATLAB's built-in function polyfit to determine μ and N0 from given data. The results show that μ is approximately 1.884 and N0 is around 38.547 for both methods.

Uploaded by

josa.palaganas.up
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
5 views

Chapter6_Prob28

The document presents a method to model bacteria growth using the equation ln Nt - ln N0 = μ (t - t0), where μ is the growth rate constant. It provides a solution using both a user-defined function LinReg and MATLAB's built-in function polyfit to determine μ and N0 from given data. The results show that μ is approximately 1.884 and N0 is around 38.547 for both methods.

Uploaded by

josa.palaganas.up
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

1

6.28 Bacteria growth rate can be modeled with the equation ln N t – ln N 0 = μ ( t – t 0 ) , where μ is the
growth rate constant, and N t and N 0 are the numbers of bacteria at times t and t 0 , respectively. Determine
μ and N 0 such that the equation will best fit the following data. Use t 0 = 0 .

t (h) 0 2 4 6 8
N (cells/ml) 35 1990 70,800 2,810,000 141,250,000

(a) Use the user-defined function LinReg developed in Problem 6.19.


(b) Use MATLAB’s built-in function polyfit.
Solution
Rewriting the equation: ln N t = μt + ln N 0 – μt 0
Linear regression form: y = a 1 t + a 0
Thus: y = ln N t a 1 = μ and a 0 = ln N 0 – μt 0 .
Once a 1 and a 0 are calculated, a 0 and α are determined by μ = a 1 and N 0 = e a0 + μt0 .

(a)
A script file that solves Part (a):

clear, clc
t0=0;
t=0:2:8;
Nt=[35 1990 70800 2810000 141250000];
y=log(Nt);
[a,Er] = LinReg(t, y);
mu=a(1)
N0=exp(a(2)+mu*t0)

When the script is executed the following results are displayed in the Command Window:
mu =
1.883709642829710
N0 =
38.54703848340930

(b)

Excerpts from this work may be reproduced by instructors for distribution on a not-for-profit basis
for testing or instructional purposes only to students enrolled in courses for which the textbook
has been adopted. Any other reproduction or translation of this work beyond that permitted by
Sections 107 or 108 of the 1976 United States Copyright Act without the permission of the
copyright owner is unlawful.
2

A script file that solves Part (b):

clear, clc
t0=0;
t=0:2:8;
Nt=[35 1990 70800 2810000 141250000];
y=log(Nt);
p=polyfit(t,y,1);
mu=p(1)
N0=exp(p(2)+mu*t0)

When the script is executed the following results are displayed in the Command Window:
mu =
1.883709642829710
N0 =
38.547038483409182

Excerpts from this work may be reproduced by instructors for distribution on a not-for-profit basis
for testing or instructional purposes only to students enrolled in courses for which the textbook
has been adopted. Any other reproduction or translation of this work beyond that permitted by
Sections 107 or 108 of the 1976 United States Copyright Act without the permission of the
copyright owner is unlawful.

You might also like