A flat histogram stochastic growth algorithm

Size: px
Start display at page:

Download "A flat histogram stochastic growth algorithm"

Transcription

1 A flat histogram stochastic growth algorithm (with applicatio to polymers ad proteis) Thomas Prellberg TU Clausthal Joit work with Jaroslaw Krawczyk, TU Clausthal Thomas Prellberg, Göttige, Ja 4, 24 p./24

2 Itroductio Thomas Prellberg, Göttige, Ja 4, 24 p.2/24

3 Modellig of Polymers i Solutio Polymers: log chais of moomers Coarse-Graiig : beads o a chai Excluded Volume : miimal distace Cotact with solvet: effective short-rage iteractio Good/bad solvet: repellig/attractig iteractio Cosequece: chais clump together Eight polymers i a bad solvet (Grassberger, FZ Jülich) Thomas Prellberg, Göttige, Ja 4, 24 p.3/24

4 Lattice Model Self-Avoidig Walks with Iteractios Physical space lattice Z 3 (or Z d ) Polymer i solutio radom walk with self-avoidace Quality of solvet short-rage iteractio ɛ Properties of moomers i, j iteractio ɛ = ɛ i,j Thomas Prellberg, Göttige, Ja 4, 24 p.4/24

5 Lattice Model Self-Avoidig Walks with Iteractios Physical space lattice Z 3 (or Z d ) Polymer i solutio radom walk with self-avoidace Quality of solvet short-rage iteractio ɛ Properties of moomers i, j iteractio ɛ = ɛ i,j Two examples: ISAW model: iteractio ɛ i,j = of iterest: thermodyamic limit (V = ad ) HP model: two types of moomers H ad P ɛ HH =, ɛ HP = ɛ P H = ɛ P P = of iterest: fixed fiite sequece, desity of states Thomas Prellberg, Göttige, Ja 4, 24 p.4/24

6 Lattice Model: ISAW = 26 m = 7 2 e βɛ i,j Boltzma weight ISAW model: iteractio ɛ i,j =, ω = e β Partitio fuctio: Z (ω) = m C,mω m C,m umber of SAW with legth ad m iteractios Thomas Prellberg, Göttige, Ja 4, 24 p.5/24

7 Lattice Model: HP model HP model: iteractio ɛhh =, ɛhp = ɛp H = ɛp P = Hydrophobic Polar Groudstate of sequece with 85 moomers (d = 2) Thomas Prellberg, Göttige, Ja 4, 24 p.6/24

8 Why Simulatios? ISAW model: Tricritical phase trasitio, d u = 3 I priciple uderstood, however surprisig details e.g. pseudo-first-order trasitio for d > 3 No good uderstadig of collapsed regime Thomas Prellberg, Göttige, Ja 4, 24 p.7/24

9 Why Simulatios? ISAW model: Tricritical phase trasitio, d u = 3 I priciple uderstood, however surprisig details e.g. pseudo-first-order trasitio for d > 3 No good uderstadig of collapsed regime HP model: Toy model for proteis Desig of sequeces with specific groud state structure Desity of states foldig dyamics Thomas Prellberg, Göttige, Ja 4, 24 p.7/24

10 Why Simulatios? ISAW model: Tricritical phase trasitio, d u = 3 I priciple uderstood, however surprisig details e.g. pseudo-first-order trasitio for d > 3 No good uderstadig of collapsed regime HP model: Toy model for proteis Desig of sequeces with specific groud state structure Desity of states foldig dyamics Most iterestig ope questios i collapsed regime Collapsed regime is otoriously difficult to simulate Thomas Prellberg, Göttige, Ja 4, 24 p.7/24

11 Stochastic Growth Algorithms Thomas Prellberg, Göttige, Ja 4, 24 p.8/24

12 PERM: Go With The Wiers PERM = Prued ad Eriched Rosebluth Method Rosebluth Method: kietic growth Grassberger, Phys Rev E 56 (997) 3682 /3 /2 trapped Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

13 PERM: Go With The Wiers PERM = Prued ad Eriched Rosebluth Method Rosebluth Method: kietic growth Grassberger, Phys Rev E 56 (997) 3682 /3 /2 trapped Erichmet: weight too large make copies of cofiguratio Pruig: weight too small remove cofiguratio occasioally Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

14 PERM: Go With The Wiers PERM = Prued ad Eriched Rosebluth Method Rosebluth Method: kietic growth Grassberger, Phys Rev E 56 (997) 3682 /3 /2 trapped Erichmet: weight too large make copies of cofiguratio Pruig: weight too small remove cofiguratio occasioally Observatio: kietic growth weights ad iteractios balace each other at suitable temperatures (i the collapse regio) Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

15 PERM cotiued PERM well suited for simulatio of collapsig polymers Erichmet creates correlated tours = = max Error estimatio a priori impossible, a posteriori difficult Thomas Prellberg, Göttige, Ja 4, 24 p./24

16 PERM cotiued PERM well suited for simulatio of collapsig polymers Erichmet creates correlated tours = = max Error estimatio a priori impossible, a posteriori difficult PERM is blid - it beefits from visual aids, such as Markovia aticipatio (learig from experiece) Thomas Prellberg, Göttige, Ja 4, 24 p./24

17 PERM further developmets A sigificat improvemet: PERM = ew PERM Eforce distict erichmet steps Crossover to exact eumeratio Hsu et al, J Chem Phys 8 (23) 444 Thomas Prellberg, Göttige, Ja 4, 24 p./24

18 PERM further developmets A sigificat improvemet: PERM = ew PERM Eforce distict erichmet steps Crossover to exact eumeratio Curret work: flatperm = flat histogram PERM Hsu et al, J Chem Phys 8 (23) 444 TP ad JK, cod-mat/32253, PRL (24) flatperm samples a geeralised multicaoical esemble Covers the whole temperature rage i oe simulatio! Related: multicaoical chai growth algorithm Bachma ad Jake, PRL (23) Thomas Prellberg, Göttige, Ja 4, 24 p./24

19 PERM further developmets A sigificat improvemet: PERM = ew PERM Eforce distict erichmet steps Crossover to exact eumeratio Curret work: flatperm = flat histogram PERM Hsu et al, J Chem Phys 8 (23) 444 TP ad JK, cod-mat/32253, PRL (24) flatperm samples a geeralised multicaoical esemble Covers the whole temperature rage i oe simulatio! Related: multicaoical chai growth algorithm Bachma ad Jake, PRL (23) Applicatios: liear ad brached polymers, proteis, percolatio,... Thomas Prellberg, Göttige, Ja 4, 24 p./24

