Derivations for maximum likelihood estimation of particle size distribution using in situ video imaging
|
|
- Moses Mosley
- 5 years ago
- Views:
Transcription
1 2 TWMCC Texs-Wisconsin Modeling nd Control Consortium 1 Technicl report numer 27-1 Derivtions for mximum likelihood estimtion of prticle size distriution using in situ video imging Pul A. Lrsen nd Jmes B. Rwlings Mrch 16, 27 Astrct This technicl report contins supplementl mteril for the mnuscript Mximum likelihood estimtion of prticle size distriution for high-spect-rtio prticles using in situ video imging [2]. This report gives the derivtion of the proility densities used to formulte the likelihood function of the MLE. 1 Mximum likelihood estimtion of PSD Let X k = (X 1k,..., X T k ) e T -dimensionl rndom vector in which X ik gives the numer of non-order prticles of size clss i oserved in imge k. A non-order prticle is prticle tht is completely enclosed within the imging volume. A order prticle, on the other hnd, is only prtilly enclosed within the imging volume such tht only portion of the prticle is oservle. For order prticles, only the oserved length (i.e. the length of the portion of the prticle tht is inside the imging volume) cn e mesured. Accordingly, we let Y k = (Y 1k,..., Y T k ) e T -dimensionl rndom vector in which Y jk gives the numer of order prticles with oserved lengths in size clss j tht re oserved in imge k. We denote the oserved dt, or the reliztions of the rndom vectors X k nd Y k, s x k nd y k, respectively. The prticle popultion is represented completely y the vectors ρ = (ρ i,..., ρ T ) nd S = (S 1,..., S T +1 ) in which ρ i represents the numer of prticles of size clss i per unit volume nd S i is the lower ound of size clss i. Given the dt x nd y (the suscript k denoting the imge index is removed here for simplicity), the mximum likelihood estimtor of ρ is defined s ˆρ = rg mx p XY (x 1, y 1, x 2, y 2,..., x T, y T ρ) (1) ρ 1
2 TWMCC Technicl Report in which in which the suscript indictes the use of order prticle mesurements nd p XY is the joint proility density for X nd Y. In other words, we wnt to determine the vlue of ρ tht mximizes the proility of oserving exctly x 1 non-order prticles of size clss 1, y 1 order prticles of size clss 1, x 2 non-order prticles of size clss 2, y 2 order prticles of size clss 2, nd so on. A simplified expression for p XY cn e otined y noting tht, t lest t low solids concentrtions, the oservtions X 1, Y 1,... X T, Y T cn e ssumed to e independent. This ssumption mens tht the oserved numer of prticles of given size clss depends only on the density of prticles in tht sme size clss. At high solids concentrtions, this ssumption seems unresonle ecuse the numer of prticle oservtions in given size clss is reduced due to occlusions y prticles in other size clsses. At low concentrtions, however, the likelihood of occlusion is low. The independence ssumption does not imply tht the oservtions re not correlted. Rther, the ssumption implies tht ny correltion etween oservtions is due to their dependence on common set of prmeters. As n exmple, if we oserve lrge numer of non-order prticles, we would expect to lso oserve lrge numer of order prticles. This correltion cn e explined y noting tht the proility densities for oth order nd non-order oservtions depend on common prmeter, nmely, the density of prticles. Given the independence ssumption, we express p XY s T T p XY = p Xi (x i ρ) p Yj (y j ρ) (2) i=1 in which p Xi nd p Yj re the proility densities for the rndom vriles X i nd Y j. Using Eqution (2), the estimtor in Eqution (1) cn e reformulted s ˆρ = rg min ρ j=1 T log p Xi (x i ρ) i=1 T log p Yj (y j ρ) (3) The proility densities p Xi nd p Yj re derived in the following sections. These derivtions show tht X i P(m Xi ), or tht X i hs Poisson distriution with prmeter m Xi = ρ i α i, in which α i is function of the field of view, depth of field, nd the lower nd upper ounds of size clss i. Furthermore, Y j P(m Yj ), in which m Yj = T i=1 ρ iβ ij 2 Derivtion of proility densities The proility densities p Xi nd p Yi in Eqution (3) cn e derived given the prticle geometry nd the sptil nd orienttionl proility distriutions. Here, we derive p Xi nd p Yi for needle-like prticles ssuming the prticles re rndomly, uniformly distriuted, oth in their 3-dimensionl sptil loction nd in their orienttion in the plne perpendiculr to the opticl xis. To simplify the discussion, we initilly present the derivtion ssuming 2-dimensionl system with monodisperse, verticlly-oriented prticles. Lter, we relx these ssumptions nd present the derivtion for rndomly-oriented, polydisperse prticles in 3-dimensionl spce. j=1
3 TWMCC Technicl Report l l/2 A n l l A l/2 I I () S () S Figure 1: Depiction of hypotheticl system of verticlly-oriented prticles rndomly nd uniformly distriuted in spce. 2.1 Non-order prticles Let S e squre domin in R 2 with dimension B nd re A S = B 2. Let I e rectngulr domin in R 2 with horizontl nd verticl dimensions nd, respectively, nd re A I =. Assume A I A S nd I S. Let n tot e the totl numer of verticllyoriented prticles with midpoints rndomly nd uniformly distriuted in S, nd define ρ = n tot /A S s the density of prticles per unit re. Let the length of ll prticles e l, with l < min (, ), nd define A n s the re of the domin in which prticles re inside I ut do not touch the order of I. Becuse the prticles re oriented verticlly, it is esy to show tht A n = ( l), s depicted in Figure 1(). Finlly, let X e rndom vrile denoting the numer of non-order prticles ppering in I. Assuming the loction of ech prticle in S is independent of the remining prticles loctions, the proility tht specific prticle will pper entirely within I is given y p = A n /A S. Given the ove ssumptions, this proility is constnt for ll prticles. The proility of oserving x non-order prticles in I is nlogous to the proility of oserving x successes in n tot Bernoulli trils in which the proility of success in ech tril is p. Thus, X is inomil rndom vrile with proility distriution p X (x) = ( ntot x ) p x (1 p) ntot x Now, ssume B while keeping ρ constnt. Then n tot nd p = A n /A S = A n ρ/n tot while Np = ρa n remins constnt. The limiting distriution of X is therefore Poisson p X (x) = e m X m x X, m X = ρa n x! To extend the nlysis to polydisperse, rndomly-oriented needles, we discretize the length scle into T size clsses nd let X = (X 1,..., X T ) e T -dimensionl rndom vector
4 TWMCC Technicl Report l l/2 sin θ A n (l, θ) θ l/2 sin θ l/2 cos θ l/2 cos θ Figure 2: Depiction of geometricl properties used to derive the non-order re function A n (l, θ). in which X i gives the numer of non-order prticles of size clss i oserved in single imge. An orienttion Θ nd length L re ssigned to ech prticle, where Θ 1, Θ 2,..., Θ ntot re i.i.d. with density function p Θ (θ), θ [ π/2, π/2) nd L 1, L 2,..., L ntot re i.i.d. with density function p L (l), l (, inf). Θ nd L re independent of ech other nd independent of the prticle s sptil loction. We define S s the T+1-dimensionl vector of reks etween size clsses. A prticle of length l elongs to size clss i if S i l < S i+1. Let i = S i+1 S i. Our gol is to determine the proility tht prticle of size clss i will pper entirely inside the imge I, given its density ρ i. Following the pproch used to solve the Buffon-Lplce needle prolem [3, p. 4], Figure 2 shows geometriclly tht for given orienttion θ nd length l, A n (l, θ) cn e clculted s A n (l, θ) = { ( l cos θ)( l sin θ) θ π/2 ( l cos θ)( + l sin θ) π/2 θ (4) The proility tht given prticle in size clss i will pper entirely within I is given y p i = Si+1 π 2 S i π 2 Si+1 π 2 S i A n (l, θ)p Θ (θ)p L (l)dθdl π 2 A S p Θ (θ)p L (l)dθdl (5) Thus, the proility tht specific prticle of size clss i will pper entirely within the imge is given y p i = α i /A S, where α i is the numertor in Eqution (5). Following the sme rguments s ove, we cn show tht for n infinitely lrge system, X i is Poisson rndom vrile with prmeter m Xi = ρ i α i. Extending the nlysis to three-dimensionl spce is trivil ecuse we ssume the prticles re oriented in the plne perpendiculr to the cmer s opticl xis nd ssume no interction etween prticles. Thus, for three-dimensionl system, X i is Poisson rndom vrile with prmeter m Xi = ρ i α i, with α i = α i d, in which d is the depth of field. Assuming Θ is distriuted uniformly nd L is distriuted uniformly cross ech size clss, α i cn e clculted s follows. Let S i = S i+1 S i, S mx = 2 + 2, S i,mx =
5 TWMCC Technicl Report S mx S i, nd ssume >. For S i+1, α i = d ( S 3 i+1 S 3 ) i + ( S 2 i+1 S 2 ) i + 3π S i π S i For < S i, S i+1, α i = d ( ( [S i+1 Si+1 2 π S sin i S i+1 ( ) ) S i ( S 2i sin + Si + 2 log S i+1 + Si ( S S i + Si 2 i+1 2 Si 2 2 ) ) For S i, S i+1 S mx, α i = d [ ( ( ) ( {S i+1 Si+1 2 π S 2 + Si sin cos i S i+1 ( ) ( ))] S i [ Si Si (sin cos Si Si + 2 log S i+1 + Si log S i+1 + Si S i + Si 2 2 S i + Si 2 2 ( ) S i 1 ( S 3 3 i+1 Si 3 ) } ) S i+1 For S i S mx nd S i+1 > S mx, [ d α i = {S mx Smx π S ( ( ) ( Smx sin cos i,mx S mx ( ) ( ))] S i [ Si Si (sin cos Si Si + 2 log S mx + Smx log S mx + Smx 2 2 S i + Si 2 2 S i + Si 2 2 ( ) S i,mx 1 ( S 3 3 mx Si 3 ) } ))] S mx ))] 2.2 Border prticles As efore, we simplify the discussion y first presenting the derivtion of p Yi for monodisperse, verticlly-oriented prticles. Let Y e rndom vrile denoting the totl numer of order prticles ppering in I. Define A s the re of the domin in which prticles
6 TWMCC Technicl Report l l A I l A l I A A A S () () (c) S S Figure 3: Depiction of hypotheticl system of verticlly-oriented prticles rndomly nd uniformly distriuted in spce. touch the order of I, s depicted in Figure 3(). For the present system, A = 2l. The proility tht specific prticle will touch the order of I is given y p = A /A S. Following the sme rguments s ove, we cn show tht for n infinitely lrge system, Y is Poisson rndom vrile with prmeter m Y = ρa. Now, ssume we would like to incorporte dditionl informtion into our estimtion y tking into ccount not only the numer of order prticles, ut lso their oserved lengths. For monodisperse popultion, these oserved lengths cn tke on vlues nywhere etween nd l. We therefore discretize the length scle on [ l] nd let j denote the size clss corresponding to the oserved length. We define Y j s rndom vrile denoting the numer of order prticles ppering in I with oserved length in size clss j. Figure 3() illustrtes this pproch for two size clsses. In this figure, A 1 is the re of the region in which prticles produce oserved lengths from to l/2, corresponding to size clss 1, while A 2 is the re of the region in which prticles produce oserved lengths from l/2 to l, corresponding to size clss 2. The proility tht specific prticle will touch the order of I nd produce n oserved length in size clss j is p = A j /A S. Thus, Y j is Poisson rndom vrile with prmeter m Yj = ρa j. In Figure 3(), A 1 = A 2. This equlity etween the res of different oserved length size clsses does not hold in generl, however, s illustrted in Figure 3(c). In this figure, we ssume ll prticles re oriented digonlly, t 45 degrees from the horizontl, nd the figure illustrtes tht A 1 > A 2. Hence, in generl, order prticles re more likely to result in oserved lengths in the lower size clsses. To extend the nlysis to polydisperse systems with rndom orienttion, we define new rndom vrile Y ij tht gives the numer of prticles in size clss i tht intersect the imge order, producing n oserved length in size clss j. Given tht the size clss of ech order prticle is unknown, we define the rndom vrile Y j s the totl numer of order prticles producing oserved lengths in size clss j, noting tht Y j = i Y ij. Our pproch is to determine the proility density for Y ij for ll i nd to use these densities to derive the proility density for Y j.
7 TWMCC Technicl Report l cos θ lθ l l (l l) sin θ θ (l l) cos θ l sin θ l θ Figure 4: Depiction of non-order re for ritrry length nd orienttion. We define the function A j (l, θ) s the re of the region in which prticle of length l nd orienttion θ produces n oserved length corresponding to size clss j. To clculte A j (l, θ), it is convenient to define n re function Ã(l, θ, l) s the re of the region in which prticles of length l nd orienttion θ either intersect or re enclosed within the imge oundry nd produce n oserved length greter thn or equl to l. Ã(l, θ, l) cn e clculted using the geometric reltionships shown in Figure 4: In this figure, the thick-lined, outer rectngle is the imge region, nd the inner rectngle is the region inside which prticle with length l nd orienttion θ will e entirely enclosed within the imge oundry, thus producing n oserved length of exctly l. A prticle with its midpoint long the perimeter of the outermost hexgon would touch the imge oundry ut give n oserved length of. A prticle with its midpoint nywhere inside the innermost hexgon will produce n oserved length greter thn or equl to l. Using the reltionships indicted in this figure, nd ssuming l l, Ã(l, θ, l) cn e clculted s { Ã(l, θ, l) ( + (l 2 l) cos θ)( + (l 2 l) sin θ) (l l) = 2 sin θ cos θ θ π/2 ( + (l 2 l) cos θ)( (l 2 l) sin θ) + (l l) 2 sin θ cos θ π/2 θ (6) If l <, Eqution (6) is vlid only for θ on ( sin (/ l), sin (/ l)). If, l < ( ) 1/2, Eqution (6) is vlid only for θ on ( sin (/ l), cos (/ l)) nd (cos (/ l), sin (/ l)). A j (l, θ) is given y Ã(l, θ, S j ) Ã(l, θ, S j+1) l S j+1 A j (l, θ) = Ã(l, θ, S j ) A n (l, θ) S j l < S j+1 (7) l < S j
8 TWMCC Technicl Report The proility tht given prticle in size clss i will pper within I nd produce n oserved length in size clss j is given y p ij = Si+1 π 2 S i π 2 Si+1 π 2 S i A j (l, θ)p Θ (θ)p L (l)dθdl π 2 A S p Θ (θ)p L (l)dθdl (8) The proility tht specific prticle in size clss i will touch the order of I nd produce n oserved length in size clss j is p ij = β ij /A S, with β ij eing the numertor in Eqution (8). Thus, for n infinitely lrge system, Y ij is Poisson rndom vrile with prmeter m Yij = ρ i βij. Assuming Y 1j, Y 2j,..., Y T j re independent, then Y j = i Y ij is lso Poisson rndom vrile with prmeter m Yj = i ρ i β ij [1, p.44]. As in the nonorder cse, the nlysis is extended to three-dimensionl spce ssuming the prticles re oriented in the plne perpendiculr to the cmer s opticl xis nd tht the prticles do not interct. Thus, for three-dimensionl system, Y j is Poisson rndom vrile with prmeter m Yj = i ρ iβ ij, with β ij = β ij d. Assuming Θ is distriuted uniformly nd L is distriuted uniformly cross ech size clss, β ij is clculted s follows. Let the length scle discretiztion e the sme for oth order nd non-order prticles. As efore, let S i = S i+1 S i, S mx = 2 + 2, S i,mx = S mx S i, nd ssume >. Then β ij is given y Ā(i, S j ) Ā(i, S j+1) i > j β ij = Ā(i, S j ) α i i = j (9) i < j in which Ā(i, S) is clculted s Ā(i, S) = d [ ( Si 2γ1 4Sγ 2 4Sγ 3 + 3S 2 ) ( γ 4 + S 2 π S i+1 S 2 ) i (γ2 + γ 3 Sγ 4 ) ] i π/2 S < sin (/S) < S < γ 1 = sin (/S) cos (/S) < S < S mx sin (/S mx ) cos (/S mx ) S > S mx 1 S < /S < S < γ 2 = ( ) S 2 2 /S < S < S mx ( ) Smx 2 2 /S mx S > S mx 1 S < 1 S 2 2 /S < S < γ 3 = ( ) S 2 2 /S < S < S mx ( ) Smx 2 2 /S mx S > S mx
9 TWMCC Technicl Report S < γ 4 = 2 /S 2 < S < ( S 2 )/S 2 < S < S mx ( Smx)/S 2 mx 2 S > S mx 3 Vlidtion of Mrginl Densities To ensure the correctness of the proility densities derived in the previous section, four different Monte Crlo simultions were crried out in which rtificil imges of prticulte popultions were generted. Figure 5 shows exmple imges generted for ech simultion. Ech of these imges hs horizontl imge dimension of =48 pixels nd verticl dimension of =48 pixels. The first row displys four simulted imges for monodisperse prticles of length.5 with N c =25 crystls per imge. The second row shows imges of prticles uniformly distriuted on [.1.9] with N c =25. The third row shows imges of prticles normlly-distriuted with µ =.5 nd σ =.4/3 with N c =25, nd the fourth row shows exmple imges for simultions of prticles uniformly-distriuted on [.1 2.] with N c =15. For ech simultion, 2, rtificil imges were generted. Bsed on the oservtions in these 2, imges, histogrm ws generted for ech size clss giving the frequency of oservtions for oth order nd non-order prticles. These histogrms re compred with the theoreticl mrginl densities in Figures 6 13.
10 TWMCC Technicl Report Figure 5: Exmple imges for simultions of vrious prticle popultions. Row 1: monodisperse prticles of length.5, Nc =25. Row 2: prticles uniformly distriuted on [.1.9]. Row 3: prticles normlly-distriuted with µ =.5 nd σ =.4/3 with Nc =25. Row 4: prticles uniformly-distriuted on [.1 2.], Nc =15.
11 TWMCC Technicl Report Figure 6: Comprison of theoreticl nd simulted mrginl densities for rndomlyoriented, monodisperse prticles of length.5 nd mesured y prtitioning [.1.9] into ten ins. Results re for non-order prticles.
12 TWMCC Technicl Report () () (c) (d) Figure 7: Comprison of theoreticl nd simulted mrginl densities for rndomlyoriented, monodisperse prticles of length.5 nd mesured y prtitioning [.1.9] into ten ins. Results re for order prticles.
13 TWMCC Technicl Report Figure 8: Comprison of theoreticl nd simulted mrginl densities for rndomly-oriented prticles distriuted uniformly on [.1.9] nd mesured y prtitioning [.1.9] into ten ins. Results re for non-order prticles.
14 TWMCC Technicl Report Figure 9: Comprison of theoreticl nd simulted mrginl densities for rndomly-oriented prticles distriuted uniformly on [.1.9] nd mesured y prtitioning [.1.9] into ten ins. Results re for order prticles.
15 TWMCC Technicl Report Figure 1: Comprison of theoreticl nd simulted mrginl densities for rndomlyoriented prticles distriuted normlly nd mesured y prtitioning [.1.9] into 1 ins. Results re for non-order prticles.
16 TWMCC Technicl Report Figure 11: Comprison of theoreticl nd simulted mrginl densities for rndomlyoriented prticles distriuted normlly nd mesured y prtitioning [.1.9] into 1 ins. Results re for order prticles.
17 TWMCC Technicl Report Figure 12: Comprison of theoreticl nd simulted mrginl densities for rndomlyoriented prticles distriuted uniformly on [.4 2.] nd mesured y prtitioning [.4 1.] into 9 ins with 1th in spnning [1. 2]. Results re for non-order prticles.
18 TWMCC Technicl Report Figure 13: Comprison of theoreticl nd simulted mrginl densities for rndomlyoriented prticles distriuted uniformly on [.4 2.] nd mesured y prtitioning [.4 1.] into 9 ins with 1th in spnning [1. 2]. Results re for order prticles.
19 TWMCC Technicl Report References [1] Yvonne M. M. Bishop, Stephen E. Fienerg, nd Pul W. Hollnd. Discrete Multivrite Anlysis: nd Prctice. The MIT Press, Cmridge, Msschusetts, [2] Pul A. Lrsen nd Jmes B. Rwlings. Mximum likelihood estimtion of prticle size distriution for high-spect-rtio prticles using in situ video imging. In preprtion. [3] Herert Solomon. Geometric Proility. SIAM Pulictions, Phildelphi, PA, 1978.
Discrete Mathematics and Probability Theory Spring 2013 Anant Sahai Lecture 17
EECS 70 Discrete Mthemtics nd Proility Theory Spring 2013 Annt Shi Lecture 17 I.I.D. Rndom Vriles Estimting the is of coin Question: We wnt to estimte the proportion p of Democrts in the US popultion,
More informationDiscrete Mathematics and Probability Theory Summer 2014 James Cook Note 17
CS 70 Discrete Mthemtics nd Proility Theory Summer 2014 Jmes Cook Note 17 I.I.D. Rndom Vriles Estimting the is of coin Question: We wnt to estimte the proportion p of Democrts in the US popultion, y tking
More informationContinuous Random Variables Class 5, Jeremy Orloff and Jonathan Bloom
Lerning Gols Continuous Rndom Vriles Clss 5, 8.05 Jeremy Orloff nd Jonthn Bloom. Know the definition of continuous rndom vrile. 2. Know the definition of the proility density function (pdf) nd cumultive
More informationMath 426: Probability Final Exam Practice
Mth 46: Probbility Finl Exm Prctice. Computtionl problems 4. Let T k (n) denote the number of prtitions of the set {,..., n} into k nonempty subsets, where k n. Argue tht T k (n) kt k (n ) + T k (n ) by
More informationReview of Gaussian Quadrature method
Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge
More information2. VECTORS AND MATRICES IN 3 DIMENSIONS
2 VECTORS AND MATRICES IN 3 DIMENSIONS 21 Extending the Theory of 2-dimensionl Vectors x A point in 3-dimensionl spce cn e represented y column vector of the form y z z-xis y-xis z x y x-xis Most of the
More informationHomework Assignment 3 Solution Set
Homework Assignment 3 Solution Set PHYCS 44 6 Ferury, 4 Prolem 1 (Griffiths.5(c The potentil due to ny continuous chrge distriution is the sum of the contriutions from ech infinitesiml chrge in the distriution.
More information4.1. Probability Density Functions
STT 1 4.1-4. 4.1. Proility Density Functions Ojectives. Continuous rndom vrile - vers - discrete rndom vrile. Proility density function. Uniform distriution nd its properties. Expected vlue nd vrince of
More information7.1 Integral as Net Change and 7.2 Areas in the Plane Calculus
7.1 Integrl s Net Chnge nd 7. Ares in the Plne Clculus 7.1 INTEGRAL AS NET CHANGE Notecrds from 7.1: Displcement vs Totl Distnce, Integrl s Net Chnge We hve lredy seen how the position of n oject cn e
More informationFarey Fractions. Rickard Fernström. U.U.D.M. Project Report 2017:24. Department of Mathematics Uppsala University
U.U.D.M. Project Report 07:4 Frey Frctions Rickrd Fernström Exmensrete i mtemtik, 5 hp Hledre: Andres Strömergsson Exmintor: Jörgen Östensson Juni 07 Deprtment of Mthemtics Uppsl University Frey Frctions
More information10. AREAS BETWEEN CURVES
. AREAS BETWEEN CURVES.. Ares etween curves So res ove the x-xis re positive nd res elow re negtive, right? Wrong! We lied! Well, when you first lern out integrtion it s convenient fiction tht s true in
More information1B40 Practical Skills
B40 Prcticl Skills Comining uncertinties from severl quntities error propgtion We usully encounter situtions where the result of n experiment is given in terms of two (or more) quntities. We then need
More informationSatellite Retrieval Data Assimilation
tellite etrievl Dt Assimiltion odgers C. D. Inverse Methods for Atmospheric ounding: Theor nd Prctice World cientific Pu. Co. Hckensck N.J. 2000 Chpter 3 nd Chpter 8 Dve uhl Artist depiction of NAA terr
More informationCalculus Module C21. Areas by Integration. Copyright This publication The Northern Alberta Institute of Technology All Rights Reserved.
Clculus Module C Ares Integrtion Copright This puliction The Northern Alert Institute of Technolog 7. All Rights Reserved. LAST REVISED Mrch, 9 Introduction to Ares Integrtion Sttement of Prerequisite
More informationPolynomials and Division Theory
Higher Checklist (Unit ) Higher Checklist (Unit ) Polynomils nd Division Theory Skill Achieved? Know tht polynomil (expression) is of the form: n x + n x n + n x n + + n x + x + 0 where the i R re the
More informationLINEAR ALGEBRA APPLIED
5.5 Applictions of Inner Product Spces 5.5 Applictions of Inner Product Spces 7 Find the cross product of two vectors in R. Find the liner or qudrtic lest squres pproimtion of function. Find the nth-order
More informationSection 6: Area, Volume, and Average Value
Chpter The Integrl Applied Clculus Section 6: Are, Volume, nd Averge Vlue Are We hve lredy used integrls to find the re etween the grph of function nd the horizontl xis. Integrls cn lso e used to find
More information( ) Same as above but m = f x = f x - symmetric to y-axis. find where f ( x) Relative: Find where f ( x) x a + lim exists ( lim f exists.
AP Clculus Finl Review Sheet solutions When you see the words This is wht you think of doing Find the zeros Set function =, fctor or use qudrtic eqution if qudrtic, grph to find zeros on clcultor Find
More information7.1 Integral as Net Change Calculus. What is the total distance traveled? What is the total displacement?
7.1 Integrl s Net Chnge Clculus 7.1 INTEGRAL AS NET CHANGE Distnce versus Displcement We hve lredy seen how the position of n oject cn e found y finding the integrl of the velocity function. The chnge
More informationContinuous Random Variables
CPSC 53 Systems Modeling nd Simultion Continuous Rndom Vriles Dr. Anirn Mhnti Deprtment of Computer Science University of Clgry mhnti@cpsc.uclgry.c Definitions A rndom vrile is sid to e continuous if there
More informationHomework Assignment 6 Solution Set
Homework Assignment 6 Solution Set PHYCS 440 Mrch, 004 Prolem (Griffiths 4.6 One wy to find the energy is to find the E nd D fields everywhere nd then integrte the energy density for those fields. We know
More information( ) as a fraction. Determine location of the highest
AB Clculus Exm Review Sheet - Solutions A. Preclculus Type prolems A1 A2 A3 A4 A5 A6 A7 This is wht you think of doing Find the zeros of f ( x). Set function equl to 0. Fctor or use qudrtic eqution if
More information( ) where f ( x ) is a. AB Calculus Exam Review Sheet. A. Precalculus Type problems. Find the zeros of f ( x).
AB Clculus Exm Review Sheet A. Preclculus Type prolems A1 Find the zeros of f ( x). This is wht you think of doing A2 A3 Find the intersection of f ( x) nd g( x). Show tht f ( x) is even. A4 Show tht f
More informationChapter 6 Continuous Random Variables and Distributions
Chpter 6 Continuous Rndom Vriles nd Distriutions Mny economic nd usiness mesures such s sles investment consumption nd cost cn hve the continuous numericl vlues so tht they cn not e represented y discrete
More information5.7 Improper Integrals
458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the
More information2 b. , a. area is S= 2π xds. Again, understand where these formulas came from (pages ).
AP Clculus BC Review Chpter 8 Prt nd Chpter 9 Things to Know nd Be Ale to Do Know everything from the first prt of Chpter 8 Given n integrnd figure out how to ntidifferentite it using ny of the following
More informationThe Dirichlet Problem in a Two Dimensional Rectangle. Section 13.5
The Dirichlet Prolem in Two Dimensionl Rectngle Section 13.5 1 Dirichlet Prolem in Rectngle In these notes we will pply the method of seprtion of vriles to otin solutions to elliptic prolems in rectngle
More informationLecture 3: Curves in Calculus. Table of contents
Mth 348 Fll 7 Lecture 3: Curves in Clculus Disclimer. As we hve textook, this lecture note is for guidnce nd supplement only. It should not e relied on when prepring for exms. In this lecture we set up
More informationRiemann is the Mann! (But Lebesgue may besgue to differ.)
Riemnn is the Mnn! (But Lebesgue my besgue to differ.) Leo Livshits My 2, 2008 1 For finite intervls in R We hve seen in clss tht every continuous function f : [, b] R hs the property tht for every ɛ >
More informationAB Calculus Review Sheet
AB Clculus Review Sheet Legend: A Preclculus, B Limits, C Differentil Clculus, D Applictions of Differentil Clculus, E Integrl Clculus, F Applictions of Integrl Clculus, G Prticle Motion nd Rtes This is
More informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More information( ) where f ( x ) is a. AB/BC Calculus Exam Review Sheet. A. Precalculus Type problems. Find the zeros of f ( x).
AB/ Clculus Exm Review Sheet A. Preclculus Type prolems A1 Find the zeros of f ( x). This is wht you think of doing A2 Find the intersection of f ( x) nd g( x). A3 Show tht f ( x) is even. A4 Show tht
More informationSeptember 13 Homework Solutions
College of Engineering nd Computer Science Mechnicl Engineering Deprtment Mechnicl Engineering 5A Seminr in Engineering Anlysis Fll Ticket: 5966 Instructor: Lrry Cretto Septemer Homework Solutions. Are
More informationConsequently, the temperature must be the same at each point in the cross section at x. Let:
HW 2 Comments: L1-3. Derive the het eqution for n inhomogeneous rod where the therml coefficients used in the derivtion of the het eqution for homogeneous rod now become functions of position x in the
More informationP 1 (x 1, y 1 ) is given by,.
MA00 Clculus nd Bsic Liner Alger I Chpter Coordinte Geometr nd Conic Sections Review In the rectngulr/crtesin coordintes sstem, we descrie the loction of points using coordintes. P (, ) P(, ) O The distnce
More information[ ( ) ( )] Section 6.1 Area of Regions between two Curves. Goals: 1. To find the area between two curves
Gols: 1. To find the re etween two curves Section 6.1 Are of Regions etween two Curves I. Are of Region Between Two Curves A. Grphicl Represention = _ B. Integrl Represention [ ( ) ( )] f x g x dx = C.
More informationHeavy tail and stable distributions
Hevy til nd stle distriutions J.K. Misiewicz Deprtment of Mthemtics nd Informtion Science Technicl University of Wrsw X or its distriution hs hevy til if E(X ) =. X or its distriution hs hevy til of order
More informationReview of Probability Distributions. CS1538: Introduction to Simulations
Review of Proility Distriutions CS1538: Introduction to Simultions Some Well-Known Proility Distriutions Bernoulli Binomil Geometric Negtive Binomil Poisson Uniform Exponentil Gmm Erlng Gussin/Norml Relevnce
More informationFully Kinetic Simulations of Ion Beam Neutralization
Fully Kinetic Simultions of Ion Bem Neutrliztion Joseph Wng University of Southern Cliforni Hideyuki Usui Kyoto University E-mil: josephjw@usc.edu; usui@rish.kyoto-u.c.jp 1. Introduction Ion em emission/neutrliztion
More informationGenetic Programming. Outline. Evolutionary Strategies. Evolutionary strategies Genetic programming Summary
Outline Genetic Progrmming Evolutionry strtegies Genetic progrmming Summry Bsed on the mteril provided y Professor Michel Negnevitsky Evolutionry Strtegies An pproch simulting nturl evolution ws proposed
More information5: The Definite Integral
5: The Definite Integrl 5.: Estimting with Finite Sums Consider moving oject its velocity (meters per second) t ny time (seconds) is given y v t = t+. Cn we use this informtion to determine the distnce
More informationQuadratic Forms. Quadratic Forms
Qudrtic Forms Recll the Simon & Blume excerpt from n erlier lecture which sid tht the min tsk of clculus is to pproximte nonliner functions with liner functions. It s ctully more ccurte to sy tht we pproximte
More informationMTH 505: Number Theory Spring 2017
MTH 505: Numer Theory Spring 207 Homework 2 Drew Armstrong The Froenius Coin Prolem. Consider the eqution x ` y c where,, c, x, y re nturl numers. We cn think of $ nd $ s two denomintions of coins nd $c
More informationDesigning Information Devices and Systems I Spring 2018 Homework 7
EECS 16A Designing Informtion Devices nd Systems I Spring 2018 omework 7 This homework is due Mrch 12, 2018, t 23:59. Self-grdes re due Mrch 15, 2018, t 23:59. Sumission Formt Your homework sumission should
More informationProperties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives
Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn
More informationCoimisiún na Scrúduithe Stáit State Examinations Commission
M 30 Coimisiún n Scrúduithe Stáit Stte Exmintions Commission LEAVING CERTIFICATE EXAMINATION, 005 MATHEMATICS HIGHER LEVEL PAPER ( 300 mrks ) MONDAY, 3 JUNE MORNING, 9:30 to :00 Attempt FIVE questions
More informationI1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3
2 The Prllel Circuit Electric Circuits: Figure 2- elow show ttery nd multiple resistors rrnged in prllel. Ech resistor receives portion of the current from the ttery sed on its resistnce. The split is
More informationContinuous Joint Distributions Chris Piech CS109, Stanford University
Continuous Joint Distriutions Chris Piech CS09, Stnford University CS09 Flow Tody Discrete Joint Distriutions: Generl Cse Multinomil: A prmetric Discrete Joint Cont. Joint Distriutions: Generl Cse Lerning
More informationTests for the Ratio of Two Poisson Rates
Chpter 437 Tests for the Rtio of Two Poisson Rtes Introduction The Poisson probbility lw gives the probbility distribution of the number of events occurring in specified intervl of time or spce. The Poisson
More informationLecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.
Lecture 3 3 Solving liner equtions In this lecture we will discuss lgorithms for solving systems of liner equtions Multiplictive identity Let us restrict ourselves to considering squre mtrices since one
More informationChapter 4 Contravariance, Covariance, and Spacetime Diagrams
Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz
More information10 Vector Integral Calculus
Vector Integrl lculus Vector integrl clculus extends integrls s known from clculus to integrls over curves ("line integrls"), surfces ("surfce integrls") nd solids ("volume integrls"). These integrls hve
More informationThe graphs of Rational Functions
Lecture 4 5A: The its of Rtionl Functions s x nd s x + The grphs of Rtionl Functions The grphs of rtionl functions hve severl differences compred to power functions. One of the differences is the behvior
More informationBases for Vector Spaces
Bses for Vector Spces 2-26-25 A set is independent if, roughly speking, there is no redundncy in the set: You cn t uild ny vector in the set s liner comintion of the others A set spns if you cn uild everything
More informationSection 4: Integration ECO4112F 2011
Reding: Ching Chpter Section : Integrtion ECOF Note: These notes do not fully cover the mteril in Ching, ut re ment to supplement your reding in Ching. Thus fr the optimistion you hve covered hs een sttic
More informationBME 207 Introduction to Biomechanics Spring 2018
April 6, 28 UNIVERSITY O RHODE ISAND Deprtment of Electricl, Computer nd Biomedicl Engineering BME 27 Introduction to Biomechnics Spring 28 Homework 8 Prolem 14.6 in the textook. In ddition to prts -e,
More informationSection 6.1 Definite Integral
Section 6.1 Definite Integrl Suppose we wnt to find the re of region tht is not so nicely shped. For exmple, consider the function shown elow. The re elow the curve nd ove the x xis cnnot e determined
More informationGolden Section Search Method - Theory
Numericl Methods Golden Section Serch Method - Theory http://nm.mthforcollege.com For more detils on this topic Go to http://nm.mthforcollege.com Click on Keyword Click on Golden Section Serch Method You
More informationStudent Activity 3: Single Factor ANOVA
MATH 40 Student Activity 3: Single Fctor ANOVA Some Bsic Concepts In designed experiment, two or more tretments, or combintions of tretments, is pplied to experimentl units The number of tretments, whether
More informationSolution for Assignment 1 : Intro to Probability and Statistics, PAC learning
Solution for Assignment 1 : Intro to Probbility nd Sttistics, PAC lerning 10-701/15-781: Mchine Lerning (Fll 004) Due: Sept. 30th 004, Thursdy, Strt of clss Question 1. Bsic Probbility ( 18 pts) 1.1 (
More informationSuppose we want to find the area under the parabola and above the x axis, between the lines x = 2 and x = -2.
Mth 43 Section 6. Section 6.: Definite Integrl Suppose we wnt to find the re of region tht is not so nicely shped. For exmple, consider the function shown elow. The re elow the curve nd ove the x xis cnnot
More informationThe steps of the hypothesis test
ttisticl Methods I (EXT 7005) Pge 78 Mosquito species Time of dy A B C Mid morning 0.0088 5.4900 5.5000 Mid Afternoon.3400 0.0300 0.8700 Dusk 0.600 5.400 3.000 The Chi squre test sttistic is the sum of
More informationDefinite Integrals. The area under a curve can be approximated by adding up the areas of rectangles = 1 1 +
Definite Integrls --5 The re under curve cn e pproximted y dding up the res of rectngles. Exmple. Approximte the re under y = from x = to x = using equl suintervls nd + x evluting the function t the left-hnd
More informationChapter 2. Random Variables and Probability Distributions
Rndom Vriles nd Proilit Distriutions- 6 Chpter. Rndom Vriles nd Proilit Distriutions.. Introduction In the previous chpter, we introduced common topics of proilit. In this chpter, we trnslte those concepts
More informationNon-Linear & Logistic Regression
Non-Liner & Logistic Regression If the sttistics re boring, then you've got the wrong numbers. Edwrd R. Tufte (Sttistics Professor, Yle University) Regression Anlyses When do we use these? PART 1: find
More informationSection 3.2 Maximum Principle and Uniqueness
Section 3. Mximum Principle nd Uniqueness Let u (x; y) e smooth solution in. Then the mximum vlue exists nd is nite. (x ; y ) ; i.e., M mx fu (x; y) j (x; y) in g Furthermore, this vlue cn e otined y point
More informationMATH 144: Business Calculus Final Review
MATH 144: Business Clculus Finl Review 1 Skills 1. Clculte severl limits. 2. Find verticl nd horizontl symptotes for given rtionl function. 3. Clculte derivtive by definition. 4. Clculte severl derivtives
More informationM344 - ADVANCED ENGINEERING MATHEMATICS
M3 - ADVANCED ENGINEERING MATHEMATICS Lecture 18: Lplce s Eqution, Anltic nd Numericl Solution Our emple of n elliptic prtil differentil eqution is Lplce s eqution, lso clled the Diffusion Eqution. If
More informationragsdale (zdr82) HW2 ditmire (58335) 1
rgsdle (zdr82) HW2 ditmire (58335) This print-out should hve 22 questions. Multiple-choice questions my continue on the next column or pge find ll choices before nswering. 00 0.0 points A chrge of 8. µc
More informationThe Wave Equation I. MA 436 Kurt Bryan
1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string
More information4.4 Areas, Integrals and Antiderivatives
. res, integrls nd ntiderivtives 333. Ares, Integrls nd Antiderivtives This section explores properties of functions defined s res nd exmines some connections mong res, integrls nd ntiderivtives. In order
More informationThings to Memorize: A Partial List. January 27, 2017
Things to Memorize: A Prtil List Jnury 27, 2017 Chpter 2 Vectors - Bsic Fcts A vector hs mgnitude (lso clled size/length/norm) nd direction. It does not hve fixed position, so the sme vector cn e moved
More informationMath 0230 Calculus 2 Lectures
Mth Clculus Lectures Chpter 7 Applictions of Integrtion Numertion of sections corresponds to the text Jmes Stewrt, Essentil Clculus, Erly Trnscendentls, Second edition. Section 7. Ares Between Curves Two
More informationLecture Solution of a System of Linear Equation
ChE Lecture Notes, Dept. of Chemicl Engineering, Univ. of TN, Knoville - D. Keffer, 5/9/98 (updted /) Lecture 8- - Solution of System of Liner Eqution 8. Why is it importnt to e le to solve system of liner
More informationTImath.com Algebra 2. Constructing an Ellipse
TImth.com Algebr Constructing n Ellipse ID: 9980 Time required 60 minutes Activity Overview This ctivity introduces ellipses from geometric perspective. Two different methods for constructing n ellipse
More information4 VECTORS. 4.0 Introduction. Objectives. Activity 1
4 VECTRS Chpter 4 Vectors jectives fter studying this chpter you should understnd the difference etween vectors nd sclrs; e le to find the mgnitude nd direction of vector; e le to dd vectors, nd multiply
More informationp-adic Egyptian Fractions
p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction
More informationGoals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite
Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite
More informationKEY CONCEPTS. satisfies the differential equation da. = 0. Note : If F (x) is any integral of f (x) then, x a
KEY CONCEPTS THINGS TO REMEMBER :. The re ounded y the curve y = f(), the -is nd the ordintes t = & = is given y, A = f () d = y d.. If the re is elow the is then A is negtive. The convention is to consider
More informationChapter 0. What is the Lebesgue integral about?
Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous
More informationThe Shortest Confidence Interval for the Mean of a Normal Distribution
Interntionl Journl of Sttistics nd Proility; Vol. 7, No. 2; Mrch 208 ISSN 927-7032 E-ISSN 927-7040 Pulished y Cndin Center of Science nd Eduction The Shortest Confidence Intervl for the Men of Norml Distriution
More informationCalculus of Variations
Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function
More informationSurface maps into free groups
Surfce mps into free groups lden Wlker Novemer 10, 2014 Free groups wedge X of two circles: Set F = π 1 (X ) =,. We write cpitl letters for inverse, so = 1. e.g. () 1 = Commuttors Let x nd y e loops. The
More information38 Riemann sums and existence of the definite integral.
38 Riemnn sums nd existence of the definite integrl. In the clcultion of the re of the region X bounded by the grph of g(x) = x 2, the x-xis nd 0 x b, two sums ppered: ( n (k 1) 2) b 3 n 3 re(x) ( n These
More informationKINEMATICS OF RIGID BODIES
KINEMTICS OF RIGID ODIES Introduction In rigid body kinemtics, e use the reltionships governing the displcement, velocity nd ccelertion, but must lso ccount for the rottionl motion of the body. Description
More informationHow do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?
XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk out solving systems of liner equtions. These re prolems tht give couple of equtions with couple of unknowns, like: 6= x + x 7=
More informationPartial Derivatives. Limits. For a single variable function f (x), the limit lim
Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the right-hnd side limit equls to the left-hnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles
More informationMORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.)
MORE FUNCTION GRAPHING; OPTIMIZATION FRI, OCT 25, 203 (Lst edited October 28, 203 t :09pm.) Exercise. Let n be n rbitrry positive integer. Give n exmple of function with exctly n verticl symptotes. Give
More informationBest Approximation. Chapter The General Case
Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given
More informationDesigning Information Devices and Systems I Discussion 8B
Lst Updted: 2018-10-17 19:40 1 EECS 16A Fll 2018 Designing Informtion Devices nd Systems I Discussion 8B 1. Why Bother With Thévenin Anywy? () Find Thévenin eqiuvlent for the circuit shown elow. 2kΩ 5V
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More information5.1 How do we Measure Distance Traveled given Velocity? Student Notes
. How do we Mesure Distnce Trveled given Velocity? Student Notes EX ) The tle contins velocities of moving cr in ft/sec for time t in seconds: time (sec) 3 velocity (ft/sec) 3 A) Lel the x-xis & y-xis
More informationThe area under the graph of f and above the x-axis between a and b is denoted by. f(x) dx. π O
1 Section 5. The Definite Integrl Suppose tht function f is continuous nd positive over n intervl [, ]. y = f(x) x The re under the grph of f nd ove the x-xis etween nd is denoted y f(x) dx nd clled the
More informationRiemann Integrals and the Fundamental Theorem of Calculus
Riemnn Integrls nd the Fundmentl Theorem of Clculus Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University September 16, 2013 Outline Grphing Riemnn Sums
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More informationLecture 4 Notes, Electromagnetic Theory I Dr. Christopher S. Baird University of Massachusetts Lowell
Lecture 4 Notes, Electromgnetic Theory I Dr. Christopher S. Bird University of Msschusetts Lowell 1. Orthogonl Functions nd Expnsions - In the intervl (, ) of the vrile x, set of rel or complex functions
More informationList all of the possible rational roots of each equation. Then find all solutions (both real and imaginary) of the equation. 1.
Mth Anlysis CP WS 4.X- Section 4.-4.4 Review Complete ech question without the use of grphing clcultor.. Compre the mening of the words: roots, zeros nd fctors.. Determine whether - is root of 0. Show
More informationHomework Solution - Set 5 Due: Friday 10/03/08
CE 96 Introduction to the Theory of Computtion ll 2008 Homework olution - et 5 Due: ridy 10/0/08 1. Textook, Pge 86, Exercise 1.21. () 1 2 Add new strt stte nd finl stte. Mke originl finl stte non-finl.
More informationNumerical integration
2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter
More information