A Random Number Based Method for Monte Carlo Integration
|
|
- Mabel Pearson
- 6 years ago
- Views:
Transcription
1 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 integration This method is more efficient with wider coverage, including improper integrals, while the classical Monte Carlo integration can only handle bounded domain integrals To implement this method in computer programming, you only need a random number generator Unlike the deterministic numerical integration methods, the expected error of this method is independent of the integral dimensionality This method is powerful and dominates other numerical integral methods for the higherdimensional integrals Keywords: Monte Carlo integration; Output analysis; Trapezoid rule; Simpson rule; Error analysis; Improper integral 1 Introduction In single variable calculus ([8]), the general way to evaluate a definite integral is to find a formula for one of the antiderivatives of such that However some antiderivatives are not easy to find, and others, like and have been proved that no such closed-form formulas exist [8] There is a need to evaluate the definite integral of such functions Two different types of numerical methods are involved in approximating The first numerical method is called the deterministic integration, such as the trapezoidal rule and the Simpson rule The second numerical method is called the Monte Carlo integration In Section 2, we introduce the Monte Carlo integration in a general -dimensional integral setting Selecting a probability density function, having the same domain as the integral domain, is an important part in the Monte Carlo integration The efficiency of the Monte Carlo integration depends on the degree of difficulty for generating the random sample according to the selected density function The most popular method in practice is the simple Monte Carlo integration This method only works for bounded integral domain It does not work for any improper integrals with unbounded domains In Section 3, we propose a new method for Monte Carlo integration It is simple and efficient It can be used to approximate integrals with unbounded domains for any dimensional integrals The key idea is to convert the original integral domain to a new integral domain a unit hyper-cube The implementation of this new method is much easier A random number generator is the only thing you need In Section 4, we conduct the variance analysis It helps us to decide when the Monte Carlo simulation run should stop, and what the random sample size should be For the new proposed method, the expected value of error is ( ), which is independent of the integral dimensionality Comparing the new method with Trapezoid rule and the Simpson Rule, the new method is getting better and better when the integral dimension number is large than 8 It is particularly attractive and powerful for higher dimensional integrals 2 Monte Carlo Integration Monte Carlo ([2], [4]) integration is a simple and powerful method for approximating complicated integrals In general, we are interested in evaluating the integral of a function over a domain, ( ) (1) Here we use vector notation to indicate that is a - dimensional function In fact, Monte Carlo techniques are more attractive for estimating higher-dimensional integrals We assume that there is a probability density function (pdf) defined on the same domain Therefore (1) equals to, [ ( ) ] ( ) (2) ( ) The integral (2) is equivalent to the expectation of with respect to a random variable distributed following to the cdf [ ( ) ( ) ] (3) According to the cdf we generate independent random samples The Monte Carlo integration estimator is defined as following, (4)
2 Based on the Strong Law of Large Numbers ([1]), we have, with probability 1, This means, with probability 1, The Monte Carlo approximation estimator of [ ( ) ( ) ] is a consistent The efficiency of this method depends on how to generate the random samples The important sampling technique ([2], [7]) can be used here to reduce the Monte Carlo integration error (the sample variance) How to pick a cdf defined on is a major concern in using important sampling techniques Another important efficient factor is how to generate the random sample Sometime it may be very expensive depending on the selected cdf In practice, the most commonly used method is called the simple Monte Carlo simulation We define the volume of the integration domain in (1), The cdf selected for the simple Monte Carlo integration ( ) { In the implementation of simple Monte Carlo integration, we only need to generate random samples uniformly over the integral domain The efficiency is depending on the shape of integral domain 3 Using Random Numbers to Approximate Integrals As we discussed in the previous section, the shape of the integral domain is an issue here for the simple Monte Carlo integration It is only feasible when the shape of is a hyper-rectangle For an improper integral with unbounded domain, the simple Monte Carlo technique does not work It is impossible to generate random samples uniformly over an unbounded integral region For the general Monte Carlo integration technique, the issue here is how to select a cdf The necessary condition for is having the same domain Can we generate random samples from this density function efficiently and how? In order to solve these issues discussed above, we propose a new random number base method for the Monte Carlo integration The key idea is that we use the change variable (substitution) method to transfer the original integral domain (including unbounded case) into the new domain: a unit hyper-cube Therefore, we can easily generate random samples uniformly over a unite hypercube We introduce the new method for one-dimensional integration first This method can be easily generalized into multiple -dimensional situations What is a random number? By definition, a random number is a continuous random variable following a uniform distribution over the interval Its pdf is given by, { (5) The cumulative distribution function (cdf) is given by { (6) In connecting the random numbers, we begin with a special case Let be a one-dimension function and we want to approximate the integral value of on (7) If is a random number having cdf (5), then the expectation of is given by Therefore If are independent random numbers, we have that are independent and identically distributed random variables with a mean, which is equivalent to Therefore, by the Strong Law of Large Numbers, we have, with probability 1, We can use a large number of random numbers to approximate the definite integral (7) The natural approximation of based on random numbers is defined by, (8) We define this approximation technique as a random number based method for Monte Carlo integration In general, we want to approximate integral, (9) Here the integral lower and upper bounds and can be any real numbers, or positive and negative infinities For any arbitrary combination of and in (9), we will provide detailed steps to convert (9) as a basic integral (7) over the interval The key idea is to find a suitable variable substitution That we substitute variable in terms of This substitution will transfer the original
3 integral interval into the new integral interval The transformation process can be described as following, where the new integrand, (10) ( ) here is the Jacobian of the variable substitution If the transformation (10) exists, we can use the method (8) to approximate in (9), (11) Here will converge to the true integral value in (9) with probability 1 Let s classify all possible combinations of the integral lower and upper bounds and in (9) into different cases, and then we discuss each individual one Case A: For Case B: For here and are both finite real numbers and, we define the following substitution, implies implies, and ( ) ( ) We have done the transformation, here is a finite real number, we define a nonlinear transformation, implies, and implies Case C: For ( ) Using the negative sign to switch the integral lower and upper bounds, we have, This implies, Here may be unbounded when However this improper integral is convergent if is finite There is no problem when we perform the Monte Carlo simulation, since all pseudo random numbers are from open interval and can t be zero here is a finite real number, we define a nonlinear transformation, and implies, implies This implies, ( ) A typical application of Case C is the Normal cdf, For this special integral, the closed-form of its antiderivative does not exist
4 Case D: For we rewrite this integral into a sum of two integrals, We apply Case C method to the first integral and Case B method to the second integral, where and respectively Therefore, we have, We have discussed all possible situations of the integral (9) so far In other words, we can use the random number based method to approximate the integral (9) for any arbitrary lower and upper bounds This idea can be generalized to approximate the multiple integrals For a general -dimensional integral, (12) we may transfer (12) into the following integral, (13) Integral (13) has a unit hyper-cube domain and is equivalent to, where are independent random numbers To approximate the value of the -dimensional integral (13), we only need to generate independent random numbers over the unit hyper-cube : Therefore, we can estimate the integral value of (14) by, By the Strong Law of Large Numbers, with probability 1, We only need a random number generator ([5], [6]) to implement our method in computer programming 4 Error Analysis The implementation the Monte Carlo integration is a computer simulation The Monte Carlo integration estimate convergences to the true integral value with probability 1 In practice, we need to consider when to stop the Monte Carlo simulation, which essentially is deciding the appropriate sample size The quality of the Monte Carlo integration estimator is depending on this sample size Typical error analysis in Monte Carlo simulation is called the output analysis using sample variance to evaluate the quality of Monte Carlo estimator ([2], [7]) In this section, we are assuming that the variance is finite so that the Central Limit Theory can be applied The mean value of in (14) [ ] This implies that is an unbiased estimator variance of is given by, ( ) The where, is the variance ( ) and is the unit uniform random vector over the hyper-cube ([2], [7]) In practice, we use the following sample variance to estimate the about true variance ( ) This variance decreases asymptotically to zero as The expected value of the error is proportional to the standard deviation, which is ( ) From the Central Limit Theory, within one- error, the confidence level is about 68% Within two- error, the confidence level is about 95% Within three- error, the confidence level is about 997% We can decide the sample size based on the desired confidence level The swing digit ideal in ([9]) is a good method for reporting simulation outputs To increase accuracy in one more digit, you need to increase the sample size 100 times The expected value of the error for this new method is always ( ), which does not depend on the integral dimensionality However, for the deterministic numerical integration, it does depend on the integral dimensionality For a -dimensional integral, the Trapezoid rule provides an error of of ([3]), and the Simpson rule provides an error In comparing error orders, for lower dimension integrals, the Trapezoid rule and Simpson rule are better than the Monte Carlo integration As the dimensionality increases, for, the Monte Carlo integration is better than the Trapezoid rule, and for, the Monte Carlo integration is better than the Simpson Rule Particularly, the Monte Carlo method is powerful for high dimension integrals When the dimensional number is higher than 20, the only working numerical method is the Monte Carlo integration in many existing software packages
5 5 Conclusions In reality, an analytical form for a given integral is sometimes difficult to find We turn this situation into numerical integration The popular deterministic numerical methods are the Trapezoid rule and the Simpson rule The statistical numerical method is the Monte Carlo integration In general, you need to select a probability density function, which has the same domain as the given integral domain Generating a sample according to the selected density function is sometimes very expensive or inefficient For the most popular method - simple Monte Carlo integration, the integral domain has to be a bounded hyperrectangle This method does not work for any integral with an unbounded integral domain The new proposed method in this paper is more efficient and easier to implement in computer programming We have covered all possible integral situations, including improper integrals Detailed steps are provided for converting the given integral domain into a new domain a unit hyper-cube You only need to generate random samples uniformly over a unit hype-cube When the integral dimensional number is lower, there is no significant in errors for different numerical integration methods When the integral dimensional number is higher, say larger than 8, the new proposed method for Monte Carlo integration dominates all other numerical integration methods Transcendentals Addison Wesley, 10 th Edition, New York, NY, 2001 [9] A L Wang, A L and Kicey, C J On the Accuracy of Buffon's Needle: A Simulation Output Analysis Proceedings of the 49th Annual ACM SE Conference, ACM Press, New York, NY, , References [1] Billingsley, P Convergence of Probability Measures Wiley, 2 nd Edition, New York, NY, 1999 [2] Bratley, P, Fox, B L, and Schrage, L E A Guide to Simulation Springer-Verlag, New York, NY, 1987 [3] Burden, R L, Faires, J D, and Reynolds, A C Numerical Analysis PWS Publishers, 2 nd Edition, Boston, MA, 1979 [4] Hammersley, J M and Handscomb, D C 1964 Monte Carlo Methods, Methuen, London, 1964 [5] D H Lehmer, D H Mathematical Methods in Large- Scale Computing Units Annals of the Computation, Laboratory of Harvard University, 26, , 1951 [6] L Ecuyer, P Uniform Random Number Generation, Annals of Operations Research, 53, , 1994 [7] Ross, S M A Course in Simulation MacMillan, New York, NY, 1990 [8] Thomas Jr, G B, Weir, M D, Hass, J R, and Giordano, F R Thomas' Calculus Early
VARIANCE REDUCTION TECHNIQUES IN MONTE CARLO SIMULATIONS K. Ming Leung
POLYTECHNIC UNIVERSITY Department of Computer and Information Science VARIANCE REDUCTION TECHNIQUES IN MONTE CARLO SIMULATIONS K. Ming Leung Abstract: Techniques for reducing the variance in Monte Carlo
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 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 informationIntegration. Volume Estimation
Monte Carlo Integration Lab Objective: Many important integrals cannot be evaluated symbolically because the integrand has no antiderivative. Traditional numerical integration techniques like Newton-Cotes
More informationMonte Carlo Integration
Lab 18 Monte Carlo Integration Lab Objective: Implement Monte Carlo integration to estimate integrals. Use Monte Carlo Integration to calculate the integral of the joint normal distribution. Some multivariable
More informationProceedings of the 2014 Winter Simulation Conference A. Tolk, S. Y. Diallo, I. O. Ryzhov, L. Yilmaz, S. Buckley, and J. A. Miller, eds.
Proceedings of the 14 Winter Simulation Conference A. Tolk, S. Y. Diallo, I. O. Ryzhov, L. Yilmaz, S. Buckley, and J. A. Miller, eds. REMOVING THE INHERENT PARADOX OF THE BUFFON'S NEEDLE MONTE CARLO SIMULATION
More informationMETROPOLIS MONTE CARLO SIMULATION
POLYTECHNIC UNIVERSITY Department of Computer and Information Science METROPOLIS MONTE CARLO SIMULATION K. Ming Leung Abstract: The Metropolis method is another method of generating random deviates that
More informationComputational Methods. Randomness and Monte Carlo Methods
Computational Methods Randomness and Monte Carlo Methods Manfred Huber 2010 1 Randomness and Monte Carlo Methods Introducing randomness in an algorithm can lead to improved efficiencies Random sampling
More informationSimulation. Monte Carlo
Simulation Monte Carlo 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
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 informationLearning Objectives. Continuous Random Variables & The Normal Probability Distribution. Continuous Random Variable
Learning Objectives Continuous Random Variables & The Normal Probability Distribution 1. Understand characteristics about continuous random variables and probability distributions 2. Understand the uniform
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 informationProbabilistic (Randomized) algorithms
Probabilistic (Randomized) algorithms Idea: Build algorithms using a random element so as gain improved performance. For some cases, improved performance is very dramatic, moving from intractable to tractable.
More informationOn Kernel Density Estimation with Univariate Application. SILOKO, Israel Uzuazor
On Kernel Density Estimation with Univariate Application BY SILOKO, Israel Uzuazor Department of Mathematics/ICT, Edo University Iyamho, Edo State, Nigeria. A Seminar Presented at Faculty of Science, Edo
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 informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 13 Random Numbers and Stochastic Simulation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright
More informationSOME NOTES ON MULTIPLICATIVE CONGRUENTIAL RANDOM NUMBER GENERATORS WITH MERSENNE PRIME MODULUS Dr. James Harris*
JournaCof the South Carolina JLcademy of Science l(l):28-32 Fall 2003 SOME NOTES ON MULTIPLICATIVE CONGRUENTIAL RANDOM NUMBER GENERATORS WITH MERSENNE PRIME MODULUS 2 61-1 Dr. James Harris* *Department
More informationHalftoning and quasi-monte Carlo
Halftoning and quasi-monte Carlo Ken Hanson CCS-2, Methods for Advanced Scientific Simulations Los Alamos National Laboratory This presentation available at http://www.lanl.gov/home/kmh/ LA-UR-04-1854
More informationis bad when combined with a LCG with a small multiplier. In section 2 this observation is examined. Section 3 gives portable implementations for LCGs
Version: 25.11.92 (To appear in ACM TOMS) A Portable Uniform Random Number Generator Well Suited for the Rejection Method W. Hormann and G. Deringer University of Economics and Business Administration
More informationRounding Theorem the Possibility of Applying the Cryptosystems on the Decimal Numbers
Journal of Mathematics and Statistics 4 (): 5-20, 2008 ISSN 549-3644 2008 Science Publications Rounding Theorem the Possibility of Applying the Cryptosystems on the Decimal Numbers Rand Alfaris, Muhamad
More informationGenerating random samples from user-defined distributions
The Stata Journal (2011) 11, Number 2, pp. 299 304 Generating random samples from user-defined distributions Katarína Lukácsy Central European University Budapest, Hungary lukacsy katarina@phd.ceu.hu Abstract.
More informationProbability Models.S4 Simulating Random Variables
Operations Research Models and Methods Paul A. Jensen and Jonathan F. Bard Probability Models.S4 Simulating Random Variables In the fashion of the last several sections, we will often create probability
More informationCS 563 Advanced Topics in Computer Graphics Monte Carlo Integration: Basic Concepts. by Emmanuel Agu
CS 563 Advanced Topics in Computer Graphics Monte Carlo Integration: Basic Concepts by Emmanuel Agu Introduction The integral equations generally don t have analytic solutions, so we must turn to numerical
More informationSTAT 725 Notes Monte Carlo Integration
STAT 725 Notes Monte Carlo Integration Two major classes of numerical problems arise in statistical inference: optimization and integration. We have already spent some time discussing different optimization
More informationStochastic Function Norm Regularization of DNNs
Stochastic Function Norm Regularization of DNNs Amal Rannen Triki Dept. of Computational Science and Engineering Yonsei University Seoul, South Korea amal.rannen@yonsei.ac.kr Matthew B. Blaschko Center
More informationSampling and Monte-Carlo Integration
Sampling and Monte-Carlo Integration Sampling and Monte-Carlo Integration Last Time Pixels are samples Sampling theorem Convolution & multiplication Aliasing: spectrum replication Ideal filter And its
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 informationA Study of Cross-Validation and Bootstrap for Accuracy Estimation and Model Selection (Kohavi, 1995)
A Study of Cross-Validation and Bootstrap for Accuracy Estimation and Model Selection (Kohavi, 1995) Department of Information, Operations and Management Sciences Stern School of Business, NYU padamopo@stern.nyu.edu
More informationComputer Experiments. Designs
Computer Experiments Designs Differences between physical and computer Recall experiments 1. The code is deterministic. There is no random error (measurement error). As a result, no replication is needed.
More informationX On record with the USOE.
Textbook Alignment to the Utah Core 5th Grade Mathematics This alignment has been completed using an Independent Alignment Vendor from the USOE approved list (www.schools.utah.gov/curr/imc/indvendor.html.)
More informationSimultaneous Perturbation Stochastic Approximation Algorithm Combined with Neural Network and Fuzzy Simulation
.--- Simultaneous Perturbation Stochastic Approximation Algorithm Combined with Neural Networ and Fuzzy Simulation Abstract - - - - Keywords: Many optimization problems contain fuzzy information. Possibility
More informationMonte Carlo Simulations
Monte Carlo Simulations Lecture 2 December 8, 2014 Outline 1 Random Number Generation 2 Uniform Random Variables 3 Normal Random Variables 4 Correlated Random Variables Random Number Generation Monte Carlo
More informationChapter 3. Bootstrap. 3.1 Introduction. 3.2 The general idea
Chapter 3 Bootstrap 3.1 Introduction The estimation of parameters in probability distributions is a basic problem in statistics that one tends to encounter already during the very first course on the subject.
More informationMonte Carlo Simula/on and Copula Func/on. by Gerardo Ferrara
Monte Carlo Simula/on and Copula Func/on by Gerardo Ferrara Introduc)on A Monte Carlo method is a computational algorithm that relies on repeated random sampling to compute its results. In a nutshell,
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 informationTextbook Alignment to the Utah Core 5th Grade Mathematics
Textbook Alignment to the Utah Core 5th Grade Mathematics This alignment has been completed using an Independent Alignment Vendor from the USOE approved list (www.schools.utah.gov/curr/imc/indvendor.html.)
More informationImportance Sampling of Area Lights in Participating Media
Importance Sampling of Area Lights in Participating Media Christopher Kulla Marcos Fajardo Outline Previous Work Single Scattering Equation Importance Sampling for Point Lights Importance Sampling for
More informationMonte Carlo Integration
Chapter 2 Monte Carlo Integration This chapter gives an introduction to Monte Carlo integration. The main goals are to review some basic concepts of probability theory, to define the notation and terminology
More informationSIMULATION AND MONTE CARLO
JHU course no. 550.790, Modeling, Simulation, and Monte Carlo SIMULATION AND MONTE CARLO Some General Principles James C. Spall Johns Hopkins University Applied Physics Laboratory September 2004 Overview
More informationInclusion of Aleatory and Epistemic Uncertainty in Design Optimization
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Inclusion of Aleatory and Epistemic Uncertainty in Design Optimization Sirisha Rangavajhala
More informationModeling with Uncertainty Interval Computations Using Fuzzy Sets
Modeling with Uncertainty Interval Computations Using Fuzzy Sets J. Honda, R. Tankelevich Department of Mathematical and Computer Sciences, Colorado School of Mines, Golden, CO, U.S.A. Abstract A new method
More informationSimulation. Outline. Common Mistakes in Simulation (3 of 4) Common Mistakes in Simulation (2 of 4) Performance Modeling Lecture #8
Introduction (1 of 3) The best advice to those about to embark on a very large simulation is often the same as Punch s famous advice to those about to marry: Don t! Bratley, Fox and Schrage (1986) Simulation
More informationand F is an antiderivative of f
THE EVALUATION OF DEFINITE INTEGRALS Prepared by Ingrid Stewart, Ph.D., College of Southern Nevada Please Send Questions Comments to ingrid.stewart@csn.edu. Thank you! We have finally reached a point,
More informationTriple Integrals in Rectangular Coordinates
Triple Integrals in Rectangular Coordinates P. Sam Johnson April 10, 2017 P. Sam Johnson (NIT Karnataka) Triple Integrals in Rectangular Coordinates April 10, 2017 1 / 28 Overview We use triple integrals
More informationAssessing the Quality of the Natural Cubic Spline Approximation
Assessing the Quality of the Natural Cubic Spline Approximation AHMET SEZER ANADOLU UNIVERSITY Department of Statisticss Yunus Emre Kampusu Eskisehir TURKEY ahsst12@yahoo.com Abstract: In large samples,
More informationCombinatorial Search; Monte Carlo Methods
Combinatorial Search; Monte Carlo Methods Parallel and Distributed Computing Department of Computer Science and Engineering (DEI) Instituto Superior Técnico May 02, 2016 CPD (DEI / IST) Parallel and Distributed
More informationNon-Parametric Modeling
Non-Parametric Modeling CE-725: Statistical Pattern Recognition Sharif University of Technology Spring 2013 Soleymani Outline Introduction Non-Parametric Density Estimation Parzen Windows Kn-Nearest Neighbor
More informationAP Calculus BC Course Description
AP Calculus BC Course Description COURSE OUTLINE: The following topics define the AP Calculus BC course as it is taught over three trimesters, each consisting of twelve week grading periods. Limits and
More informationPRACTICE FINAL - MATH 1210, Spring 2012 CHAPTER 1
PRACTICE FINAL - MATH 2, Spring 22 The Final will have more material from Chapter 4 than other chapters. To study for chapters -3 you should review the old practice eams IN ADDITION TO what appears here.
More informationClassical Gradient Methods
Classical Gradient Methods Note simultaneous course at AMSI (math) summer school: Nonlin. Optimization Methods (see http://wwwmaths.anu.edu.au/events/amsiss05/) Recommended textbook (Springer Verlag, 1999):
More informationShort-Cut MCMC: An Alternative to Adaptation
Short-Cut MCMC: An Alternative to Adaptation Radford M. Neal Dept. of Statistics and Dept. of Computer Science University of Toronto http://www.cs.utoronto.ca/ radford/ Third Workshop on Monte Carlo Methods,
More informationLevel-set MCMC Curve Sampling and Geometric Conditional Simulation
Level-set MCMC Curve Sampling and Geometric Conditional Simulation Ayres Fan John W. Fisher III Alan S. Willsky February 16, 2007 Outline 1. Overview 2. Curve evolution 3. Markov chain Monte Carlo 4. Curve
More informationAn Average-Case Analysis of the k-nearest Neighbor Classifier for Noisy Domains
An Average-Case Analysis of the k-nearest Neighbor Classifier for Noisy Domains Seishi Okamoto Fujitsu Laboratories Ltd. 2-2-1 Momochihama, Sawara-ku Fukuoka 814, Japan seishi@flab.fujitsu.co.jp Yugami
More informationQuasi-Monte Carlo Methods Combating Complexity in Cost Risk Analysis
Quasi-Monte Carlo Methods Combating Complexity in Cost Risk Analysis Blake Boswell Booz Allen Hamilton ISPA / SCEA Conference Albuquerque, NM June 2011 1 Table Of Contents Introduction Monte Carlo Methods
More informationTesting unrolling optimization technique for quasi random numbers
Testing unrolling optimization technique for quasi random numbers Romain Reuillon David R.C. Hill LIMOS, UMR CNRS 6158 LIMOS, UMR CNRS 6158 Blaise Pascal University Blaise Pascal University ISIMA, Campus
More informationSimulation. Chapter Copyright 2010 Pearson Education, Inc. Publishing as Prentice Hall
Simulation Chapter 14 14-1 Chapter Topics The Monte Carlo Process Computer Simulation with Excel Spreadsheets Simulation of a Queuing System Continuous Probability Distributions Statistical Analysis of
More informationReproducibility in Stochastic Simulation
Reproducibility in Stochastic Simulation Prof. Michael Mascagni Department of Computer Science Department of Mathematics Department of Scientific Computing Graduate Program in Molecular Biophysics Florida
More informationUse of Extreme Value Statistics in Modeling Biometric Systems
Use of Extreme Value Statistics in Modeling Biometric Systems Similarity Scores Two types of matching: Genuine sample Imposter sample Matching scores Enrolled sample 0.95 0.32 Probability Density Decision
More informationUsing a Painting Metaphor to Rate Large Numbers
1 Using a Painting Metaphor to Rate Large Numbers of Objects Patrick Baudisch Integrated Publication and Information Systems Institute (IPSI) German National Research Center for Information Technology
More informationLECTURE NOTES Non-Linear Programming
CEE 6110 David Rosenberg p. 1 Learning Objectives LECTURE NOTES Non-Linear Programming 1. Write out the non-linear model formulation 2. Describe the difficulties of solving a non-linear programming model
More informationADAPTIVE METROPOLIS-HASTINGS SAMPLING, OR MONTE CARLO KERNEL ESTIMATION
ADAPTIVE METROPOLIS-HASTINGS SAMPLING, OR MONTE CARLO KERNEL ESTIMATION CHRISTOPHER A. SIMS Abstract. A new algorithm for sampling from an arbitrary pdf. 1. Introduction Consider the standard problem of
More information10703 Deep Reinforcement Learning and Control
10703 Deep Reinforcement Learning and Control Russ Salakhutdinov Machine Learning Department rsalakhu@cs.cmu.edu Policy Gradient I Used Materials Disclaimer: Much of the material and slides for this lecture
More informationSOME EXPLICIT RESULTS FOR THE DISTRIBUTION PROBLEM OF STOCHASTIC LINEAR PROGRAMMING
SOME EXPLICIT RESULTS FOR THE DISTRIBUTION PROBLEM OF STOCHASTIC LINEAR PROGRAMMING Afrooz Ansaripour Hillel J. Kumin School of Industrial and Systems Engineering University of Oklahoma, Norman, Oklahoma
More informationChapter 7 Transcendental Functions
Chapter 7 Transcendental Functions. Trapezoidal Rule and Simpson s Rule. Other Numerical Integration 3. Logarithmic and Eponential Functions 4. Inverse Trigonometric Functions 5. Comparing Symbolic Integration
More informationVISUALIZING NP-COMPLETENESS THROUGH CIRCUIT-BASED WIDGETS
University of Portland Pilot Scholars Engineering Faculty Publications and Presentations Shiley School of Engineering 2016 VISUALIZING NP-COMPLETENESS THROUGH CIRCUIT-BASED WIDGETS Steven R. Vegdahl University
More informationRecurrent Neural Network Models for improved (Pseudo) Random Number Generation in computer security applications
Recurrent Neural Network Models for improved (Pseudo) Random Number Generation in computer security applications D.A. Karras 1 and V. Zorkadis 2 1 University of Piraeus, Dept. of Business Administration,
More informationMonte Carlo Methods; Combinatorial Search
Monte Carlo Methods; Combinatorial Search Parallel and Distributed Computing Department of Computer Science and Engineering (DEI) Instituto Superior Técnico November 22, 2012 CPD (DEI / IST) Parallel and
More informationStandard 1 Students will expand number sense to include integers and perform operations with whole numbers, simple fractions, and decimals.
Stretch Standard 1 Students will expand number sense to include integers and perform operations with whole numbers, simple fractions, and decimals. Objective 1: Represent whole numbers and decimals from
More informationInterval Algorithms for Coin Flipping
IJCSNS International Journal of Computer Science and Network Security, VOL.10 No.2, February 2010 55 Interval Algorithms for Coin Flipping Sung-il Pae, Hongik University, Seoul, Korea Summary We discuss
More informationDynamic Programming. Other Topics
Dynamic Programming Other Topics 1 Objectives To explain the difference between discrete and continuous dynamic programming To discuss about multiple state variables To discuss the curse of dimensionality
More informationMONTE CARLO EVALUATION OF DEFINITE INTEGRALS
MISN-0-355 MONTE CARLO EVALUATION OF DEFINITE INTEGRALS by Robert Ehrlich MONTE CARLO EVALUATION OF DEFINITE INTEGRALS y 1. Integrals, 1-2 Dimensions a. Introduction.............................................1
More informationA Low Level Introduction to High Dimensional Sparse Grids
A Low Level Introduction to High Dimensional Sparse Grids http://people.sc.fsu.edu/ jburkardt/presentations/sandia 2007.pdf... John 1 Clayton Webster 2 1 Virginia Tech 2 Sandia National Laboratory. 21
More informationAn Excel Add-In for Capturing Simulation Statistics
2001 Joint Statistical Meetings Atlanta, GA An Excel Add-In for Capturing Simulation Statistics David C. Trindade david.trindade@sun.com Sun Microsystems, Inc. Cupertino, CA David Meade david.meade@amd.com
More informationApproximation of 3D-Parametric Functions by Bicubic B-spline Functions
International Journal of Mathematical Modelling & Computations Vol. 02, No. 03, 2012, 211-220 Approximation of 3D-Parametric Functions by Bicubic B-spline Functions M. Amirfakhrian a, a Department of Mathematics,
More informationKutztown University Kutztown, Pennsylvania. MAT 550: Foundations of Geometry
Kutztown University Kutztown, Pennsylvania I. Three semester hours; three clock hours; this is an elective course in the M.Ed. or M.A. degree programs in mathematics. II. Catalog Description: 3 s.h. Foundational
More informationCalculus I Review Handout 1.3 Introduction to Calculus - Limits. by Kevin M. Chevalier
Calculus I Review Handout 1.3 Introduction to Calculus - Limits by Kevin M. Chevalier We are now going to dive into Calculus I as we take a look at the it process. While precalculus covered more static
More informationDECISION SCIENCES INSTITUTE. Exponentially Derived Antithetic Random Numbers. (Full paper submission)
DECISION SCIENCES INSTITUTE (Full paper submission) Dennis Ridley, Ph.D. SBI, Florida A&M University and Scientific Computing, Florida State University dridley@fsu.edu Pierre Ngnepieba, Ph.D. Department
More informationUsing Modified Euler Method (MEM) for the Solution of some First Order Differential Equations with Initial Value Problems (IVPs).
Using Modified Euler Method (MEM) for the Solution of some First Order Differential Equations with Initial Value Problems (IVPs). D.I. Lanlege, Ph.D. * ; U.M. Garba, B.Sc.; and A. Aluebho, B.Sc. Department
More informationGENERATING RANDOM NUMBERS IN APL
GENERATING RANDOM NUMBERS IN APL by Thomas N. Herzog Loyola College of Baltimore 99 . ~.....-:- 100 APL is a vector-matrix based computer programming language which is a powerful tool for solving actuarial
More informationEstimating and Sampling Graphs with Multidimensional Random Walks. Group 2: Mingyan Zhao, Chengle Zhang, Biao Yin, Yuchen Liu
Estimating and Sampling Graphs with Multidimensional Random Walks Group 2: Mingyan Zhao, Chengle Zhang, Biao Yin, Yuchen Liu Complex Network Motivation Social Network Biological Network erence: www.forbe.com
More informationExperimental Study on Bound Handling Techniques for Multi-Objective Particle Swarm Optimization
Experimental Study on Bound Handling Techniques for Multi-Objective Particle Swarm Optimization adfa, p. 1, 2011. Springer-Verlag Berlin Heidelberg 2011 Devang Agarwal and Deepak Sharma Department of Mechanical
More informationThree Different Algorithms for Generating Uniformly Distributed Random Points on the N-Sphere
Three Different Algorithms for Generating Uniformly Distributed Random Points on the N-Sphere Jan Poland Oct 4, 000 Abstract We present and compare three different approaches to generate random points
More information7.6 Simple Monte Carlo Integration. 304 Chapter 7. Random Numbers
304 Chapter 7. Random Numbers psdes(&lword,&irword); Pseudo-DES encode the words. itemp=jflone (jflmsk & irword); Mask to a floating number between 1 and ++(*idum); 2. return (*(float *)&itemp)-1.0; Subtraction
More informationA noninformative Bayesian approach to small area estimation
A noninformative Bayesian approach to small area estimation Glen Meeden School of Statistics University of Minnesota Minneapolis, MN 55455 glen@stat.umn.edu September 2001 Revised May 2002 Research supported
More informationFunctions of Several Variables
Chapter 3 Functions of Several Variables 3.1 Definitions and Examples of Functions of two or More Variables In this section, we extend the definition of a function of one variable to functions of two or
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 informationEvaluation of Parallel Programs by Measurement of Its Granularity
Evaluation of Parallel Programs by Measurement of Its Granularity Jan Kwiatkowski Computer Science Department, Wroclaw University of Technology 50-370 Wroclaw, Wybrzeze Wyspianskiego 27, Poland kwiatkowski@ci-1.ci.pwr.wroc.pl
More informationDiscrete Optimization. Lecture Notes 2
Discrete Optimization. Lecture Notes 2 Disjunctive Constraints Defining variables and formulating linear constraints can be straightforward or more sophisticated, depending on the problem structure. The
More informationAn Evolutionary Algorithm for Minimizing Multimodal Functions
An Evolutionary Algorithm for Minimizing Multimodal Functions D.G. Sotiropoulos, V.P. Plagianakos and M.N. Vrahatis University of Patras, Department of Mamatics, Division of Computational Mamatics & Informatics,
More informationNondeterministic Query Algorithms
Journal of Universal Computer Science, vol. 17, no. 6 (2011), 859-873 submitted: 30/7/10, accepted: 17/2/11, appeared: 28/3/11 J.UCS Nondeterministic Query Algorithms Alina Vasilieva (Faculty of Computing,
More informationEast Penn School District Secondary Curriculum
East Penn School District Secondary Curriculum A Planned Course Statement for Analytic Geometry and Calculus (BC) AP Course # 360 Grade(s) 12 Department: Math ength of Period (mins.) 41 Total Clock Hours:
More informationA Superconvergent Monte Carlo Method for Multiple Integrals on the Grid
A Superconvergent Monte Carlo Method for Multiple Integrals on the Grid Sofiya Ivanovska, Emanouil Atanassov, and Aneta Karaivanova Institute for Parallel Processing - Bulgarian Academy of Sciences, Acad.
More informationA Fast Algorithm for Optimal Alignment between Similar Ordered Trees
Fundamenta Informaticae 56 (2003) 105 120 105 IOS Press A Fast Algorithm for Optimal Alignment between Similar Ordered Trees Jesper Jansson Department of Computer Science Lund University, Box 118 SE-221
More informationThe Limit Concept. Introduction to Limits. Definition of Limit. Example 1. Example 2. Example 3 4/7/2015
4/7/015 The Limit Concept Introduction to Limits Precalculus 1.1 The notion of a it is a fundamental concept of calculus. We will learn how to evaluate its and how they are used in the two basic problems
More informationFinite Element Convergence for Time-Dependent PDEs with a Point Source in COMSOL 4.2
Finite Element Convergence for Time-Dependent PDEs with a Point Source in COMSOL 4.2 David W. Trott and Matthias K. Gobbert Department of Mathematics and Statistics, University of Maryland, Baltimore County,
More informationepub WU Institutional Repository
epub WU Institutional Repository Josef Leydold and Wolfgang Hörmann The Automatic Generation of One- and Multi-dimensional Distributions with Transformed Density Rejection Paper Original Citation: Leydold,
More informationOVERVIEW OF MAPLETS FOR CALCULUS 1.3
OVERVIEW OF MAPLETS FOR CALCULUS 1.3 Douglas B. Meade Dept of Mathematics Univ of South Carolina Columbia, SC 29208 meade@math.sc.edu (803) 777-6183 Philip B. Yasskin Dept of Mathematics Texas A&M University
More informationNICOLAS BOURBAKI ELEMENTS OF MATHEMATICS. General Topology. Chapters 1-4. Springer-Verlag Berlin Heidelberg New York London Paris Tokyo
NICOLAS BOURBAKI ELEMENTS OF MATHEMATICS General Topology Chapters 1-4 Springer-Verlag Berlin Heidelberg New York London Paris Tokyo ADVICE TO THE READER v CONTENTS OF THE ELEMENTS OF MATHEMATICS SERIES
More informationYou ve already read basics of simulation now I will be taking up method of simulation, that is Random Number Generation
Unit 5 SIMULATION THEORY Lesson 39 Learning objective: To learn random number generation. Methods of simulation. Monte Carlo method of simulation You ve already read basics of simulation now I will be
More information