20 Algorithm details - kietic growth View kietic growth as approximate eumeratio Thomas Prellberg, Göttige, Ja 4, 24 p.2/24

21 Algorithm details - kietic growth View kietic growth as approximate eumeratio Exact eumeratio: choose all a cotiuatios with equal weight Kietic growth: chose oe cotiuatio with a-fold weight (a may be zero). Thomas Prellberg, Göttige, Ja 4, 24 p.2/24

22 Algorithm details - kietic growth View kietic growth as approximate eumeratio Exact eumeratio: choose all a cotiuatios with equal weight Kietic growth: chose oe cotiuatio with a-fold weight (a may be zero). A step cofiguratio gets assiged a weight W = k= S growth chais with weights W (i) a k give estimate C est = W = S i W (i) Thomas Prellberg, Göttige, Ja 4, 24 p.2/24

23 Algorithm details - pruig/erichmet W (i) is estimate of C Cosider ratio r = W (i) /C est r > erichmet step: make c = mi( r, a ) distict copies with weight c W (i) r < pruig step: cotiue growig with probability r ad weight C est Thomas Prellberg, Göttige, Ja 4, 24 p.3/24

24 Algorithm details - pruig/erichmet W (i) is estimate of C Cosider ratio r = W (i) /C est r > erichmet step: make c = mi( r, a ) distict copies with weight c W (i) r < pruig step: cotiue growig with probability r ad weight C est Cosequeces Number of samples geerated for each is roughly costat Ideally, ubiased radom walk i cofiguratio size We have a flat histogram algorithm Thomas Prellberg, Göttige, Ja 4, 24 p.3/24

25 From PERM to flatperm PERM: estimate umber of cofiguratios C C est = W r = W (i) /C est Thomas Prellberg, Göttige, Ja 4, 24 p.4/24

26 From PERM to flatperm PERM: estimate umber of cofiguratios C C est = W r = W (i) /C est Cosider eergy E m = ɛm, temperature β = /k B T thermal PERM: estimate partitio fuctio Z (β) Z est (β) = W exp( βe) r = W (i) exp( βe m (i) )/Z est (β) Thomas Prellberg, Göttige, Ja 4, 24 p.4/24

27 From PERM to flatperm PERM: estimate umber of cofiguratios C C est = W r = W (i) /C est Cosider eergy E m = ɛm, temperature β = /k B T thermal PERM: estimate partitio fuctio Z (β) Z est (β) = W exp( βe) r = W (i) exp( βe m (i) )/Z est (β) flatperm: estimate desity of states C,m C est,m = W,m r = W (i),m/c est,m Thomas Prellberg, Göttige, Ja 4, 24 p.4/24

28 From PERM to flatperm PERM: estimate umber of cofiguratios C C est = W r = W (i) /C est Cosider eergy E m = ɛm, temperature β = /k B T thermal PERM: estimate partitio fuctio Z (β) Z est (β) = W exp( βe) r = W (i) exp( βe m (i) )/Z est (β) flatperm: estimate desity of states C,m C est,m = W,m r = W (i),m/c est,m Geeralizatio to more microcaoical parameters possible Thomas Prellberg, Göttige, Ja 4, 24 p.4/24

29 flatperm details Implemetatio details Parameter free implemetatio Equilibrate by slowly icreasig maximal size of cofiguratios Reduce correlatios by cosiderig effective sample size Thomas Prellberg, Göttige, Ja 4, 24 p.5/24

30 flatperm details Implemetatio details Parameter free implemetatio Equilibrate by slowly icreasig maximal size of cofiguratios Reduce correlatios by cosiderig effective sample size What if we do t eed the whole histogram? Restrictio of histogram possible e.g. impose iteral eergy restrictio m/ < u cutoff Thomas Prellberg, Göttige, Ja 4, 24 p.5/24

31 flatperm details Implemetatio details Parameter free implemetatio Equilibrate by slowly icreasig maximal size of cofiguratios Reduce correlatios by cosiderig effective sample size What if we do t eed the whole histogram? Restrictio of histogram possible e.g. impose iteral eergy restrictio m/ < u cutoff What if we do t wat a flat histogram? Sample profile ca be tued Selectively perform pruig/erichmet with respect to profile shape f,m Thomas Prellberg, Göttige, Ja 4, 24 p.5/24

32 Simulatio Results Thomas Prellberg, Göttige, Ja 4, 24 p.6/24

33 ) Simulatio results: HP model Sequece I (4 Moomers, HPHPHHPHPHHPPH, d = 3):! ' ( %& # $ "! Pedagogical example, egieered for ative groud state Perfect agreemet with exact eumeratio Thomas Prellberg, Göttige, Ja 4, 24 p.7/24

34 # Simulatio results: HP model Sequece II (85 Moomers, d = 2):! " Ivestigated several other sequeces i d = 2 ad d = 3 Collapsed regime accessible Reproduced kow groud state eergies Obtaied C,m over large rage Thomas Prellberg, Göttige, Ja 4, 24 p.8/24

35 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Equilibratio with delay.: after t tours growth up to legth t Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

36 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size:,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

37 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size:,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

38 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 2,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

39 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 3,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

40 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 4,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

41 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 5,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

42 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 6,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

43 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 7,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

44 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 8,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

45 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 9,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

46 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size:,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

47 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size:,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

48 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 2,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

49 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 3,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

50 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 4,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

51 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 5,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

52 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 6,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

53 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 7,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

54 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 8,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

55 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 9,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

56 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 2,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

57 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 2,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

58 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 22,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

59 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 23,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

60 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 24,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

61 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 25,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

62 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 26,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

63 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 27,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

64 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 28,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

65 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 29,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

66 Simulatio: Equilibratio for 2d ISAW 2d ISAW simulatio up to = 24 Total sample size: 3,, log (C m ) S m m/ m/ Thomas Prellberg, Göttige, Ja 4, 24 p.9/24

67 " ( Simulatio results: 2d ISAW e+8 e+7 e d ISAW up to = 24 Oe simulatio suffices 3 orders of magitude! & '! #%$ Thomas Prellberg, Göttige, Ja 4, 24 p.2/24

