Stochastic Modelling
|
|
- Patrick Fitzgerald
- 6 years ago
- Views:
Transcription
1 Stochastic Modelling Simulating Random Walks and Markov Chains This lab sheets is available for downloading from as is the spreadsheet mentioned in section 1. The lab sheet was developed using Microsoft Excel for Office 97, but should work on subsequent versions as well. 1 The Pseudo-Random Numbers You can download a spreadsheet from the web page mentioned above, or you can start a new spreadsheet. If you choose to start a new one, you will need a column of pseudo-random numbers for the simulation. Name one of the worksheets Rand and fill the range A1:A200 with a single formula =RAND(). Use Insert Name Define (or equivalent) to call this range Random. [Note: There are various problems with using Excel's pseudo-random number generator: Excel insists on recalculating all its random numbers each time you do anything; the generator is not very sophisticated and does not rate highly in comparison with other generators; it is hard to set seeds in such a way that a sequence of pseudo-random numbers can be reproduced. That is why the spreadsheet supplied on the web page uses a different approach to generate the pseudorandom numbers.] 2 The Simple Random Walk * 1 Start a new worksheet and call it "SRW". The simple random walk has a single parameter, p, so set aside a cell to hold the value, and name the cell p. To start off with, use the value 0.5. Column A will hold the values of the random walk, column B the increments (jumps). So, in cell B5, enter the formula =IF(Random<p, 1, -1), and copy this down into the next 50 cells. Remember that the IF statement takes the form IF(Condition, value if condition true, value if condition false). In A4 enter the starting value (0 will do), then in A5 the formula =A4+B5, which can again be copied down for 50 rows. While this column is still selected, click on Chart Wizard and go through the steps required to obtain a Line chart of the SRW. Change the value of p by a small amount such as ±0.01 and observe what happens to the SRW as shown on the chart. Since we are using the same seed all the time, changes are relatively small. Now change the seed on the Generator sheet (choose any number in the range 0 to 1) and see how the chart changes. 3 Barriers *We shall use columns C and D to construct another SRW, but this time we shall impose some barriers on it. We shall need to store the location and type of each barrier. Type labels UpperLocation, LowerLocation, UpperType and LowerType into cells D2, D3, F2 and F3, and use Insert Name Define to give the adjacent cells (E2, E3, G2, G3) these names. 1 * indicates that this has been done for you on the downloadable spreadsheet. Stochastic Modelling 2003 Session 1, page 1
2 Now decide where you want the barriers to be 0 and 5 would be simple to begin with, so enter these in E3 and E2 respectively. The barrier type will be a numerical value: 1 for a downward reflecting barrier, 0 for an absorbing barrier and 1 for an upward reflecting barrier. Choose whatever types you like and enter them in G2 and G3. The increments, which will be in column D starting from D5, are now generated using a nested IF statement: =IF(C4=UpperLocation,UpperType,IF(C4=LowerLocation,LowerType,IF(Random<p, 1, -1))). Check that you understand how this works, then enter the formula =C4+D5 in cell C5 and copy both formulas down for 50 rows. Select the values in the C column, press Copy, then click on the chart and choose Paste Special from the Edit menu to get the new data series added to the chart along with the original random walk. Try changing the barrier locations to see how much difference they make, both in the case where p = 0.5 and in the case where p is different from 0.5. See if you can predict what you are going to observe before you actually observe it. 4 General discrete distributions for increments Although Excel can be used to simulate from the Binomial and Poisson distributions, we are going to show how to use it to simulate from a more general discrete distribution. *To begin with, we need to set up a table listing the possible values and their probabilities. Start a new worksheet and name it "Discrete". At the top left corner of the new worksheet construct a region like this: Prob CumProb 0 =B1+B2 =C1+C2 =D1+D2 =E1+E2 =F1+F2 Value and give the name "ProbTable" to the range B2:F3. This is called a Lookup Table: we are going to produce random numbers between 0 and 1, look them up in the Cumulative Probability row and return the corresponding quantity from the Value column. The worksheet function to use is HLOOKUP. The syntax of the HLOOKUP function is =HLOOKUP(Lookup value, Lookup table, Row number), where the row number refers to the row of the table to look in, starting from the first row of the table which is counted as row 1. In this case the lookup values are in the first row of ProbTable (column B) and the corresponding values of the variable are in the second row (column C). Therefore a suitable increment for the random walk will be generated by the formula =HLOOKUP(Random, ProbTable, 2). Simulate this random walk for 50 steps and add the resulting trajectory to the diagram. 5 The transition matrix A general discrete-time Markov chain has transition probabilities which are likely to vary depending on the current state. It is not too hard to write a macro in Visual Basic to carry out the transitions, but it is more challenging to try to arrange a simulation mechanism using only the functions which are built into Excel. We shall use look-up tables to accomplish our goals. *Start a new worksheet and call it Markov. Enter a title for the sheet somewhere on row 1. We are going to simulate a no-claims discount scheme with 4 states: No discount, small discount, medium discount and large discount. The transition mechanism is: if there are no claims during a year, the Stochastic Modelling 2003 Session 1, page 2
3 policyholder moves to the next higher level of discount (unless already receiving Large discount); if there is one or more claim during a year the policyholder moves to the next lower discount level (unless receiving no discount). *We need a parameter: type Claim frequency into cell A3, with some value like 0.4 in B3: this represents the average number of claims per year. The Poisson distribution will be used, so the probability of no claims, q, is given by the formula =EXP(-B3): enter this formula in B4, with the label q in A4, and name the cell q (using Insert Name Define). Now enter the transition matrix into a blank area on the worksheet, starting from column E. The diagram below contains labels in column D: they are not necessary, but serve to remind us which state is which. The procedure outlined here only works if the states are numbered with integers starting from 0. Don t forget that every entry must be non-negative and that the row sums must be 1: it s probably best to make the last column contain 1 (sum of previous columns). D E F G H I J 3 No discount 0 =1-q =q 0 =1-SUM(F3:H3) 4 Small disc 1 =1-q 0 =q =1-SUM(F4:H4) 5 Med disc 2 0 =1-q 0 =1-SUM(F5:H5) 6 Large disc =1-q =1-SUM(F6:H6) The lookup function needs a cumulative transition matrix: the table shown below will calculate these in the form required. Use Insert Name Define to name the cumulative matrix (in this example, F8:J11) CumMatrix, and the vector of states (F12:J12) StateVector. D E F G H I J 7 8 No discount 0 0 =F3 =G8+G3 =H8+H3 =I8+I3 9 Small disc 1 0 =F4 =G9+G4 =H9+H4 =I9+I4 10 Med disc 2 0 =F5 =G10+G5 =H10+H5 =I10+I5 11 Large disc 3 0 =F6 =G11+G6 =H11+H6 =I11+I6 12 =E8 =E9 =E10 =E The LOOKUP function The LOOKUP function has syntax LOOKUP(lookup_value, lookup_vector, result_vector). Here the lookup_value is going to be a random number between 0 and 1. We tell Excel which row of the cumulative transition matrix to consult (the row corresponding to the current state, i). When it finds the largest entry in that row which does not exceed the lookup_value it notes how many columns it had to search, then sends back the corresponding entry from the result_vector: this will be the state j to which the chain jumps next. If this is hard to follow, type in the formulas below and see if the practical application makes it clearer. Type column labels n in A6, X n in B6, then 0 in A7 and some initial state (eg., 0) in B7. Row 8 should contain =1+A7 in A8 and =LOOKUP(Random, OFFSET(CumMatrix,B7,0,1,5), StateVector) in B8. It is this last entry which is doing all the work. The OFFSET function directs the search to the appropriate row (B7 contains the row number, where the top row is always numbered 0) of the matrix CumMatrix; the parameters 1 and 5 refer to the number of rows (1) and the number of columns (5) which are available to be consulted. Once the function has found the appropriate entry it returns the corresponding value from StateVector. Stochastic Modelling 2003 Session 1, page 3
4 A8:B8 can be copied and pasted into as many rows as you like to simulate a sequence of consecutive values of the chain. Obtain a plot of the path of the Markov chain using the Line option of Chart Wizard. Give the column of observed values of the Markov chain (B7:Bwhatever) the name Chain. 7 Observed average discount We are going to find the observed proportion of time spent in states 0, 1, 2 and 3. Type 0 in F14, 1 in G14, 2 in H14, 3 in I14. Then in F15 enter the formula =COUNTIF(Chain,"=" & F14)/COUNT(Chain) to find the observed proportion of time spent in state 0. Copy the formula into G15:I15. Now think up some figures for the percentage discount at each discount level, starting from discount 0 in state 0, and enter them in F16:I16. To calculate the average discount achieved by the simulated policyholder over the period of the simulation, all we need is the SUMPRODUCT function, which essentially does a dot product of two vectors. Thus =SUMPRODUCT(F15:I15,F16:I16) gives the required answer for the observed average discount. 8 Array functions Excel has three built-in matrix functions, which can speed up a matrix-type calculation immensely but are a little tricky to type in. The functions are MDETERM, MINVERSE and MMULT. We are going to start off this section by using MMULT to produce successive calculations of the probability that the chain is in state i at time n starting from state 0. Theoretically, if v n is the vector of probabilities given by v n,i = P(X n =i), then v n satisfies the recurrence relation v n =v n-1 P. We are going to need to multiply a row vector v by the transition matrix P. Select the whole transition matrix (F3:I6) and give it a name, such as Transition_matrix. Now type in the initial vector of probabilities. Assuming that the motorist starts with no discount, the vector will be ( ): enter these four values into cells F18:I18. We have to multiply this by the matrix P to get the probabilities of being in the various states at the next time point. So select the cells F19:I19 to store the results of the calculation, type in the formula =MMULT(F18:I18,Transition_matrix) but DO NOT PRESS ENTER press Ctrl-Shift-Enter instead. All cells of the new array (F19:I19) will be filled with the same formula, shown inside braces (curly brackets). The Ctrl-Shift-Enter method is the way all arrays are entered. If you make a mistake in entering the formula you will have to clear the whole array at once, then reselect the whole area to have another go. Copy this formula down into several more blank rows by dragging the drag-handle at the bottom right corner of the selected area. Watch the probabilities converge to their long-run limits π. 9 Matrix inversion Let us try to automate the procedure used to calculate π. The first thing to bear in mind is that we need to solve πp = π, or in other words π(i P) = 0. To start with, then, let us get the identity matrix into the worksheet. Stochastic Modelling 2003 Session 1, page 4
5 Choose a 5 by 5 area such as L3:P7. Enter the names of the states (0, 1, 2, 3 in this case) into the first column (L3:L6) and the last row (M7:P7). We can fill the matrix using a single array function: select the matrix area M3:P6, type in the formula =IF(L3:L6=M7:P7,1,0) and press Ctrl-Shift-Enter. (This function returns a value of 1 if the row is equal to the column, 0 otherwise.) Immediately below this, in M8:P11, we can now enter the matrix I P: just use a simple subtraction formula like =M3-F3 and copy it across the whole of M8:P11. (Don't use Ctrl-Shift-Enter here as we will want to change the last column later on and Excel will not let you change just part of an array.) We know, however, that the matrix I P is not invertible, that there are multiple solutions to the equation π(i P) = 0, that the last equation in the set is redundant and that the solution we want has Σπ i = 1. This means that we need to replace the last column of the matrix I P by a column of 1's, resulting in a matrix M, say, then solving the equation πm = ( ), which is to say π = ( ) M 1. So you need to replace the last column of the matrix I P (P8:P11) by a column of 1's. Now select the region M13:P16, type in the formula =MINVERSE(M8:P11) and press Ctrl-Shift-Enter: as if by magic, the inverse of the matrix M appears in the selected area. You are probably expecting that we will now have to type in ( ) and perform a matrix multiplication to obtain the equilibrium probability vector π. In fact, though, a little thought will reveal that ( ) M 1 is simply the bottom row of the matrix M 1, so in fact the equilibrium distribution is now sitting in M16:P Multiple runs Suppose we are interested in a quantity whose distribution cannot easily be calculated, such as the position at time 50 of the random walk whose increments follow the general discrete distribution introduced in section 5. Even if we can't find the exact distribution by calculation, we can simulate the process a few times and calculate the sample mean, or draw a histogram of the results, or something similar. The Excel feature which we shall use is called a Data Table. Its purpose is to allow the user to perform a long calculation several times, using several different values for an input variable, without having to make a huge number of copies of the worksheet. In our case the input variable we want to change is the seed, and the output variable we want to study is the position of the RW at time 50. The one drawback of using the Data Table is that it requires both the input and output variables to be located on the same worksheet as the Data Table itself. So let us start by inserting a new worksheet to hold the data table, name it "DataTable", and write labels Input and Output near the top. The value in the Output cell must be equal to the position of the RW at time 50, so a formula like =SRW!E54 will do the job OK. The Input cell works the other way round: we need to ensure that a value entered in the Input cell is then copied across into the Seed cell on the Generator sheet. This involves replacing the Seed value with the formula =DataTable!B2, making sure that the Input cell contains a suitable value for the seed. Now set up the data table on the DataTable worksheet: in cells A5 to A20, to hold the different seeds, enter the formula =RAND(). Then in B4 use the formula =B1, so that the output value appears there. Next, select the range from A4 to B20, choose Data Table; leaving the Row Input Cell blank, enter $B$2 in the Column Input Cell and press Enter. Stochastic Modelling 2003 Session 1, page 5
6 The right-hand column of the table now fills with the values which the output cell would contain if the input cell contained the value shown in the left-hand column. Since the input cell contains only the seed for the random number generator, these are the output values for a number of independent runs of the simulation. You can construct a histogram or bar chart of the results, or investigate sample mean and variance. It goes without saying that a 'proper' simulation requires a good many more runs than 16, but naturally the data table can be made as large as you like, within the limits imposed by Excel (maximum is 65,000 or so). Once you have a collection of independently-generated output values you can analyse it in the way you would analyse any data set: find the sample mean and variance, plot a histogram of the output values, test for normality, etc. Stochastic Modelling 2003 Session 1, page 6
COMPUTING AND DATA ANALYSIS WITH EXCEL. Matrix manipulation and systems of linear equations
COMPUTING AND DATA ANALYSIS WITH EXCEL Matrix manipulation and systems of linear equations Outline 1 Matrices Addition Subtraction Excel functions that return more than one cell Solving systems of linear
More informationUsing Microsoft Excel
Using Microsoft Excel Objective: Students will gain familiarity with using Excel to record data, display data properly, use built-in formulae to do calculations, and plot and fit data with linear functions.
More informationSome hints for the Radioactive Decay lab
Some hints for the Radioactive Decay lab Edward Stokan, March 7, 2011 Plotting a histogram using Microsoft Excel The way I make histograms in Excel is to put the bounds of the bin on the top row beside
More informationPhysics 201 Lab 2 Air Drag Simulation
Physics 201 Lab 2 Air Drag Simulation Jan 28, 2013 Equipment Initial Set Up Type the data from Table 1 into the appropriate cells. By preceding the content of the cell with an equal sign (as in cell A6)
More informationDISCRETE RANDOM VARIABLES EXCEL LAB #3
DISCRETE RANDOM VARIABLES EXCEL LAB #3 ECON/BUSN 180: Quantitative Methods for Economics and Business Department of Economics and Business Lake Forest College Lake Forest, IL 60045 Copyright, 2011 Overview
More informationSome Excel Problems ELEMENTARY MATHEMATICS FOR BIOLOGISTS 2013
ELEMENTARY MATHEMATICS FOR BIOLOGISTS 2013 Some Excel Problems It should be possible to solve all the problems presented in this handout by using the information given in the various sessions which are
More informationExperiment: Oscillations of a Mass on a Spring
Physics NYC F17 Objective: Theory: Experiment: Oscillations of a Mass on a Spring A: to verify Hooke s law for a spring and measure its elasticity constant. B: to check the relationship between the period
More information9. Using Excel matrices functions to calculate partial autocorrelations
9 Using Excel matrices functions to calculate partial autocorrelations In order to use a more elegant way to calculate the partial autocorrelations we need to borrow some equations from stochastic modelling,
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 informationBoyle s Law and Charles Law Activity
Boyle s Law and Charles Law Activity Introduction: This simulation helps you to help you fully understand 2 Gas Laws: Boyle s Law and Charles Law. These laws are very simple to understand, but are also
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 information8: Statistical Distributions
: Statistical Distributions The Uniform Distribution 1 The Normal Distribution The Student Distribution Sample Calculations The Central Limit Theory Calculations with Samples Histograms & the Normal Distribution
More informationHow many states. Record high temperature
Record high temperature How many states Class Midpoint Label 94.5 99.5 94.5-99.5 0 97 99.5 104.5 99.5-104.5 2 102 102 104.5 109.5 104.5-109.5 8 107 107 109.5 114.5 109.5-114.5 18 112 112 114.5 119.5 114.5-119.5
More informationFigure 1. Distillation Train. Table 1. Stream compositions.
CM 3450 Drills 2 9/8/2010 1. A stream containing compounds,, and are fed to a series of distillation columns as shown in Figure 1 with corresponding stream compositions given in Table 1. Figure 1. Distillation
More informationLinear Algebra, Vectors and Matrices
Linear Algebra, Vectors and Matrices Prof. Manuela Pedio 20550 Quantitative Methods for Finance August 2018 Outline of the Course Lectures 1 and 2 (3 hours, in class): Linear and non-linear functions on
More informationIEOR 6711: Professor Whitt. Introduction to Markov Chains
IEOR 6711: Professor Whitt Introduction to Markov Chains 1. Markov Mouse: The Closed Maze We start by considering how to model a mouse moving around in a maze. The maze is a closed space containing nine
More informationContents. 13. Graphs of Trigonometric Functions 2 Example Example
Contents 13. Graphs of Trigonometric Functions 2 Example 13.19............................... 2 Example 13.22............................... 5 1 Peterson, Technical Mathematics, 3rd edition 2 Example 13.19
More informationConformational Analysis of n-butane
Conformational Analysis of n-butane In this exercise you will calculate the Molecular Mechanics (MM) single point energy of butane in various conformations with respect to internal rotation around the
More information1. The Basic X-Y Scatter Plot
1. The Basic X-Y Scatter Plot EXCEL offers a wide range of plots; however, this discussion will be restricted to generating XY scatter plots in various formats. The easiest way to begin is to highlight
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 informationA Scientific Model for Free Fall.
A Scientific Model for Free Fall. I. Overview. This lab explores the framework of the scientific method. The phenomenon studied is the free fall of an object released from rest at a height H from the ground.
More informationSimple Linear Regression
CHAPTER 13 Simple Linear Regression CHAPTER OUTLINE 13.1 Simple Linear Regression Analysis 13.2 Using Excel s built-in Regression tool 13.3 Linear Correlation 13.4 Hypothesis Tests about the Linear Correlation
More informationHow to Make or Plot a Graph or Chart in Excel
This is a complete video tutorial on How to Make or Plot a Graph or Chart in Excel. To make complex chart like Gantt Chart, you have know the basic principles of making a chart. Though I have used Excel
More informationModule 17. Diffusion in solids III. Lecture 17. Diffusion in solids III
Module 17 Diffusion in solids III Lecture 17 Diffusion in solids III 1 NPTEL Phase II : IIT Kharagpur : Prof. R. N. Ghosh, Dept of Metallurgical and Materials Engineering Keywords: Numerical problems in
More informationEXPERIMENT 2 Reaction Time Objectives Theory
EXPERIMENT Reaction Time Objectives to make a series of measurements of your reaction time to make a histogram, or distribution curve, of your measured reaction times to calculate the "average" or mean
More informationLab 2 Worksheet. Problems. Problem 1: Geometry and Linear Equations
Lab 2 Worksheet Problems Problem : Geometry and Linear Equations Linear algebra is, first and foremost, the study of systems of linear equations. You are going to encounter linear systems frequently in
More informationM61 1 M61.1 PC COMPUTER ASSISTED DETERMINATION OF ANGULAR ACCELERATION USING TORQUE AND MOMENT OF INERTIA
M61 1 M61.1 PC COMPUTER ASSISTED DETERMINATION OF ANGULAR ACCELERATION USING TORQUE AND MOMENT OF INERTIA PRELAB: Before coming to the lab, you must write the Object and Theory sections of your lab report
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Experiment 03: Work and Energy
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.01 Fall Term 2010 Experiment 03: Work and Energy Purpose of the Experiment: In this experiment you allow a cart to roll down an inclined
More informationLesson 24: True and False Number Sentences
NYS COMMON CE MATHEMATICS CURRICULUM Lesson 24 6 4 Student Outcomes Students identify values for the variables in equations and inequalities that result in true number sentences. Students identify values
More informationGravity Modelling Forward Modelling Of Synthetic Data
Gravity Modelling Forward Modelling Of Synthetic Data After completing this practical you should be able to: The aim of this practical is to become familiar with the concept of forward modelling as a tool
More information1. AN INTRODUCTION TO DESCRIPTIVE STATISTICS. No great deed, private or public, has ever been undertaken in a bliss of certainty.
CIVL 3103 Approximation and Uncertainty J.W. Hurley, R.W. Meier 1. AN INTRODUCTION TO DESCRIPTIVE STATISTICS No great deed, private or public, has ever been undertaken in a bliss of certainty. - Leon Wieseltier
More informationExperiment 0 ~ Introduction to Statistics and Excel Tutorial. Introduction to Statistics, Error and Measurement
Experiment 0 ~ Introduction to Statistics and Excel Tutorial Many of you already went through the introduction to laboratory practice and excel tutorial in Physics 1011. For that reason, we aren t going
More informationIn this activity, students will compare weather data from to determine if there is a warming trend in their community.
Overview: In this activity, students will compare weather data from 1910-2000 to determine if there is a warming trend in their community. Objectives: The student will: use the Internet to locate scientific
More informationNCSS Statistical Software. Harmonic Regression. This section provides the technical details of the model that is fit by this procedure.
Chapter 460 Introduction This program calculates the harmonic regression of a time series. That is, it fits designated harmonics (sinusoidal terms of different wavelengths) using our nonlinear regression
More informationExperiment 4 Free Fall
PHY9 Experiment 4: Free Fall 8/0/007 Page Experiment 4 Free Fall Suggested Reading for this Lab Bauer&Westfall Ch (as needed) Taylor, Section.6, and standard deviation rule ( t < ) rule in the uncertainty
More informationAn area chart emphasizes the trend of each value over time. An area chart also shows the relationship of parts to a whole.
Excel 2003 Creating a Chart Introduction Page 1 By the end of this lesson, learners should be able to: Identify the parts of a chart Identify different types of charts Create an Embedded Chart Create a
More informationCE-300 Mathcad Tutorial
CE-00 Mathcad Tutorial The purposes of this tutorial are: (1) to get you acquainted with the Mathcad procedures and synta you will use to solve typical problems in CE-00 and (2) to demonstrate how to set
More informationExperiment 13. Dilutions and Data Handling in a Spreadsheet rev 1/2013
Absorbance Experiment 13 Dilutions and Data Handling in a Spreadsheet rev 1/2013 GOAL: This lab experiment will provide practice in making dilutions using pipets and introduce basic spreadsheet skills
More informationChapter 1 Review of Equations and Inequalities
Chapter 1 Review of Equations and Inequalities Part I Review of Basic Equations Recall that an equation is an expression with an equal sign in the middle. Also recall that, if a question asks you to solve
More informationA (Mostly) Correctly Formatted Sample Lab Report. Brett A. McGuire Lab Partner: Microsoft Windows Section AB2
A (Mostly) Correctly Formatted Sample Lab Report Brett A. McGuire Lab Partner: Microsoft Windows Section AB2 August 26, 2008 Abstract Your abstract should not be indented and be single-spaced. Abstracts
More informationStar Cluster Photometry and the H-R Diagram
Star Cluster Photometry and the H-R Diagram Contents Introduction Star Cluster Photometry... 1 Downloads... 1 Part 1: Measuring Star Magnitudes... 2 Part 2: Plotting the Stars on a Colour-Magnitude (H-R)
More informationChem 7040 Statistical Thermodynamics Problem Set #6 - Computational Due 1 Oct at beginning of class
Chem 7040 Statistical Thermodynamics Problem Set #6 - Computational Due 1 Oct at beginning of class Good morning, gang! I was coerced into giving a last-minute Chem Dept recruiting talk in balmy Minnesota
More informationRadiological Control Technician Training Fundamental Academic Training Study Guide Phase I
Module 1.01 Basic Mathematics and Algebra Part 4 of 9 Radiological Control Technician Training Fundamental Academic Training Phase I Coordinated and Conducted for the Office of Health, Safety and Security
More informationEXCELLING WITH BIOLOGICAL MODELS FROM THE CLASSROOM T0 RESEARCH
EXCELLING WITH BIOLOGICAL MODELS FROM THE CLASSROOM T0 RESEARCH Timothy D. Comar Benedictine University Department of Mathematics 5700 College Road Lisle, IL 60532 tcomar@ben.edu Introduction Computer
More informationChapter 9 Ingredients of Multivariable Change: Models, Graphs, Rates
Chapter 9 Ingredients of Multivariable Change: Models, Graphs, Rates 9.1 Multivariable Functions and Contour Graphs Although Excel can easily draw 3-dimensional surfaces, they are often difficult to mathematically
More informationGeology Geomath Computer Lab Quadratics and Settling Velocities
Geology 351 - Geomath Computer Lab Quadratics and Settling Velocities In Chapter 3 of Mathematics: A simple tool for geologists, Waltham takes us through a brief review of quadratic equations and their
More informationPhotometry of Supernovae with Makali i
Photometry of Supernovae with Makali i How to perform photometry specifically on supernovae targets using the free image processing software, Makali i This worksheet describes how to use photometry to
More informationMeasurement: The Basics
I. Introduction Measurement: The Basics Physics is first and foremost an experimental science, meaning that its accumulated body of knowledge is due to the meticulous experiments performed by teams of
More informationFalling Bodies (last
Dr. Larry Bortner Purpose Falling Bodies (last edited ) To investigate the motion of a body under constant acceleration, specifically the motion of a mass falling freely to Earth. To verify the parabolic
More informationACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER / Lines and Their Equations
ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER 1 017/018 DR. ANTHONY BROWN. Lines and Their Equations.1. Slope of a Line and its y-intercept. In Euclidean geometry (where
More informationAppendix B Microsoft Office Specialist exam objectives maps
B 1 Appendix B Microsoft Office Specialist exam objectives maps This appendix covers these additional topics: A Excel 2003 Specialist exam objectives with references to corresponding material in Course
More informationα m ! m or v T v T v T α m mass
FALLING OBJECTS (WHAT TO TURN IN AND HOW TO DO SO) In the real world, because of air resistance, objects do not fall indefinitely with constant acceleration. One way to see this is by comparing the fall
More informationChapter 14: Finding the Equilibrium Solution and Exploring the Nature of the Equilibration Process
Chapter 14: Finding the Equilibrium Solution and Exploring the Nature of the Equilibration Process Taking Stock: In the last chapter, we learned that equilibrium problems have an interesting dimension
More informationBoyle s Law: A Multivariable Model and Interactive Animated Simulation
Boyle s Law: A Multivariable Model and Interactive Animated Simulation Using tools available in Excel, we will turn a multivariable model into an interactive animated simulation. Projectile motion, Boyle's
More informationMATH 118 FINAL EXAM STUDY GUIDE
MATH 118 FINAL EXAM STUDY GUIDE Recommendations: 1. Take the Final Practice Exam and take note of questions 2. Use this study guide as you take the tests and cross off what you know well 3. Take the Practice
More informationElectric Fields and Equipotentials
OBJECTIVE Electric Fields and Equipotentials To study and describe the two-dimensional electric field. To map the location of the equipotential surfaces around charged electrodes. To study the relationship
More informationSME 864 Mark Urban-Lurain
SME 864 Mark Urban-Lurain 1 Import data from non-excel sources Probe Software Web sites Other sources Organize data Structure file for analysis Clean values Analyze Summarize Statistics Graph 2 Get files
More informationx y = 2 x + 2y = 14 x = 2, y = 0 x = 3, y = 1 x = 4, y = 2 x = 5, y = 3 x = 6, y = 4 x = 7, y = 5 x = 0, y = 7 x = 2, y = 6 x = 4, y = 5
List six positive integer solutions for each of these equations and comment on your results. Two have been done for you. x y = x + y = 4 x =, y = 0 x = 3, y = x = 4, y = x = 5, y = 3 x = 6, y = 4 x = 7,
More informationWEATHER AND CLIMATE COMPLETING THE WEATHER OBSERVATION PROJECT CAMERON DOUGLAS CRAIG
WEATHER AND CLIMATE COMPLETING THE WEATHER OBSERVATION PROJECT CAMERON DOUGLAS CRAIG Introduction The Weather Observation Project is an important component of this course that gets you to look at real
More informationStochastic Modelling Unit 1: Markov chain models
Stochastic Modelling Unit 1: Markov chain models Russell Gerrard and Douglas Wright Cass Business School, City University, London June 2004 Contents of Unit 1 1 Stochastic Processes 2 Markov Chains 3 Poisson
More informationLAB 5 INSTRUCTIONS LINEAR REGRESSION AND CORRELATION
LAB 5 INSTRUCTIONS LINEAR REGRESSION AND CORRELATION In this lab you will learn how to use Excel to display the relationship between two quantitative variables, measure the strength and direction of the
More informationRemember that C is a constant and ë and n are variables. This equation now fits the template of a straight line:
CONVERTING NON-LINEAR GRAPHS INTO LINEAR GRAPHS Linear graphs have several important attributes. First, it is easy to recognize a graph that is linear. It is much more difficult to identify if a curved
More informationThe Boundary Problem: Markov Chain Solution
MATH 529 The Boundary Problem: Markov Chain Solution Consider a random walk X that starts at positive height j, and on each independent step, moves upward a units with probability p, moves downward b units
More informationExperiment IV. To find the velocity of waves on a string by measuring the wavelength and frequency of standing waves.
Experiment IV The Vibrating String I. Purpose: To find the velocity of waves on a string by measuring the wavelength and frequency of standing waves. II. References: Serway and Jewett, 6th Ed., Vol., Chap.
More informationMATH 56A: STOCHASTIC PROCESSES CHAPTER 1
MATH 56A: STOCHASTIC PROCESSES CHAPTER. Finite Markov chains For the sake of completeness of these notes I decided to write a summary of the basic concepts of finite Markov chains. The topics in this chapter
More informationEureka Lessons for 6th Grade Unit FIVE ~ Equations & Inequalities
Eureka Lessons for 6th Grade Unit FIVE ~ Equations & Inequalities These 2 lessons can easily be taught in 2 class periods. If you like these lessons, please consider using other Eureka lessons as well.
More informationMEASUREMENT OF THE CHARGE TO MASS RATIO (e/m e ) OF AN ELECTRON
MEASUREMENT OF THE CHARGE TO MASS RATIO (e/m e ) OF AN ELECTRON Object This experiment will allow you to observe and understand the motion of a charged particle in a magnetic field and to measure the ratio
More informationStatistical Analysis of Data
Statistical Analysis of Data Goals and Introduction In this experiment, we will work to calculate, understand, and evaluate statistical measures of a set of data. In most cases, when taking data, there
More informationChapter 5 Simplifying Formulas and Solving Equations
Chapter 5 Simplifying Formulas and Solving Equations Look at the geometry formula for Perimeter of a rectangle P = L W L W. Can this formula be written in a simpler way? If it is true, that we can simplify
More informationIntroduction to Measurement Physics 114 Eyres
1 Introduction to Measurement Physics 114 Eyres 6/5/2016 Module 1: Measurement 1 2 Significant Figures Count all non-zero digits Count zeros between non-zero digits Count zeros after the decimal if also
More information17-Nov-2015 PHYS MAXWELL WHEEL. To test the conservation of energy in a system with gravitational, translational and rotational energies.
Objective MAXWELL WHEEL To test the conservation of energy in a system with gravitational, translational and rotational energies. Introduction A wheel is suspended by two cords wrapped on its axis. After
More informationMoving into the information age: From records to Google Earth
Moving into the information age: From records to Google Earth David R. R. Smith Psychology, School of Life Sciences, University of Hull e-mail: davidsmith.butterflies@gmail.com Introduction Many of us
More informationMath Circle at FAU 10/27/2018 SOLUTIONS
Math Circle at FAU 10/27/2018 SOLUTIONS 1. At the grocery store last week, small boxes of facial tissue were priced at 4 boxes for $5. This week they are on sale at 5 boxes for $4. Find the percent decrease
More informationASTR Astrometry of Asteroids Lab Exercise Due: March 4, 2011
Student 1 Name: Student 2 Name: Student 3 Name: NetID: NetID: NetID: ASTR 150 - Astrometry of Asteroids Lab Exercise Due: March 4, 2011 This lab exercise will be graded out of 100 points (point values
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 informationData Structures & Database Queries in GIS
Data Structures & Database Queries in GIS Objective In this lab we will show you how to use ArcGIS for analysis of digital elevation models (DEM s), in relationship to Rocky Mountain bighorn sheep (Ovis
More informationMEASUREMENT OF THE CHARGE TO MASS RATIO (e/m e ) OF AN ELECTRON
MEASUREMENT OF THE CHARGE TO MASS RATIO (e/m e ) OF AN ELECTRON Object This experiment will allow you to observe and understand the motion of a charged particle in a magnetic field and to measure the ratio
More information[Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty.]
Math 43 Review Notes [Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty Dot Product If v (v, v, v 3 and w (w, w, w 3, then the
More informationTo factor an expression means to write it as a product of factors instead of a sum of terms. The expression 3x
Factoring trinomials In general, we are factoring ax + bx + c where a, b, and c are real numbers. To factor an expression means to write it as a product of factors instead of a sum of terms. The expression
More informationPhysics Lab 202P-3. Electric Fields and Superposition: A Virtual Lab NAME: LAB PARTNERS:
Physics Lab 202P-3 Electric Fields and Superposition: A Virtual Lab NAME: LAB PARTNERS: LAB SECTION: LAB INSTRUCTOR: DATE: EMAIL ADDRESS: Penn State University Created by nitin samarth Physics Lab 202P-3
More informationLarge Scale Structure of the Universe Lab
Large Scale Structure of the Universe Lab Introduction: Since the mid-1980 s astronomers have gathered data allowing, for the first time, a view of the structure of the Universe in three-dimensions. You
More informationMath 132. Population Growth: Raleigh and Wake County
Math 132 Population Growth: Raleigh and Wake County S. R. Lubkin Application Ask anyone who s been living in Raleigh more than a couple of years what the biggest issue is here, and if the answer has nothing
More informationLooking hard at algebraic identities.
Looking hard at algebraic identities. Written by Alastair Lupton and Anthony Harradine. Seeing Double Version 1.00 April 007. Written by Anthony Harradine and Alastair Lupton. Copyright Harradine and Lupton
More informationPhysics Lab 202P-9. Magnetic Fields & Electric Current NAME: LAB PARTNERS:
Physics Lab 202P-9 Magnetic Fields & Electric Current NAME: LAB PARTNERS: LAB SECTION: LAB INSTRUCTOR: DATE: EMAIL ADDRESS: Penn State University Created by nitin samarth Physics Lab 202P-9 Page 1 of 22
More informationEXPERIMENT: REACTION TIME
EXPERIMENT: REACTION TIME OBJECTIVES to make a series of measurements of your reaction time to make a histogram, or distribution curve, of your measured reaction times to calculate the "average" or "mean"
More informationConservation of Momentum
Learning Goals Conservation of Momentum After you finish this lab, you will be able to: 1. Use Logger Pro to analyze video and calculate position, velocity, and acceleration. 2. Use the equations for 2-dimensional
More informationSeptember 12, Math Analysis Ch 1 Review Solutions. #1. 8x + 10 = 4x 30 4x 4x 4x + 10 = x = x = 10.
#1. 8x + 10 = 4x 30 4x 4x 4x + 10 = 30 10 10 4x = 40 4 4 x = 10 Sep 5 7:00 AM 1 #. 4 3(x + ) = 5x 7(4 x) 4 3x 6 = 5x 8 + 7x CLT 3x = 1x 8 +3x +3x = 15x 8 +8 +8 6 = 15x 15 15 x = 6 15 Sep 5 7:00 AM #3.
More informationOUTLINE. Deterministic and Stochastic With spreadsheet program : Integrated Mathematics 2
COMPUTER SIMULATION OUTLINE In this module, we will focus on the act simulation, taking mathematical models and implement them on computer systems. Simulation & Computer Simulations Mathematical (Simulation)
More informationWISE Regression/Correlation Interactive Lab. Introduction to the WISE Correlation/Regression Applet
WISE Regression/Correlation Interactive Lab Introduction to the WISE Correlation/Regression Applet This tutorial focuses on the logic of regression analysis with special attention given to variance components.
More informationLAB Exercise #4 - Answers The Traction Vector and Stress Tensor. Introduction. Format of lab. Preparation reading
LAB Exercise #4 - Answers The Traction Vector and Stress Tensor Due: Thursday, 26 February 2009 (Special Thanks to D.D. Pollard who pioneered this exercise in 1991) Introduction Stress concentrations in
More informationRESERVE DESIGN INTRODUCTION. Objectives. In collaboration with Wendy K. Gram. Set up a spreadsheet model of a nature reserve with two different
RESERVE DESIGN In collaboration with Wendy K. Gram Objectives Set up a spreadsheet model of a nature reserve with two different habitats. Calculate and compare abundances of species with different habitat
More informationStudying Topography, Orographic Rainfall, and Ecosystems (STORE)
Studying Topography, Orographic Rainfall, and Ecosystems (STORE) Introduction Basic Lesson 3: Using Microsoft Excel to Analyze Weather Data: Topography and Temperature This lesson uses NCDC data to compare
More information4.7.1 Computing a stationary distribution
At a high-level our interest in the rest of this section will be to understand the limiting distribution, when it exists and how to compute it To compute it, we will try to reason about when the limiting
More informationLAB 2 - ONE DIMENSIONAL MOTION
Name Date Partners L02-1 LAB 2 - ONE DIMENSIONAL MOTION OBJECTIVES Slow and steady wins the race. Aesop s fable: The Hare and the Tortoise To learn how to use a motion detector and gain more familiarity
More informationCOUNCIL ROCK HIGH SCHOOL MATHEMATICS. A Note Guideline of Algebraic Concepts. Designed to assist students in A Summer Review of Algebra
COUNCIL ROCK HIGH SCHOOL MATHEMATICS A Note Guideline of Algebraic Concepts Designed to assist students in A Summer Review of Algebra [A teacher prepared compilation of the 7 Algebraic concepts deemed
More informationThe Celsius temperature scale is based on the freezing point and the boiling point of water. 12 degrees Celsius below zero would be written as
Prealgebra, Chapter 2 - Integers, Introductory Algebra 2.1 Integers In the real world, numbers are used to represent real things, such as the height of a building, the cost of a car, the temperature of
More informationMITOCW ocw f99-lec30_300k
MITOCW ocw-18.06-f99-lec30_300k OK, this is the lecture on linear transformations. Actually, linear algebra courses used to begin with this lecture, so you could say I'm beginning this course again by
More information25. Strassen s Fast Multiplication of Matrices Algorithm and Spreadsheet Matrix Multiplications
25.1 Introduction 25. Strassen s Fast Multiplication of Matrices Algorithm and Spreadsheet Matrix Multiplications We will use the notation A ij to indicate the element in the i-th row and j-th column of
More informationFree Fall. v gt (Eq. 4) Goals and Introduction
Free Fall Goals and Introduction When an object is subjected to only a gravitational force, the object is said to be in free fall. This is a special case of a constant-acceleration motion, and one that
More informationEXERCISE 13: SINGLE-SPECIES, MULTIPLE-SEASON OCCUPANCY MODELS
EXERCISE 13: SINGLE-SPECIES, MULTIPLE-SEASON OCCUPANCY MODELS Please cite this work as: Donovan, T. M. and J. Hines. 2007. Exercises in occupancy modeling and estimation.
More information