Monte Carlo Methods:
|
|
- Alexia Jennings
- 5 years ago
- Views:
Transcription
1 Short Course on Computational Monte Carlo Methods: Fundamentals Shuang Zhao Assistant Professor Computer Science Department University of California, Irvine Shuang Zhao 1
2 Teaching Objective Introducing Monte Carlo integration As a stochastic quadrature method Advantages/disadvantages Introducing the stochastics model of light transport in tissue Show how these simple ideas can be extended to the simulation of light transport VP Short Course 2017 Shuang Zhao 2
3 Outline 1. Monte Carlo integration A powerful numerical tool for estimating complex integrals 2. Radiative transfer The physical framework governing light propagation in translucent materials (including human tissue) 3. Monte Carlo solution to the radiative transfer equation Basically, Shuang Zhao 3
4 Shuang Zhao 4
5 Monte Carlo Integration VP Short Course 2017 Shuang Zhao 5
6 Why Monte Carlo? A powerful tool for numerically estimating complex integrals The gold standard approach to simulate light transport in tissue Advantages: Provides estimates with quantifiable uncertainty Adaptable to systems with complex geometries Shuang Zhao 6
7 Monte Carlo Integration A powerful framework for computing integrals Numerical Nondeterministic (i.e., using randomness) Scalable to high-dimensional problems Shuang Zhao 7
8 Random Variables (Discrete) random variable X Possible outcomes: x 1, x 2,, x n with probability masses p 1, p 2,, p n such that E.g., fair coin Outcomes: x 1 = head, x 2 = tail Probabilities: p 1 = p 2 = ½ Shuang Zhao 8
9 Random Variables (Continuous) random variable X Possible outcomes: with probability density function (PDF) p(x) satisfying Shuang Zhao 9
10 Strong Law of Large Numbers Let x 1, x 2,, x n be n independent observations (aka. samples) of X Sample mean Actual mean Shuang Zhao 10
11 Example: Evaluating π Unit circle S Let X be a point uniformly distributed in the square Let, then Circle area = π Square area = 4 Shuang Zhao 11
12 Example: Evaluating π Unit circle S Simple solution for computing π: Generate n samples x 1,, x n independently Circle area = π Square area = 4 Compute Live demo Shuang Zhao 12
13 Example: Evaluating π Shuang Zhao 13
14 Integral One-dimensional High-dimensional Shuang Zhao 14
15 Deterministic Integration Quadrature rule: Scales poorly with high dimensionality: Needs n m bins for a m-dimensional problem We have a high-dimensional problem! Shuang Zhao 15
16 Monte Carlo Integration: Overview Goal: Estimating Idea: Constructing random variable Such that is called an unbiased estimator of But how? Shuang Zhao 16
17 Monte Carlo Integration Let p() be any probability density function over Γ and X be a random variable with density p Let, then: To estimate : strong law of large numbers Shuang Zhao 17
18 Monte Carlo Integration Goal: to estimate Pick a probability density function p(x) Generate n independent samples: How? Evaluate for j = 1, 2,, n Return sample mean: Shuang Zhao 18
19 PDF Sampling A universal method for geniting i.i.d samples from (almost) arbitrary one-dimensional PDFs Given a 1D distribution with PDF f(x), its cumulative density function (CDF) is given by Let, then follows the distribution given by f VP Short Course 2017 Shuang Zhao 19
20 How to pick density function p()? In theory (Almost) anything In practice: Uniform distributions (almost) always work As long as the domain is bounded Choice of p() greatly affects the effectiveness (i.e., convergence rate) of the resulting estimator Shuang Zhao 20
21 Monte Carlo Integration Hello, World! Estimating Algorithm: Draw i.i.d x 1, x 2,, x n uniformly from [0, 1) Return Live demo Shuang Zhao 21
22 Monte Carlo Integration Hello, World! Estimating Shuang Zhao 22
23 Monte Carlo Integration Goal: to estimate Pick a probability density function p(x) Generate n independent samples: Evaluate for j = 1, 2,, n Return sample mean: Shuang Zhao 23
24 Central Limit Theorem Let variables with Let be a sequence of i.i.d. random for all i, then Sample mean Estimation error In case of Monte Carlo integration, we have random variable <I> with. It follows that Confidence interval The size of the confidence interval shrinks at a rate of Shuang Zhao 24
25 Convergence Rate of MC Methods Monte Carlo methods have convergence rates of Evaluating π Evaluating definite integral Shuang Zhao 25
26 Radiative Transfer VP Short Course 2017 Shuang Zhao 26
27 Recap: Radiance In Lecture 1, we learnt that radiance L(r, Ω, t) is the central quantity of interest for many biomedical problems For example: defines the fluence at r and time t Shuang Zhao 27
28 Recap: Radiance Radiance generally depends on tissue geometry & optical properties as well as the source Example: Collimated beam Tissue The rest of this lecture focuses on estimating the radiance field given these information Shuang Zhao 28
29 Radiative Transfer A physical model describing light propagation in translucent materials based on geometric optics Shuang Zhao 29
30 Translucent Materials (Materials allowing light to scatter within) Human skin Cloud Milk Shuang Zhao 30
31 Geometric Optics Describe light as rays Light travels in straight lines in homogeneous media Accurate where the wave characteristics (e.g., diffraction, interference) of light propagation are not prominent Shuang Zhao 31
32 Stochastic Model Light enters a material and scatters around before eventually leaving or being absorbed X Absorbed Scattered Translucent material (e.g., human tissue) Shuang Zhao 32
33 Radiative Transfer Equation (RTE) Governs the radiance field L inside volume Boundary condition: radiance L on the boundary Shuang Zhao 33
34 Radiative Transfer Equation (RTE) Governs the radiance field L inside volume Boundary condition: radiance L on the boundary Steady-state version (i.e., assuming ) First-order integro-differential equation Shuang Zhao 34
35 Radiative Transfer Equation (RTE) In-scattering Out-scattering & absorption Emission Differential radiance In-scattering Out-scattering & absorption Emission Shuang Zhao 35
36 Radiative Transfer Equation (RTE) Differential radiance In-scattering Out-scattering & absorption Emission Scattering coefficient: Phase function: Extinction coefficient: Source term:, a probability density for Ω Shuang Zhao 36
37 Integral Form of the RTE (IRTE) Integro-differential equation Integral equation It is desirable to rewrite the RTE as an integral equation Better suited to the stochastic model Easier to solve numerically Shuang Zhao 37
38 Integral Form of the RTE (IRTE) Transmittance In-scattering Emission Transmittance Boundary cond. Transmittance: Shuang Zhao 38
39 Monte Carlo Solution to the Radiative Transfer Equation VP Short Course 2017 Shuang Zhao 39
40 General Measurements Many biomedical problems requires numerically evaluating integrals of the form f is called the measurement function and is given by the detector Example 1: An internal detector at at this point: measuring the fluence Shuang Zhao 40
41 General Measurements Example 2: A surface detector measuring irradiance: where θ is the angle between -Ω and the surface normal at r Shuang Zhao 41
42 General Measurements General measurements can again be estimated using Monte Carlo integration By picking some density p(r, Ω), we have an unbiased estimator Next, we focus on estimating the radiance L, which is the solution of the RTE Shuang Zhao 42
43 Recap: the IRTE where and Transmittance In-scattering Emission For simpler derivations, we assume the medium to be infinite, causing the second term to vanish Shuang Zhao 43
44 Solving the IRTE using Monte Carlo where and Consider estimating L(r 0, Ω 0 ) for some fixed r 0 and Ω 0 : Transmittance In-scattering Emission MC integration Shuang Zhao 44
45 Solving the IRTE using Monte Carlo where and Transmittance In-scattering Emission F 0 contains another integral, so we apply MC integration again: Transmittance In-scattering Emission Shuang Zhao 45
46 Solving the IRTE using Monte Carlo where and Putting everything together, we get Transmittance In-scattering Emission To obtain L(r 1, Ω 1 ), we simply repeat this whole process starting with r 1 and Ω 1 Shuang Zhao 46
47 Solving the IRTE using Monte Carlo Transmittance In-scattering Emission Pseudocode: Radiance(r,Ω): draw τ from some pdf p(τ) r 1 = r τ*ω draw Ω 1 from some pdf p(ω 1 ) return T(r 1 r)/p(τ)* ( μ s (r 1 )*p(r 1,Ω 1 Ω)*Radiance(r 1,Ω 1 ) + Q(r 1,Ω) ) Shuang Zhao 47
48 Solving the IRTE using Monte Carlo Transmittance In-scattering Emission This process effectively constructs a photon path or biography: (r 0, Ω 0 ), (r 1, Ω 1 ), (r 2, Ω 2 ), This construction starts from r 0 on the detector. There exists an adjoint version that starts from the source. Shuang Zhao 48
49 Solving the IRTE using Monte Carlo Transmittance In-scattering Emission Pseudocode: Radiance(r,Ω): draw τ from some pdf p(τ) Free distance sampling r 1 = r τ*ω draw Ω 1 from some pdf p(ω 1 ) Incident dir. sampling return T(r 1 r 0 )/p(τ)* ( μ s (r 1 )*p(r 1,Ω 1 Ω)*Radiance(r 1,Ω 1 ) + Q(r 1,Ω) ) Shuang Zhao 49
50 Sampling Incident Direction Transmittance In-scattering Emission Transmittance In-scattering Emission To draw the incident direction, we simply follows the stochastic model by setting p(ω 1 ) based on the phase function at r 1 and Ω 0 This cancels out the phase function term in the nominator Shuang Zhao 50
51 Sampling Free Distance Transmittance In-scattering Emission To sample the free-flight distance τ, it is desirable to cancel out the transmittance term This can be done by setting, yielding In-scattering Emission Shuang Zhao 51
52 Sampling Free Distance For homogeneous media, we have and We can draw τ from this exponential distribution using the PDF sampling method: The CDF is Let, then Shuang Zhao 52
53 Stochastic Model vs. Monte Carlo The Monte Carlo solution to the RTE is equivalent to the stochastic model captures the likelihood for light to collide with some particle in the medium follows the phase function Scattered Shuang Zhao 53 This Biophotonics is NOT 2017 the only solution: other probabilities can be used if necessary
54 Absorption Weighting In-scattering Emission During the construction of the photon biography, the weight of a photon decreases at each collision due to the scaling factor of is usually called the single-scattering albedo Two interpretations: at each collision, The photon gets absorbed with probability (1 α) The photon s weight scales down with a factor of α Shuang Zhao 54
55 Analog In-scattering Emission Pseudocode: Radiance(r,Ω): draw τ from r 1 = r τ*ω draw Ω 1 from rad = Q(r 1,Ω)/μ t (r 1 ) if rand() < μ s (r 1 )/μ t (r 1 ): rad += Radiance(r 1,Ω 1 ) return rad Shuang Zhao 55
56 Discrete Absorption Weighting In-scattering Emission Pseudocode: Radiance(r,Ω): draw τ from r 1 = r τ*ω draw Ω 1 from rad = Q(r 1,Ω)/μ t (r 1 ) rad += μ s (r 1 )/μ t (r 1 )*Radiance(r 1,Ω 1 ) return rad Shuang Zhao 56
57 Analog vs. DAW Analog is faster DAW offers cleaner results Shuang Zhao 57
58 Monte Carlo for Inverse Problem Monte Carlo method can also be applied to solve the inverse radiative transfer problem i.e., given some measurement I, need to solve for μ t, μ s, and p The inverse problem is generally more challenging and modeled as an optimization where I 0 is the provided measurement and I is that given by scattering parameters μ t, μ s, and p Also need to determine which norm ( ) to use Shuang Zhao 58
59 Monte Carlo integration is also used to render photorealistic images Shuang Zhao 59
60 Colored Explosion [Kutz et al. 2017] Shuang Zhao 60
61 Soap, Olive Oil, Curacao, and Milk [Gkioulekas et al. 2013] Shuang Zhao 61
62 Take Home Message Radiative transfer is an accurate model for light propagation in translucent materials when wave effects and not prominent Monte Carlo solution to the RTE provides a gold standard for simulating light propagation in tissue The flexibility of the Monte Carlo framework allows future research toward more efficient estimators Shuang Zhao 62
Monte Carlo Integration II & Sampling from PDFs
Monte Carlo Integration II & Sampling from PDFs CS295, Spring 2017 Shuang Zhao Computer Science Department University of California, Irvine CS295, Spring 2017 Shuang Zhao 1 Last Lecture Direct illumination
More informationSTOCHASTIC & DETERMINISTIC SOLVERS
STOCHASTIC & DETERMINISTIC SOLVERS Outline Spatial Scales of Optical Technologies & Mathematical Models Prototype RTE Problems 1-D Transport in Slab Geometry: Exact Solution Stochastic Models: Monte Carlo
More informationMonte Carlo Radiation Transfer I
Monte Carlo Radiation Transfer I Monte Carlo Photons and interactions Sampling from probability distributions Optical depths, isotropic emission, scattering Monte Carlo Basics Emit energy packet, hereafter
More informationHigh-Order Similarity Relations in Radiative Transfer: Supplementary Document
High-Order Similarity Relations in Radiative Transfer: Supplementary Document Shuang hao Ravi Ramamoorthi Kavita Bala Cornell University University of California, Berkeley Cornell University Overview This
More informationMonte Carlo Integration. Computer Graphics CMU /15-662, Fall 2016
Monte Carlo Integration Computer Graphics CMU 15-462/15-662, Fall 2016 Talk Announcement Jovan Popovic, Senior Principal Scientist at Adobe Research will be giving a seminar on Character Animator -- Monday
More informationMCRT: L4 A Monte Carlo Scattering Code
MCRT: L4 A Monte Carlo Scattering Code Plane parallel scattering slab Optical depths & physical distances Emergent flux & intensity Internal intensity moments Constant density slab, vertical optical depth
More informationDetectors in Nuclear Physics: Monte Carlo Methods. Dr. Andrea Mairani. Lectures I-II
Detectors in Nuclear Physics: Monte Carlo Methods Dr. Andrea Mairani Lectures I-II INTRODUCTION Sampling from a probability distribution Sampling from a probability distribution X λ Sampling from a probability
More informationEnrico Fermi and the FERMIAC. Mechanical device that plots 2D random walks of slow and fast neutrons in fissile material
Monte Carlo History Statistical sampling Buffon s needles and estimates of π 1940s: neutron transport in fissile material Origin of name People: Ulam, von Neuman, Metropolis, Teller Other areas of use:
More informationAnalysis of Scattering of Radiation in a Plane-Parallel Atmosphere. Stephanie M. Carney ES 299r May 23, 2007
Analysis of Scattering of Radiation in a Plane-Parallel Atmosphere Stephanie M. Carney ES 299r May 23, 27 TABLE OF CONTENTS. INTRODUCTION... 2. DEFINITION OF PHYSICAL QUANTITIES... 3. DERIVATION OF EQUATION
More informationWeighting Functions and Atmospheric Soundings: Part I
Weighting Functions and Atmospheric Soundings: Part I Ralf Bennartz Cooperative Institute for Meteorological Satellite Studies University of Wisconsin Madison Outline What we want to know and why we need
More informationMonte Carlo Radiation Transport Kenny Wood
MCRT: L0 Some background, what previous courses students should look over, gentle introduction/recap of probabilities Get an idea of computer programming experience of the class Overview of course structure
More informationMCRT L10: Scattering and clarification of astronomy/medical terminology
MCRT L10: Scattering and clarification of astronomy/medical terminology What does the scattering? Shape of scattering Sampling from scattering phase functions Co-ordinate frames Refractive index changes
More informationChapter V: Cavity theories
Chapter V: Cavity theories 1 Introduction Goal of radiation dosimetry: measure of the dose absorbed inside a medium (often assimilated to water in calculations) A detector (dosimeter) never measures directly
More informationMultiple-source optical diffusion approximation for a multilayer scattering medium
Multiple-source optical diffusion approximation for a multilayer scattering medium Joseph. Hollmann 1 and ihong V. Wang 1,2, * 1 Optical Imaging aboratory, Department of Biomedical Engineering, Texas A&M
More informationULTRAFAST LASER PULSE TRAIN RADIATION TRANSFER IN A SCATTERING-ABSORBING 3D MEDIUM WITH AN INHOMOGENEITY
Heat Transfer Research 46(9), 861 879 (2015) ULTRAFAST LASER PULSE TRAIN RADIATION TRANSFER IN A SCATTERING-ABSORBING 3D MEDIUM WITH AN INHOMOGENEITY Masato Akamatsu 1,* & Zhixiong Guo 2 1 Graduate School
More informationScientific Computing: Monte Carlo
Scientific Computing: Monte Carlo Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Spring 2012 April 5th and 12th, 2012 A. Donev (Courant Institute)
More informationMonte Carlo Radiation Transport Kenny Wood
Monte Carlo Radiation Transport Kenny Wood kw25@st-andrews.ac.uk A practical approach to the numerical simulation of radiation transport Develop programs for the random walks of photons and neutrons using
More informationStatistics 100A Homework 1 Solutions
Problem Statistics 00A Homework Solutions Ryan Rosario Suppose we flip a fair coin 4 times independently. () What is the sample space? By definition, the sample space, denoted as Ω, is the set of all possible
More informationRadiative Transfer Multiple scattering: two stream approach 2
Radiative Transfer Multiple scattering: two stream approach 2 N. Kämpfer non Institute of Applied Physics University of Bern 28. Oct. 24 Outline non non Interpretation of some specific cases Semi-infinite
More informationI ν. di ν. = α ν. = (ndads) σ ν da α ν. = nσ ν = ρκ ν
Absorption Consider a beam passing through an absorbing medium. Define the absorption coefficient, α ν, by ie the fractional loss in intensity in travelling a distance ds is α ν ds (convention: positive
More information2 The Radiative Transfer Equation
9 The Radiative Transfer Equation. Radiative transfer without absorption and scattering Free space or homogeneous space I (r,,) I (r,,) r -r d da da Figure.: Following a pencil of radiation in free space
More informationFundametals of Rendering - Radiometry / Photometry
Fundametals of Rendering - Radiometry / Photometry Physically Based Rendering by Pharr & Humphreys Chapter 5: Color and Radiometry Chapter 6: Camera Models - we won t cover this in class Realistic Rendering
More informationLichtausbreitung in streuenden Medien: Prinzip und Anwendungsbeispiele
Lichtausbreitung in streuenden Medien: Prinzip und Anwendungsbeispiele Alwin Kienle 06.12.2013 Institut für Lasertechnologien in der Medizin und Meßtechnik an der Universität Ulm Overview 1) Theory of
More informationInfluence of the grain shape on the albedo and light extinction in snow
Influence of the grain shape on the albedo and light extinction in snow Q. Libois 1, G. Picard 1, L. Arnaud 1, M. Dumont 2, J. France 3, C. Carmagnola 2, S. Morin 2, and M. King 3 1 Laboratoire de Glaciologie
More informationNeutron Interactions Part I. Rebecca M. Howell, Ph.D. Radiation Physics Y2.5321
Neutron Interactions Part I Rebecca M. Howell, Ph.D. Radiation Physics rhowell@mdanderson.org Y2.5321 Why do we as Medical Physicists care about neutrons? Neutrons in Radiation Therapy Neutron Therapy
More informationLecture 4* Inherent optical properties, IOP Theory. Loss due to absorption. IOP Theory 12/2/2008
Lecture 4* Inherent optical properties, part I IOP Theory What properties of a medium affect the radiance field as it propagate through it? 1. Sinks of photons (absorbers) 2. Sources of photons (internal
More informationCMPSCI 240: Reasoning about Uncertainty
CMPSCI 240: Reasoning about Uncertainty Lecture 2: Sets and Events Andrew McGregor University of Massachusetts Last Compiled: January 27, 2017 Outline 1 Recap 2 Experiments and Events 3 Probabilistic Models
More informationClassical and Improved Diffusion Theory for Subsurface Scattering
Classical and Improved Diffusion Theory for Subsurface Scattering Ralf Habel 1 Per H. Christensen 2 Wojciech Jarosz 1 1 Disney Research Zürich 2 Pixar Animation Studios Technical report, June 2013 Abstract
More informationStatistical Methods in Particle Physics. Lecture 2
Statistical Methods in Particle Physics Lecture 2 October 17, 2011 Silvia Masciocchi, GSI Darmstadt s.masciocchi@gsi.de Winter Semester 2011 / 12 Outline Probability Definition and interpretation Kolmogorov's
More informationMCSHAPE: A MONTE CARLO CODE FOR SIMULATION OF POLARIZED PHOTON TRANSPORT
Copyright JCPDS - International Centre for Diffraction Data 2003, Advances in X-ray Analysis, Volume 46. 363 MCSHAPE: A MONTE CARLO CODE FOR SIMULATION OF POLARIZED PHOTON TRANSPORT J.E. Fernández, V.
More informationFast simulation of finite-beam optical coherence tomography of inhomogeneous turbid media
Fast simulation of finite-beam optical coherence tomography of inhomogeneous turbid media by Micheal Sobhy A thesis submitted to the Faculty of Graduate Studies of the University of Manitoba, in partial
More informationRandom variables. DS GA 1002 Probability and Statistics for Data Science.
Random variables DS GA 1002 Probability and Statistics for Data Science http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall17 Carlos Fernandez-Granda Motivation Random variables model numerical quantities
More informationModern Methods of Data Analysis - WS 07/08
Modern Methods of Data Analysis Lecture VII (26.11.07) Contents: Maximum Likelihood (II) Exercise: Quality of Estimators Assume hight of students is Gaussian distributed. You measure the size of N students.
More informationThe mathematics of scattering and absorption and emission
The mathematics of scattering and absorption and emission The transmittance of an layer depends on its optical depth, which in turn depends on how much of the substance the radiation has to pass through,
More informationPhysics for Scientists & Engineers 2
Light as Waves Physics for Scientists & Engineers 2 Spring Semester 2005 Lecture 41! In the previous chapter we discussed light as rays! These rays traveled in a straight line except when they were reflected
More informationMAT 271E Probability and Statistics
MAT 71E Probability and Statistics Spring 013 Instructor : Class Meets : Office Hours : Textbook : Supp. Text : İlker Bayram EEB 1103 ibayram@itu.edu.tr 13.30 1.30, Wednesday EEB 5303 10.00 1.00, Wednesday
More informationMath 105 Course Outline
Math 105 Course Outline Week 9 Overview This week we give a very brief introduction to random variables and probability theory. Most observable phenomena have at least some element of randomness associated
More information6. MC data analysis. Doing experiments or complex simulations can be quite time-consuming, expensive, or both.
6. MC data analysis [Numerical Recipes 15.6, own derivation] Monte Carlo simulations can be a highly useful tool in the analysis of experimental data, or data produced by other computer simulations. Doing
More informationWaves & Oscillations
Physics 42200 Waves & Oscillations Lecture 32 Electromagnetic Waves Spring 2016 Semester Matthew Jones Electromagnetism Geometric optics overlooks the wave nature of light. Light inconsistent with longitudinal
More information1: PROBABILITY REVIEW
1: PROBABILITY REVIEW Marek Rutkowski School of Mathematics and Statistics University of Sydney Semester 2, 2016 M. Rutkowski (USydney) Slides 1: Probability Review 1 / 56 Outline We will review the following
More information1. The most important aspects of the quantum theory.
Lecture 5. Radiation and energy. Objectives: 1. The most important aspects of the quantum theory: atom, subatomic particles, atomic number, mass number, atomic mass, isotopes, simplified atomic diagrams,
More informationElectricity & Optics
Physics 24100 Electricity & Optics Lecture 26 Chapter 33 sec. 1-4 Fall 2017 Semester Professor Koltick Interference of Light Interference phenomena are a consequence of the wave-like nature of light Electric
More information6. MC data analysis. [Numerical Recipes 15.6, own derivation] Basics of Monte Carlo simulations, Kai Nordlund
6. MC data analysis [Numerical Recipes 15.6, own derivation] Monte Carlo simulations can be a highly useful tool in the analysis of experimental data, or data produced by other computer simulations. The
More information16 : Markov Chain Monte Carlo (MCMC)
10-708: Probabilistic Graphical Models 10-708, Spring 2014 16 : Markov Chain Monte Carlo MCMC Lecturer: Matthew Gormley Scribes: Yining Wang, Renato Negrinho 1 Sampling from low-dimensional distributions
More information1.225J J (ESD 205) Transportation Flow Systems
1.225J J (ESD 25) Transportation Flow Systems Lecture 9 Simulation Models Prof. Ismail Chabini and Prof. Amedeo R. Odoni Lecture 9 Outline About this lecture: It is based on R16. Only material covered
More informationCHARGED PARTICLE INTERACTIONS
CHARGED PARTICLE INTERACTIONS Background Charged Particles Heavy charged particles Charged particles with Mass > m e α, proton, deuteron, heavy ion (e.g., C +, Fe + ), fission fragment, muon, etc. α is
More informationMidterm Exam 1 Solution
EECS 126 Probability and Random Processes University of California, Berkeley: Fall 2015 Kannan Ramchandran September 22, 2015 Midterm Exam 1 Solution Last name First name SID Name of student on your left:
More informationLecture 2: From Classical to Quantum Model of Computation
CS 880: Quantum Information Processing 9/7/10 Lecture : From Classical to Quantum Model of Computation Instructor: Dieter van Melkebeek Scribe: Tyson Williams Last class we introduced two models for deterministic
More informationDISCRETE RANDOM VARIABLES: PMF s & CDF s [DEVORE 3.2]
DISCRETE RANDOM VARIABLES: PMF s & CDF s [DEVORE 3.2] PROBABILITY MASS FUNCTION (PMF) DEFINITION): Let X be a discrete random variable. Then, its pmf, denoted as p X(k), is defined as follows: p X(k) :=
More informationRandom Processes. DS GA 1002 Probability and Statistics for Data Science.
Random Processes DS GA 1002 Probability and Statistics for Data Science http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall17 Carlos Fernandez-Granda Aim Modeling quantities that evolve in time (or space)
More informationProbability Theory Review
Cogsci 118A: Natural Computation I Lecture 2 (01/07/10) Lecturer: Angela Yu Probability Theory Review Scribe: Joseph Schilz Lecture Summary 1. Set theory: terms and operators In this section, we provide
More informationTutorial 1 : Probabilities
Lund University ETSN01 Advanced Telecommunication Tutorial 1 : Probabilities Author: Antonio Franco Emma Fitzgerald Tutor: Farnaz Moradi January 11, 2016 Contents I Before you start 3 II Exercises 3 1
More informationMath 151. Rumbos Fall Solutions to Review Problems for Exam 2. Pr(X = 1) = ) = Pr(X = 2) = Pr(X = 3) = p X. (k) =
Math 5. Rumbos Fall 07 Solutions to Review Problems for Exam. A bowl contains 5 chips of the same size and shape. Two chips are red and the other three are blue. Draw three chips from the bowl at random,
More informationNumerical Methods I Monte Carlo Methods
Numerical Methods I Monte Carlo Methods Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Dec. 9th, 2010 A. Donev (Courant Institute) Lecture
More informationRobert Collins CSE586, PSU Intro to Sampling Methods
Robert Collins Intro to Sampling Methods CSE586 Computer Vision II Penn State Univ Robert Collins A Brief Overview of Sampling Monte Carlo Integration Sampling and Expected Values Inverse Transform Sampling
More informationFrequency domain photon migration in the -P 1 approximation: Analysis of ballistic, transport, and diffuse regimes
Frequency domain photon migration in the -P 1 approximation: Analysis of ballistic, transport, and diffuse regimes J. S. You, 1,2 C. K. Hayakawa, 2 and V. Venugopalan 1,2,3, * 1 Department of Biomedical
More informationE X A M. Probability Theory and Stochastic Processes Date: December 13, 2016 Duration: 4 hours. Number of pages incl.
E X A M Course code: Course name: Number of pages incl. front page: 6 MA430-G Probability Theory and Stochastic Processes Date: December 13, 2016 Duration: 4 hours Resources allowed: Notes: Pocket calculator,
More informationParticle Interactions in Detectors
Particle Interactions in Detectors Dr Peter R Hobson C.Phys M.Inst.P. Department of Electronic and Computer Engineering Brunel University, Uxbridge Peter.Hobson@brunel.ac.uk http://www.brunel.ac.uk/~eestprh/
More informationBrandon C. Kelly (Harvard Smithsonian Center for Astrophysics)
Brandon C. Kelly (Harvard Smithsonian Center for Astrophysics) Probability quantifies randomness and uncertainty How do I estimate the normalization and logarithmic slope of a X ray continuum, assuming
More informationLecture 2 Solutions to the Transport Equation
Lecture 2 Solutions to the Transport Equation Equation along a ray I In general we can solve the static transfer equation along a ray in some particular direction. Since photons move in straight lines
More informationLecture 26. Regional radiative effects due to anthropogenic aerosols. Part 2. Haze and visibility.
Lecture 26. Regional radiative effects due to anthropogenic aerosols. Part 2. Haze and visibility. Objectives: 1. Attenuation of atmospheric radiation by particulates. 2. Haze and Visibility. Readings:
More informationLecture 3. Discrete Random Variables
Math 408 - Mathematical Statistics Lecture 3. Discrete Random Variables January 23, 2013 Konstantin Zuev (USC) Math 408, Lecture 3 January 23, 2013 1 / 14 Agenda Random Variable: Motivation and Definition
More informationIEOR 3106: Introduction to Operations Research: Stochastic Models. Professor Whitt. SOLUTIONS to Homework Assignment 1
IEOR 3106: Introduction to Operations Research: Stochastic Models Professor Whitt SOLUTIONS to Homework Assignment 1 Probability Review: Read Chapters 1 and 2 in the textbook, Introduction to Probability
More informationStatistics: Learning models from data
DS-GA 1002 Lecture notes 5 October 19, 2015 Statistics: Learning models from data Learning models from data that are assumed to be generated probabilistically from a certain unknown distribution is a crucial
More information= nm. = nm. = nm
Chemistry 60 Analytical Spectroscopy PROBLEM SET 5: Due 03/0/08 1. At a recent birthday party, a young friend (elementary school) noticed that multicolored rings form across the surface of soap bubbles.
More informationRobert Collins CSE586, PSU Intro to Sampling Methods
Intro to Sampling Methods CSE586 Computer Vision II Penn State Univ Topics to be Covered Monte Carlo Integration Sampling and Expected Values Inverse Transform Sampling (CDF) Ancestral Sampling Rejection
More informationProperties of Continuous Probability Distributions The graph of a continuous probability distribution is a curve. Probability is represented by area
Properties of Continuous Probability Distributions The graph of a continuous probability distribution is a curve. Probability is represented by area under the curve. The curve is called the probability
More informationThe atom cont. +Investigating EM radiation
The atom cont. +Investigating EM radiation Announcements: First midterm is 7:30pm on Sept 26, 2013 Will post a past midterm exam from 2011 today. We are covering Chapter 3 today. (Started on Wednesday)
More informationCourse: ESO-209 Home Work: 1 Instructor: Debasis Kundu
Home Work: 1 1. Describe the sample space when a coin is tossed (a) once, (b) three times, (c) n times, (d) an infinite number of times. 2. A coin is tossed until for the first time the same result appear
More informationElectromagnetic spectra
Properties of Light Waves, particles and EM spectrum Interaction with matter Absorption Reflection, refraction and scattering Polarization and diffraction Reading foci: pp 175-185, 191-199 not responsible
More informationMODERN OPTICS. P47 Optics: Unit 9
MODERN OPTICS P47 Optics: Unit 9 Course Outline Unit 1: Electromagnetic Waves Unit 2: Interaction with Matter Unit 3: Geometric Optics Unit 4: Superposition of Waves Unit 5: Polarization Unit 6: Interference
More informationIGD-TP Exchange Forum n 5 WG1 Safety Case: Handling of uncertainties October th 2014, Kalmar, Sweden
IGD-TP Exchange Forum n 5 WG1 Safety Case: Handling of uncertainties October 28-30 th 2014, Kalmar, Sweden Comparison of probabilistic and alternative evidence theoretical methods for the handling of parameter
More informationThe Radiative Transfer Equation
The Radiative Transfer Equation R. Wordsworth April 11, 215 1 Objectives Derive the general RTE equation Derive the atmospheric 1D horizontally homogenous RTE equation Look at heating/cooling rates in
More informationnoise = function whose amplitude is is derived from a random or a stochastic process (i.e., not deterministic)
SIMG-716 Linear Imaging Mathematics I, Handout 05 1 1-D STOCHASTIC FUCTIOS OISE noise = function whose amplitude is is derived from a random or a stochastic process (i.e., not deterministic) Deterministic:
More informationMAT 271E Probability and Statistics
MAT 7E Probability and Statistics Spring 6 Instructor : Class Meets : Office Hours : Textbook : İlker Bayram EEB 3 ibayram@itu.edu.tr 3.3 6.3, Wednesday EEB 6.., Monday D. B. Bertsekas, J. N. Tsitsiklis,
More informationNorthwestern University Department of Electrical Engineering and Computer Science
Northwestern University Department of Electrical Engineering and Computer Science EECS 454: Modeling and Analysis of Communication Networks Spring 2008 Probability Review As discussed in Lecture 1, probability
More informationik () uk () Today s menu Last lecture Some definitions Repeatability of sensing elements
Last lecture Overview of the elements of measurement systems. Sensing elements. Signal conditioning elements. Signal processing elements. Data presentation elements. Static characteristics of measurement
More informationSuper-exponential extinction of radiation in a negatively-correlated random medium
Super-exponential extinction of radiation in a negatively-correlated random medium Raymond A. Shaw 1 Alexander B. Kostinski Daniel D. Lanterman 2 Department of Physics, Michigan Technological University,
More informationDiffuse Optical Spectroscopy: Analysis of BiophotonicsSignals Short Course in Computational Biophotonics Albert Cerussi Day 4/Lecture #2
Diffuse Optical Spectroscopy: Analysis of BiophotonicsSignals 2013 Short Course in Computational Biophotonics Albert Cerussi Day 4/Lecture #2 Disclosures No financial conflicts to disclose Learning Objectives
More informationRobert Collins CSE586, PSU Intro to Sampling Methods
Intro to Sampling Methods CSE586 Computer Vision II Penn State Univ Topics to be Covered Monte Carlo Integration Sampling and Expected Values Inverse Transform Sampling (CDF) Ancestral Sampling Rejection
More informationParametric Models. Dr. Shuang LIANG. School of Software Engineering TongJi University Fall, 2012
Parametric Models Dr. Shuang LIANG School of Software Engineering TongJi University Fall, 2012 Today s Topics Maximum Likelihood Estimation Bayesian Density Estimation Today s Topics Maximum Likelihood
More informationPhysics 2D Lecture Slides Week of May 11,2009. Sunil Sinha UCSD Physics
Physics 2D Lecture Slides Week of May 11,2009 Sunil Sinha UCSD Physics Recap!! Wave Packet : Localization Finite # of diff. Monochromatic waves always produce INFINTE sequence of repeating wave groups
More informationThe Interaction of Light and Matter: α and n
The Interaction of Light and Matter: α and n The interaction of light and matter is what makes life interesting. Everything we see is the result of this interaction. Why is light absorbed or transmitted
More informationLecture 2: Repetition of probability theory and statistics
Algorithms for Uncertainty Quantification SS8, IN2345 Tobias Neckel Scientific Computing in Computer Science TUM Lecture 2: Repetition of probability theory and statistics Concept of Building Block: Prerequisites:
More informationBrief Introduction to: Transport Theory/Monte Carlo Techniques/MCNP
22:54 Neutron Interactions and Applications Brief Introduction to: Transport Theory/Monte Carlo Techniques/MCNP The behavior of individual neutrons and nuclei cannot be predicted. However, the average
More informationME 476 Solar Energy UNIT TWO THERMAL RADIATION
ME 476 Solar Energy UNIT TWO THERMAL RADIATION Unit Outline 2 Electromagnetic radiation Thermal radiation Blackbody radiation Radiation emitted from a real surface Irradiance Kirchhoff s Law Diffuse and
More informationName: Firas Rassoul-Agha
Midterm 1 - Math 5010 - Spring 016 Name: Firas Rassoul-Agha Solve the following 4 problems. You have to clearly explain your solution. The answer carries no points. Only the work does. CALCULATORS ARE
More informationMonte Carlo Methods in High Energy Physics I
Helmholtz International Workshop -- CALC 2009, July 10--20, Dubna Monte Carlo Methods in High Energy Physics CALC2009 - July 20 10, Dubna 2 Contents 3 Introduction Simple definition: A Monte Carlo technique
More informationStatistical Methods for Particle Physics Lecture 1: parameter estimation, statistical tests
Statistical Methods for Particle Physics Lecture 1: parameter estimation, statistical tests http://benasque.org/2018tae/cgi-bin/talks/allprint.pl TAE 2018 Benasque, Spain 3-15 Sept 2018 Glen Cowan Physics
More informationUniversity of Cyprus. Reflectance and Diffuse Spectroscopy
University of Cyprus Biomedical Imaging and Applied Optics Reflectance and Diffuse Spectroscopy Spectroscopy What is it? from the Greek: spectro = color + scope = look at or observe = measuring/recording
More information1. Waves and Particles 2. Interference of Waves 3. Wave Nature of Light
1. Waves and Particles 2. Interference of Waves 3. Wave Nature of Light 1. Double-Slit Eperiment reading: Chapter 22 2. Single-Slit Diffraction reading: Chapter 22 3. Diffraction Grating reading: Chapter
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 1 (2/3/04) Overview -- Interactions, Distributions, Cross Sections, Applications
.54 Neutron Interactions and Applications (Spring 004) Chapter 1 (/3/04) Overview -- Interactions, Distributions, Cross Sections, Applications There are many references in the vast literature on nuclear
More informationModeling Uncertainty in the Earth Sciences Jef Caers Stanford University
Probability theory and statistical analysis: a review Modeling Uncertainty in the Earth Sciences Jef Caers Stanford University Concepts assumed known Histograms, mean, median, spread, quantiles Probability,
More informationLecture 2. October 21, Department of Biostatistics Johns Hopkins Bloomberg School of Public Health Johns Hopkins University.
Lecture 2 Department of Biostatistics Johns Hopkins Bloomberg School of Public Health Johns Hopkins University October 21, 2007 1 2 3 4 5 6 Define probability calculus Basic axioms of probability Define
More informationModule 5 : Plane Waves at Media Interface. Lecture 36 : Reflection & Refraction from Dielectric Interface (Contd.) Objectives
Objectives In this course you will learn the following Reflection and Refraction with Parallel Polarization. Reflection and Refraction for Normal Incidence. Lossy Media Interface. Reflection and Refraction
More informationProbability measures A probability measure, P, is a real valued function from the collection of possible events so that the following
This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike License. Your use of this material constitutes acceptance of that license and the conditions of use of materials on this
More informationOptics.
Optics www.optics.rochester.edu/classes/opt100/opt100page.html Course outline Light is a Ray (Geometrical Optics) 1. Nature of light 2. Production and measurement of light 3. Geometrical optics 4. Matrix
More informationRecursive Estimation
Recursive Estimation Raffaello D Andrea Spring 08 Problem Set : Bayes Theorem and Bayesian Tracking Last updated: March, 08 Notes: Notation: Unless otherwise noted, x, y, and z denote random variables,
More informationDeep Learning for Computer Vision
Deep Learning for Computer Vision Lecture 3: Probability, Bayes Theorem, and Bayes Classification Peter Belhumeur Computer Science Columbia University Probability Should you play this game? Game: A fair
More informationMAT 271E Probability and Statistics
MAT 271E Probability and Statistics Spring 2011 Instructor : Class Meets : Office Hours : Textbook : Supp. Text : İlker Bayram EEB 1103 ibayram@itu.edu.tr 13.30 16.30, Wednesday EEB? 10.00 12.00, Wednesday
More information