Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede
|
|
- Preston Hardy
- 6 years ago
- Views:
Transcription
1 Fll Anlss of Epermentl Mesurements B. Esensten/rev. S. Errede Monte Crlo Methods/Technques: These re mong the most powerful tools for dt nlss nd smulton of eperments. The growth of ther mportnce s closel ted to dvnces n computer memor sze nd speed. We wll defne Monte Crlo technque s one whch mes use of rndom numbers n order to solve problem. The computtonl technques tht re most useful to phscsts re usull smulton nd ntegrton. (There s lttle dstncton between them!) These technques m be ppled to problems of both sttstcl nd determnstc nture. Hstorcll, the frst lrge-scle use of MC smulton ws ssocted wth studes of neutron sctterng nd bsorpton, whch re rndom/stochstc processes. Rndom numbers were used to generte populton,.e. smple of neutron sctterng events, whch ws ntended to mmc rel events. Then phscs prmeters such s bsorpton lengths were estmted from the MC events just s f the hd been collected n rel eperment. Of course, n rel eperment one cnnot see n detl the behvor of ndvdul neutrons n medum, whle n MC eperment, the mcroscopc behvor s wht s beng smulted. The ngredents of such MC method re neutron-nucleus (or nucleon) cross sectons, mcroscopc nd mcroscopc propertes of the medum, etc. (n.b. In prncple, gven these sme nput prmeters, the results could lso hve been obtned nltcll b ver comple mult-dmensonl ntegrls,.e. determnstcll.) The Bss for Monte Crlo Integrton: Consder contnuous, rel, sngle-vlued functon f ( ) defned (t lest) on the ntervl. Let us smple the functon b obtnng the vlues f ( ) f ( ) f ( ),,,,,, where re rndoml drwn from the unform dstrbuton on tht ntervl,.e. (,) = We determne the smple men: f f ( ) U. The Lw of Lrge umbers ss tht under certn condtons : lm f E[ f ( ) ] = fˆ on the ntervl. Ths cn be wrtten s: lm = ( ) ( ) f f d The epresson on the LHS of ths relton s clled the Monte Crlo estmte of the Integrl. It converges to the correct (nltc) vlue s provded tht the certn condtons hold: f ( ) must be ntegrble. f ( ) must be everwhere fnte. f ( ) must be t lest pecewse contnuous. Of course, n rel MC clculton wll remn fnte, ts mmum sze s determned b the vlble computer tme nd CPU speed. P598AEM Lecture otes 6
2 Fll Anlss of Epermentl Mesurements B. Esensten/rev. S. Errede We often csull cll f the Monte Crlo Integrl. Techncll/formll-speng, t s the MC ntegrl estmte of the ntegrl f ( ) d. The convergence of the Monte Crlo estmte lm f ( ) of the ntegrl ( ) = m be fst or slow, dependng on the nture of the functon f ( ). The convergence s, f d however, sttstcl,.e. for ver lrge (but fnte) we cn occsonll hve lrge fluctutons of f bout ˆf. Ths mens tht f ˆ f < δ onl wth some probblt. Fnll, the Centrl Lmt Theorem tell us tht the smple men f wll be (ppromtel) Gussn/normll-dstrbuted for lrge. Generlzton: Suppose nsted of the ntervl (, ) we re nterested n the ntervl (, ) where the b lm f ( ) f ( ) d b = re rndoml drwn from the unform dstrbuton (, ) U b. b. Then: If we lred hve vlble source of rndom numbers r whch re unforml dstrbuted on the ntervl (, ), we cn smpl me chnge vrbles: = + ( b ) r to obtn rndom numbers whch re unforml dstrbuted on the ntervl ( b, ), s shown n the fgure below: b Clerl, ths s lner -to- mppng of the ntervl (, ) onto the ntervl (, ) More formll, let f ( r ) be the P.D.F. of r (.e. U (,) ). Let ( ) Then: g ( ) d = f ( r) dr nd: g( ) = f ( r) Snce = + ( b ) r, then: g( ) b r r b. g be the P.D.F. of. d dr =. Ths s ( ) b U, b. It s constnt, nd ( ) g d=. P598AEM Lecture otes 6
3 Fll Anlss of Epermentl Mesurements B. Esensten/rev. S. Errede Propertes of the Monte Crlo Integrl estmte f : If the vrnce f ( ) = of f ( ) s fnte, then the MC ntegrl estmte f f ( ) s consstent,.e. t converges to the true vlue E[ f ( ) ] fˆ = = for ver lrge. The MC ntegrl estmte f f ( ) s unbsed for ll,.e. E[ f ( ) ] fˆ We hve lred seen ths propert of smple mens. The MC ntegrl estmte f f ( ) The stndrd devton f = =. s (smptotcll) Gussn/normll-dstrbuted. = of the MC ntegrl estmte f f ( ) f s: = f Ths s true for ll, but t s useful onl when the smple sze s lrge enough tht the MC ntegrl estmte s Gussn/normll-dstrbuted. An Emple of the MC Integrl Estmte - Buffon s eedle: The erlest (?) documented MC clculton ws reported b Georges-Lous LeClerc, Comte de Buffon n Ess d Arthmetque Morle (777) nd ws used to clculte the vlue of. The prncple of Buffon s method ws: -- Devse n eperment where the sttstcl probblt p of the eperment s outcome nvolved. #of successes -- Do the eperment, nd estmte the verge probblt s p =. totl #of trls -- The conceptul eperment Buffon cme up wth conssted of rndoml (n poston nd ngle) tossng needles of length d onto grd of prllel lnes tht re lso seprted b dstnce d. The verge per-trl probblt of ht needle ntersectng one of the prllel lnes (see fgure below) cn be shown to be p = : ( ) d θ mss ht All needles hve length d: Consder needle wth gven θ (= ngle between the needle nd the norml to grd lne, s shown n the bove fgure). The component of the length d of the needle perpendculr to the grd lne s d cosθ, thus the probblt tht lne wll be crossed (.e. ntersected) b needle s: p( θ ) d cosθ d cosθ of p ( θ ) over mn trls. ote tht there s no dependence on the poston of the needle, just ts ngulr orentton. = =. The verge per-trl probblt p wll be the verge P598AEM Lecture otes 6 3
4 Fll Anlss of Epermentl Mesurements B. Esensten/rev. S. Errede Snce θ s unforml (.e. rndoml) dstrbuted over the ngulr ntervl (, ), nd snce p ( θ ) = cosθ, t s therefore suffcent to verge over just the reduced ntervl (, ). Then: cosθ dθ cosθdθ θ ( ) p cosθ = = = sn = = ( ) When we estmte probblt b the rto of hts to tres usng MC methods, we cll tht the Ht-Or-Mss Monte Crlo. Ths method s generll the lest effcent mens of estmton,.e t hs lrger vrnce for gven thn n other MC method. It lso s tpcll much more CPU tme ntensve, computtonll. Let us clculte the vrnce of our estmte of s obtned b ths method: For our needle-throwng eperment, the frequenc of success s defned s: p = = #successes (n.b. = mens s estmted b ). totl# trls Snce the verge per-trl probblt of needle httng lne on the grd s p = =.637 (n.b. p s not smll!) then the verge per-trl probblt of needle mssng the prllel lnes on the grd s q p = =.637 =.363. Thus, we relze tht sttstcll, the needle ht-or-mss eperment obes bnoml probblt dstrbuton:!! P( ;, p) = p q = p p!!!! ( ) ( ) wth epectton vlue for successes n trls of: ( ) ˆ! E[ ] = P( ;, p) = p ( p) = p = =! ( )! The vrnce of s: ( ) = E ( E ) = pq= p( p) vr [ ] [ ] Hence the -sgm stndrd devton of s: ( ) pq p( p) = vr = = However, here n ths stuton we re nterested n the frequenc of success whch we clled f. (n.b. do not confuse ths wth the smple men of the MC ntegrton process!!!) It ws shown bove tht f = frequenc of success ( ) ( ) ( ), hence for fed totl # of trls : f = = pq = p p = = {here} P598AEM Lecture otes 6 4
5 Fll Anlss of Epermentl Mesurements B. Esensten/rev. S. Errede Of course, f we pretend tht we don t pror now to better thn 3.4, or /7 {or the Indn Stte Legslture s vlue of ( 3) }, or... we get slghtl dfferent vlue for f. However, the conclusons won t chnge. Fnll, snce = f, the vrnce of our MC ntegrl estmte of s: or: 4 f = f = f = 4 f f f ( ) ( ) ( ) p p = = = = = f.37 f f If we clculted v Buffon s method, we d get -stndrd devton uncertntes of,,, e.g. fter totl of,, needle throws, we d hve determned = 3.46 ±.4. Clerl, ths s not n es/ccurte w to clculte. And even f we could me gzllons of throws, we would need to worr bout phscl bses n the epermentl method e.g. lc of perfect rndomness n the throw of ech needle... We cn estmte usng nother MC technque: the Crude Monte Crlo,.e. the estmton of n ntegrl b sum: Snce: cosθdθ =, then: = cosθdθ cosθ cosθ where the verge of cosθ s ten over the ntervl (, ) We cn estmte cosθ cosθ b the smple men. If we drw rndom vlues of θ from the unform dstrbuton (, ) U nd defne: z cosθ, then: = lm z lm cosθ = (n.b. Ths s ctull rell sll, snce we need to now n order to get the dstrbuton, but... ) Ths s n unbsed estmte, snce: z cosθ = = = = = P598AEM Lecture otes 6 5
6 Fll Anlss of Epermentl Mesurements B. Esensten/rev. S. Errede et, we clculte the vrnce of z from z z : z z = cosθ cosθ j = cos cos cos θ + θ θ j = j= = j= j= j = z = + = + cos θ cos cos θ θ j = j= j= j = = j= j= j = How mn terms do we hve n ech of these sums? Consder the cse for = 4: j We see tht there re terms wth = j nd = ( ) terms wth j z = + ( ) = + Thus: z z = + = + z = z = = 4.38 z Fnll:.5 = Z compred to:.37 for the Ht-or-Mss MC. Ths s generl result: the Crude MC s lws more effcent thn the Ht-or-Mss MC. P598AEM Lecture otes 6 6
7 Fll Anlss of Epermentl Mesurements B. Esensten/rev. S. Errede We cn now determne the vrnce on n MC ntegrl estmte: Let: I f ( ). Then: I f ( ) = f ( ) =. = Also: I f f = f + f j = j= ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) f ( ) f f I I ( ) ( ) ( ) I f f f = + = Thus: I = f ( ) or: I = f ( ). We see tht we cn lws decrese ncrese b fctor of n order to decrese the stndrd devton estmte b fctor of. I b ncresng, but t s slow process. must I of the MC ntegrl However, there re some trcs tht enble us to decrese I even more, whch nvolve modfng the functon f ( ). For emple, suppose f ( ) = g( ) + h( ), nd we lso now nltcll tht: ( ) g d (= number). Then f h( ) s smll,.e. hs smller vrnce thn f ( ) cn be comprble to g( ) s) then the estmte of the MC ntegrl wll converge much fster! = G (whose vrnce ote tht there s n entre culture devoted to fndng trcs nd developng technques for mprovng the convergence of the estmtes of MC ntegrls. We wll not pursue them... If we need to crr out -dmensonl ntegrton, there s no problem n dong t b MC methods. However, there re usull (dependng on f ( ) ) more effcent computtonl/ numercl methods thn MC. On the other hnd, for problems nvolvng d-dmensonl ntegrtons, t s not so cler tht the MC pproch s less effcent. For emple, f we tbulte the -sgm uncertnt n the estmte of the MC ntegrl s functon of the number of ponts, (.e. evlutons of the functon tht one must perform), comprng s for -dmensonl vs. d-dmensonl ntegrls: uncertnt of ntegrl for: -dmenson: d-dmensons: MC: Trpezodl rule: Smpson s rule: 4 Guss rule: m d 4 d ( m ) (n.b. n the lst row, m s the order, e.g. 6 ponts, 3 ponts, etc.) d P598AEM Lecture otes 6 7
8 Fll Anlss of Epermentl Mesurements B. Esensten/rev. S. Errede Sooner or lter, for d lrge, the MC ntegrton method wll converge fster thn the other f s not well-enough behved to numercl methods. In rel lfe, oftentmes the functon ( ) relbl use the other numercl ntegrton technques, whle the MC ntegrton method demnds f {other thn stted on p. of these Lect. otes}. reltvel ver lttle from the functon ( ) How do we crr out obtnng estmtes of MC ntegrls n d-dmensons? In -dmenson, we hve lred seen tht: b b ( ) f d = (s estmted b) f ( ) = where s drwn from U(, b ) In -dmensons, the procedure s ver smple f the regon of ntegrton s rectngulr: d b d c d b (, ) f dd = f (, ) c = c b where nd numbers; re pr of ndependent rndom s drwn from (, ) drwn from U( c, d ). U b nd s It s strghtforwrd to generlze ths method for the cse of rectngulr regon n d-dmensons. In -dmensons, the regon of ntegrton often m be more comple, for emple bounded g, = C, s shown n the fgure below: b contour ( ) m g(, ) = C mn m mn Here, we wll hve to scrfce some effcenc n the MC ntegrton method. We fnd rectngle bounded b mn, m, mn, m tht entrel (but just brel) encloses the contour g(, ) = C. P598AEM Lecture otes 6 8
9 Fll Anlss of Epermentl Mesurements B. Esensten/rev. S. Errede We then unforml populte the rectngle wth rndom pont prs (, ), wth drwn from (, ) drwn from U(, ). We then test ech rndom pont pr (, ) U mn m nd mn to see f t les wthn the contour C. If t does, we eep t. If not, we go on to the net pr. Clerl the ept pont prs (, ) bounded b the contour g(, ) = C. m wll be rndoml unforml dstrbuted over the re The penlt n effcentl genertng rndom pont prs (, ) choosng the smllest possble rectngle contnng the contour g(, ) Agn, ths method cn be generlzed to d-dmensonl MC ntegrton. A generl comment: wll clerl be mnmzed b = C. In mn-dmensonl MC ntegrl t s etremel es to set up the MC ntegrton procedure. However, there s n dvntge to nltcll dong s much of the ntegrton s possble, snce then the ponts wll go further n the remnng dmensons, thus helpng to reduce the overll computton tme For emple, f we need to use = 6 ponts n 6-dmensonl MC ntegrl, we re effectvel spredng (onl!) ponts long ech of the s es. If, however, we could ntegrte out three of the vrbles, then we would get comprble ccurc wth onl 3 ponts nd the clculton would requre / th of the computer tme! P598AEM Lecture otes 6 9
Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede. with respect to λ. 1. χ λ χ λ ( ) λ, and thus:
More on χ nd errors : uppose tht we re fttng for sngle -prmeter, mnmzng: If we epnd The vlue χ ( ( ( ; ( wth respect to. χ n Tlor seres n the vcnt of ts mnmum vlue χ ( mn χ χ χ χ + + + mn mnmzes χ, nd
More information4. Eccentric axial loading, cross-section core
. Eccentrc xl lodng, cross-secton core Introducton We re strtng to consder more generl cse when the xl force nd bxl bendng ct smultneousl n the cross-secton of the br. B vrtue of Snt-Vennt s prncple we
More informationDefinition of Tracking
Trckng Defnton of Trckng Trckng: Generte some conclusons bout the moton of the scene, objects, or the cmer, gven sequence of mges. Knowng ths moton, predct where thngs re gong to project n the net mge,
More informationRemember: Project Proposals are due April 11.
Bonformtcs ecture Notes Announcements Remember: Project Proposls re due Aprl. Clss 22 Aprl 4, 2002 A. Hdden Mrov Models. Defntons Emple - Consder the emple we tled bout n clss lst tme wth the cons. However,
More informationApplied Statistics Qualifier Examination
Appled Sttstcs Qulfer Exmnton Qul_june_8 Fll 8 Instructons: () The exmnton contns 4 Questons. You re to nswer 3 out of 4 of them. () You my use ny books nd clss notes tht you mght fnd helpful n solvng
More informationUNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS. M.Sc. in Economics MICROECONOMIC THEORY I. Problem Set II
Mcroeconomc Theory I UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS MSc n Economcs MICROECONOMIC THEORY I Techng: A Lptns (Note: The number of ndctes exercse s dffculty level) ()True or flse? If V( y )
More informationRank One Update And the Google Matrix by Al Bernstein Signal Science, LLC
Introducton Rnk One Updte And the Google Mtrx y Al Bernsten Sgnl Scence, LLC www.sgnlscence.net here re two dfferent wys to perform mtrx multplctons. he frst uses dot product formulton nd the second uses
More informationPrinciple Component Analysis
Prncple Component Anlyss Jng Go SUNY Bufflo Why Dmensonlty Reducton? We hve too mny dmensons o reson bout or obtn nsghts from o vsulze oo much nose n the dt Need to reduce them to smller set of fctors
More informationCENTROID (AĞIRLIK MERKEZİ )
CENTOD (ĞLK MEKEZİ ) centrod s geometrcl concept rsng from prllel forces. Tus, onl prllel forces possess centrod. Centrod s tougt of s te pont were te wole wegt of pscl od or sstem of prtcles s lumped.
More informationInternational Journal of Pure and Applied Sciences and Technology
Int. J. Pure Appl. Sc. Technol., () (), pp. 44-49 Interntonl Journl of Pure nd Appled Scences nd Technolog ISSN 9-67 Avlle onlne t www.jopst.n Reserch Pper Numercl Soluton for Non-Lner Fredholm Integrl
More informationChapter Newton-Raphson Method of Solving a Nonlinear Equation
Chpter.4 Newton-Rphson Method of Solvng Nonlner Equton After redng ths chpter, you should be ble to:. derve the Newton-Rphson method formul,. develop the lgorthm of the Newton-Rphson method,. use the Newton-Rphson
More informationIntroduction to Numerical Integration Part II
Introducton to umercl Integrton Prt II CS 75/Mth 75 Brn T. Smth, UM, CS Dept. Sprng, 998 4/9/998 qud_ Intro to Gussn Qudrture s eore, the generl tretment chnges the ntegrton prolem to ndng the ntegrl w
More informationVariable time amplitude amplification and quantum algorithms for linear algebra. Andris Ambainis University of Latvia
Vrble tme mpltude mplfcton nd quntum lgorthms for lner lgebr Andrs Ambns Unversty of Ltv Tlk outlne. ew verson of mpltude mplfcton;. Quntum lgorthm for testng f A s sngulr; 3. Quntum lgorthm for solvng
More informationReview of linear algebra. Nuno Vasconcelos UCSD
Revew of lner lgebr Nuno Vsconcelos UCSD Vector spces Defnton: vector spce s set H where ddton nd sclr multplcton re defned nd stsf: ) +( + ) (+ )+ 5) λ H 2) + + H 6) 3) H, + 7) λ(λ ) (λλ ) 4) H, - + 8)
More informationQuiz: Experimental Physics Lab-I
Mxmum Mrks: 18 Totl tme llowed: 35 mn Quz: Expermentl Physcs Lb-I Nme: Roll no: Attempt ll questons. 1. In n experment, bll of mss 100 g s dropped from heght of 65 cm nto the snd contner, the mpct s clled
More informationLecture 36. Finite Element Methods
CE 60: Numercl Methods Lecture 36 Fnte Element Methods Course Coordntor: Dr. Suresh A. Krth, Assocte Professor, Deprtment of Cvl Engneerng, IIT Guwht. In the lst clss, we dscussed on the ppromte methods
More informationDCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x)
DCDM BUSINESS SCHOOL NUMEICAL METHODS (COS -8) Solutons to Assgnment Queston Consder the followng dt: 5 f() 8 7 5 () Set up dfference tble through fourth dfferences. (b) Wht s the mnmum degree tht n nterpoltng
More information6 Roots of Equations: Open Methods
HK Km Slghtly modfed 3//9, /8/6 Frstly wrtten t Mrch 5 6 Roots of Equtons: Open Methods Smple Fed-Pont Iterton Newton-Rphson Secnt Methods MATLAB Functon: fzero Polynomls Cse Study: Ppe Frcton Brcketng
More informationLecture 4: Piecewise Cubic Interpolation
Lecture notes on Vrtonl nd Approxmte Methods n Appled Mthemtcs - A Perce UBC Lecture 4: Pecewse Cubc Interpolton Compled 6 August 7 In ths lecture we consder pecewse cubc nterpolton n whch cubc polynoml
More informationDemand. Demand and Comparative Statics. Graphically. Marshallian Demand. ECON 370: Microeconomic Theory Summer 2004 Rice University Stanley Gilbert
Demnd Demnd nd Comrtve Sttcs ECON 370: Mcroeconomc Theory Summer 004 Rce Unversty Stnley Glbert Usng the tools we hve develoed u to ths ont, we cn now determne demnd for n ndvdul consumer We seek demnd
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 informationES.182A Topic 32 Notes Jeremy Orloff
ES.8A Topic 3 Notes Jerem Orloff 3 Polr coordintes nd double integrls 3. Polr Coordintes (, ) = (r cos(θ), r sin(θ)) r θ Stndrd,, r, θ tringle Polr coordintes re just stndrd trigonometric reltions. In
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 informationESCI 342 Atmospheric Dynamics I Lesson 1 Vectors and Vector Calculus
ESI 34 tmospherc Dnmcs I Lesson 1 Vectors nd Vector lculus Reference: Schum s Outlne Seres: Mthemtcl Hndbook of Formuls nd Tbles Suggested Redng: Mrtn Secton 1 OORDINTE SYSTEMS n orthonorml coordnte sstem
More information8. INVERSE Z-TRANSFORM
8. INVERSE Z-TRANSFORM The proce by whch Z-trnform of tme ere, nmely X(), returned to the tme domn clled the nvere Z-trnform. The nvere Z-trnform defned by: Computer tudy Z X M-fle trn.m ued to fnd nvere
More informationSUMMER KNOWHOW STUDY AND LEARNING CENTRE
SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18
More informationMath 497C Sep 17, Curves and Surfaces Fall 2004, PSU
Mth 497C Sep 17, 004 1 Curves nd Surfces Fll 004, PSU Lecture Notes 3 1.8 The generl defnton of curvture; Fox-Mlnor s Theorem Let α: [, b] R n be curve nd P = {t 0,...,t n } be prtton of [, b], then the
More informationChapter Newton-Raphson Method of Solving a Nonlinear Equation
Chpter 0.04 Newton-Rphson Method o Solvng Nonlner Equton Ater redng ths chpter, you should be ble to:. derve the Newton-Rphson method ormul,. develop the lgorthm o the Newton-Rphson method,. use the Newton-Rphson
More informationLecture 5 Single factor design and analysis
Lectue 5 Sngle fcto desgn nd nlss Completel ndomzed desgn (CRD Completel ndomzed desgn In the desgn of expements, completel ndomzed desgns e fo studng the effects of one pm fcto wthout the need to tke
More informationWork and Energy (Work Done by a Varying Force)
Lecture 1 Chpter 7 Physcs I 3.5.14 ork nd Energy (ork Done y Vryng Force) Course weste: http://fculty.uml.edu/andry_dnylov/techng/physcsi Lecture Cpture: http://echo36.uml.edu/dnylov13/physcs1fll.html
More informationLecture 3 Gaussian Probability Distribution
Introduction Lecture 3 Gussin Probbility Distribution Gussin probbility distribution is perhps the most used distribution in ll of science. lso clled bell shped curve or norml distribution Unlike the binomil
More informationM/G/1/GD/ / System. ! Pollaczek-Khinchin (PK) Equation. ! Steady-state probabilities. ! Finding L, W q, W. ! π 0 = 1 ρ
M/G//GD/ / System! Pollcze-Khnchn (PK) Equton L q 2 2 λ σ s 2( + ρ ρ! Stedy-stte probbltes! π 0 ρ! Fndng L, q, ) 2 2 M/M/R/GD/K/K System! Drw the trnston dgrm! Derve the stedy-stte probbltes:! Fnd L,L
More information1 Part II: Numerical Integration
Mth 4 Lb 1 Prt II: Numericl Integrtion This section includes severl techniques for getting pproimte numericl vlues for definite integrls without using ntiderivtives. Mthemticll, ect nswers re preferble
More informationCISE 301: Numerical Methods Lecture 5, Topic 4 Least Squares, Curve Fitting
CISE 3: umercl Methods Lecture 5 Topc 4 Lest Squres Curve Fttng Dr. Amr Khouh Term Red Chpter 7 of the tetoo c Khouh CISE3_Topc4_Lest Squre Motvton Gven set of epermentl dt 3 5. 5.9 6.3 The reltonshp etween
More informationGAUSS ELIMINATION. Consider the following system of algebraic linear equations
Numercl Anlyss for Engneers Germn Jordnn Unversty GAUSS ELIMINATION Consder the followng system of lgebrc lner equtons To solve the bove system usng clsscl methods, equton () s subtrcted from equton ()
More informationModel Fitting and Robust Regression Methods
Dertment o Comuter Engneerng Unverst o Clorn t Snt Cruz Model Fttng nd Robust Regresson Methods CMPE 64: Imge Anlss nd Comuter Vson H o Fttng lnes nd ellses to mge dt Dertment o Comuter Engneerng Unverst
More informationEffects of polarization on the reflected wave
Lecture Notes. L Ros PPLIED OPTICS Effects of polrzton on the reflected wve Ref: The Feynmn Lectures on Physcs, Vol-I, Secton 33-6 Plne of ncdence Z Plne of nterfce Fg. 1 Y Y r 1 Glss r 1 Glss Fg. Reflecton
More informationStatistics 423 Midterm Examination Winter 2009
Sttstcs 43 Mdterm Exmnton Wnter 009 Nme: e-ml: 1. Plese prnt your nme nd e-ml ddress n the bove spces.. Do not turn ths pge untl nstructed to do so. 3. Ths s closed book exmnton. You my hve your hnd clcultor
More informationON SIMPSON S INEQUALITY AND APPLICATIONS. 1. Introduction The following inequality is well known in the literature as Simpson s inequality : 2 1 f (4)
ON SIMPSON S INEQUALITY AND APPLICATIONS SS DRAGOMIR, RP AGARWAL, AND P CERONE Abstrct New neultes of Smpson type nd ther pplcton to udrture formule n Numercl Anlyss re gven Introducton The followng neulty
More informationFig. 1. Open-Loop and Closed-Loop Systems with Plant Variations
ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses
More informationThe Schur-Cohn Algorithm
Modelng, Estmton nd Otml Flterng n Sgnl Processng Mohmed Njm Coyrght 8, ISTE Ltd. Aendx F The Schur-Cohn Algorthm In ths endx, our m s to resent the Schur-Cohn lgorthm [] whch s often used s crteron for
More informationChapter 5 Supplemental Text Material R S T. ij i j ij ijk
Chpter 5 Supplementl Text Mterl 5-. Expected Men Squres n the Two-fctor Fctorl Consder the two-fctor fxed effects model y = µ + τ + β + ( τβ) + ε k R S T =,,, =,,, k =,,, n gven s Equton (5-) n the textook.
More informationPhysics 121 Sample Common Exam 2 Rev2 NOTE: ANSWERS ARE ON PAGE 7. Instructions:
Physcs 121 Smple Common Exm 2 Rev2 NOTE: ANSWERS ARE ON PAGE 7 Nme (Prnt): 4 Dgt ID: Secton: Instructons: Answer ll 27 multple choce questons. You my need to do some clculton. Answer ech queston on the
More information90 S.S. Drgomr nd (t b)du(t) =u()(b ) u(t)dt: If we dd the bove two equltes, we get (.) u()(b ) u(t)dt = p(; t)du(t) where p(; t) := for ll ; t [; b]:
RGMIA Reserch Report Collecton, Vol., No. 1, 1999 http://sc.vu.edu.u/οrgm ON THE OSTROWSKI INTEGRAL INEQUALITY FOR LIPSCHITZIAN MAPPINGS AND APPLICATIONS S.S. Drgomr Abstrct. A generlzton of Ostrowsk's
More information523 P a g e. is measured through p. should be slower for lesser values of p and faster for greater values of p. If we set p*
R. Smpth Kumr, R. Kruthk, R. Rdhkrshnn / Interntonl Journl of Engneerng Reserch nd Applctons (IJERA) ISSN: 48-96 www.jer.com Vol., Issue 4, July-August 0, pp.5-58 Constructon Of Mxed Smplng Plns Indexed
More informationLOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN FRACTAL HEAT TRANSFER
Yn, S.-P.: Locl Frctonl Lplce Seres Expnson Method for Dffuson THERMAL SCIENCE, Yer 25, Vol. 9, Suppl., pp. S3-S35 S3 LOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN
More informationPROPERTIES OF AREAS In general, and for an irregular shape, the definition of the centroid at position ( x, y) is given by
PROPERTES OF RES Centroid The concept of the centroid is prol lred fmilir to ou For plne shpe with n ovious geometric centre, (rectngle, circle) the centroid is t the centre f n re hs n is of smmetr, the
More informationElectrochemical Thermodynamics. Interfaces and Energy Conversion
CHE465/865, 2006-3, Lecture 6, 18 th Sep., 2006 Electrochemcl Thermodynmcs Interfces nd Energy Converson Where does the energy contrbuton F zϕ dn come from? Frst lw of thermodynmcs (conservton of energy):
More informationCS667 Lecture 6: Monte Carlo Integration 02/10/05
CS667 Lecture 6: Monte Crlo Integrtion 02/10/05 Venkt Krishnrj Lecturer: Steve Mrschner 1 Ide The min ide of Monte Crlo Integrtion is tht we cn estimte the vlue of n integrl by looking t lrge number of
More informationChapter Runge-Kutta 2nd Order Method for Ordinary Differential Equations
Cter. Runge-Kutt nd Order Metod or Ordnr Derentl Eutons Ater redng ts cter ou sould be ble to:. understnd te Runge-Kutt nd order metod or ordnr derentl eutons nd ow to use t to solve roblems. Wt s te Runge-Kutt
More information7.2 Volume. A cross section is the shape we get when cutting straight through an object.
7. Volume Let s revew the volume of smple sold, cylnder frst. Cylnder s volume=se re heght. As llustrted n Fgure (). Fgure ( nd (c) re specl cylnders. Fgure () s rght crculr cylnder. Fgure (c) s ox. A
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More informationES.182A Topic 30 Notes Jeremy Orloff
ES82A opic 3 Notes Jerem Orloff 3 Non-independent vribles: chin rule Recll the chin rule: If w = f, ; nd = r, t, = r, t then = + r t r t r t = + t t t r nfortuntel, sometimes there re more complicted reltions
More informationOn the Statistical Uncertainties of Time-domain-based Assessment of Stability Failures: Confidence Interval for the Mean and Variance of a Time Series
Interntonl Shp Stblt Workshop 3 Proceedngs of the 3 th Interntonl Shp Stblt Workshop, Brest 3-6 September On the Sttstcl Uncertntes of Tme-domn-bsed Assessment of Stblt Flures: Confdence Intervl for the
More informationPhysics for Scientists and Engineers I
Phscs for Scentsts nd Engneers I PHY 48, Secton 4 Dr. Betr Roldán Cuen Unverst of Centrl Flord, Phscs Deprtment, Orlndo, FL Chpter - Introducton I. Generl II. Interntonl Sstem of Unts III. Converson of
More informationThe Number of Rows which Equal Certain Row
Interntonl Journl of Algebr, Vol 5, 011, no 30, 1481-1488 he Number of Rows whch Equl Certn Row Ahmd Hbl Deprtment of mthemtcs Fcult of Scences Dmscus unverst Dmscus, Sr hblhmd1@gmlcom Abstrct Let be X
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 informationDennis Bricker, 2001 Dept of Industrial Engineering The University of Iowa. MDP: Taxi page 1
Denns Brcker, 2001 Dept of Industrl Engneerng The Unversty of Iow MDP: Tx pge 1 A tx serves three djcent towns: A, B, nd C. Ech tme the tx dschrges pssenger, the drver must choose from three possble ctons:
More information, MATHS H.O.D.: SUHAG R.KARIYA, BHOPAL, CONIC SECTION PART 8 OF
DOWNLOAD FREE FROM www.tekoclsses.com, PH.: 0 903 903 7779, 98930 5888 Some questions (Assertion Reson tpe) re given elow. Ech question contins Sttement (Assertion) nd Sttement (Reson). Ech question hs
More informationFor the percentage of full time students at RCC the symbols would be:
Mth 17/171 Chpter 7- ypothesis Testing with One Smple This chpter is s simple s the previous one, except it is more interesting In this chpter we will test clims concerning the sme prmeters tht we worked
More informationSubstitution Matrices and Alignment Statistics. Substitution Matrices
Susttuton Mtrces nd Algnment Sttstcs BMI/CS 776 www.ostt.wsc.edu/~crven/776.html Mrk Crven crven@ostt.wsc.edu Ferur 2002 Susttuton Mtrces two oulr sets of mtrces for roten seuences PAM mtrces [Dhoff et
More informationCALIBRATION OF SMALL AREA ESTIMATES IN BUSINESS SURVEYS
CALIBRATION OF SMALL AREA ESTIMATES IN BUSINESS SURVES Rodolphe Prm, Ntle Shlomo Southmpton Sttstcl Scences Reserch Insttute Unverst of Southmpton Unted Kngdom SAE, August 20 The BLUE-ETS Project s fnnced
More informationWe partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.
Mth 255 - Vector lculus II Notes 4.2 Pth nd Line Integrls We begin with discussion of pth integrls (the book clls them sclr line integrls). We will do this for function of two vribles, but these ides cn
More informationOperations with Polynomials
38 Chpter P Prerequisites P.4 Opertions with Polynomils Wht you should lern: How to identify the leding coefficients nd degrees of polynomils How to dd nd subtrct polynomils How to multiply polynomils
More informationINTRODUCTION TO COMPLEX NUMBERS
INTRODUCTION TO COMPLEX NUMBERS The numers -4, -3, -, -1, 0, 1,, 3, 4 represent the negtve nd postve rel numers termed ntegers. As one frst lerns n mddle school they cn e thought of s unt dstnce spced
More informationTopic 1 Notes Jeremy Orloff
Topic 1 Notes Jerem Orloff 1 Introduction to differentil equtions 1.1 Gols 1. Know the definition of differentil eqution. 2. Know our first nd second most importnt equtions nd their solutions. 3. Be ble
More informationdu = C dy = 1 dy = dy W is invertible with inverse U, so that y = W(t) is exactly the same thing as t = U(y),
29. Differentil equtions. The conceptul bsis of llometr Did it occur to ou in Lecture 3 wh Fiboncci would even cre how rpidl rbbit popultion grows? Mbe he wnted to et the rbbits. If so, then he would be
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 informationUniform Circular Motion
Unfom Ccul Moton Unfom ccul Moton An object mong t constnt sped n ccle The ntude of the eloct emns constnt The decton of the eloct chnges contnuousl!!!! Snce cceleton s te of chnge of eloct:!! Δ Δt The
More information( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that
Arc Length of Curves in Three Dimensionl Spce If the vector function r(t) f(t) i + g(t) j + h(t) k trces out the curve C s t vries, we cn mesure distnces long C using formul nerly identicl to one tht we
More informationStratified Extreme Ranked Set Sample With Application To Ratio Estimators
Journl of Modern Appled Sttstcl Metods Volume 3 Issue Artcle 5--004 Strtfed Extreme Rned Set Smple Wt Applcton To Rto Estmtors Hn M. Smw Sultn Qboos Unversty, smw@squ.edu.om t J. Sed Sultn Qboos Unversty
More informationTwo Coefficients of the Dyson Product
Two Coeffcents of the Dyson Product rxv:07.460v mth.co 7 Nov 007 Lun Lv, Guoce Xn, nd Yue Zhou 3,,3 Center for Combntorcs, LPMC TJKLC Nnk Unversty, Tnjn 30007, P.R. Chn lvlun@cfc.nnk.edu.cn gn@nnk.edu.cn
More informationA Family of Multivariate Abel Series Distributions. of Order k
Appled Mthemtcl Scences, Vol. 2, 2008, no. 45, 2239-2246 A Fmly of Multvrte Abel Seres Dstrbutons of Order k Rupk Gupt & Kshore K. Ds 2 Fculty of Scence & Technology, The Icf Unversty, Agrtl, Trpur, Ind
More informationPES 1120 Spring 2014, Spendier Lecture 6/Page 1
PES 110 Sprng 014, Spender Lecture 6/Page 1 Lecture today: Chapter 1) Electrc feld due to charge dstrbutons -> charged rod -> charged rng We ntroduced the electrc feld, E. I defned t as an nvsble aura
More informationPartially Observable Systems. 1 Partially Observable Markov Decision Process (POMDP) Formalism
CS294-40 Lernng for Rootcs nd Control Lecture 10-9/30/2008 Lecturer: Peter Aeel Prtlly Oservle Systems Scre: Dvd Nchum Lecture outlne POMDP formlsm Pont-sed vlue terton Glol methods: polytree, enumerton,
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 informationCHI-SQUARE DIVERGENCE AND MINIMIZATION PROBLEM
CHI-SQUARE DIVERGENCE AND MINIMIZATION PROBLEM PRANESH KUMAR AND INDER JEET TANEJA Abstrct The mnmum dcrmnton nformton prncple for the Kullbck-Lebler cross-entropy well known n the lterture In th pper
More information( ) ( )()4 x 10-6 C) ( ) = 3.6 N ( ) = "0.9 N. ( )ˆ i ' ( ) 2 ( ) 2. q 1 = 4 µc q 2 = -4 µc q 3 = 4 µc. q 1 q 2 q 3
3 Emple : Three chrges re fed long strght lne s shown n the fgure boe wth 4 µc, -4 µc, nd 3 4 µc. The dstnce between nd s. m nd the dstnce between nd 3 s lso. m. Fnd the net force on ech chrge due to the
More informationOperations with Matrices
Section. Equlit of Mtrices Opertions with Mtrices There re three ws to represent mtri.. A mtri cn be denoted b n uppercse letter, such s A, B, or C.. A mtri cn be denoted b representtive element enclosed
More informationMath 131. Numerical Integration Larson Section 4.6
Mth. Numericl Integrtion Lrson Section. This section looks t couple of methods for pproimting definite integrls numericlly. The gol is to get good pproimtion of the definite integrl in problems where n
More informationChapter 3 Single Random Variables and Probability Distributions (Part 2)
Chpter 3 Single Rndom Vriles nd Proilit Distriutions (Prt ) Contents Wht is Rndom Vrile? Proilit Distriution Functions Cumultive Distriution Function Proilit Densit Function Common Rndom Vriles nd their
More informationMultiple view geometry
EECS 442 Computer vson Multple vew geometry Perspectve Structure from Moton - Perspectve structure from moton prolem - mgutes - lgerc methods - Fctorzton methods - Bundle djustment - Self-clrton Redng:
More informationAnnouncements. Image Formation: Outline. The course. Image Formation and Cameras (cont.)
nnouncements Imge Formton nd Cmers (cont.) ssgnment : Cmer & Lenses, gd Trnsformtons, nd Homogrph wll be posted lter tod. CSE 5 Lecture 5 CS5, Fll CS5, Fll CS5, Fll The course rt : The phscs of mgng rt
More information1. Extend QR downwards to meet the x-axis at U(6, 0). y
In the digrm, two stright lines re to be drwn through so tht the lines divide the figure OPQRST into pieces of equl re Find the sum of the slopes of the lines R(6, ) S(, ) T(, 0) Determine ll liner functions
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 informationAP Calculus Multiple Choice: BC Edition Solutions
AP Clculus Multiple Choice: BC Edition Solutions J. Slon Mrch 8, 04 ) 0 dx ( x) is A) B) C) D) E) Divergent This function inside the integrl hs verticl symptotes t x =, nd the integrl bounds contin this
More informationPLEASE SCROLL DOWN FOR ARTICLE
Ths rtcle ws downloded by:ntonl Cheng Kung Unversty] On: 1 September 7 Access Detls: subscrpton number 7765748] Publsher: Tylor & Frncs Inform Ltd Regstered n Englnd nd Wles Regstered Number: 17954 Regstered
More informationAbhilasha Classes Class- XII Date: SOLUTION (Chap - 9,10,12) MM 50 Mob no
hlsh Clsses Clss- XII Dte: 0- - SOLUTION Chp - 9,0, MM 50 Mo no-996 If nd re poston vets of nd B respetvel, fnd the poston vet of pont C n B produed suh tht C B vet r C B = where = hs length nd dreton
More informationW. We shall do so one by one, starting with I 1, and we shall do it greedily, trying
Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)
More information13: Diffusion in 2 Energy Groups
3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups
More informationTorsion, Thermal Effects and Indeterminacy
ENDS Note Set 7 F007bn orson, herml Effects nd Indetermncy Deformton n orsonlly Loded Members Ax-symmetrc cross sectons subjected to xl moment or torque wll remn plne nd undstorted. At secton, nternl torque
More information4. More general extremum principles and thermodynamic potentials
4. More generl etremum prncples nd thermodynmc potentls We hve seen tht mn{u(s, X )} nd m{s(u, X)} mply one nother. Under certn condtons, these prncples re very convenent. For emple, ds = 1 T du T dv +
More information13.3. The Area Bounded by a Curve. Introduction. Prerequisites. Learning Outcomes
The Are Bounded b Curve 3.3 Introduction One of the importnt pplictions of integrtion is to find the re bounded b curve. Often such n re cn hve phsicl significnce like the work done b motor, or the distnce
More informationfractions Let s Learn to
5 simple lgebric frctions corne lens pupil retin Norml vision light focused on the retin concve lens Shortsightedness (myopi) light focused in front of the retin Corrected myopi light focused on the retin
More informationNUMERICAL MODELLING OF A CILIUM USING AN INTEGRAL EQUATION
NUEICAL ODELLING OF A CILIU USING AN INTEGAL EQUATION IHAI EBICAN, DANIEL IOAN Key words: Cl, Numercl nlyss, Electromgnetc feld, gnetton. The pper presents fst nd ccurte method to model the mgnetc behvour
More information5.1 Estimating with Finite Sums Calculus
5.1 ESTIMATING WITH FINITE SUMS Emple: Suppose from the nd to 4 th hour of our rod trip, ou trvel with the cruise control set to ectl 70 miles per hour for tht two hour stretch. How fr hve ou trveled during
More informationPHYS 2421 Fields and Waves
PHYS 242 Felds nd Wves Instucto: Joge A. López Offce: PSCI 29 A, Phone: 747-7528 Textook: Unvesty Physcs e, Young nd Feedmn 23. Electc potentl enegy 23.2 Electc potentl 23.3 Clcultng electc potentl 23.4
More informationMath Sequences and Series RETest Worksheet. Short Answer
Mth 0- Nme: Sequences nd Series RETest Worksheet Short Answer Use n infinite geometric series to express 353 s frction [ mrk, ll steps must be shown] The popultion of community ws 3 000 t the beginning
More informationBeam based calibration for beam position monitor
TBLA0 Bem bsed clbrton for bem poston montor 5-SEP-05 IBIC05 Melbourne M. Tejm KEK Contents Clbrton t the bennn Includn mppn of BPM hed lnment nd n clbrton of electrc crcut. Bem bsed lnment BBA n net er
More informationIn this Chapter. Chap. 3 Markov chains and hidden Markov models. Probabilistic Models. Example: CpG Islands
In ths Chpter Chp. 3 Mrov chns nd hdden Mrov models Bontellgence bortory School of Computer Sc. & Eng. Seoul Ntonl Unversty Seoul 5-74, Kore The probblstc model for sequence nlyss HMM (hdden Mrov model)
More information