68 % Simulatio results: 3d ISAW e+7 e d ISAW up to = 24 Oe simulatio suffices 4 orders of magitude # $ "! Thomas Prellberg, Göttige, Ja 4, 24 p.2/24

69 Summary Thomas Prellberg, Göttige, Ja 4, 24 p.22/24

70 Coclusio: A Promisig New Algorithm Reviewed Stochastic Growth Algorithms for Polymers Preseted flat histogram versio of PERM Oe simulatio for complete desity of states! (the rage ca also be selectively restricted) Applicatios: HP model, ISAW Thomas Prellberg, Göttige, Ja 4, 24 p.23/24

71 Coclusio: A Promisig New Algorithm Reviewed Stochastic Growth Algorithms for Polymers Preseted flat histogram versio of PERM Oe simulatio for complete desity of states! (the rage ca also be selectively restricted) Applicatios: HP model, ISAW Outlook: applicatios to further models, e.g. Square lattice trees / brached polymers: with A. Rechitzer Coective costat estimate 5.435(5) (compare to 5.434(7) with Pivot algorithm) Two-dimesioal desity of states with oe simulatio Absorbig collapsig polymers: with A. Owczarek Exteded Domb-Joyce model: with J. Krawczyk Thomas Prellberg, Göttige, Ja 4, 24 p.23/24

72 The Ed Thomas Prellberg, Göttige, Ja 4, 24 p.24/24

Rare event sampling with stochastic growth algorithms

Rare event sampling with stochastic growth algorithms Rare event sampling with stochastic growth algorithms School of Mathematical Sciences Queen Mary, University of London CAIMS 2012 June 24-28, 2012 Topic Outline 1 Blind 2 Pruned and Enriched Flat Histogram

More information

Physica A 392 (2013) Contents lists available at SciVerse ScienceDirect. Physica A. journal homepage:

Physica A 392 (2013) Contents lists available at SciVerse ScienceDirect. Physica A. journal homepage: Physica A 392 (2013) 1602 1610 Cotets lists available at SciVerse ScieceDirect Physica A joural homepage: www.elsevier.com/locate/physa Semi-flexible iteractig self-avoidig trails o the square lattice

More information

Anomalous critical behavior in the polymer collapse transition of three-dimensional lattice trails

Anomalous critical behavior in the polymer collapse transition of three-dimensional lattice trails PHYSICAL REVIEW E 86, 023 (202) Aomalous critical behavior i the polymer collapse trasitio of three-dimesioal lattice trails Adrea Bedii * ad Aleksader L. Owczarek Departmet of Mathematics ad Statistics,

More information

Available online at ScienceDirect. Physics Procedia 57 (2014 ) 48 52

Available online at  ScienceDirect. Physics Procedia 57 (2014 ) 48 52 Available olie at www.sciecedirect.com ScieceDirect Physics Procedia 57 (214 ) 48 52 27th Aual CSP Workshops o Recet Developmets i Computer Simulatio Studies i Codesed Matter Physics, CSP 214 Applicatio

More information

Markov Decision Processes

Markov Decision Processes Markov Decisio Processes Defiitios; Statioary policies; Value improvemet algorithm, Policy improvemet algorithm, ad liear programmig for discouted cost ad average cost criteria. Markov Decisio Processes

More information

The Random Walk For Dummies

The Random Walk For Dummies The Radom Walk For Dummies Richard A Mote Abstract We look at the priciples goverig the oe-dimesioal discrete radom walk First we review five basic cocepts of probability theory The we cosider the Beroulli

More information

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B.

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B. Amme 3500 : System Dyamics ad Cotrol Root Locus Course Outlie Week Date Cotet Assigmet Notes Mar Itroductio 8 Mar Frequecy Domai Modellig 3 5 Mar Trasiet Performace ad the s-plae 4 Mar Block Diagrams Assig

More information

Formation of A Supergain Array and Its Application in Radar

