Markov processes, lab 2
|
|
- Lora Allen
- 6 years ago
- Views:
Transcription
1 Lunds Universitet Matematikcentrum Matematisk statistik FMSF15/MASC03 Markovprocesser Markov processes, lab 2 The first part of the lab is about simple Poisson processes. You will simulate and analyse Poisson processes for various intensities, on the line as well as in the plane. In addition, the intensity of a real process will be estimated. In the final part, you will do both simulation and estimation for a non-homogeneous processes. 1 Preparations Read through the instructions and answer the following questions. It is required that you bring your written answers with you to the computer lab. In the beginning of the computer session the lab supervisors will check your answers before you carry out the rest of the assignment. 1. What is meant by independent increments of a stochastic process? 2. What is the distribution of the number of events N(t) in the interval (0, t] for a Poisson process? What is the mean and standard deviation of N(t)? 3. What is the distribution of the time between two consecutive events in a Poisson process? 4. We have observed a process that may be modelled as a Poisson process in the interval (0, t]. (a) Derive the ML estimator of the intensity λ. (b) How can one compute a confidence interval for the estimator in (a)? 5. If one has observed two independent Poisson processes with different intensities λ 1 and λ 2, during different time intervals, how can one compute a confidence interval for the difference λ 1 λ 2? 6. How can one simulate a non-homogeneous Poisson process? 1
2 2 The Poisson process Log in at one of the PCs in the computer room MH:230 or MH:231, using the user name and password which you have obtained for this course. Click on the icon MClogin on the desktop and login again with the same user name and password. This will attach the hard drive L: where your working directory will be saved. Note that you need to do this before you start up the software package Matlab. Choose the latest version of Matlab from the Start menu. If you have problems either logging in or starting Matlab, ask for help. After starting Matlab write mh init( fmsf15 ). This will add the directory containing the m-files which you will use in this assignment. 2.1 Simulation of exponentially distributed holding times The distance between two events in a Poisson process is exponentially distributed with intensity λ, where λ is the intensity of the Poisson process. To simulate a Poisson process we need to simulate these distances, but standard MATLAB can only generate uniform and normal random variables (the Statistics Toolbox can do a lot more, however). We can generate exponentially distributed random variables through a transformation of uniform ones. Assume that we want to simulate random variables with the distribution function F (y), and that F (y) is invertible. Denote its inverse by F 1 (y) and let U be a random variable that has a uniform distribution on (0, 1). If we put Y = F 1 (U), then P (Y y) = P (F 1 (U) y) = P (F (F 1 (U)) F (y)) = P (U F (y)) = F (y). This shows that the transformed variable Y has distribution function F, and we can obtain exponentially distributed variables from a uniform U as Y = ln (1 U)/λ, where the intensity of the exponential distribution is λ. In MATLAB one can generate a vector of 100 exponentially distributed random variables with intensity λ through the commands: >> u=rand(1,100); >> y=-log(1-u)/lambda; 2.2 The counting process We now have a vector y of inter-event times, but our interest is often in the absolute time points of the events. We can obtain these by summing the interevent times. MATLAB provides a function cumsum, doing just that. Through 2
3 T=cumsum(y) we obtain a vector containing the time points of the simulated events. The value of the process is the number of events so far, hence 1 at the time point of the first event, 2 at the time point of the second one, and so on. Such a vector may be created as N=1:100, if we want 100 events. The process can be plotted by stairs(t,n). Choose an intensity λ and make a few simulations using different random sequences u. Plot them in the same figure using hold. 2.3 Estimation of a Poisson intensity The intensity λ of a Poisson process is the expected number of events per time unit. The parameter λ may be estimated by dividing the number of events during a time interval by the length of that time interval. Estimate the intensity for some realisations and compute confidence intervals. Do the cover they the true λ? Make a somewhat longer realisation and compute a confidence interval. How does the width of the confidence interval change with the number of events? 2.4 Estimation of the intensity of a real process The file coal.dat contains data about coal mine catastrophes in the United Kingdom from 1851 to The first column contains day of month, the second column contains the month, the third column contains the year, the fourth column contains the number of days elapsed thst year, the fifth column is the number of days since the last catastrophe and the sixth column contains the number of casualties in the catastrophe. The number of catastrophes since the start time may be viewed as a stochastic process that we may plot through summing the inter-event times. Load the data using the command load coal.dat and sum the inter-event times by T=cumsum(coal(:,5)). Put N=1:size(coal,1) and plot the process. Does is look like a Poisson process? 3
4 Estimate the intensity λ and compute a confidence interval. Give the answer in the unit catastrophes/year. The intensity does not appear constant throughout the whole period. Split the data into two time periods with two different intensities λ 1 and λ 2. Estimate the intensities and compute a confidence interval for the difference λ 1 λ 2. Plot the process with years on the x-axis by, for example, >> stairs(coal(:,3)+coal(:,4)/365,n) When do the changes occur? What can the underlying reasons be? Also look at the accumulative casualties process by >> stairs(coal(:,3)+coal(:,4)/365,cumsum(coal(:,6))) 2.5 A different method of simulation on the line and in the plane When we studied the Poisson process in one dimension, we used a method based on inter-event times. The lecture notes also describe a different method. This method is based on the result that, given n events until t, these n events are uniformly and independently distributed on [0, t]. Hence we can simulate a Poisson process on [0, t] by first drawing a number from Po(λt) and then drawing n random variables T i from a uniform distribution on the interval. Make a few simulations on [0, 100] for some different intensities. Use the function porand to simulate the Poisson variable. Use sort to sort the event times. Plot Poisson events as marks using e.g. plot(t,0*t, x ). The method can also be extended to simulate a Poisson process on the plane. The intensity is the the number of points per area unit. Simulate a Poisson process in the area 0 x 1, 0 y 1, and plot the result for some different intensities. Use, for example plot(x,y, x ). 4
5 Simulate a Poisson process on the triangle with corners in (0, 0), (0, 1) and (1, 0). Is there more that one way of doing that? 3 Non-homogeneous Poisson processes 3.1 Simulation Find a way to simulate a non-homogeneous process over the interval [0, t] with t = 10, having intensity λ(u) = 4(1 cos(u 2 /4)). Plot the Poisson points T i and the intensity function in the same graph, using e.g. >> u=0:0.01:10; >> plot(u,4*(1-cos(u.^2/4)),t,0*t, x ) Is the result as you expect it? 3.2 Simulation and estimation We will now consider an inhomogeneous Poisson process, N(t), on an interval t [0, S] with intensity function λ(t) = θ 1 + θ 2 cos(2πt/s) + θ 3 sin(2πt/s) where θ 1 θ θ 2 3 to ensure λ(t) 0 for all 0 t S. Given event times T 1,..., T N(S), the negative log-likelihood function is l(θ) = Λ(S) + ln(n(s)!) N(S) k=1 ln(λ(t k )) where Λ(S) = S 0 λ(t)dt = θ 1S. The following pre-written functions will be used in this part: inhom_poisson_simulate: Simulates a random sample from the model. inhom_poisson_est: Calculates estimates of the model parameters, the estimation covariances, and approximate 95% confidence intervals. 5
6 Use help inhom_poisson_simulate, etc for more information and also read the Matlab code. Start with parameter values S = 10 and θ = [20, 15, 0]. Use inhom_poisson_simulate and simulate a sample with these parameters. As before, plot the intensity function and simulated event points and convince yourself that the results are as expected. The pre-written function inhom_poisson_est mentioned above uses a robustified Newton optimisation on l(θ) to find the Maximum Likelihood estimate of the parameters for a given process sample. Find the estimates of the parameters using the sample you generated above. Do the confidence interval cover the true parameter values? The accuracy of the estimates depends on the true parameter values. Simulate new samples with different parameter values and run the estimation function for each sample. How low can you make the intensity before the estimates deteriorate significantly? Also increase the intensity; are the three parameters equally easy to estimate accurately? Why? 6
Markov processes, lab 1
Lunds Universitet Matematikcentrum Matematisk statistik FMSF15/MASC03 Markovprocesser Markov processes, lab 1 The aim of the lab is to demonstrate how Markov chains work and how one can use MATLAB as a
More informationSenior astrophysics Lab 2: Evolution of a 1 M star
Senior astrophysics Lab 2: Evolution of a 1 M star Name: Checkpoints due: Friday 13 April 2018 1 Introduction This is the rst of two computer labs using existing software to investigate the internal structure
More informationPolarSync Quick Start
PolarSync Quick Start Installation and Use In this Quick Start guide, we will cover installing the PolarSync program and using it as a teacher, student or guest. I. Installing PolarSync... 1 II. Teacher
More informationManual for SOA Exam MLC.
Chapter 10. Poisson processes. Section 10.5. Nonhomogenous Poisson processes Extract from: Arcones Fall 2009 Edition, available at http://www.actexmadriver.com/ 1/14 Nonhomogenous Poisson processes Definition
More informationLecture 20. Poisson Processes. Text: A Course in Probability by Weiss STAT 225 Introduction to Probability Models March 26, 2014
Lecture 20 Text: A Course in Probability by Weiss 12.1 STAT 225 Introduction to Probability Models March 26, 2014 Whitney Huang Purdue University 20.1 Agenda 1 2 20.2 For a specified event that occurs
More informationLinear Motion with Constant Acceleration
Linear Motion 1 Linear Motion with Constant Acceleration Overview: First you will attempt to walk backward with a constant acceleration, monitoring your motion with the ultrasonic motion detector. Then
More informationLaboratory handouts, ME 340
Laboratory handouts, ME 340 This document contains summary theory, solved exercises, prelab assignments, lab instructions, and report assignments for Lab 4. 2014-2016 Harry Dankowicz, unless otherwise
More informationThe story of the film so far... Mathematics for Informatics 4a. Continuous-time Markov processes. Counting processes
The story of the film so far... Mathematics for Informatics 4a José Figueroa-O Farrill Lecture 19 28 March 2012 We have been studying stochastic processes; i.e., systems whose time evolution has an element
More informationUniversity of Minnesota Nano Center Standard Operating Procedure
University of Minnesota Nano Center Standard Operating Procedure Equipment Name: Zeta Potential Analyzer Model: Stabino Location: PAN 185 Badger Name: Not on Badger Revision Number: 0-Inital release Revisionist:
More informationLab 1: Numerical Solution of Laplace s Equation
Lab 1: Numerical Solution of Laplace s Equation ELEC 3105 last modified August 27, 2012 1 Before You Start This lab and all relevant files can be found at the course website. You will need to obtain an
More informationSteelSmart System Cold Formed Steel Design Software Download & Installation Instructions
Step 1 - Login or Create an Account at the ASI Portal: Login: https://portal.appliedscienceint.com/account/login Create Account: https://portal.appliedscienceint.com/account/register 2 0 1 7 A p p l i
More informationLaboratory handout 5 Mode shapes and resonance
laboratory handouts, me 34 82 Laboratory handout 5 Mode shapes and resonance In this handout, material and assignments marked as optional can be skipped when preparing for the lab, but may provide a useful
More informationThe exponential distribution and the Poisson process
The exponential distribution and the Poisson process 1-1 Exponential Distribution: Basic Facts PDF f(t) = { λe λt, t 0 0, t < 0 CDF Pr{T t) = 0 t λe λu du = 1 e λt (t 0) Mean E[T] = 1 λ Variance Var[T]
More information1 Overview of Simulink. 2 State-space equations
Modelling and simulation of engineering systems Simulink Exercise 1 - translational mechanical systems Dr. M. Turner (mct6@sun.engg.le.ac.uk 1 Overview of Simulink Simulink is a package which runs in the
More informationPre-Lab Exercise Full Name:
L07 Rotational Motion and the Moment of Inertia 1 Pre-Lab Exercise Full Name: Lab Section: Hand this in at the beginning of the lab period. The grade for these exercises will be included in your lab grade
More informationRelated topics Velocity, acceleration, force, gravitational acceleration, kinetic energy, and potential energy
Related topics Velocity, acceleration, force, gravitational acceleration, kinetic energy, and potential energy Principle A mass, which is connected to a cart via a silk thread, drops to the floor. The
More information15EE103L ELECTRIC CIRCUITS LAB RECORD
15EE103L ELECTRIC CIRCUITS LAB RECORD REGISTER NO: NAME OF THE STUDENT: SEMESTER: DEPARTMENT: INDEX SHEET S.No. Date of Experiment Name of the Experiment Date of submission Marks Staff Sign 1 Verification
More informationComputational Chemistry Lab Module: Conformational Analysis of Alkanes
Introduction Computational Chemistry Lab Module: Conformational Analysis of Alkanes In this experiment, we will use CAChe software package to model the conformations of butane, 2-methylbutane, and substituted
More informationApplied Probability. School of Mathematics and Statistics, University of Sheffield. (University of Sheffield) Applied Probability / 8
Applied Probability School of Mathematics and Statistics, University of Sheffield 2018 19 (University of Sheffield) Applied Probability 2018 19 1 / 8 Introduction You will have seen probability models.
More informationEnvironmental Systems Research Institute
Introduction to ArcGIS ESRI Environmental Systems Research Institute Redlands, California 2 ESRI GIS Development Arc/Info (coverage model) Versions 1-7 from 1980 1999 Arc Macro Language (AML) ArcView (shapefile
More informationSRV02-Series Rotary Experiment # 1. Position Control. Student Handout
SRV02-Series Rotary Experiment # 1 Position Control Student Handout SRV02-Series Rotary Experiment # 1 Position Control Student Handout 1. Objectives The objective in this experiment is to introduce the
More informationMahdi karbasian* & Zoubi Ibrahim
International Journal of Industrial Engineering & Production Research (010) pp. 105-110 September 010, Volume 1, Number International Journal of Industrial Engineering & Production Research ISSN: 008-4889
More informationES205 Analysis and Design of Engineering Systems: Lab 1: An Introductory Tutorial: Getting Started with SIMULINK
ES205 Analysis and Design of Engineering Systems: Lab 1: An Introductory Tutorial: Getting Started with SIMULINK What is SIMULINK? SIMULINK is a software package for modeling, simulating, and analyzing
More informationZetasizer Nano-ZS User Instructions
Zetasizer Nano-ZS User Instructions 1. Activate the instrument computer by logging in to CORAL. If needed, log in to the local instrument computer Username: zetasizer. Password: zetasizer. 2. Instrument
More informationNewton's 2 nd Law. . Your end results should only be interms of m
Newton's nd Law Introduction: In today's lab you will demonstrate the validity of Newton's Laws in predicting the motion of a simple mechanical system. The system that you will investigate consists of
More informationComputer simulation of radioactive decay
Computer simulation of radioactive decay y now you should have worked your way through the introduction to Maple, as well as the introduction to data analysis using Excel Now we will explore radioactive
More informationMolecular Modeling and Conformational Analysis with PC Spartan
Molecular Modeling and Conformational Analysis with PC Spartan Introduction Molecular modeling can be done in a variety of ways, from using simple hand-held models to doing sophisticated calculations on
More informationIIE5. Modul Electricity II. Earth s magnetic field
IIE5 Modul Electricity II Earth s magnetic field In the present experiment, the earth s magnetic field is measured for different axes of rotation of induced voltage. From the amplitude and the frequency
More informationPhotoelectric Photometry of the Pleiades Student Manual
Name: Lab Partner: Photoelectric Photometry of the Pleiades Student Manual A Manual to Accompany Software for the Introductory Astronomy Lab Exercise Edited by Lucy Kulbago, John Carroll University 11/24/2008
More informationIE 581 Introduction to Stochastic Simulation
1. List criteria for choosing the majorizing density r (x) when creating an acceptance/rejection random-variate generator for a specified density function f (x). 2. Suppose the rate function of a nonhomogeneous
More informationAssignment 2: Conformation Searching (50 points)
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard September 16, 2015 Assignment 2: Conformation Searching (50 points) In this assignment, you will use the Spartan software package to investigate some conformation
More informationEVA Tutorial #2 PEAKS OVER THRESHOLD APPROACH. Rick Katz
1 EVA Tutorial #2 PEAKS OVER THRESHOLD APPROACH Rick Katz Institute for Mathematics Applied to Geosciences National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu Home page: www.isse.ucar.edu/staff/katz/
More informationSuperCELL Data Programmer and ACTiSys IR Programmer User s Guide
SuperCELL Data Programmer and ACTiSys IR Programmer User s Guide This page is intentionally left blank. SuperCELL Data Programmer and ACTiSys IR Programmer User s Guide The ACTiSys IR Programmer and SuperCELL
More informationUSE OF MATLAB TO UNDERSTAND BASIC MATHEMATICS
USE OF MATLAB TO UNDERSTAND BASIC MATHEMATICS Sanjay Gupta P. G. Department of Mathematics, Dev Samaj College For Women, Punjab ( India ) ABSTRACT In this paper, we talk about the ways in which computer
More informationCentre for Mathematical Sciences HT 2018 Mathematical Statistics
Lund University Stationary stochastic processes Centre for Mathematical Sciences HT 2018 Mathematical Statistics Computer exercise 1 in Stationary stochastic processes, HT 18. The purpose of this computer
More informationLab VI Light Emitting Diodes ECE 476
Lab VI Light Emitting Diodes ECE 476 I. Purpose This experiment examines the properties of light emitting diodes. The use of the monochromator for studying the spectrum of light sources is also examined.
More informationDOAS measurements of Atmospheric Species
Practical Environmental Measurement Techniques: DOAS measurements of Atmospheric Species Last change of document: April 14, 2014 Supervisor: Dr. Andreas Richter, room U2090, tel 62103, and Dr. Folkard
More informationCAS Exam MAS-1. Howard Mahler. Stochastic Models
CAS Exam MAS-1 Howard Mahler Stochastic Models 2019 Stochastic Models, HCM 12/31/18, Page 1 The slides are in the same order as the sections of my study guide. Section # Section Name 1 Introduction 2 Exponential
More informationCHAPTER 3 MATHEMATICAL AND SIMULATION TOOLS FOR MANET ANALYSIS
44 CHAPTER 3 MATHEMATICAL AND SIMULATION TOOLS FOR MANET ANALYSIS 3.1 INTRODUCTION MANET analysis is a multidimensional affair. Many tools of mathematics are used in the analysis. Among them, the prime
More informationLAB 2 1. Measurement of 2. Binomial Distribution
LAB 2 Gan Phys 3700 1. Measurement of π In this exercise we will determine a value for π by throwing darts: a) Determine π by throwing a dart 100 or more times. Use an 8x11 inch sheet of paper with a circle
More informationNewton's Laws and Atwood's Machine
Newton's Laws and Atwood's Machine Purpose: In this lab we will verify Newton's Second Law of Motion within estimated uncertainty and propose an explanation if verification is not within estimated uncertainty.
More informationExamination paper for TMA4275 Lifetime Analysis
Department of Mathematical Sciences Examination paper for TMA4275 Lifetime Analysis Academic contact during examination: Ioannis Vardaxis Phone: 95 36 00 26 Examination date: Saturday May 30 2015 Examination
More informationIntroduction to Computer Tools and Uncertainties
Experiment 1 Introduction to Computer Tools and Uncertainties 1.1 Objectives To become familiar with the computer programs and utilities that will be used throughout the semester. To become familiar with
More informationON SITE SYSTEMS Chemical Safety Assistant
ON SITE SYSTEMS Chemical Safety Assistant CS ASSISTANT WEB USERS MANUAL On Site Systems 23 N. Gore Ave. Suite 200 St. Louis, MO 63119 Phone 314-963-9934 Fax 314-963-9281 Table of Contents INTRODUCTION
More informationLab 5: Calculating an equilibrium constant
Chemistry 162 The following write-up is inaccurate for the particular chemicals we are using. Please have all sections up through and including the data tables ready before class on Wednesday, February
More informationMATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f
MATRIAL MCHANICS, S226 COMPUTR LAB 4 MICRO MCHANICS 2 2 2 f m f f m T m f m f f m v v + + = + PART A SPHRICAL PARTICL INCLUSION Consider a solid granular material, a so called particle composite, shown
More informationDetecting Galactic HI line using 4-m SRT
Detecting Galactic HI line using 4-m SRT 1 Goal of the experiment The final goal of the experiment is to detect the galactic HI line emission and to understand the physics behind it. In this experiment,
More informationLaboration 8a. Relaxation, T 1 -measurement with inversion recovery
, T 1 -measurement with inversion recovery KR Theory The way the magnetizations returns to equilibrium, relaxation, is a very important concept in NMR, for example, due to the fact that the rate of relaxation
More informationRewrite logarithmic equations 2 3 = = = 12
EXAMPLE 1 Rewrite logarithmic equations Logarithmic Form a. log 2 8 = 3 Exponential Form 2 3 = 8 b. log 4 1 = 0 4 0 = 1 log 12 = 1 c. 12 12 1 = 12 log 4 = 1 d. 1/4 1 4 1 = 4 GUIDED PRACTICE for Example
More informationSolutions For Stochastic Process Final Exam
Solutions For Stochastic Process Final Exam (a) λ BMW = 20 0% = 2 X BMW Poisson(2) Let N t be the number of BMWs which have passes during [0, t] Then the probability in question is P (N ) = P (N = 0) =
More informationNINE CHOICE SERIAL REACTION TIME TASK
instrumentation and software for research NINE CHOICE SERIAL REACTION TIME TASK MED-STATE NOTATION PROCEDURE SOF-700RA-8 USER S MANUAL DOC-025 Rev. 1.3 Copyright 2013 All Rights Reserved MED Associates
More informationD.T.M: TRANSFER TEXTBOOKS FROM ONE SCHOOL TO ANOTHER
Destiny Textbook Manager allows users with full access to transfer Textbooks from one school site to another and receive transfers from the warehouse In this tutorial you will learn how to: Requirements:
More informationFinite Element Modules for Enhancing Undergraduate Transport Courses: Application to Fuel Cell Fundamentals
Finite Element Modules for Enhancing Undergraduate Transport Courses: Application to Fuel Cell Fundamentals Originally published in 007 American Society for Engineering Education Conference Proceedings
More informationSession-Based Queueing Systems
Session-Based Queueing Systems Modelling, Simulation, and Approximation Jeroen Horters Supervisor VU: Sandjai Bhulai Executive Summary Companies often offer services that require multiple steps on the
More informationPhotoelectric Photometry of the Pleiades Student Manual
Photoelectric Photometry of the Pleiades Student Manual A Manual to Accompany Software for the Introductory Astronomy Lab Exercise Document SM 2: Version 1.1.1 lab Department of Physics Gettysburg College
More informationMarkov decision processes and interval Markov chains: exploiting the connection
Markov decision processes and interval Markov chains: exploiting the connection Mingmei Teo Supervisors: Prof. Nigel Bean, Dr Joshua Ross University of Adelaide July 10, 2013 Intervals and interval arithmetic
More informationTHE GEIGER-MULLER TUBE AND THE STATISTICS OF RADIOACTIVITY
GMstats. THE GEIGER-MULLER TUBE AN THE STATISTICS OF RAIOACTIVITY This experiment examines the Geiger-Muller counter, a device commonly used for detecting and counting ionizing radiation. Various properties
More informationCentre for Mathematical Sciences HT 2017 Mathematical Statistics. Study chapters 6.1, 6.2 and in the course book.
Lund University Stationary stochastic processes Centre for Mathematical Sciences HT 2017 Mathematical Statistics Computer exercise 2 in Stationary stochastic processes, HT 17. The purpose with this computer
More informationAugust 7, 2007 NUMERICAL SOLUTION OF LAPLACE'S EQUATION
August 7, 007 NUMERICAL SOLUTION OF LAPLACE'S EQUATION PURPOSE: This experiment illustrates the numerical solution of Laplace's Equation using a relaxation method. The results of the relaxation method
More informationExponential Distribution and Poisson Process
Exponential Distribution and Poisson Process Stochastic Processes - Lecture Notes Fatih Cavdur to accompany Introduction to Probability Models by Sheldon M. Ross Fall 215 Outline Introduction Exponential
More informationComposite FEM Lab-work
Composite FEM Lab-work You may perform these exercises in groups of max 2 persons. You may also between exercise 5 and 6. Be critical on the results obtained! Exercise 1. Open the file exercise1.inp in
More informationNewton s Second Law. Sample
Newton s Second Law Experiment 4 INTRODUCTION In your discussion of Newton s first law, you learned that when the sum of the forces acting on an object is zero, its velocity does not change. However, when
More informationMATLAB BASICS. Instructor: Prof. Shahrouk Ahmadi. TA: Kartik Bulusu
MATLAB BASICS Instructor: Prof. Shahrouk Ahmadi 1. What are M-files TA: Kartik Bulusu M-files are files that contain a collection of MATLAB commands or are used to define new MATLAB functions. For the
More informationHW5 computer problem: modelling a dye molecule Objectives By the end of this lab you should: Know how to manipulate complex numbers and matrices. Prog
HW5 computer problem: modelling a dye molecule Objectives By the end of this lab you should: Know how to manipulate complex numbers and matrices. Programming environment: Spyder. To implement this, do
More informationLecture 1. Introduction to stochastic processes: concepts, definitions and applications.
Lecture 1. Introduction to stochastic processes: concepts, definitions and applications. Jesper Rydén Department of Mathematics, Uppsala University jesper@math.uu.se Stationary Stochastic Processes Fall
More informationReview. DS GA 1002 Statistical and Mathematical Models. Carlos Fernandez-Granda
Review DS GA 1002 Statistical and Mathematical Models http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall16 Carlos Fernandez-Granda Probability and statistics Probability: Framework for dealing with
More informationChem 253. Tutorial for Materials Studio
Chem 253 Tutorial for Materials Studio This tutorial is designed to introduce Materials Studio 7.0, which is a program used for modeling and simulating materials for predicting and rationalizing structure
More informationfirst name (print) last name (print) brock id (ab17cd) (lab date)
(ta initials) first name (print) last name (print) brock id (ab17cd) (lab date) Experiment 1 Archimedes principle Prelab preparation Print a copy of this experiment to bring to your scheduled lab session.
More informationMatlab code for stochastic model updating
Matlab code for stochastic model updating The Matlab code provided performs the stochastic model updating methods described in [1] and [2] using the data obtained from experiments. This document contains
More informationRobustness and Distribution Assumptions
Chapter 1 Robustness and Distribution Assumptions 1.1 Introduction In statistics, one often works with model assumptions, i.e., one assumes that data follow a certain model. Then one makes use of methodology
More information1-D Convection-Diffusion Lab
Computational Fluid Dynamics -D Convection-Diffusion Lab The lab. uses scientificworkplace symbolic calculus and maths editor software (SWP) This file Concevtion-Diffusion-Lab is available from Blackboard
More informationPoisson Processes for Neuroscientists
Poisson Processes for Neuroscientists Thibaud Taillefumier This note is an introduction to the key properties of Poisson processes, which are extensively used to simulate spike trains. For being mathematical
More informationpka AND MOLAR MASS OF A WEAK ACID
Experiment 10 pka AND MOLAR MASS OF A WEAK ACID Adapted by the Chemistry Faculty of Eastern Michigan University from EQUL 305,written by Richard C. Bell, Lebanon Valley College, published by Chemical Education
More informationLab 2: Photon Counting with a Photomultiplier Tube
Lab 2: Photon Counting with a Photomultiplier Tube 1 Introduction 1.1 Goals In this lab, you will investigate properties of light using a photomultiplier tube (PMT). You will assess the quantitative measurements
More informationQualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama
Qualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama Instructions This exam has 7 pages in total, numbered 1 to 7. Make sure your exam has all the pages. This exam will be 2 hours
More informationShown below is a sample titration curve for a diprotic acid. Note the two equivalence points.
EXPERIMENT 9 Titration Curve for a Polyprotic Acid INTRODUCTION Other than by strength and concentration, another way of classifying acids involves the number of H + ions an acid can donate. A monoprotic
More informationParameter Estimation and Fitting to Data
Parameter Estimation and Fitting to Data Parameter estimation Maximum likelihood Least squares Goodness-of-fit Examples Elton S. Smith, Jefferson Lab 1 Parameter estimation Properties of estimators 3 An
More informationDetermination of Density 1
Introduction Determination of Density 1 Authors: B. D. Lamp, D. L. McCurdy, V. M. Pultz and J. M. McCormick* Last Update: February 1, 2013 Not so long ago a statistical data analysis of any data set larger
More informationSkyGlobe Planetarium
SkyGlobe Planetarium Introduction: This exercise will simulate the night sky and demonstrate a number of principles of the celestial sphere and the motions of the Earth and planets. Getting Started: 1.
More informationSRV02-Series Rotary Experiment # 7. Rotary Inverted Pendulum. Student Handout
SRV02-Series Rotary Experiment # 7 Rotary Inverted Pendulum Student Handout SRV02-Series Rotary Experiment # 7 Rotary Inverted Pendulum Student Handout 1. Objectives The objective in this experiment is
More informationDetermination of the Equivalent Weight and the K a or K b for a Weak Acid or Base
INTRODUCTION Determination of the Equivalent Weight and the K a or K b for a Weak Acid or Base Chemists frequently make use of the equivalent weight (eq. wt.) as the basis for volumetric calculations.
More informationBayesian course - problem set 5 (lecture 6)
Bayesian course - problem set 5 (lecture 6) Ben Lambert November 30, 2016 1 Stan entry level: discoveries data The file prob5 discoveries.csv contains data on the numbers of great inventions and scientific
More informationExercises Stochastic Performance Modelling. Hamilton Institute, Summer 2010
Exercises Stochastic Performance Modelling Hamilton Institute, Summer Instruction Exercise Let X be a non-negative random variable with E[X ]
More informationThe Coupled Pendulum Experiment
The Coupled Pendulum Experiment In this lab you will briefly study the motion of a simple pendulum, after which you will couple two pendulums and study the properties of this system. 1. Introduction to
More informationAcid-Base Titration Curves Using a ph Meter
Acid-Base Titration Curves Using a ph Meter Introduction: In this experiment you will use a ph sensor to collect volume and ph data as you titrate two acids with sodium hydroxide. You will obtain titration
More informationSolution: The process is a compound Poisson Process with E[N (t)] = λt/p by Wald's equation.
Solutions Stochastic Processes and Simulation II, May 18, 217 Problem 1: Poisson Processes Let {N(t), t } be a homogeneous Poisson Process on (, ) with rate λ. Let {S i, i = 1, 2, } be the points of the
More informationEEE161 Applied Electromagnetics Laboratory 2
Dr. Milica Marković Applied Electromagnetics Laboratory page 1 EEE161 Applied Electromagnetics Laboratory 2 Instructor: Dr. Milica Marković Office: Riverside Hall 5026 Email: milica@csus.edu Web:http://gaia.ecs.csus.edu/
More informationRADIATION PROCEDURES MANUAL Procedure Cover Sheet
RADIATION PROCEDURES MANUAL Procedure Cover Sheet Procedure Title: Radioactive Material Inventory Procedure Number: TSO-09-16-REV 0 Effective Date: June 18, 2009 Approved By: Date: 11 August, 2009 Technical
More information171S1.1 Graphs, Functions, Models August 23, 2012
MAT 171 D10 8:00 9:15 MAT 171 D11 9:30 10:45 Aug 21 7:11 AM Aug 21 7:14 AM MAT 171 D14 2:00 3:15 MAT 171 D15 3:30 4:45 Aug 21 7:14 AM Aug 21 7:14 AM 1 MAT 171 Dr. Claude Moore, CFCC CHAPTER 1: Graphs,
More informationQ = (c) Assuming that Ricoh has been working continuously for 7 days, what is the probability that it will remain working at least 8 more days?
IEOR 4106: Introduction to Operations Research: Stochastic Models Spring 2005, Professor Whitt, Second Midterm Exam Chapters 5-6 in Ross, Thursday, March 31, 11:00am-1:00pm Open Book: but only the Ross
More informationFleXScan User Guide. for version 3.1. Kunihiko Takahashi Tetsuji Yokoyama Toshiro Tango. National Institute of Public Health
FleXScan User Guide for version 3.1 Kunihiko Takahashi Tetsuji Yokoyama Toshiro Tango National Institute of Public Health October 2010 http://www.niph.go.jp/soshiki/gijutsu/index_e.html User Guide version
More informationFigure 2.1 The Inclined Plane
PHYS-101 LAB-02 One and Two Dimensional Motion 1. Objectives The objectives of this experiment are: to measure the acceleration due to gravity using one-dimensional motion, i.e. the motion of an object
More informationIntroduction to Reliability Theory (part 2)
Introduction to Reliability Theory (part 2) Frank Coolen UTOPIAE Training School II, Durham University 3 July 2018 (UTOPIAE) Introduction to Reliability Theory 1 / 21 Outline Statistical issues Software
More informationThe Phase Change Lab: Freezing and Melting of Water
The Phase Change Lab: Freezing and Melting of Water Experiment 3 Freezing temperature is the temperature at which a substance turns from a liquid to a solid. Melting temperature is the temperature at which
More informationEmpirical Models Interpolation Polynomial Models
Mathematical Modeling Lia Vas Empirical Models Interpolation Polynomial Models Lagrange Polynomial. Recall that two points (x 1, y 1 ) and (x 2, y 2 ) determine a unique line y = ax + b passing them (obtained
More information1. What two kinds of patterns/variation in carbon dioxide concentrations ([CO 2 ]) can you observe?
Case Study 1: Atmospheric carbon dioxide increase The Keeling curve describes the patterns of carbon dioxide concentrations in the atmosphere from the 1950s onwards. The original curve is shown below (left).
More informationGAS CHROMATOGRAPHY (GC)
GAS CHROMATOGRAPHY (GC) Pre-Lab Questions Questions are to be answered before the beginning of the laboratory. The answers are due at the beginning of each experiment (the questions are for credit and
More informationCOMPOSITION OF CONCURRENT FORCES
COMPOSITION OF CONCURRENT FORCES OBJECTIVE: To see if the result of applying three forces on an object can be determined by ADDING the three forces as VECTORS. GENERAL PROCEDURE: In this experiment your
More informationRENEWAL PROCESSES AND POISSON PROCESSES
1 RENEWAL PROCESSES AND POISSON PROCESSES Andrea Bobbio Anno Accademico 1997-1998 Renewal and Poisson Processes 2 Renewal Processes A renewal process is a point process characterized by the fact that the
More informationTitrator 3.0 Tutorial: Calcite precipitation
Titrator 3.0 Tutorial: Calcite precipitation November 2008 Steve Cabaniss A. Introduction This brief tutorial is intended to acquaint you with some of the features of the program Titrator. It assumes that
More information