Simulation. Monte Carlo
|
|
- Gyles Marsh
- 5 years ago
- Views:
Transcription
1 Simulation Monte Carlo
2 Monte Carlo simulation Outcome of a single stochastic simulation run is always random A single instance of a random variable Goal of a simulation experiment is to get knowledge about the distribution of the random variable (mean, variance) Right value is deterministic but it can not be determined explicitly.
3 Buffons needle Classical example of simulation experiment where exact result is known. Count of Buffon presented a method to define the value for p in Throw a needle of length l on a plane that has parallel lines with distance d. Count how often the needle crosses a line. Compute experimental probability for hits P= #hits/#trials
4 Buffons needle Needle hits a line if The distance from the center of the needle to closest line is less than l sin a, where a is the angle between the needle and the line d l a Angle ~ Unif(0, p/2) Center ~ Unif (0,d/2)
5 Buffons needle Probability of hit can be computed using the volume of the area bounded by the sinusoidal curve. p= 2l/(pd) Hence estimate is p = 2l/(pd) for experimental value for p. d/2 l/2 p/2
6 Buffons needle Result of a single throw is random So is the average of N throws. What do we know after N throws? Can we define the distribution of the average P after N throws. Or at least expectation and variance P is an average of N independent random variables Single attempts obey binomial distribution with mean p (=2l/(pd)) E(P)=p.
7 Buffons needle Variance of the result in one throw is p(1-p) (result is a Bin(p) variable) Variance of the average of N independent trials is p(1-p)/n I.e. Var(P) = p(1-p)/n So now we have observation of a random variable with known variance. We can estimate the relationship between the sample average and the expectation.
8 Confidence interval Assume we know a sample average of a random variable Where is the true expectation with, say 99% probability. Define sc. confidence interval for which P( P-d < p< P+d) >0.99. Can be defined if distribution of P is known. P is average of N independent Bin variables. For large N, P is approximately normally distributed. d is of form c(p)n^(-1/2).
9 Monte Carlo -integration Buffon s needle was sampling a variable expectation has a formula as definite integral. The approach can be used generally to approximate integrals. Integrate f on [a,b] given that 0<f<c If x is Unif(a,b) and y is Unif(0,c), determine experimentally the probability p for y< f(x). The sought integral is p(b-a)c.
10 Monte Carlo -integration More experiments lead to more accurate estimate for p. Confidence interval (error) is proportional to N^(-1/2). Not efficient for one dimensional integrals Length of confidence interval and asymptotic behavior depend only on p, not the dimension of the integral. Efficient way to get rough estimates for high dimensional integrals.
11 Monte Carlo Previous Monte Carlo does not apply directly to all cases Unbounded interval or function Possible to give up y variable Compute only E(f(x)) Cheaper but error analysis is more demanding Flat upper bound c can be replaced Find a pdf g such that f(x)< cg(x) Draw x:s from distribution g Aim for success rate p ~ 1
12 Monte Carlo applications Typical Monte Carlo case is (very) high dimensional integral arising from modelling ray propagation in material. Each collision is modelled with multidimensional integral (probabilities for absorption, scattering as functions of incident angle, energy, particle shapes, surface properties, adsorption in free path, etc) For single ray the complexity grows only linearily with number of collisions.
13 M C Example Consider scattering of laser beam from a material layer MSc thesis of Jukka Räbinä 2005 Goal is to simulate different statistics of the scattered image using Monte-Carlo
14 Experimental set up
15 Simulate propagation of ray in cloud of particles Basically ray tracing Positions and scattering directions of particles are random Scattering
16 Compute the intensity, center of mass etc of the scatter image captured by camera Goal of simulation I.e. an integral of a function involving the intensity distribution.
17 Simulation experiment Send parallel rays with normally distributed intensity Collect the (few) rays scattered to the camera
18 Simulation results Three different implementations of M-C Each with 100M simulated rays Differences in execution times and confidence intervals Differences can be explained after learning about variance reduction methods
Simulation. Variance reduction Example
Simulation Variance reduction Example M C Example Consider scattering of laser beam from a material layer MSc thesis of Jukka Räbinä 2005 Goal is to simulate different statistics of the scattered image
More informationA Random Number Based Method for Monte Carlo Integration
A Random Number Based Method for Monte Carlo Integration J Wang and G Harrell Department Math and CS, Valdosta State University, Valdosta, Georgia, USA Abstract - A new method is proposed for Monte Carlo
More informationLecture 34: Curves defined by Parametric equations
Curves defined by Parametric equations When the path of a particle moving in the plane is not the graph of a function, we cannot describe it using a formula that express y directly in terms of x, or x
More informationGAMES Webinar: Rendering Tutorial 2. Monte Carlo Methods. Shuang Zhao
GAMES Webinar: Rendering Tutorial 2 Monte Carlo Methods Shuang Zhao Assistant Professor Computer Science Department University of California, Irvine GAMES Webinar Shuang Zhao 1 Outline 1. Monte Carlo integration
More informationAdvanced Graphics. Path Tracing and Photon Mapping Part 2. Path Tracing and Photon Mapping
Advanced Graphics Path Tracing and Photon Mapping Part 2 Path Tracing and Photon Mapping Importance Sampling Combine importance sampling techniques Reflectance function (diffuse + specular) Light source
More informationAnti-aliasing and Monte Carlo Path Tracing. Brian Curless CSE 557 Autumn 2017
Anti-aliasing and Monte Carlo Path Tracing Brian Curless CSE 557 Autumn 2017 1 Reading Required: Marschner and Shirley, Section 13.4 (online handout) Pharr, Jakob, and Humphreys, Physically Based Ray Tracing:
More informationAnti-aliasing and Monte Carlo Path Tracing
Reading Required: Anti-aliasing and Monte Carlo Path Tracing Brian Curless CSE 557 Autumn 2017 Marschner and Shirley, Section 13.4 (online handout) Pharr, Jakob, and Humphreys, Physically Based Ray Tracing:
More informationGlobal Illumination The Game of Light Transport. Jian Huang
Global Illumination The Game of Light Transport Jian Huang Looking Back Ray-tracing and radiosity both computes global illumination Is there a more general methodology? It s a game of light transport.
More informationSPH3UW Unit 7.2 Reflection Page 1 of 7
SPH3UW Unit 7.2 Reflection Page 1 of 7 Notes Physics Tool box Law of Reflection On flat surfaces, the angle of incidence equals the angle of reflection. Diffuse Reflection when light is incident on a rough
More informationBasics of treatment planning II
Basics of treatment planning II Sastry Vedam PhD DABR Introduction to Medical Physics III: Therapy Spring 2015 Monte Carlo Methods 1 Monte Carlo! Most accurate at predicting dose distributions! Based on
More informationPractical Course WS12/13 Introduction to Monte Carlo Localization
Practical Course WS12/13 Introduction to Monte Carlo Localization Cyrill Stachniss and Luciano Spinello 1 State Estimation Estimate the state of a system given observations and controls Goal: 2 Bayes Filter
More informationNumerical Integration
Lecture 12: Numerical Integration (with a focus on Monte Carlo integration) Computer Graphics CMU 15-462/15-662, Fall 2015 Review: fundamental theorem of calculus Z b f(x)dx = F (b) F (a) a f(x) = d dx
More informationPropagation and Reflection of Light
Al-Saudia Virtual Academy Online tuition Pakistan Online Tutor Pakistan Propagation and Reflection of Light Q1. Define reflection of light. State the laws of reflection. Ans: REFLECTION OF LIGHT: When
More informationPhysics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. - Statistics and Error Analysis -
Physics 736 Experimental Methods in Nuclear-, Particle-, and Astrophysics - Statistics and Error Analysis - Karsten Heeger heeger@wisc.edu Feldman&Cousin what are the issues they deal with? what processes
More informationLecture Ray Model of Light. Physics Help Q&A: tutor.leiacademy.org
Lecture 1201 Ray Model of Light Physics Help Q&A: tutor.leiacademy.org Reflection of Light A ray of light, the incident ray, travels in a medium. When it encounters a boundary with a second medium, part
More informationIn the real world, light sources emit light particles, which travel in space, reflect at objects or scatter in volumetric media (potentially multiple
1 In the real world, light sources emit light particles, which travel in space, reflect at objects or scatter in volumetric media (potentially multiple times) until they are absorbed. On their way, they
More informationExperiment 8 Wave Optics
Physics 263 Experiment 8 Wave Optics In this laboratory, we will perform two experiments on wave optics. 1 Double Slit Interference In two-slit interference, light falls on an opaque screen with two closely
More informationReflection and Refraction of Light
PC1222 Fundamentals of Physics II Reflection and Refraction of Light 1 Objectives Investigate for reflection of rays from a plane surface, the dependence of the angle of reflection on the angle of incidence.
More information13 Distribution Ray Tracing
13 In (hereafter abbreviated as DRT ), our goal is to render a scene as accurately as possible. Whereas Basic Ray Tracing computed a very crude approximation to radiance at a point, in DRT we will attempt
More informationChapter 3. Physical phenomena: plane parallel plate. This chapter provides an explanation about how rays of light physically behave when
Chapter 3 Physical phenomena: plane parallel plate This chapter provides an explanation about how rays of light physically behave when propagating through different medium (different index of refraction).
More informationSynrad3D Photon propagation and scattering simulation
Synrad3D Photon propagation and scattering simulation G. Dugan, D. Sagan CLASSE, Cornell University, Ithaca, NY 14853 USA Abstract As part of the Bmad software library, a program called Synrad3D has been
More informationTracePro s Monte Carlo Raytracing Methods, reducing statistical noise, memory usage and raytrace times
TracePro s Monte Carlo Raytracing Methods, reducing statistical noise, memory usage and raytrace times Presented by : Lambda Research Corporation 25 Porter Rd. Littleton, MA 01460 www.lambdares.com Moderator:
More informationPractice problems from old exams for math 233
Practice problems from old exams for math 233 William H. Meeks III October 26, 2012 Disclaimer: Your instructor covers far more materials that we can possibly fit into a four/five questions exams. These
More informationAt the interface between two materials, where light can be reflected or refracted. Within a material, where the light can be scattered or absorbed.
At the interface between two materials, where light can be reflected or refracted. Within a material, where the light can be scattered or absorbed. The eye sees by focusing a diverging bundle of rays from
More information2017 Summer Course on Optical Oceanography and Ocean Color Remote Sensing. Monte Carlo Simulation
2017 Summer Course on Optical Oceanography and Ocean Color Remote Sensing Curtis Mobley Monte Carlo Simulation Delivered at the Darling Marine Center, University of Maine July 2017 Copyright 2017 by Curtis
More informationMonte Carlo Integration
Lecture 11: Monte Carlo Integration Computer Graphics and Imaging UC Berkeley CS184/284A, Spring 2016 Reminder: Quadrature-Based Numerical Integration f(x) Z b a f(x)dx x 0 = a x 1 x 2 x 3 x 4 = b E.g.
More informationMonte Carlo Ray Tracing. Computer Graphics CMU /15-662
Monte Carlo Ray Tracing Computer Graphics CMU 15-462/15-662 TODAY: Monte Carlo Ray Tracing How do we render a photorealistic image? Put together many of the ideas we ve studied: - color - materials - radiometry
More informationOptics. a- Before the beginning of the nineteenth century, light was considered to be a stream of particles.
Optics 1- Light Nature: a- Before the beginning of the nineteenth century, light was considered to be a stream of particles. The particles were either emitted by the object being viewed or emanated from
More informationLight and the Properties of Reflection & Refraction
Light and the Properties of Reflection & Refraction OBJECTIVE To study the imaging properties of a plane mirror. To prove the law of reflection from the previous imaging study. To study the refraction
More informationA 10-minute introduction to. SynRad+ A test-particle Monte Carlo simulator for synchrotron radiation
A 10-minute introduction to SynRad+ A test-particle Monte Carlo simulator for synchrotron radiation 1 The basics First, let s learn the SynRad+ terminology and the interface in a few slides. Or, if you
More informationCOMPUTER SIMULATION TECHNIQUES FOR ACOUSTICAL DESIGN OF ROOMS - HOW TO TREAT REFLECTIONS IN SOUND FIELD SIMULATION
J.H. Rindel, Computer simulation techniques for the acoustical design of rooms - how to treat reflections in sound field simulation. ASVA 97, Tokyo, 2-4 April 1997. Proceedings p. 201-208. COMPUTER SIMULATION
More informationPhoton Maps. The photon map stores the lighting information on points or photons in 3D space ( on /near 2D surfaces)
Photon Mapping 1/36 Photon Maps The photon map stores the lighting information on points or photons in 3D space ( on /near 2D surfaces) As opposed to the radiosity method that stores information on surface
More informationMATH 31A HOMEWORK 9 (DUE 12/6) PARTS (A) AND (B) SECTION 5.4. f(x) = x + 1 x 2 + 9, F (7) = 0
FROM ROGAWSKI S CALCULUS (2ND ED.) SECTION 5.4 18.) Express the antiderivative F (x) of f(x) satisfying the given initial condition as an integral. f(x) = x + 1 x 2 + 9, F (7) = 28.) Find G (1), where
More information1) Find. a) b) c) d) e) 2) The function g is defined by the formula. Find the slope of the tangent line at x = 1. a) b) c) e) 3) Find.
1 of 7 1) Find 2) The function g is defined by the formula Find the slope of the tangent line at x = 1. 3) Find 5 1 The limit does not exist. 4) The given function f has a removable discontinuity at x
More informationModeling Custom Surface Roughness with LucidShape 2D Scatter Curve BSDF Material
WHITE PAPER Modeling Custom Surface Roughness with LucidShape 2D Scatter Curve BSDF Material Author Andreas Bielawny, Ph.D. CAE Synopsys, Inc. Abstract LucidShape accurately simulates how light interacts
More informationReflection AB5 Concave Mirror. Teacher s Notes
Reflection: Concave Mirror Teacher s Notes Main Topic Subtopic Learning Level Technology Level Activity Type Required Equipment Optional Equipment & Color Reflection Middle Low Student and Optical Set
More informationMonte Carlo Techniques. Professor Stephen Sekula Guest Lecture PHY 4321/7305 Sep. 3, 2014
Monte Carlo Techniques Professor Stephen Sekula Guest Lecture PHY 431/7305 Sep. 3, 014 What are Monte Carlo Techniques? Computational algorithms that rely on repeated random sampling in order to obtain
More informationPhysics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. - Statistical Methods -
Physics 736 Experimental Methods in Nuclear-, Particle-, and Astrophysics - Statistical Methods - Karsten Heeger heeger@wisc.edu Course Schedule and Reading course website http://neutrino.physics.wisc.edu/teaching/phys736/
More informationREFLECTION & REFRACTION
REFLECTION & REFRACTION OBJECTIVE: To study and verify the laws of reflection and refraction using a plane mirror and a glass block. To see the virtual images that can be formed by the reflection and refraction
More informationMULTI-DIMENSIONAL MONTE CARLO INTEGRATION
CS580: Computer Graphics KAIST School of Computing Chapter 3 MULTI-DIMENSIONAL MONTE CARLO INTEGRATION 2 1 Monte Carlo Integration This describes a simple technique for the numerical evaluation of integrals
More informationISyE 6416: Computational Statistics Spring Lecture 13: Monte Carlo Methods
ISyE 6416: Computational Statistics Spring 2017 Lecture 13: Monte Carlo Methods Prof. Yao Xie H. Milton Stewart School of Industrial and Systems Engineering Georgia Institute of Technology Determine area
More informationChapter 33 The Nature and Propagation of Light by C.-R. Hu
Chapter 33 The Nature and Propagation of Light by C.-R. Hu Light is a transverse wave of the electromagnetic field. In 1873, James C. Maxwell predicted it from the Maxwell equations. The speed of all electromagnetic
More informationUNIVERSITY OF CAMBRIDGE INTERNATIONAL EXAMINATIONS General Certificate of Education Ordinary Level
UNIVERSITY OF CAMBRIDGE INTERNATIONAL EXAMINATIONS General Certificate of Education Ordinary Level * 9 4 95570362* ADDITIONAL MATHEMATICS 4037/12 Paper 1 May/June 2010 Additional Materials: Answer Booklet/Paper
More informationPhysics 1C. Lecture 22A. "There are two ways of spreading light: to be the candle or the mirror that reflects it." --Edith Wharton
Physics 1C Lecture 22A "There are two ways of spreading light: to be the candle or the mirror that reflects it." --Edith Wharton The Nature of Light An interesting question developed as to the nature of
More informationRadiography and Tomography Lab 2: Computing Radiographs
Radiography and Tomography Lab 2: Computing Radiographs An Abbreviated View of Radiography Transmission Radiography is the process of measuring and recording changes in a high-energy particle beam (X-rays,
More informationProbability and Statistics for Final Year Engineering Students
Probability and Statistics for Final Year Engineering Students By Yoni Nazarathy, Last Updated: April 11, 2011. Lecture 1: Introduction and Basic Terms Welcome to the course, time table, assessment, etc..
More informationThis part of the course addresses the idea of path guiding and how it can be used to reduce variance of Monte Carlo volumetric light transport
This part of the course addresses the idea of path guiding and how it can be used to reduce variance of Monte Carlo volumetric light transport simulation. We will also show that the theory of zero-variance
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
1 DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. A Direct Simulation-Based Study of Radiance in a Dynamic Ocean LONG-TERM GOALS Dick K.P. Yue Center for Ocean Engineering
More informationACCURACY AND EFFICIENCY OF MONTE CARLO METHOD. Julius Goodman. Bechtel Power Corporation E. Imperial Hwy. Norwalk, CA 90650, U.S.A.
- 430 - ACCURACY AND EFFICIENCY OF MONTE CARLO METHOD Julius Goodman Bechtel Power Corporation 12400 E. Imperial Hwy. Norwalk, CA 90650, U.S.A. ABSTRACT The accuracy of Monte Carlo method of simulating
More informationFigure 1: Derivation of Bragg s Law
What is Bragg s Law and why is it Important? Bragg s law refers to a simple equation derived by English physicists Sir W. H. Bragg and his son Sir W. L. Bragg in 1913. This equation explains why the faces
More informationPrecalculus 2 Section 10.6 Parametric Equations
Precalculus 2 Section 10.6 Parametric Equations Parametric Equations Write parametric equations. Graph parametric equations. Determine an equivalent rectangular equation for parametric equations. Determine
More informationPhysics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. Lecture 14
Physics 736 Experimental Methods in Nuclear-, Particle-, and Astrophysics Lecture 14 Karsten Heeger heeger@wisc.edu Course Schedule and Reading course website http://neutrino.physics.wisc.edu/teaching/phys736/
More informationRay Optics I. Last time, finished EM theory Looked at complex boundary problems TIR: Snell s law complex Metal mirrors: index complex
Phys 531 Lecture 8 20 September 2005 Ray Optics I Last time, finished EM theory Looked at complex boundary problems TIR: Snell s law complex Metal mirrors: index complex Today shift gears, start applying
More informationChapter 35. The Nature of Light and the Laws of Geometric Optics
Chapter 35 The Nature of Light and the Laws of Geometric Optics Introduction to Light Light is basic to almost all life on Earth. Light is a form of electromagnetic radiation. Light represents energy transfer
More informationdq dt I = Irradiance or Light Intensity is Flux Φ per area A (W/m 2 ) Φ =
Radiometry (From Intro to Optics, Pedrotti -4) Radiometry is measurement of Emag radiation (light) Consider a small spherical source Total energy radiating from the body over some time is Q total Radiant
More informationspecular diffuse reflection.
Lesson 8 Light and Optics The Nature of Light Properties of Light: Reflection Refraction Interference Diffraction Polarization Dispersion and Prisms Total Internal Reflection Huygens s Principle The Nature
More informationLecture 17: Solid Modeling.... a cubit on the one side, and a cubit on the other side Exodus 26:13
Lecture 17: Solid Modeling... a cubit on the one side, and a cubit on the other side Exodus 26:13 Who is on the LORD's side? Exodus 32:26 1. Solid Representations A solid is a 3-dimensional shape with
More informationPHY 171 Lecture 6 (January 18, 2012)
PHY 171 Lecture 6 (January 18, 2012) Light Throughout most of the next 2 weeks, we will be concerned with the wave properties of light, and phenomena based on them (interference & diffraction). Light also
More informationMotivation. Advanced Computer Graphics (Fall 2009) CS 283, Lecture 11: Monte Carlo Integration Ravi Ramamoorthi
Advanced Computer Graphics (Fall 2009) CS 283, Lecture 11: Monte Carlo Integration Ravi Ramamoorthi http://inst.eecs.berkeley.edu/~cs283 Acknowledgements and many slides courtesy: Thomas Funkhouser, Szymon
More informationImperfections and Errors (Chapter 6)
Imperfections and Errors (Chapter 6) EC4630 Radar and Laser Cross Section Fall 011 Prof. D. Jenn jenn@nps.navy.mil www.nps.navy.mil/jenn AY011 1 Imperfections and Errors Imperfections and errors are deviations
More informationLECTURE 15 REFLECTION & REFRACTION. Instructor: Kazumi Tolich
LECTURE 15 REFLECTION & REFRACTION Instructor: Kazumi Tolich Lecture 15 2 18.1 The ray model of light Source of light rays Ray diagrams Seeing objects Shadows 18.2 Reflection Diffuse reflection The plane
More informationGG450 4/5/2010. Today s material comes from p and in the text book. Please read and understand all of this material!
GG450 April 6, 2010 Seismic Reflection I Today s material comes from p. 32-33 and 81-116 in the text book. Please read and understand all of this material! Back to seismic waves Last week we talked about
More informationRealization of Improved Solar UV Diffusers
Realization of Improved Solar UV Diffusers Tomi Pulli 1, Josef Schreder 2, Allard Partosoebroto 3, Gregor Huelsen 4, Janne Askola 1, Joop Mes 3, Petri Kärhä 1,5, and Julian Gröbner 4 1. Metrology Research
More information4X Beam Expander P/N SPZ17022 and Optional Beam Splitter P/N SPZ17027 User Notes
4X Beam Expander P/N SPZ17022 and Optional Beam Splitter P/N SPZ17027 User Notes Ophir-Spiricon Inc. 60 West 1000 North Logan, UT 84321 For Sales, Service or Technical Support Phone (435) 753-3729 Fax
More informationHome Lab 7 Refraction, Ray Tracing, and Snell s Law
Home Lab Week 7 Refraction, Ray Tracing, and Snell s Law Home Lab 7 Refraction, Ray Tracing, and Snell s Law Activity 7-1: Snell s Law Objective: Verify Snell s law Materials Included: Laser pointer Cylindrical
More informationMonte Carlo Integration
Lecture 15: Monte Carlo Integration Computer Graphics and Imaging UC Berkeley Reminder: Quadrature-Based Numerical Integration f(x) Z b a f(x)dx x 0 = a x 1 x 2 x 3 x 4 = b E.g. trapezoidal rule - estimate
More informationEngineered Diffusers Intensity vs Irradiance
Engineered Diffusers Intensity vs Irradiance Engineered Diffusers are specified by their divergence angle and intensity profile. The divergence angle usually is given as the width of the intensity distribution
More informationLight and refractive index
17 Fig. 7.1 shows a ray of light incident on a rectangular glass block at point X. W P X air glass Q R S Fig. 7.1 The ray of light is refracted at X. On Fig. 7.1, (a) draw the normal at X, [1] (b) draw
More informationDirection Fields; Euler s Method
Direction Fields; Euler s Method It frequently happens that we cannot solve first order systems dy (, ) dx = f xy or corresponding initial value problems in terms of formulas. Remarkably, however, this
More informationAnti-aliasing and Monte Carlo Path Tracing. Brian Curless CSE 457 Autumn 2017
Anti-aliasing and Monte Carlo Path Tracing Brian Curless CSE 457 Autumn 2017 1 Reading Required: Marschner and Shirley, Section 13.4 (online handout) Further reading: Pharr, Jakob, and Humphreys, Physically
More informationRay Optics. Lecture 23. Chapter 34. Physics II. Course website:
Lecture 23 Chapter 34 Physics II Ray Optics Course website: http://faculty.uml.edu/andriy_danylov/teaching/physicsii Today we are going to discuss: Chapter 34: Section 34.1-3 Ray Optics Ray Optics Wave
More informationCh. 25 The Reflection of Light
Ch. 25 The Reflection of Light 25. Wave fronts and rays We are all familiar with mirrors. We see images because some light is reflected off the surface of the mirror and into our eyes. In order to describe
More informationRandom Numbers Random Walk
Random Numbers Random Walk Computational Physics Random Numbers Random Walk Outline Random Systems Random Numbers Monte Carlo Integration Example Random Walk Exercise 7 Introduction Random Systems Deterministic
More informationChapter 33 cont. The Nature of Light and Propagation of Light (lecture 2) Dr. Armen Kocharian
Chapter 33 cont The Nature of Light and Propagation of Light (lecture 2) Dr. Armen Kocharian Polarization of Light Waves The direction of polarization of each individual wave is defined to be the direction
More informationLight. Form of Electromagnetic Energy Only part of Electromagnetic Spectrum that we can really see
Light Form of Electromagnetic Energy Only part of Electromagnetic Spectrum that we can really see Facts About Light The speed of light, c, is constant in a vacuum. Light can be: REFLECTED ABSORBED REFRACTED
More informationPhysics 11 Chapter 18: Ray Optics
Physics 11 Chapter 18: Ray Optics "... Everything can be taken from a man but one thing; the last of the human freedoms to choose one s attitude in any given set of circumstances, to choose one s own way.
More informationModeling the ORTEC EX-100 Detector using MCNP
Modeling the ORTEC EX-100 Detector using MCNP MCNP is a general-purpose Monte Carlo radiation transport code for modeling the interaction of radiation with materials based on composition and density. MCNP
More informationNested Sampling: Introduction and Implementation
UNIVERSITY OF TEXAS AT SAN ANTONIO Nested Sampling: Introduction and Implementation Liang Jing May 2009 1 1 ABSTRACT Nested Sampling is a new technique to calculate the evidence, Z = P(D M) = p(d θ, M)p(θ
More informationProbabilistic Robotics
Probabilistic Robotics Discrete Filters and Particle Filters Models Some slides adopted from: Wolfram Burgard, Cyrill Stachniss, Maren Bennewitz, Kai Arras and Probabilistic Robotics Book SA-1 Probabilistic
More informationSIMULATION AND VISUALIZATION IN THE EDUCATION OF COHERENT OPTICS
SIMULATION AND VISUALIZATION IN THE EDUCATION OF COHERENT OPTICS J. KORNIS, P. PACHER Department of Physics Technical University of Budapest H-1111 Budafoki út 8., Hungary e-mail: kornis@phy.bme.hu, pacher@phy.bme.hu
More information12. Reflection and Refraction
12. Reflection and Refraction Wave fronts Spherical waves are solutions to Maxwell s equations Molecules scatter spherical waves Spherical waves can add up to plane waves Reflection and refraction from
More informationdq dt I = Irradiance or Light Intensity is Flux Φ per area A (W/m 2 ) Φ =
Radiometry (From Intro to Optics, Pedrotti -4) Radiometry is measurement of Emag radiation (light) Consider a small spherical source Total energy radiating from the body over some time is Q total Radiant
More informationdistribution ray-tracing
distribution ray-tracing 1 distribution ray-tracing use many rays to compute average values over pixel areas, time, area lights, reflected directions,... 2 antialiasing origin compute average color subtended
More informationThe Rendering Equation. Computer Graphics CMU /15-662, Fall 2016
The Rendering Equation Computer Graphics CMU 15-462/15-662, Fall 2016 Review: What is radiance? Radiance at point p in direction N is radiant energy ( #hits ) per unit time, per solid angle, per unit area
More informationOrthogonal Projection Matrices. Angel and Shreiner: Interactive Computer Graphics 7E Addison-Wesley 2015
Orthogonal Projection Matrices 1 Objectives Derive the projection matrices used for standard orthogonal projections Introduce oblique projections Introduce projection normalization 2 Normalization Rather
More informationIntroduction to Mobile Robotics Bayes Filter Particle Filter and Monte Carlo Localization. Wolfram Burgard
Introduction to Mobile Robotics Bayes Filter Particle Filter and Monte Carlo Localization Wolfram Burgard 1 Motivation Recall: Discrete filter Discretize the continuous state space High memory complexity
More informationA 10-minute introduction to. Molfow+ A test-particle Monte Carlo simulator for UHV systems
A 10-minute introduction to Molfow+ A test-particle Monte Carlo simulator for UHV systems 1 The basics First, let s learn the Molflow terminology and the interface in a few slides. Or, if you prefer learning
More informationMATHEMATICS SYLLABUS SECONDARY 5th YEAR
European Schools Office of the Secretary-General Pedagogical Development Unit Ref.: 011-01-D-7-en- Orig.: EN MATHEMATICS SYLLABUS SECONDARY 5th YEAR 4 period/week course APPROVED BY THE JOINT TEACHING
More informationUncertainty analysis in spatial environmental modelling. Kasia Sawicka Athens, 29 Mar 2017
Uncertainty analysis in spatial environmental modelling Kasia Sawicka Athens, 29 Mar 2017 Uncertainty propagation overview Input ± U Model (± U) Output ± U Input data Model parameters Model structure Solution
More informationdistribution ray tracing
distribution ray tracing computer graphics distribution ray tracing 2006 fabio pellacini 1 distribution ray tracing use many rays to compute average values over pixel areas, time, area lights, reflected
More informationWhat is Monte Carlo Modeling*?
What is Monte Carlo Modeling*? Monte Carlo Modeling is a statisitcal method used here to simulate radiative transfer by simulating photon (or more exactly light rays/beams) interaction with a medium. MC
More informationRules for Deviation of Light Rays During Refraction
REFLECTION OF LIGHT Refraction of light is the phenomenon due to which a ray of light deviates from its path, at the surface of separation of two media, when the ray of light is travelling from one optical
More informationCS580: Ray Tracing. Sung-Eui Yoon ( 윤성의 ) Course URL:
CS580: Ray Tracing Sung-Eui Yoon ( 윤성의 ) Course URL: http://sglab.kaist.ac.kr/~sungeui/gcg/ Recursive Ray Casting Gained popularity in when Turner Whitted (1980) recognized that recursive ray casting could
More informationOutline. Monte Carlo Radiation Transport Modeling Overview (MCNP5/6) Monte Carlo technique: Example. Monte Carlo technique: Introduction
Monte Carlo Radiation Transport Modeling Overview () Lecture 7 Special Topics: Device Modeling Outline Principles of Monte Carlo modeling Radiation transport modeling with Utilizing Visual Editor (VisEd)
More informationKlaus Dehmelt EIC Detector R&D Weekly Meeting November 28, 2011 GEM SIMULATION FRAMEWORK
Klaus Dehmelt EIC Detector R&D Weekly Meeting November 28, 2011 GEM SIMULATION FRAMEWORK Overview GEM Simulation Framework in the context of Simulation Studies for a High Resolution Time Projection Chamber
More informationRadiometry (From Intro to Optics, Pedrotti 1-4) Radiometry is measurement of Emag radiation (light) Consider a small spherical source Assume a black
Radiometry (From Intro to Optics, Pedrotti -4) Radiometry is measurement of Emag radiation (light) Consider a small spherical source Assume a black body type emitter: uniform emission Total energy radiating
More informationGhost and Stray Light Analysis using TracePro. February 2012 Webinar
Ghost and Stray Light Analysis using TracePro February 2012 Webinar Moderator: Andy Knight Technical Sales Manager Lambda Research Corporation Presenter: Michael Gauvin Vice President of Sales Lambda Research
More informationRendering Smoke & Clouds
Rendering Smoke & Clouds Game Design Seminar 2007 Jürgen Treml Talk Overview 1. Introduction to Clouds 2. Virtual Clouds based on physical Models 1. Generating Clouds 2. Rendering Clouds using Volume Rendering
More informationChapter 1. Introduction
Chapter 1 Introduction A Monte Carlo method is a compuational method that uses random numbers to compute (estimate) some quantity of interest. Very often the quantity we want to compute is the mean of
More informationAn interesting related problem is Buffon s Needle which was first proposed in the mid-1700 s.
Using Monte Carlo to Estimate π using Buffon s Needle Problem An interesting related problem is Buffon s Needle which was first proposed in the mid-1700 s. Here s the problem (in a simplified form). Suppose
More information