Formation of A Supergain Array and Its Application in Radar Formatio of A Supergai Array ad ts Applicatio i Radar Tra Cao Quye, Do Trug Kie ad Bach Gia Duog. Research Ceter for Electroic ad Telecommuicatios, College of Techology (Coltech, Vietam atioal Uiversity,

More information

Exercises and Problems

Exercises and Problems HW Chapter 4: Oe-Dimesioal Quatum Mechaics Coceptual Questios 4.. Five. 4.4.. is idepedet of. a b c mu ( E). a b m( ev 5 ev) c m(6 ev ev) Exercises ad Problems 4.. Model: Model the electro as a particle

More information

1. pn junction under bias 2. I-Vcharacteristics

1. pn junction under bias 2. I-Vcharacteristics Lecture 10 The p Juctio (II) 1 Cotets 1. p juctio uder bias 2. I-Vcharacteristics 2 Key questios Why does the p juctio diode exhibit curret rectificatio? Why does the juctio curret i forward bias icrease

More information

Studying Traffic Flow Using Simulation Cellular Automata Model

Studying Traffic Flow Using Simulation Cellular Automata Model Iter atioal Joural of Pure ad Applied Mathematics Volume 113 No. 6 2017, 413 421 ISSN: 1311-8080 (prited versio); ISSN: 1314-3395 (o-lie versio) url: http://www.ijpam.eu ijpam.eu Studyig Traffic Flow Usig

More information

Sequential Monte Carlo Methods - A Review. Arnaud Doucet. Engineering Department, Cambridge University, UK

Sequential Monte Carlo Methods - A Review. Arnaud Doucet. Engineering Department, Cambridge University, UK Sequetial Mote Carlo Methods - A Review Araud Doucet Egieerig Departmet, Cambridge Uiversity, UK http://www-sigproc.eg.cam.ac.uk/ ad2/araud doucet.html ad2@eg.cam.ac.uk Istitut Heri Poicaré - Paris - 2

More information

ACO Comprehensive Exam 9 October 2007 Student code A. 1. Graph Theory

ACO Comprehensive Exam 9 October 2007 Student code A. 1. Graph Theory 1. Graph Theory Prove that there exist o simple plaar triagulatio T ad two distict adjacet vertices x, y V (T ) such that x ad y are the oly vertices of T of odd degree. Do ot use the Four-Color Theorem.

More information

RAINFALL PREDICTION BY WAVELET DECOMPOSITION

RAINFALL PREDICTION BY WAVELET DECOMPOSITION RAIFALL PREDICTIO BY WAVELET DECOMPOSITIO A. W. JAYAWARDEA Departmet of Civil Egieerig, The Uiversit of Hog Kog, Hog Kog, Chia P. C. XU Academ of Mathematics ad Sstem Scieces, Chiese Academ of Scieces,

More information

Statistical Noise Models and Diagnostics

Statistical Noise Models and Diagnostics L. Yaroslavsky: Advaced Image Processig Lab: A Tutorial, EUSIPCO2 LECTURE 2 Statistical oise Models ad Diagostics 2. Statistical models of radom iterfereces: (i) Additive sigal idepedet oise model: r =

More information

Lecture 11: Decision Trees

Lecture 11: Decision Trees ECE9 Sprig 7 Statistical Learig Theory Istructor: R. Nowak Lecture : Decisio Trees Miimum Complexity Pealized Fuctio Recall the basic results of the last lectures: let X ad Y deote the iput ad output spaces

More information

Simulation of Discrete Event Systems

Simulation of Discrete Event Systems Simulatio of Discrete Evet Systems Uit 9 Queueig Models Fall Witer 2014/2015 Uiv.-Prof. Dr.-Ig. Dipl.-Wirt.-Ig. Christopher M. Schlick Chair ad Istitute of Idustrial Egieerig ad Ergoomics RWTH Aache Uiversity

More information

BIRKHOFF ERGODIC THEOREM

BIRKHOFF ERGODIC THEOREM BIRKHOFF ERGODIC THEOREM Abstract. We will give a proof of the poitwise ergodic theorem, which was first proved by Birkhoff. May improvemets have bee made sice Birkhoff s orgial proof. The versio we give

More information

Random Signals and Noise Winter Semester 2017 Problem Set 12 Wiener Filter Continuation

Random Signals and Noise Winter Semester 2017 Problem Set 12 Wiener Filter Continuation Radom Sigals ad Noise Witer Semester 7 Problem Set Wieer Filter Cotiuatio Problem (Sprig, Exam A) Give is the sigal W t, which is a Gaussia white oise with expectatio zero ad power spectral desity fuctio

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science. BACKGROUND EXAM September 30, 2004.

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science. BACKGROUND EXAM September 30, 2004. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Departmet of Electrical Egieerig ad Computer Sciece 6.34 Discrete Time Sigal Processig Fall 24 BACKGROUND EXAM September 3, 24. Full Name: Note: This exam is closed

More information

Solution to problem set 2, Phys 210A

Solution to problem set 2, Phys 210A Solutio to problem set 2, Phys 20A Zhiyua Su Dated: April 27, 206 I. INFORMATION ENTROPY AND BOLTZMANN DISTRIBUTION a Usig Lagrage multipliers, we wat to maximize f p lp λe α From p f 0, we get p p lp

More information

Overview of Estimation

Overview of Estimation Topic Iferece is the problem of turig data ito kowledge, where kowledge ofte is expressed i terms of etities that are ot preset i the data per se but are preset i models that oe uses to iterpret the data.

More information

THE SI AND SIR EPIDEMICS ON GENERAL NETWORKS DAVID A L D O U S (BERKELEY)

THE SI AND SIR EPIDEMICS ON GENERAL NETWORKS DAVID A L D O U S (BERKELEY) PROBABILITY AND MATHEMATICAL STATISTICS Vol. 37, Fasc. 2 (2017), pp. 000 000 doi:10.19195/0208-4147.37.2.2 THE SI AND SIR EPIDEMICS ON GENERAL NETWORKS BY DAVID A L D O U S (BERKELEY) Dedicated to Tomasz

More information

HOMEWORK 2 SOLUTIONS

HOMEWORK 2 SOLUTIONS HOMEWORK SOLUTIONS CSE 55 RANDOMIZED AND APPROXIMATION ALGORITHMS 1. Questio 1. a) The larger the value of k is, the smaller the expected umber of days util we get all the coupos we eed. I fact if = k

More information

Basics of Probability Theory (for Theory of Computation courses)

Basics of Probability Theory (for Theory of Computation courses) Basics of Probability Theory (for Theory of Computatio courses) Oded Goldreich Departmet of Computer Sciece Weizma Istitute of Sciece Rehovot, Israel. oded.goldreich@weizma.ac.il November 24, 2008 Preface.

More information

EECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1

EECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1 EECS564 Estimatio, Filterig, ad Detectio Hwk 2 Sols. Witer 25 4. Let Z be a sigle observatio havig desity fuctio where. p (z) = (2z + ), z (a) Assumig that is a oradom parameter, fid ad plot the maximum

More information

Lecture 2: Monte Carlo Simulation

Lecture 2: Monte Carlo Simulation STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?

More information

EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY

EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA, 016 MODULE : Statistical Iferece Time allowed: Three hours Cadidates should aswer FIVE questios. All questios carry equal marks. The umber

More information

The continuous configurational Boltzmann biased direct Monte Carlo method for free energy properties of polymer chains

The continuous configurational Boltzmann biased direct Monte Carlo method for free energy properties of polymer chains The cotiuous cofiguratioal Boltzma biased direct Mote Carlo method for free eergy properties of polymer chais Jiro Sadaobu Polymer ad Materials Research Laboratories, Teiji Limited, 2- Hiode-cho, Iwakui-shi,

More information

Optimal Design of Accelerated Life Tests with Multiple Stresses

Optimal Design of Accelerated Life Tests with Multiple Stresses Optimal Desig of Accelerated Life Tests with Multiple Stresses Y. Zhu ad E. A. Elsayed Departmet of Idustrial ad Systems Egieerig Rutgers Uiversity 2009 Quality & Productivity Research Coferece IBM T.

More information

MECHANICAL INTEGRITY DESIGN FOR FLARE SYSTEM ON FLOW INDUCED VIBRATION

MECHANICAL INTEGRITY DESIGN FOR FLARE SYSTEM ON FLOW INDUCED VIBRATION MECHANICAL INTEGRITY DESIGN FOR FLARE SYSTEM ON FLOW INDUCED VIBRATION Hisao Izuchi, Pricipal Egieerig Cosultat, Egieerig Solutio Uit, ChAS Project Operatios Masato Nishiguchi, Egieerig Solutio Uit, ChAS

