Leslie matrices and Markov chains.
|
|
- Gilbert Damian Waters
- 5 years ago
- Views:
Transcription
1 Leslie matrices and Markov chains. Example. Suppose a certain species of insect can be divided into 2 classes, eggs and adults. 10% of eggs survive for 1 week to become adults, each adult yields an average of 100 eggs per week, and 60% of adults survive an additional week. If a population of these insects currently consists of 100 eggs and 50 adults, find the state of the population one week from now, two weeks from now, and three weeks from now. Solution (by hand)..6(50) = 30 adults will survive, and.10(100) = 10 eggs will survive to adulthood, giving a total of 40 adults. The 50 adults will produce 50(100) = 5000 eggs. Better solution (using matrices). Let 0 100, > L = matrix(c(0,100,.10,.60), nrow=2, ncol=2, byrow=true) > L [,2] [1,] [2,] > x = matrix(c(100,50), nrow=2, ncol=1) > x [1,] 100 [2,] 50 > L %*% x [1,] 5000 [2,] 40 Not surprisingly, gives this information. The matrix L is a so called Leslie matrix. Since gives the state of the population after 7 days, will give its state 7 days after that. > L %*% (L %*% x) [1,] 4000 [2,] 524 After 14 days, we ll have 4000 eggs, and 524 adults.
2 Of course, this gets a little clumsy after a while. We d like to be able to calculate, but > L^2 %*% x [1,] [2,] 19 obviously is wrong. This is because L^2 computes L*L, and not the matrix product L %*% L. Fortunately, there is a package we can install that remedies this small annoyance. If you select Packages Load Package from the taskbar, load the expm (matrix exponential) package, then from the command line, type > library(expm) Then, > (L %^% 2) %*% x [1,] 4000 [2,] 524 gives the proper result. In general, M %^% n will compute the n th power of matrix M. To find the state of the population after 21 days, > (L %^% 3) %*% x [1,] [2,] does the trick. It should now be clear that this population, left unchecked, is about to explode. Example. See the previous problem. But, suppose instead each adult produces only 5 eggs per week. Find the state of the population after 7, 14, and 21 days. Solution. The calculations > L = matrix(c(0,2,.10,.60), nrow=2, ncol=2, byrow=true) > L %*% x [1,] 0 [2,] 10
3 > (L %^% 2) %*% x [1,] 20 [2,] 6 > (L %^% 3) %*% x [1,] 12.0 [2,] 5.6 > (L %^% 4) %*% x [1,] [2,] 4.56 > (L %^% 5) %*% x [1,] [2,] > (L %^% 6) %*% x [1,] [2,] > (L %^% 7) %*% x [1,] [2,] > (L %^% 8) %*% x [1,] [2,] suggests that the population is slowly going to die out. Markov chains. In the above example, there were two possible states of insect (eggs, adults), and there was a certain probability (possibly zero!) of transitioning from one state to the other. In general, given a collection of states", and a collection of probabilities that give the probabilities of transitioning to and from the states, we call this a Markov Chain (or a Markov model).
4 Example. A small downtown area consists of four street corners. Corner 1 Corner 2 BAR HOME A drunken man is wandering the area. From corner 1, it is equally likely that his next move will be to either corner 2 or the bar. From corner 2, it is equally likely that his next move will be to corner 1 or home. Once the man reaches either the bar or his home, he will remain there with probability 1. We create a matrix M where column i gives the probability of transitioning from state i each of the other states in a single move. Let states 1, 2, 3, and 4 be corner 1, corner 2, the bar, and home respectively Note! If you encounter this topic again someday, you might see these matrices defined a bit differently, where the rows (not the columns!) tell you the probability of moving from a given state to another. In most applications, it makes little difference but it s something to be aware of. In general, if you just swap rows and columns, one says you ve transposed the matrix in question. Supposing the man begins at corner 1, let Calculate and interpret the meanings of a. b. c. > M = matrix(c(0,.5,0,0,.5,0,0,0,.5,0,1,0,0,.5,0,1), nrow=4, ncol=4, byrow=true) > M [,2] [,3] [,4] [1,] [2,] [3,] [4,] > x = matrix(c(1,0,0,0), nrow=4, ncol=1) > x
5 [1,] 1 [2,] 0 [3,] 0 [4,] 0 > M %*% x [1,] 0.0 [2,] 0.5 [3,] 0.5 [4,] 0.0 Not surprisingly, if you begin in state 1, then after 1 move, there is a 50% probability of being in state 2, and a 50% probability of being in state 3. > (M %^% 2) %*% x [1,] 0.25 [2,] 0.00 [3,] 0.50 [4,] 0.25 Note that if you begin in state 1, there is a 25% probability you return to state 1 in 2 moves (move to 2 with probability.5 and back to 1 in the next move with probability.5). 0% probability of ending up in state 2 after 2 moves (check the assumptions to verify that this is impossible). 50% probability of ending up in state 3 (the bar) after 2 moves (you could go directly there on move 1 with probability.5, and if you move to state 2, it will be impossible to return to the bar on the next move). 25% probability of being home after 2 moves (move to state 2 then home is the only route). Thus it should not surprise you that will tell you the probability of ending up in each of the states after n moves. > (M %^% 100) %*% x [1,] e 31 [2,] e+00 [3,] e 01 [4,] e 01
6 Note that after 100 moves, the probability of being in state 1 is for all practical purposes 0, and the probabilities of ending up at the bar or at home are approximately 2/3 and 1/3 respectively. We d actually say the vector [0,0,2/3,1/3] is an eigenvector of the transition matrix. Example. Consider a population that contains two types of individuals, type A s and type B s. During each move, we do two things: 1. Pick one member of the population at random to die. 2. Pick one member of the population at random that is chosen to reproduce (possibly the same one as in step one!), making an identical copy of itself. This selection process is often referred to as a Moran Process. Let the states be the current number of A s in the population (so the states are 0, 1, 2 and 3). Write down the transition matrix for this process. Column 1 (from the state with 0 A s). If there are no A s, none can appear, so we ll remain in this state with probability 1. Column 2 (from the state with 1 A). Here, we start with the population looking like ABB. Three things are possible. After one move, there could be 0, 1, or 2 A s. P(0 A s after 1 move) = P(an A is selected to die) x P(a B is chosen to reproduce) = 1/3 x 2/3 = 2/9 P(2 A s after 1 move) = P(a B is selected to die) x P(an A is chosen to reproduce) = 2/3 x 1/3 = 2/9 P(1 A after 1 move) = 1 (2/9 + 2/9) = 1 4/9 = 5/9 Column 3 (from the state with 2 A s). Here, we start with the population looking like AAB. Three things are possible. After one move, there could be 1, 2, or 3 A s. P(1 A) = P(an A is selected to die) x P(a B is chosen to reproduce) = 2/3 x 1/3 = 2/9 P(3 A) = P(a B is selected to die) x P(an A is chosen to reproduce) = 1/3 x 2/3 = 2/9 P(2 A) = 5/9 Column 4 (from the state with 3 A s). Here, the only possibility is that we will continue to have 3 A s (one will die, but they ll be replaced with another A). 1 2/ /9 2/9 0 Thus, the matrix of this process is 0 2/9 5/ /9 1 In R,
7 > N = matrix(c(1,0,0,0,2/9,5/9,2/9,0,0,2/9,5/9,2/9,0,0,0,1), nrow=4, ncol=4) > N [,2] [,3] [,4] [1,] [2,] [3,] [4,] Note that we let byrow default to FALSE, and we filled the matrix by columns. As before, if you calculate a large power of N, that will tell you what the long term state of this system will be. > N %^% 100 [,2] [,3] [,4] [1,] e e 01 0 [2,] e e 12 0 [3,] e e 12 0 [4,] e e 01 1 For something a little cleaner to look at, let s round all the entries of this matrix to 3 decimal places. > round (N %^% 100, 3) [,2] [,3] [,4] [1,] [2,] [3,] [4,] If we start with a population of zero or three A s, we ll stay in that state with probability 1. Starting with 1 A, we ll eventually end up with 0 A s with probability 2/3, and 3 A s with probability 1/3. Starting with 2 A s, the A s will disappear with probability 1/3, and the B s will vanish with probability 2/3. In the homework, one of the exercises will explore what happens if either A or B has a reproductive advantage.
Eigenvalues and eigenvectors.
Eigenvalues and eigenvectors. Example. A population of a certain animal can be broken into three classes: eggs, juveniles, and adults. Eggs take one year to become juveniles, and juveniles take one year
More informationA simple dynamic model of labor market
1 A simple dynamic model of labor market Let e t be the number of employed workers, u t be the number of unemployed worker. An employed worker has probability of 0.8 to be employed in the next period,
More informationLecture 10: Powers of Matrices, Difference Equations
Lecture 10: Powers of Matrices, Difference Equations Difference Equations A difference equation, also sometimes called a recurrence equation is an equation that defines a sequence recursively, i.e. each
More informationMATH 310, REVIEW SHEET 2
MATH 310, REVIEW SHEET 2 These notes are a very short summary of the key topics in the book (and follow the book pretty closely). You should be familiar with everything on here, but it s not comprehensive,
More information8 Square matrices continued: Determinants
8 Square matrices continued: Determinants 8.1 Introduction Determinants give us important information about square matrices, and, as we ll soon see, are essential for the computation of eigenvalues. You
More informationMatrix-Vector Products and the Matrix Equation Ax = b
Matrix-Vector Products and the Matrix Equation Ax = b A. Havens Department of Mathematics University of Massachusetts, Amherst January 31, 2018 Outline 1 Matrices Acting on Vectors Linear Combinations
More informationLecture 3: Special Matrices
Lecture 3: Special Matrices Feedback of assignment1 Random matrices The magic matrix commend magic() doesn t give us random matrix. Random matrix means we will get different matrices each time when we
More informationMarkov Chains and Pandemics
Markov Chains and Pandemics Caleb Dedmore and Brad Smith December 8, 2016 Page 1 of 16 Abstract Markov Chain Theory is a powerful tool used in statistical analysis to make predictions about future events
More informationMath 2J Lecture 16-11/02/12
Math 2J Lecture 16-11/02/12 William Holmes Markov Chain Recap The population of a town is 100000. Each person is either independent, democrat, or republican. In any given year, each person can choose to
More informationAnswers in blue. If you have questions or spot an error, let me know. 1. Find all matrices that commute with A =. 4 3
Answers in blue. If you have questions or spot an error, let me know. 3 4. Find all matrices that commute with A =. 4 3 a b If we set B = and set AB = BA, we see that 3a + 4b = 3a 4c, 4a + 3b = 3b 4d,
More informationTMA 4265 Stochastic Processes Semester project, fall 2014 Student number and
TMA 4265 Stochastic Processes Semester project, fall 2014 Student number 730631 and 732038 Exercise 1 We shall study a discrete Markov chain (MC) {X n } n=0 with state space S = {0, 1, 2, 3, 4, 5, 6}.
More informationLecture 3: Sizes of Infinity
Math/CS 20: Intro. to Math Professor: Padraic Bartlett Lecture 3: Sizes of Infinity Week 2 UCSB 204 Sizes of Infinity On one hand, we know that the real numbers contain more elements than the rational
More informationTopic 15 Notes Jeremy Orloff
Topic 5 Notes Jeremy Orloff 5 Transpose, Inverse, Determinant 5. Goals. Know the definition and be able to compute the inverse of any square matrix using row operations. 2. Know the properties of inverses.
More informationHomework Set #1 Solutions
Homework Set #1 Solutions Exercises 1.1 (p. 10) Assignment: Do #33, 34, 1, 3,, 29-31, 17, 19, 21, 23, 2, 27 33. (a) True. (p. 7) (b) False. It has five rows and six columns. (c) False. The definition given
More informationMath 4310 Solutions to homework 1 Due 9/1/16
Math 0 Solutions to homework Due 9//6. An element [a] Z/nZ is idempotent if [a] 2 [a]. Find all idempotent elements in Z/0Z and in Z/Z. Solution. First note we clearly have [0] 2 [0] so [0] is idempotent
More information2 Systems of Linear Equations
2 Systems of Linear Equations A system of equations of the form or is called a system of linear equations. x + 2y = 7 2x y = 4 5p 6q + r = 4 2p + 3q 5r = 7 6p q + 4r = 2 Definition. An equation involving
More informationMatrix exponential. R study group 2010
Matrix exponential R study group 2010 Niels Richard Hansen University of Copenhagen August 11, 2010 The objective is to compute the matrix exponential e X for a square p p matrix efficiently and reliably
More informationLTCC. Exercises. (1) Two possible weather conditions on any day: {rainy, sunny} (2) Tomorrow s weather depends only on today s weather
1. Markov chain LTCC. Exercises Let X 0, X 1, X 2,... be a Markov chain with state space {1, 2, 3, 4} and transition matrix 1/2 1/2 0 0 P = 0 1/2 1/3 1/6. 0 0 0 1 (a) What happens if the chain starts in
More informationMATH Mathematics for Agriculture II
MATH 10240 Mathematics for Agriculture II Academic year 2018 2019 UCD School of Mathematics and Statistics Contents Chapter 1. Linear Algebra 1 1. Introduction to Matrices 1 2. Matrix Multiplication 3
More informationDynamic Programming: Matrix chain multiplication (CLRS 15.2)
Dynamic Programming: Matrix chain multiplication (CLRS.) The problem Given a sequence of matrices A, A, A,..., A n, find the best way (using the minimal number of multiplications) to compute their product.
More informationElementary Linear Algebra
Matrices J MUSCAT Elementary Linear Algebra Matrices Definition Dr J Muscat 2002 A matrix is a rectangular array of numbers, arranged in rows and columns a a 2 a 3 a n a 2 a 22 a 23 a 2n A = a m a mn We
More information1 GSW Sets of Systems
1 Often, we have to solve a whole series of sets of simultaneous equations of the form y Ax, all of which have the same matrix A, but each of which has a different known vector y, and a different unknown
More informationMAT 1302B Mathematical Methods II
MAT 1302B Mathematical Methods II Alistair Savage Mathematics and Statistics University of Ottawa Winter 2015 Lecture 19 Alistair Savage (uottawa) MAT 1302B Mathematical Methods II Winter 2015 Lecture
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 informationFinite-Horizon Statistics for Markov chains
Analyzing FSDT Markov chains Friday, September 30, 2011 2:03 PM Simulating FSDT Markov chains, as we have said is very straightforward, either by using probability transition matrix or stochastic update
More informationMath 166: Topics in Contemporary Mathematics II
Math 166: Topics in Contemporary Mathematics II Xin Ma Texas A&M University November 26, 2017 Xin Ma (TAMU) Math 166 November 26, 2017 1 / 14 A Review A Markov process is a finite sequence of experiments
More informationPart III: A Simplex pivot
MA 3280 Lecture 31 - More on The Simplex Method Friday, April 25, 2014. Objectives: Analyze Simplex examples. We were working on the Simplex tableau The matrix form of this system of equations is called
More informationAnd, even if it is square, we may not be able to use EROs to get to the identity matrix. Consider
.2. Echelon Form and Reduced Row Echelon Form In this section, we address what we are trying to achieve by doing EROs. We are trying to turn any linear system into a simpler one. But what does simpler
More informationMODEL ANSWERS TO THE THIRD HOMEWORK
MODEL ANSWERS TO THE THIRD HOMEWORK 1 (i) We apply Gaussian elimination to A First note that the second row is a multiple of the first row So we need to swap the second and third rows 1 3 2 1 2 6 5 7 3
More informationMATH 320, WEEK 11: Eigenvalues and Eigenvectors
MATH 30, WEEK : Eigenvalues and Eigenvectors Eigenvalues and Eigenvectors We have learned about several vector spaces which naturally arise from matrix operations In particular, we have learned about the
More informationHonors Advanced Mathematics Determinants page 1
Determinants page 1 Determinants For every square matrix A, there is a number called the determinant of the matrix, denoted as det(a) or A. Sometimes the bars are written just around the numbers of the
More informationPackage matrixnormal
Version 0.0.0 Date 2018-10-26 Type Package Title The Matrix Normal Distribution Depends R (>= 3.5.0) Imports mvtnorm (>= 1.0.8), utils (>= 3.5.0) Package matrixnormal October 30, 2018 Suggests datasets
More informationLinear Algebra Handout
Linear Algebra Handout References Some material and suggested problems are taken from Fundamentals of Matrix Algebra by Gregory Hartman, which can be found here: http://www.vmi.edu/content.aspx?id=779979.
More informationChapter 4. Solving Systems of Equations. Chapter 4
Solving Systems of Equations 3 Scenarios for Solutions There are three general situations we may find ourselves in when attempting to solve systems of equations: 1 The system could have one unique solution.
More informationMarkov Chains Handout for Stat 110
Markov Chains Handout for Stat 0 Prof. Joe Blitzstein (Harvard Statistics Department) Introduction Markov chains were first introduced in 906 by Andrey Markov, with the goal of showing that the Law of
More informationChapter 5. Matrices. 5.1 Basic matrix notation. This material is in Chapter 1 of Anton & Rorres.
Chapter 5. Matrices This material is in Chapter of Anton & Rorres. 5. Basic matrix notation We recall that a matrix is a rectangular array or table of numbers. We call the individual numbers entries of
More informationMATH 221: SOLUTIONS TO SELECTED HOMEWORK PROBLEMS
MATH 221: SOLUTIONS TO SELECTED HOMEWORK PROBLEMS 1. HW 1: Due September 4 1.1.21. Suppose v, w R n and c is a scalar. Prove that Span(v + cw, w) = Span(v, w). We must prove two things: that every element
More informationEigenvectors and Hermitian Operators
7 71 Eigenvalues and Eigenvectors Basic Definitions Let L be a linear operator on some given vector space V A scalar λ and a nonzero vector v are referred to, respectively, as an eigenvalue and corresponding
More informationDiagonalization. MATH 1502 Calculus II Notes. November 4, 2008
Diagonalization MATH 1502 Calculus II Notes November 4, 2008 We want to understand all linear transformations L : R n R m. The basic case is the one in which n = m. That is to say, the case in which the
More informationLecture 5: Latin Squares and Magic
Latin Squares Instructor: Padraic Bartlett Lecture 5: Latin Squares and Magic Week Mathcamp 0 Today s application is to magic! Not the friendship kind though ; instead we re going to talk about magic squares
More informationDot Products, Transposes, and Orthogonal Projections
Dot Products, Transposes, and Orthogonal Projections David Jekel November 13, 2015 Properties of Dot Products Recall that the dot product or standard inner product on R n is given by x y = x 1 y 1 + +
More informationVolume in n Dimensions
Volume in n Dimensions MA 305 Kurt Bryan Introduction You ve seen that if we have two vectors v and w in two dimensions then the area spanned by these vectors can be computed as v w = v 1 w 2 v 2 w 1 (where
More informationStochastic processes and Markov chains (part II)
Stochastic processes and Markov chains (part II) Wessel van Wieringen w.n.van.wieringen@vu.nl Department of Epidemiology and Biostatistics, VUmc & Department of Mathematics, VU University Amsterdam, The
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 information30.3. LU Decomposition. Introduction. Prerequisites. Learning Outcomes
LU Decomposition 30.3 Introduction In this Section we consider another direct method for obtaining the solution of systems of equations in the form AX B. Prerequisites Before starting this Section you
More informationSection 1.1 Algorithms. Key terms: Algorithm definition. Example from Trapezoidal Rule Outline of corresponding algorithm Absolute Error
Section 1.1 Algorithms Key terms: Algorithm definition Example from Trapezoidal Rule Outline of corresponding algorithm Absolute Error Approximating square roots Iterative method Diagram of a general iterative
More informationLecture 10: The Normal Distribution. So far all the random variables have been discrete.
Lecture 10: The Normal Distribution 1. Continuous Random Variables So far all the random variables have been discrete. We need a different type of model (called a probability density function) for continuous
More informationA Note on the Eigenvalues and Eigenvectors of Leslie matrices. Ralph Howard Department of Mathematics University of South Carolina
A Note on the Eigenvalues and Eigenvectors of Leslie matrices Ralph Howard Department of Mathematics University of South Carolina Vectors and Matrices A size n vector, v, is a list of n numbers put in
More informationWhen they compared their results, they had an interesting discussion:
27 2.5 Making My Point A Solidify Understanding Task Zac and Sione were working on predicting the number of quilt blocks in this pattern: CC BY Camille King https://flic.kr/p/hrfp When they compared their
More informationTrue or False? Lamarck s Theory of Evolution. Jean-Baptiste Pierre Antoine de Monet, Chevalier de Lamarck
True or False? We know what it is, we ve seen the evidence, but Aim: How does evolution happen? Charles Darwin was the 1 st scientist to offer an explanation for how Evolution happens. Jean-Baptiste Pierre
More informationMATH 320, WEEK 7: Matrices, Matrix Operations
MATH 320, WEEK 7: Matrices, Matrix Operations 1 Matrices We have introduced ourselves to the notion of the grid-like coefficient matrix as a short-hand coefficient place-keeper for performing Gaussian
More informationMath 123, Week 2: Matrix Operations, Inverses
Math 23, Week 2: Matrix Operations, Inverses Section : Matrices We have introduced ourselves to the grid-like coefficient matrix when performing Gaussian elimination We now formally define general matrices
More informationSENSITIVITY AND ELASTICITY ANALYSES
1 SENSITIVITY AND ELASTICITY ANALYSES Objectives Using the stage-based matrix model for a sea turtle population, conduct a sensitivity analysis of model parameters to determine the absolute contribution
More informationThe Derivative of a Function
The Derivative of a Function James K Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University March 1, 2017 Outline A Basic Evolutionary Model The Next Generation
More informationMAT1302F Mathematical Methods II Lecture 19
MAT302F Mathematical Methods II Lecture 9 Aaron Christie 2 April 205 Eigenvectors, Eigenvalues, and Diagonalization Now that the basic theory of eigenvalues and eigenvectors is in place most importantly
More informationHandout 1 EXAMPLES OF SOLVING SYSTEMS OF LINEAR EQUATIONS
22M:33 J. Simon page 1 of 7 22M:33 Summer 06 J. Simon Example 1. Handout 1 EXAMPLES OF SOLVING SYSTEMS OF LINEAR EQUATIONS 2x 1 + 3x 2 5x 3 = 10 3x 1 + 5x 2 + 6x 3 = 16 x 1 + 5x 2 x 3 = 10 Step 1. Write
More informationEigenspaces. (c) Find the algebraic multiplicity and the geometric multiplicity for the eigenvaules of A.
Eigenspaces 1. (a) Find all eigenvalues and eigenvectors of A = (b) Find the corresponding eigenspaces. [ ] 1 1 1 Definition. If A is an n n matrix and λ is a scalar, the λ-eigenspace of A (usually denoted
More informationMath Camp. Justin Grimmer. Associate Professor Department of Political Science Stanford University. September 9th, 2016
Math Camp Justin Grimmer Associate Professor Department of Political Science Stanford University September 9th, 2016 Justin Grimmer (Stanford University) Methodology I September 9th, 2016 1 / 61 Where
More informationOne-to-one functions and onto functions
MA 3362 Lecture 7 - One-to-one and Onto Wednesday, October 22, 2008. Objectives: Formalize definitions of one-to-one and onto One-to-one functions and onto functions At the level of set theory, there are
More informationIterated Functions. Tom Davis November 5, 2009
Iterated Functions Tom Davis tomrdavis@earthlink.net http://www.geometer.org/mathcircles November 5, 2009 Abstract In this article we will examine various properties of iterated functions. If f(x) is a
More informationDiscrete Mathematics and Probability Theory Summer 2014 James Cook Note 5
CS 70 Discrete Mathematics and Probability Theory Summer 2014 James Cook Note 5 Modular Arithmetic In several settings, such as error-correcting codes and cryptography, we sometimes wish to work over a
More informationThe minimal polynomial
The minimal polynomial Michael H Mertens October 22, 2015 Introduction In these short notes we explain some of the important features of the minimal polynomial of a square matrix A and recall some basic
More information2.5 Multivariate Curve Resolution (MCR)
2.5 Multivariate Curve Resolution (MCR) Lecturer: Dr. Lionel Blanchet The Multivariate Curve Resolution (MCR) methods are widely used in the analysis of mixtures in chemistry and biology. The main interest
More information19. TAYLOR SERIES AND TECHNIQUES
19. TAYLOR SERIES AND TECHNIQUES Taylor polynomials can be generated for a given function through a certain linear combination of its derivatives. The idea is that we can approximate a function by a polynomial,
More informationMATH10212 Linear Algebra B Homework Week 4
MATH22 Linear Algebra B Homework Week 4 Students are strongly advised to acquire a copy of the Textbook: D. C. Lay Linear Algebra and its Applications. Pearson, 26. ISBN -52-2873-4. Normally, homework
More informationExamples True or false: 3. Let A be a 3 3 matrix. Then there is a pattern in A with precisely 4 inversions.
The exam will cover Sections 6.-6.2 and 7.-7.4: True/False 30% Definitions 0% Computational 60% Skip Minors and Laplace Expansion in Section 6.2 and p. 304 (trajectories and phase portraits) in Section
More informationDeterminants and Scalar Multiplication
Properties of Determinants In the last section, we saw how determinants interact with the elementary row operations. There are other operations on matrices, though, such as scalar multiplication, matrix
More informationEigenvalues and Eigenvectors
Eigenvalues and Eigenvectors MAT 67L, Laboratory III Contents Instructions (1) Read this document. (2) The questions labeled Experiments are not graded, and should not be turned in. They are designed for
More informationPerform the same three operations as above on the values in the matrix, where some notation is given as a shorthand way to describe each operation:
SECTION 2.1: SOLVING SYSTEMS OF EQUATIONS WITH A UNIQUE SOLUTION In Chapter 1 we took a look at finding the intersection point of two lines on a graph. Chapter 2 begins with a look at a more formal approach
More informationNumerical Methods Lecture 2 Simultaneous Equations
Numerical Methods Lecture 2 Simultaneous Equations Topics: matrix operations solving systems of equations pages 58-62 are a repeat of matrix notes. New material begins on page 63. Matrix operations: Mathcad
More informationLecture 3: Latin Squares and Groups
Latin Squares Instructor: Padraic Bartlett Lecture 3: Latin Squares and Groups Week 2 Mathcamp 2012 In our last lecture, we came up with some fairly surprising connections between finite fields and Latin
More informationChapter 29 out of 37 from Discrete Mathematics for Neophytes: Number Theory, Probability, Algorithms, and Other Stuff by J. M.
29 Markov Chains Definition of a Markov Chain Markov chains are one of the most fun tools of probability; they give a lot of power for very little effort. We will restrict ourselves to finite Markov chains.
More informationVector Spaces. 9.1 Opening Remarks. Week Solvable or not solvable, that s the question. View at edx. Consider the picture
Week9 Vector Spaces 9. Opening Remarks 9.. Solvable or not solvable, that s the question Consider the picture (,) (,) p(χ) = γ + γ χ + γ χ (, ) depicting three points in R and a quadratic polynomial (polynomial
More informationRatios, Proportions, Unit Conversions, and the Factor-Label Method
Ratios, Proportions, Unit Conversions, and the Factor-Label Method Math 0, Littlefield I don t know why, but presentations about ratios and proportions are often confused and fragmented. The one in your
More informationEigenvalue and Eigenvector Homework
Eigenvalue and Eigenvector Homework Olena Bormashenko November 4, 2 For each of the matrices A below, do the following:. Find the characteristic polynomial of A, and use it to find all the eigenvalues
More informationChapter 1. Vectors, Matrices, and Linear Spaces
1.7 Applications to Population Distributions 1 Chapter 1. Vectors, Matrices, and Linear Spaces 1.7. Applications to Population Distributions Note. In this section we break a population into states and
More informationMatrix Algebra: Definitions and Basic Operations
Section 4 Matrix Algebra: Definitions and Basic Operations Definitions Analyzing economic models often involve working with large sets of linear equations. Matrix algebra provides a set of tools for dealing
More information18.06 Problem Set 8 - Solutions Due Wednesday, 14 November 2007 at 4 pm in
806 Problem Set 8 - Solutions Due Wednesday, 4 November 2007 at 4 pm in 2-06 08 03 Problem : 205+5+5+5 Consider the matrix A 02 07 a Check that A is a positive Markov matrix, and find its steady state
More information22m:033 Notes: 3.1 Introduction to Determinants
22m:033 Notes: 3. Introduction to Determinants Dennis Roseman University of Iowa Iowa City, IA http://www.math.uiowa.edu/ roseman October 27, 2009 When does a 2 2 matrix have an inverse? ( ) a a If A =
More information5.3 Other Algebraic Functions
5.3 Other Algebraic Functions 397 5.3 Other Algebraic Functions This section serves as a watershed for functions which are combinations of polynomial, and more generally, rational functions, with the operations
More informationx(t+ δt) - x(t) = slope δt t+δt
Techniques of Physics Worksheet 2 Classical Vibrations and Waves Introduction You will have encountered many different examples of wave phenomena in your courses and should be familiar with most of the
More informationNumerical Analysis Lecture Notes
Numerical Analysis Lecture Notes Peter J Olver 8 Numerical Computation of Eigenvalues In this part, we discuss some practical methods for computing eigenvalues and eigenvectors of matrices Needless to
More informationLecture 5: Special Functions and Operations
Lecture 5: Special Functions and Operations Feedback of Assignment2 Rotation Transformation To rotate by angle θ counterclockwise, set your transformation matrix A as [ ] cos θ sin θ A =. sin θ cos θ We
More informationWarm-up Simple methods Linear recurrences. Solving recurrences. Misha Lavrov. ARML Practice 2/2/2014
Solving recurrences Misha Lavrov ARML Practice 2/2/2014 Warm-up / Review 1 Compute 100 k=2 ( 1 1 ) ( = 1 1 ) ( 1 1 ) ( 1 1 ). k 2 3 100 2 Compute 100 k=2 ( 1 1 ) k 2. Homework: find and solve problem Algebra
More informationMATH 521, WEEK 2: Rational and Real Numbers, Ordered Sets, Countable Sets
MATH 521, WEEK 2: Rational and Real Numbers, Ordered Sets, Countable Sets 1 Rational and Real Numbers Recall that a number is rational if it can be written in the form a/b where a, b Z and b 0, and a number
More informationMA 1125 Lecture 15 - The Standard Normal Distribution. Friday, October 6, Objectives: Introduce the standard normal distribution and table.
MA 1125 Lecture 15 - The Standard Normal Distribution Friday, October 6, 2017. Objectives: Introduce the standard normal distribution and table. 1. The Standard Normal Distribution We ve been looking at
More informationBasic Definitions: Indexed Collections and Random Functions
Chapter 1 Basic Definitions: Indexed Collections and Random Functions Section 1.1 introduces stochastic processes as indexed collections of random variables. Section 1.2 builds the necessary machinery
More informationIntroduction to R, Part I
Introduction to R, Part I Basic math, variables, and variable types Tim Frasier Copyright Tim Frasier This work is licensed under the Creative Commons Attribution 4.0 International license. Click here
More informationMITOCW ocw f99-lec01_300k
MITOCW ocw-18.06-f99-lec01_300k Hi. This is the first lecture in MIT's course 18.06, linear algebra, and I'm Gilbert Strang. The text for the course is this book, Introduction to Linear Algebra. And the
More informationMatrices MA1S1. Tristan McLoughlin. November 9, Anton & Rorres: Ch
Matrices MA1S1 Tristan McLoughlin November 9, 2014 Anton & Rorres: Ch 1.3-1.8 Basic matrix notation We have studied matrices as a tool for solving systems of linear equations but now we want to study them
More informationMath 308 Midterm Answers and Comments July 18, Part A. Short answer questions
Math 308 Midterm Answers and Comments July 18, 2011 Part A. Short answer questions (1) Compute the determinant of the matrix a 3 3 1 1 2. 1 a 3 The determinant is 2a 2 12. Comments: Everyone seemed to
More informationMatrices related to linear transformations
Math 4326 Fall 207 Matrices related to linear transformations We have encountered several ways in which matrices relate to linear transformations. In this note, I summarize the important facts and formulas
More informationQuantitative Understanding in Biology Module III: Linear Difference Equations Lecture I: Introduction to Linear Difference Equations
Quantitative Understanding in Biology Module III: Linear Difference Equations Lecture I: Introduction to Linear Difference Equations Introductory Remarks This section of the course introduces dynamic systems;
More informationMetric Predicted Variable With One Nominal Predictor Variable
Metric Predicted Variable With One Nominal Predictor Variable Tim Frasier Copyright Tim Frasier This work is licensed under the Creative Commons Attribution 4.0 International license. Click here for more
More informationIntroduction to Digital Evolution Handout Answers
Introduction to Digital Evolution Handout Answers Note to teacher: The questions in this handout and the suggested answers (in red, below) are meant to guide discussion, not be an assessment. It is recommended
More informationof 8 28/11/ :25
Paul's Online Math Notes Home Content Chapter/Section Downloads Misc Links Site Help Contact Me Differential Equations (Notes) / First Order DE`s / Modeling with First Order DE's [Notes] Differential Equations
More informationLinear Algebra Basics
Linear Algebra Basics For the next chapter, understanding matrices and how to do computations with them will be crucial. So, a good first place to start is perhaps What is a matrix? A matrix A is an array
More informationThe Exciting Guide To Probability Distributions Part 2. Jamie Frost v1.1
The Exciting Guide To Probability Distributions Part 2 Jamie Frost v. Contents Part 2 A revisit of the multinomial distribution The Dirichlet Distribution The Beta Distribution Conjugate Priors The Gamma
More informationLTCC. Exercises solutions
1. Markov chain LTCC. Exercises solutions (a) Draw a state space diagram with the loops for the possible steps. If the chain starts in state 4, it must stay there. If the chain starts in state 1, it will
More informationLecture 14: The Rank-Nullity Theorem
Math 108a Professor: Padraic Bartlett Lecture 14: The Rank-Nullity Theorem Week 6 UCSB 2013 In today s talk, the last before we introduce the concept of matrices, we prove what is arguably the strongest
More information