More information

On forward improvement iteration for stopping problems

On forward improvement iteration for stopping problems O forward improvemet iteratio for stoppig problems Mathematical Istitute, Uiversity of Kiel, Ludewig-Mey-Str. 4, D-24098 Kiel, Germay irle@math.ui-iel.de Albrecht Irle Abstract. We cosider the optimal

More information

Asymptotic Coupling and Its Applications in Information Theory

Asymptotic Coupling and Its Applications in Information Theory Asymptotic Couplig ad Its Applicatios i Iformatio Theory Vicet Y. F. Ta Joit Work with Lei Yu Departmet of Electrical ad Computer Egieerig, Departmet of Mathematics, Natioal Uiversity of Sigapore IMS-APRM

More information

Lecture 27. Capacity of additive Gaussian noise channel and the sphere packing bound

Lecture 27. Capacity of additive Gaussian noise channel and the sphere packing bound Lecture 7 Ageda for the lecture Gaussia chael with average power costraits Capacity of additive Gaussia oise chael ad the sphere packig boud 7. Additive Gaussia oise chael Up to this poit, we have bee

More information

MATH 1910 Workshop Solution

MATH 1910 Workshop Solution MATH 90 Workshop Solutio Fractals Itroductio: Fractals are atural pheomea or mathematical sets which exhibit (amog other properties) self similarity: o matter how much we zoom i, the structure remais the

More information

Advanced Stochastic Processes.

Advanced Stochastic Processes. Advaced Stochastic Processes. David Gamarik LECTURE 2 Radom variables ad measurable fuctios. Strog Law of Large Numbers (SLLN). Scary stuff cotiued... Outlie of Lecture Radom variables ad measurable fuctios.

More information

Applied Mathematics Letters. On the properties of Lucas numbers with binomial coefficients

Applied Mathematics Letters. On the properties of Lucas numbers with binomial coefficients Applied Mathematics Letters 3 (1 68 7 Cotets lists available at ScieceDirect Applied Mathematics Letters joural homepage: wwwelseviercom/locate/aml O the properties of Lucas umbers with biomial coefficiets

More information

Machine Learning Brett Bernstein

Machine Learning Brett Bernstein Machie Learig Brett Berstei Week Lecture: Cocept Check Exercises Starred problems are optioal. Statistical Learig Theory. Suppose A = Y = R ad X is some other set. Furthermore, assume P X Y is a discrete

More information

ECEN 655: Advanced Channel Coding Spring Lecture 7 02/04/14. Belief propagation is exact on tree-structured factor graphs.

ECEN 655: Advanced Channel Coding Spring Lecture 7 02/04/14. Belief propagation is exact on tree-structured factor graphs. ECEN 655: Advaced Chael Codig Sprig 014 Prof. Hery Pfister Lecture 7 0/04/14 Scribe: Megke Lia 1 4-Cycles i Gallager s Esemble What we already kow: Belief propagatio is exact o tree-structured factor graphs.

More information

Algorithms for Clustering

Algorithms for Clustering CR2: Statistical Learig & Applicatios Algorithms for Clusterig Lecturer: J. Salmo Scribe: A. Alcolei Settig: give a data set X R p where is the umber of observatio ad p is the umber of features, we wat

More information

Estimation for Complete Data

Estimation for Complete Data Estimatio for Complete Data complete data: there is o loss of iformatio durig study. complete idividual complete data= grouped data A complete idividual data is the oe i which the complete iformatio of

More information

1 Review of Probability & Statistics

1 Review of Probability & Statistics 1 Review of Probability & Statistics a. I a group of 000 people, it has bee reported that there are: 61 smokers 670 over 5 960 people who imbibe (drik alcohol) 86 smokers who imbibe 90 imbibers over 5

More information

Renormalization and disorder : a simple toy model

Renormalization and disorder : a simple toy model Reormalizatio ad disorder : a simple toy model Berard DERRIDA Collège de Frace Chaire de Physique Statistique Laboratoire de Probabilités, Statistique et Modélisatio (LPSM) Kick-o Coferece - Jue 18th to

More information

This is the story of how I discovered

This is the story of how I discovered EXPLORATIONS BY JOHN A. MAJOR, ASA Bayesia Dragos: A Cautioary Note Here Be Dragos implied by a rage of plausible model Aciet map label for uexplored parameter values. regios My challege was to get estimates

More information

An Introduction to Randomized Algorithms

An Introduction to Randomized Algorithms A Itroductio to Radomized Algorithms The focus of this lecture is to study a radomized algorithm for quick sort, aalyze it usig probabilistic recurrece relatios, ad also provide more geeral tools for aalysis

More information

ON THE DELOCALIZED PHASE OF THE RANDOM PINNING MODEL

ON THE DELOCALIZED PHASE OF THE RANDOM PINNING MODEL O THE DELOCALIZED PHASE OF THE RADOM PIIG MODEL JEA-CHRISTOPHE MOURRAT Abstract. We cosider the model of a directed polymer pied to a lie of i.i.d. radom charges, ad focus o the iterior of the delocalized

More information

The picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled

The picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled 1 Lecture : Area Area ad distace traveled Approximatig area by rectagles Summatio The area uder a parabola 1.1 Area ad distace Suppose we have the followig iformatio about the velocity of a particle, how

More information

A Genetic Algorithm for Solving General System of Equations

A Genetic Algorithm for Solving General System of Equations A Geetic Algorithm for Solvig Geeral System of Equatios Győző Molárka, Edit Miletics Departmet of Mathematics, Szécheyi Istvá Uiversity, Győr, Hugary molarka@sze.hu, miletics@sze.hu Abstract: For solvig

More information

Bohr s Atomic Model Quantum Mechanical Model

Bohr s Atomic Model Quantum Mechanical Model September 7, 0 - Summary - Itroductio to Atomic Theory Bohr s Atomic Model Quatum Mechaical Model 3- Some Defiitio 3- Projects Temperature Pressure Website Subject Areas Plasma is a Mixture of electros,

More information

Problem Set 4 Due Oct, 12

Problem Set 4 Due Oct, 12 EE226: Radom Processes i Systems Lecturer: Jea C. Walrad Problem Set 4 Due Oct, 12 Fall 06 GSI: Assae Gueye This problem set essetially reviews detectio theory ad hypothesis testig ad some basic otios

More information

SOME THEORY AND PRACTICE OF STATISTICS by Howard G. Tucker

SOME THEORY AND PRACTICE OF STATISTICS by Howard G. Tucker SOME THEORY AND PRACTICE OF STATISTICS by Howard G. Tucker CHAPTER 9. POINT ESTIMATION 9. Covergece i Probability. The bases of poit estimatio have already bee laid out i previous chapters. I chapter 5

More information

Entropy Rates and Asymptotic Equipartition

Entropy Rates and Asymptotic Equipartition Chapter 29 Etropy Rates ad Asymptotic Equipartitio Sectio 29. itroduces the etropy rate the asymptotic etropy per time-step of a stochastic process ad shows that it is well-defied; ad similarly for iformatio,

More information

Information Theory and Statistics Lecture 4: Lempel-Ziv code

Information Theory and Statistics Lecture 4: Lempel-Ziv code Iformatio Theory ad Statistics Lecture 4: Lempel-Ziv code Łukasz Dębowski ldebowsk@ipipa.waw.pl Ph. D. Programme 203/204 Etropy rate is the limitig compressio rate Theorem For a statioary process (X i)

More information

Introduction to Computational Molecular Biology. Gibbs Sampling

Introduction to Computational Molecular Biology. Gibbs Sampling 18.417 Itroductio to Computatioal Molecular Biology Lecture 19: November 16, 2004 Scribe: Tushara C. Karuarata Lecturer: Ross Lippert Editor: Tushara C. Karuarata Gibbs Samplig Itroductio Let s first recall

More information

SOLUTIONS: ECE 606 Homework Week 7 Mark Lundstrom Purdue University (revised 3/27/13) e E i E T

SOLUTIONS: ECE 606 Homework Week 7 Mark Lundstrom Purdue University (revised 3/27/13) e E i E T SOUIONS: ECE 606 Homework Week 7 Mark udstrom Purdue Uiversity (revised 3/27/13) 1) Cosider a - type semicoductor for which the oly states i the badgap are door levels (i.e. ( E = E D ). Begi with the

More information

Molecular Mechanisms of Gas Diffusion in CO 2 Hydrates

Molecular Mechanisms of Gas Diffusion in CO 2 Hydrates Supportig Iformatio Molecular Mechaisms of Gas Diffusio i CO Hydrates Shuai Liag, * Deqig Liag, Negyou Wu,,3 Lizhi Yi, ad Gaowei Hu,3 Key Laboratory of Gas Hydrate, Guagzhou Istitute of Eergy Coversio,

More information

WHEN IS THE (CO)SINE OF A RATIONAL ANGLE EQUAL TO A RATIONAL NUMBER? 1. My Motivation Some Sort of an Introduction

WHEN IS THE (CO)SINE OF A RATIONAL ANGLE EQUAL TO A RATIONAL NUMBER? 1. My Motivation Some Sort of an Introduction WHEN IS THE (CO)SINE OF A RATIONAL ANGLE EQUAL TO A RATIONAL NUMBER? JÖRG JAHNEL 1. My Motivatio Some Sort of a Itroductio Last term I taught Topological Groups at the Göttige Georg August Uiversity. This

More information

FIR Filter Design: Part I

FIR Filter Design: Part I EEL3: Discrete-Time Sigals ad Systems FIR Filter Desig: Part I. Itroductio FIR Filter Desig: Part I I this set o otes, we cotiue our exploratio o the requecy respose o FIR ilters. First, we cosider some

More information

Axioms of Measure Theory

Axioms of Measure Theory MATH 532 Axioms of Measure Theory Dr. Neal, WKU I. The Space Throughout the course, we shall let X deote a geeric o-empty set. I geeral, we shall ot assume that ay algebraic structure exists o X so that

More information

On Equivalence of Martingale Tail Bounds and Deterministic Regret Inequalities

On Equivalence of Martingale Tail Bounds and Deterministic Regret Inequalities O Equivalece of Martigale Tail Bouds ad Determiistic Regret Iequalities Sasha Rakhli Departmet of Statistics, The Wharto School Uiversity of Pesylvaia Dec 16, 2015 Joit work with K. Sridhara arxiv:1510.03925

More information

ECONOMETRIC THEORY. MODULE XIII Lecture - 34 Asymptotic Theory and Stochastic Regressors

ECONOMETRIC THEORY. MODULE XIII Lecture - 34 Asymptotic Theory and Stochastic Regressors ECONOMETRIC THEORY MODULE XIII Lecture - 34 Asymptotic Theory ad Stochastic Regressors Dr. Shalabh Departmet of Mathematics ad Statistics Idia Istitute of Techology Kapur Asymptotic theory The asymptotic

More information

AIT. Blackbody Radiation IAAT

AIT. Blackbody Radiation IAAT 3 1 Blackbody Radiatio Itroductio 3 2 First radiatio process to look at: radiatio i thermal equilibrium with itself: blackbody radiatio Assumptios: 1. Photos are Bosos, i.e., more tha oe photo per phase

More information

The Wave Function and Quantum Reality

The Wave Function and Quantum Reality The Wave Fuctio ad Quatum Reality Sha Gao Uit for History ad Philosophy of Sciece & Cetre for Time, SOPHI Uiversity of Sydey, Sydey, NSW 006, Australia Abstract. We ivestigate the meaig of the wave fuctio

More information

Mathematical Statistics - MS

Mathematical Statistics - MS Paper Specific Istructios. The examiatio is of hours duratio. There are a total of 60 questios carryig 00 marks. The etire paper is divided ito three sectios, A, B ad C. All sectios are compulsory. Questios

More information

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015 ECE 8527: Itroductio to Machie Learig ad Patter Recogitio Midterm # 1 Vaishali Ami Fall, 2015 tue39624@temple.edu Problem No. 1: Cosider a two-class discrete distributio problem: ω 1 :{[0,0], [2,0], [2,2],

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall Midterm Solutions

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall Midterm Solutions MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.65/5.070J Fall 0 Midterm Solutios Problem Suppose a radom variable X is such that P(X > ) = 0 ad P(X > E) > 0 for every E > 0. Recall that the large deviatios rate

More information

AN ALMOST LINEAR RECURRENCE. Donald E. Knuth Calif. Institute of Technology, Pasadena, Calif.

AN ALMOST LINEAR RECURRENCE. Donald E. Knuth Calif. Institute of Technology, Pasadena, Calif. AN ALMOST LINEAR RECURRENCE Doald E. Kuth Calif. Istitute of Techology, Pasadea, Calif. form A geeral liear recurrece with costat coefficiets has the U 0 = a l* U l = a 2 " ' " U r - l = a r ; u = b, u,

More information

Information-based Feature Selection

Information-based Feature Selection Iformatio-based Feature Selectio Farza Faria, Abbas Kazeroui, Afshi Babveyh Email: {faria,abbask,afshib}@staford.edu 1 Itroductio Feature selectio is a topic of great iterest i applicatios dealig with

More information

Ch3 Discrete Time Fourier Transform

Ch3 Discrete Time Fourier Transform Ch3 Discrete Time Fourier Trasform 3. Show that the DTFT of [] is give by ( k). e k 3. Determie the DTFT of the two sided sigal y [ ],. 3.3 Determie the DTFT of the causal sequece x[ ] A cos( 0 ) [ ],

More information

Quantum Simulation: Solving Schrödinger Equation on a Quantum Computer

Quantum Simulation: Solving Schrödinger Equation on a Quantum Computer Purdue Uiversity Purdue e-pubs Birc Poster Sessios Birc Naotechology Ceter 4-14-008 Quatum Simulatio: Solvig Schrödiger Equatio o a Quatum Computer Hefeg Wag Purdue Uiversity, wag10@purdue.edu Sabre Kais

More information

LONG SNAKES IN POWERS OF THE COMPLETE GRAPH WITH AN ODD NUMBER OF VERTICES

LONG SNAKES IN POWERS OF THE COMPLETE GRAPH WITH AN ODD NUMBER OF VERTICES J Lodo Math Soc (2 50, (1994, 465 476 LONG SNAKES IN POWERS OF THE COMPLETE GRAPH WITH AN ODD NUMBER OF VERTICES Jerzy Wojciechowski Abstract I [5] Abbott ad Katchalski ask if there exists a costat c >

More information

UNIVERSITY OF NORTH CAROLINA Department of Statistics Chapel Hill, N. C. FOR MULTIVARIATE POPULATIONS. J. N. Srivastava.

UNIVERSITY OF NORTH CAROLINA Department of Statistics Chapel Hill, N. C. FOR MULTIVARIATE POPULATIONS. J. N. Srivastava. UNIVERSITY OF NORTH CAROLINA Departmet of Statistics Chapel Hill N. C. A NOTE ON THE BEST LINEAR UNBIASED ESTIMATES FOR MULTIVARIATE POPULATIONS by J. N. Srivastava November 1962 Cotract No. AF 49(638)-213

More information

18.S096: Homework Problem Set 1 (revised)

18.S096: Homework Problem Set 1 (revised) 8.S096: Homework Problem Set (revised) Topics i Mathematics of Data Sciece (Fall 05) Afoso S. Badeira Due o October 6, 05 Exteded to: October 8, 05 This homework problem set is due o October 6, at the

More information

http://www.xelca.l/articles/ufo_ladigsbaa_houte.aspx imulatio Output aalysis 3/4/06 This lecture Output: A simulatio determies the value of some performace measures, e.g. productio per hour, average queue

More information

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka) 7 Phoos ad coductio electros i solids Hiroshi Matsuoa I this chapter we will discuss a miimal microscopic model for phoos i a solid ad a miimal microscopic model for coductio electros i a simple metal.

More information

PC5215 Numerical Recipes with Applications - Review Problems

PC5215 Numerical Recipes with Applications - Review Problems PC55 Numerical Recipes with Applicatios - Review Problems Give the IEEE 754 sigle precisio bit patter (biary or he format) of the followig umbers: 0 0 05 00 0 00 Note that it has 8 bits for the epoet,

More information

MIXED REVIEW of Problem Solving

MIXED REVIEW of Problem Solving MIXED REVIEW of Problem Solvig STATE TEST PRACTICE classzoe.com Lessos 2.4 2.. MULTI-STEP PROBLEM A ball is dropped from a height of 2 feet. Each time the ball hits the groud, it bouces to 70% of its previous

More information

SDS 321: Introduction to Probability and Statistics

SDS 321: Introduction to Probability and Statistics SDS 321: Itroductio to Probability ad Statistics Lecture 23: Cotiuous radom variables- Iequalities, CLT Puramrita Sarkar Departmet of Statistics ad Data Sciece The Uiversity of Texas at Austi www.cs.cmu.edu/

More information

6. Kalman filter implementation for linear algebraic equations. Karhunen-Loeve decomposition

6. Kalman filter implementation for linear algebraic equations. Karhunen-Loeve decomposition 6. Kalma filter implemetatio for liear algebraic equatios. Karhue-Loeve decompositio 6.1. Solvable liear algebraic systems. Probabilistic iterpretatio. Let A be a quadratic matrix (ot obligatory osigular.

More information

Analysis of Algorithms -Quicksort-

Analysis of Algorithms -Quicksort- Aalysis of Algorithms -- Adreas Ermedahl MRTC (Mälardales Real-Time Research Ceter) adreas.ermedahl@mdh.se Autum 2004 Proposed by C.A.R. Hoare i 962 Worst- case ruig time: Θ( 2 ) Expected ruig time: Θ(

More information

The Boolean Ring of Intervals

The Boolean Ring of Intervals MATH 532 Lebesgue Measure Dr. Neal, WKU We ow shall apply the results obtaied about outer measure to the legth measure o the real lie. Throughout, our space X will be the set of real umbers R. Whe ecessary,

More information

Lecture 11. Applying Boltzmann Statistics. Heat capacities. Elasticity of a Polymer. C V of molecules for real!! When equipartition fails

Lecture 11. Applying Boltzmann Statistics. Heat capacities. Elasticity of a Polymer. C V of molecules for real!! When equipartition fails Lecture 11 Applyig Boltzma Statistics Elasticity of a Polymer Heat capacities C V of molecules for real!! Whe equipartitio fails Readig for this Lecture: Elemets Ch 8 Readig for Lecture 12: Elemets Ch

More information

Modified Ratio Estimators Using Known Median and Co-Efficent of Kurtosis

Modified Ratio Estimators Using Known Median and Co-Efficent of Kurtosis America Joural of Mathematics ad Statistics 01, (4): 95-100 DOI: 10.593/j.ajms.01004.05 Modified Ratio s Usig Kow Media ad Co-Efficet of Kurtosis J.Subramai *, G.Kumarapadiya Departmet of Statistics, Podicherry

More information

Disjoint set (Union-Find)

Disjoint set (Union-Find) CS124 Lecture 7 Fall 2018 Disjoit set (Uio-Fid) For Kruskal s algorithm for the miimum spaig tree problem, we foud that we eeded a data structure for maitaiig a collectio of disjoit sets. That is, we eed

More information

Lecture 3: August 31

Lecture 3: August 31 36-705: Itermediate Statistics Fall 018 Lecturer: Siva Balakrisha Lecture 3: August 31 This lecture will be mostly a summary of other useful expoetial tail bouds We will ot prove ay of these i lecture,

More information

Polynomials with Rational Roots that Differ by a Non-zero Constant. Generalities

Polynomials with Rational Roots that Differ by a Non-zero Constant. Generalities Polyomials with Ratioal Roots that Differ by a No-zero Costat Philip Gibbs The problem of fidig two polyomials P(x) ad Q(x) of a give degree i a sigle variable x that have all ratioal roots ad differ by

More information

January 25, 2017 INTRODUCTION TO MATHEMATICAL STATISTICS

January 25, 2017 INTRODUCTION TO MATHEMATICAL STATISTICS Jauary 25, 207 INTRODUCTION TO MATHEMATICAL STATISTICS Abstract. A basic itroductio to statistics assumig kowledge of probability theory.. Probability I a typical udergraduate problem i probability, we

More information

Notes on iteration and Newton s method. Iteration

Notes on iteration and Newton s method. Iteration Notes o iteratio ad Newto s method Iteratio Iteratio meas doig somethig over ad over. I our cotet, a iteratio is a sequece of umbers, vectors, fuctios, etc. geerated by a iteratio rule of the type 1 f

More information

Pattern Classification

Pattern Classification Patter Classificatio All materials i these slides were tae from Patter Classificatio (d ed) by R. O. Duda, P. E. Hart ad D. G. Stor, Joh Wiley & Sos, 000 with the permissio of the authors ad the publisher

More information

Generalized Semi- Markov Processes (GSMP)

Generalized Semi- Markov Processes (GSMP) Geeralized Semi- Markov Processes (GSMP) Summary Some Defiitios Markov ad Semi-Markov Processes The Poisso Process Properties of the Poisso Process Iterarrival times Memoryless property ad the residual

More information

On Random Line Segments in the Unit Square

On Random Line Segments in the Unit Square O Radom Lie Segmets i the Uit Square Thomas A. Courtade Departmet of Electrical Egieerig Uiversity of Califoria Los Ageles, Califoria 90095 Email: tacourta@ee.ucla.edu I. INTRODUCTION Let Q = [0, 1] [0,

More information

Application to Random Graphs

Application to Random Graphs A Applicatio to Radom Graphs Brachig processes have a umber of iterestig ad importat applicatios. We shall cosider oe of the most famous of them, the Erdős-Réyi radom graph theory. 1 Defiitio A.1. Let

More information

Olli Simula T / Chapter 1 3. Olli Simula T / Chapter 1 5

Olli Simula T / Chapter 1 3. Olli Simula T / Chapter 1 5 Sigals ad Systems Sigals ad Systems Sigals are variables that carry iformatio Systemstake sigals as iputs ad produce sigals as outputs The course deals with the passage of sigals through systems T-6.4

More information

1 Inferential Methods for Correlation and Regression Analysis

1 Inferential Methods for Correlation and Regression Analysis 1 Iferetial Methods for Correlatio ad Regressio Aalysis I the chapter o Correlatio ad Regressio Aalysis tools for describig bivariate cotiuous data were itroduced. The sample Pearso Correlatio Coefficiet

More information

Roberto s Notes on Infinite Series Chapter 1: Sequences and series Section 3. Geometric series

Roberto s Notes on Infinite Series Chapter 1: Sequences and series Section 3. Geometric series Roberto s Notes o Ifiite Series Chapter 1: Sequeces ad series Sectio Geometric series What you eed to kow already: What a ifiite series is. The divergece test. What you ca le here: Everythig there is to

More information

A Note on Effi cient Conditional Simulation of Gaussian Distributions. April 2010

A Note on Effi cient Conditional Simulation of Gaussian Distributions. April 2010 A Note o Effi ciet Coditioal Simulatio of Gaussia Distributios A D D C S S, U B C, V, BC, C April 2010 A Cosider a multivariate Gaussia radom vector which ca be partitioed ito observed ad uobserved compoetswe

More information

Confidence interval for the two-parameter exponentiated Gumbel distribution based on record values

Confidence interval for the two-parameter exponentiated Gumbel distribution based on record values Iteratioal Joural of Applied Operatioal Research Vol. 4 No. 1 pp. 61-68 Witer 2014 Joural homepage: www.ijorlu.ir Cofidece iterval for the two-parameter expoetiated Gumbel distributio based o record values

More information

Accuracy of TEXTOR He-beam diagnostics

Accuracy of TEXTOR He-beam diagnostics Accuracy of TEXTOR He-beam diagostics O. Schmitz, I.L.Beigma *, L.A. Vaishtei *, A. Pospieszczyk, B. Schweer, M. Krychoviak, U. Samm ad the TEXTOR team Forschugszetrum Jülich,, Jülich, Germay *Lebedev

More information

Stochastic Simulation

Stochastic Simulation Stochastic Simulatio 1 Itroductio Readig Assigmet: Read Chapter 1 of text. We shall itroduce may of the key issues to be discussed i this course via a couple of model problems. Model Problem 1 (Jackso

More information

On Algorithm for the Minimum Spanning Trees Problem with Diameter Bounded Below

On Algorithm for the Minimum Spanning Trees Problem with Diameter Bounded Below O Algorithm for the Miimum Spaig Trees Problem with Diameter Bouded Below Edward Kh. Gimadi 1,2, Alexey M. Istomi 1, ad Ekateria Yu. Shi 2 1 Sobolev Istitute of Mathematics, 4 Acad. Koptyug aveue, 630090

More information