Math 4 - Summary Notes

Similar documents
Linear Algebra: Lecture Notes. Dr Rachel Quinlan School of Mathematics, Statistics and Applied Mathematics NUI Galway

Linear Algebra. Matrices Operations. Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0.

[Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty.]

Linear Algebra Highlights

Math Linear Algebra Final Exam Review Sheet

Chapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form

Math 520 Exam 2 Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

MATH 240 Spring, Chapter 1: Linear Equations and Matrices

Math 4A Notes. Written by Victoria Kala Last updated June 11, 2017

Linear Algebra March 16, 2019

All of my class notes can be found at

NOTES FOR LINEAR ALGEBRA 133

Math Camp Lecture 4: Linear Algebra. Xiao Yu Wang. Aug 2010 MIT. Xiao Yu Wang (MIT) Math Camp /10 1 / 88

2. Linear algebra. matrices and vectors. linear equations. range and nullspace of matrices. function of vectors, gradient and Hessian

c c c c c c c c c c a 3x3 matrix C= has a determinant determined by

ECON2285: Mathematical Economics

Final Review Sheet. B = (1, 1 + 3x, 1 + x 2 ) then 2 + 3x + 6x 2

Linear Algebra. The analysis of many models in the social sciences reduces to the study of systems of equations.

SAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear Algebra

M. Matrices and Linear Algebra

MTH Linear Algebra. Study Guide. Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education

Digital Workbook for GRA 6035 Mathematics

3 Matrix Algebra. 3.1 Operations on matrices

MATH 310, REVIEW SHEET 2

LS.1 Review of Linear Algebra

Equality: Two matrices A and B are equal, i.e., A = B if A and B have the same order and the entries of A and B are the same.

Knowledge Discovery and Data Mining 1 (VO) ( )

STEP Support Programme. STEP 2 Matrices Topic Notes

A matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and

SAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear Algebra

MAT Linear Algebra Collection of sample exams

A matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and

Chapter 9: Systems of Equations and Inequalities

Chapter 3. Determinants and Eigenvalues

1 Matrices and Systems of Linear Equations. a 1n a 2n

Introduction to Quantitative Techniques for MSc Programmes SCHOOL OF ECONOMICS, MATHEMATICS AND STATISTICS MALET STREET LONDON WC1E 7HX

Elementary Linear Algebra

Math Camp II. Basic Linear Algebra. Yiqing Xu. Aug 26, 2014 MIT

4. Determinants.

MTH 464: Computational Linear Algebra

1. General Vector Spaces

Linear Algebra Primer

A Brief Outline of Math 355

Topic 1: Matrix diagonalization

IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET

Introduction to Matrices

Massachusetts Institute of Technology Department of Economics Statistics. Lecture Notes on Matrix Algebra

Math Camp Notes: Linear Algebra I

TOPIC III LINEAR ALGEBRA

CS 246 Review of Linear Algebra 01/17/19

1 Last time: determinants

Chapter 1: Systems of linear equations and matrices. Section 1.1: Introduction to systems of linear equations

Some Notes on Linear Algebra

MATRICES. knowledge on matrices Knowledge on matrix operations. Matrix as a tool of solving linear equations with two or three unknowns.

Lecture Notes in Linear Algebra

Chapter 2: Matrix Algebra

Linear Algebra: Matrix Eigenvalue Problems

Chapter 5 Eigenvalues and Eigenvectors

Solving a system by back-substitution, checking consistency of a system (no rows of the form

Chapter 2: Matrices and Linear Systems

Fundamentals of Engineering Analysis (650163)

IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET

Systems of Linear Equations and Matrices

Introduction to Matrix Algebra

Glossary of Linear Algebra Terms. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Math Bootcamp An p-dimensional vector is p numbers put together. Written as. x 1 x =. x p

Linear Algebra Primer

Linear Algebra: Lecture notes from Kolman and Hill 9th edition.

ENGI 9420 Lecture Notes 2 - Matrix Algebra Page Matrix operations can render the solution of a linear system much more efficient.

1 9/5 Matrices, vectors, and their applications

Systems of Linear Equations and Matrices

Sometimes the domains X and Z will be the same, so this might be written:

Finite Mathematics Chapter 2. where a, b, c, d, h, and k are real numbers and neither a and b nor c and d are both zero.

ECON0702: Mathematical Methods in Economics

Review of Basic Concepts in Linear Algebra

Math 302 Outcome Statements Winter 2013

NOTES on LINEAR ALGEBRA 1

Review of Linear Algebra

Elementary maths for GMT

1 Matrices and Systems of Linear Equations

Getting Started with Communications Engineering. Rows first, columns second. Remember that. R then C. 1

3 (Maths) Linear Algebra

Matrices and Linear Algebra

Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition

Linear Algebra: Linear Systems and Matrices - Quadratic Forms and Deniteness - Eigenvalues and Markov Chains

OR MSc Maths Revision Course

Topic 2 Quiz 2. choice C implies B and B implies C. correct-choice C implies B, but B does not imply C

1. In this problem, if the statement is always true, circle T; otherwise, circle F.

MATH 310, REVIEW SHEET

Graduate Mathematical Economics Lecture 1

Econ Slides from Lecture 7

Topic 15 Notes Jeremy Orloff

Math 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination

MATRIX ALGEBRA. or x = (x 1,..., x n ) R n. y 1 y 2. x 2. x m. y m. y = cos θ 1 = x 1 L x. sin θ 1 = x 2. cos θ 2 = y 1 L y.

1. Linear systems of equations. Chapters 7-8: Linear Algebra. Solution(s) of a linear system of equations (continued)

Basic Concepts in Linear Algebra

Things we can already do with matrices. Unit II - Matrix arithmetic. Defining the matrix product. Things that fail in matrix arithmetic

MATH 20F: LINEAR ALGEBRA LECTURE B00 (T. KEMP)

REVIEW FOR EXAM II. The exam covers sections , the part of 3.7 on Markov chains, and

ELEMENTARY LINEAR ALGEBRA

Transcription:

Math 4 - Summary Notes July, 0 Chat: what is class about? Mathematics of Economics, NOT Economics! Lecturer NOT Economist - have some training, but reluctant to discuss strict economic models because much of my knowledge comes from non-standard places (actuarial exams, etc) In particular, the IS-LM model in st section of is NOT examinable and will only be covered lightly here If you know about these things from other places, do the questions as they ll help! The first part of the class is Linear Algebra It gets everywhere and is probably the most applicable area of undergrad math Depending on yr background it can be applied to traffic networks, cell division models, physical motion, engineering problems, etc, etc Why so useful? Primarily cos so easy (relatively) Linear problems can mostly be solved exactly using algorithms Non-linear problems often cannot Sacrifice correctness in the model for easy of solution Best fitting straight line, etc Easier then complex curve Suppose have supply-demand model with two curves (quantity =f(price)) q s = f (p), q d = g(p) (f monotone down, g monotone up) Equilibrium price/quantity given where curves meet Hard to find exactly, but if the curves were straight lines NOT a course of economic concepts! Will use some to illustrate the math, but the math is independent Economics is ALL models - choosing right from start how to balance accuracy and easy of solution/analysis Math different: models approximate exactness Draw pic! Math fits between Reality and Explanation/Prediction of Reality Modeling is the first link, solution is the second Pure math is mostly concerned with the second link - economics equally with both Graph of model for economic growth versus interest rates Different models - which is best? 7 Systems of Linear Equations 7 Solving Systems of Linear Equations Definition 7 A linear equation in unknowns is an equation of the form αx + βy = γ, where α, β, γ are given constants and x, y are variables A system of two equations in two unknowns is a pair, α x + β y = γ, α x + β y = γ

Notice: 6 constants and two variables A solution to a system is any pair of numbers (x, y) which satisfies both equations in the system The solution set is the complete collection of solutions A system is consistent if it has at least one solution and inconsistent if it has no solutions Graph linear equations as straight lines: example x + y =, x y = 4, x + y = Can solve all three possible systems of two equations in two unknowns Plot and solve ((, ), (0, ), (, )) by inspection Pictures: what can happen? Two lines either intersect at a point (most likely), are parallel and distinct (next most likely) or are the same line (least likely) Pics and examples x + 4x = x + 4x = x + 4x =,, - geometry = lines x + x = 6x + 8x = 6x + 8x = 0 We have proved the following: Theorem 7 The solution set of a system of two equations in two unknowns has either none, one or infinite solutions Suppose add a third equation system of three equations in two unknowns Most likely no solutions Picture Substitution and elimination Generally difficult to solve systems graphically if you want numerical solutions rather than just knowing they exist Substitution = solve one equation for x (in terms of y), or vice versa, then substitute in Eg with above Elimination involves adding/subtracting multiples of the equations to eliminate a variable Eg More examples: Market equilibrium Examples Suppose that the supplied and demanded quantities q s, q d of a good are related to the price p by the following equations: q s = 5 + p, q d = 40 p Find the equilibrium price and quantity of goods Answer: at equilibrium, supply = demand, hence 5 + p = 40 p 5p = 45 p = 9 q s = q d = Two related goods have the following supply and demand equations: q s = 6 + 4p, q d = 50 p + p q s = 4 + p, q s = 0 p + p Find the equilibrium prices and quantities of the two goods Answer: supply = demand qi s = qi d results in two equations: 6 + 4p = 50 p + p 6p p = 56 = 4 + p = 0 p + p p + p = 4 = 9p = 90,

so that p = 0, p = 4, q = 4, q = 4 Let c be a constant and consider the system of equations x y =, cx + y = For each value of c find how many solutions the system has and calculate them (in terms of c) Ans: (x, y) = +c (7, 9 c) if c = / If so, no solutions 4 As the previous question except now x + y =, cx + cy = Ans: (x, y) = c (c +, ) if c = 0, No soln if c = 0 c = gives infinite solns ((x, y) = ( y, y)) 7 Linear Systems in n-variables Definition 7 Linear Equation in n unknowns: a x + + a n x n = b Here a i, b are constants, x i variables (unknowns) Linear system of m equations in n unknowns (m n-system) a x + + a n x n = b a x + + a n x n = b a m x + + a mn x n = b m A solution to a system is any n-tuple (x,, x n ) which satisfies every equation in the system The solution set is the complete collection of solutions Common to use x, y for m and x, y, z for m systems Definition 74 Two linear systems are equivalent if they have the same solution set (systems don t have to have same number of equations) x + 4x = x = equivalent to equivalent to x + x = 0 x = x + x = x + 5x =, etc, etc

Row Operations Three ways of changing a system to get an equivalent one: Multiply both sides of an equation by a non-zero constant Add a multiple of one equation to another Change the order of two equations Idea - use row operations to put a system in a nice form so we can read off solutions Examples Get (,, ) Get ( + z,, z) infinity of solutions Some facts: x + 4y z = 9 x + y z = 9 x y + z = x + y z = 0 x + y z = 7 x y 4z = All systems of linear equations have either no, one, or infinite numbers of solutions A system is consistent if it has at least one solution and inconsistent if it has no solutions A system with more equations than unknowns (m > n) is overdetermined If m < n the system is underdetermined A consistent underdetermined system has infinite solutions Table of under/over/square versus expected numbers of solns Matrix Arrays Definition 75 The augmented matrix of an m n system is the matrix a a n b a n a nn b n Don t need to put the line in, but it s often helpful Use row operations on these matrices and systems to solve them side by side See that using augmented matrix beams that can avoid having o write down the variables x i (saves time and makes mistakes less likely) Know how to go between a linear system and its augmented matrix Example: x + x + x = 4 x + x = x x = 0 Solution: (/,/,/) - use bottom row to eliminate x terms in others first 4

Row Echelon Form Want to do row operations until matrix is a simple as possible (roughly speaking as many zeros as possible with s going diagonally (Reduced Row Echelon Form) Every matrix has one and it doesn t matter what row ops you do, will always get to the same RREF Once got there can either read off the solns or easily tell there are none 0 Example: 0 reduce to row of zeros has no solution: rd row is 0= contradiction 0 Left part of augmented matrix isn t triangular, but is in a useful and easy to read form Definition 76 A is in Reduced Row Echelon Form if Any entirely zero rows are at the bottom First = 0 entry in each non-zero row is a If a row isn t all zero, then the next row has more leading zeros (not necessarily one place to the left - book mistake!) 4 Each column with a leading has zeros elsewhere 4 0 0 0 Examples 0 0 0 0 0 0 4 0 5 0 Homogeneous Systems 0 0 0 0 0 0 0 which yields x =, x =, x = which yields (x, x, x, x 4 ) = ( x, x, x, ) A system is homogeneous if all the equations are of the form a i x + a i x + + a in x n = 0 Clearly any homogeneous system is consistent, since (x,, x n ) = (0,, 0) is a (trivial) solution If the system is underdetermined then it automatically has infinitely many solutions Example Let c be a constant Find all the solutions to the homogeneous system x + y + z = 0 6x y + 4z = 0 cx + cy z = 0 Only trivial unless c = in which case get line (λ, 0, λ) (Ignore linear independence discussion - it s the wrong use of the phrase Ignore free-variables too) 5

8 Matrices 8 General Notation Matrices are used to summarize information and to facilitate many calculations at once Will see the latter in action later Definition 8 A matrix is a rectangular array of objects (usually numbers) in parentheses (square or round) A matrix with m rows and n columns is described as an m n matrix A (column) vector is an m matrix Write A, B, C for matrices (capitals) Entries are a ij for the ith row, jth column of A Sometimes write A = (a ij ) Write column vectors boldface (underlined) x Scalars are usually written with lower case letters Sometimes write a matrix in terms of its column vectors: A = a a m a n a nm = (a,, a n ), a i = Definition 8 Two matrices are equal if they have the same dimensions and their corresponding entries are equal A matrix is square if m = n (same number of rows and columns) A square matrix is diagonal if a ij = 0 for i = j The identity matrix I n is the diagonal n n matrix such that a ii = A square matrix all of whose entries are zero is called the zero matrix or null matrix 8 Basic Matrix Operations Scalar Multiplication (example) αa Matrix Addition (same size example) A + B Multiplication matrix by vector Ax: need, A m n and x R n Linear systems are matrix mult Examples 4 Multiplication matrix by matrix AB AB = (Ab,, Ab n ) Need A m r, B r n then AB is m n Note: AB might be defined even when BA is not Even when they are, the two are generally not equal ( ) ( ) p q Examples Let A = and anything but B = for AB = BA 4 7 q p + q Cost and profits Suppose a grocer buys and sells apples, bananas and oranges at the following cents per fruit: a i a ni 6

Fruit Cost to buy Cost to sell Apple 40 55 Banana 0 60 Orange 50 80 Suppose the grocer buys b = (b a, b b, b o ) T and sells s = (s a, s b, s o ) T of each fruit Then the profit of the grocer is 55 40 (s a, s b, s o ) 60 (b a, b b, b o ) 0 = 55s a 40b a + 60s b 0b b + 80s o 50b o 80 50 Supply and demand in many markets Suppose have a market of supplies of goods q s,, qs n and demanded quantities q d,, qd n arranged in vectors qs, q d If p is the vector of the prices of these goods, then a linear model relating quantities and prices would be q s = a + Ap, q d = b + Bp, where a, b are constant vectors of length n and A, B are square n n matrices The equilibrium pricesd would then be a solution p to the equation (A B)p = b a We shall see how to solve this for p and thus q in chapter 9 8 Matrix Transposition Definition 8 The transpose of an m n A = (a ij ) is the n m A T = (b ij ) where b ij = a ji Examples Rules - behave well with scalar mult and addition and (AB) T = B T A T (and with C) - examples Note (A T ) T = A Definition 84 A is symmetric if A = A T (note that we need A to be square) 84 Some Special Matrices Some types of matrices have special properties that are useful in studying equations Definition 85 A square matrix is idempotent if A = A Hence A n = A for all n =,,, 4, 5, Example A = (4, 6//, ) is idempotent (Ignore partitioned matrices) Definition 86 The trace of a square matrix A is the sum a + a + + a nn Theorem 87 tr(ab) = tr(ba) (regardless of whether AB = BA or not - proof in book) Example Check for two matrices (,//,-) and (,-//,) Observe tr(ab) = tr(a) tr(b) 7

9 Determinants and the Inverse Matrix 9 Defining the Inverse Matrix algebra works almost like normal algebra with the exception of commutativity of multiplication and that we have no division Inverses: =, so what is comparable idea for matrices? Only works for most but not all square matrices Definition 9 The inverse A of an n n matrix A is the n n matrix which satisfies AA = A A = I n Not every square matrix has an inverse Those that do not are known as singular Only one (at most) inverse can exist for any given matrix Why do we care? Want to solve matrix equations such as Ax = b, where A, b are given and x is unknown (of same length as b) Solution is x = A b DON T WRITE AS DIVISION!!! Generally A B = BA Examples - formula for Alternative to book: if B,C both inverse to A then B = BI = BAC = IC = C so only one possibility However ( ) a b thus = c d ad bc ( ) ( ) a b d b = c d ad bc c a ( ) d b, provided ad bc = 0 c a ( ) 0, 0 If ad bc = 0 then can t divide and so don t get an inverse (not proper argument) Examples Definition 9 ad bc is the determinant of the matrix A and is denoted A Observe that A is invertible (A exists) iff A = 0 Can write the inverse of a matrix as A = ( ) d b, if A = 0 A c a Properties of the Determinant Do example of each A T = A Swap two rows (or columns) of A changes sign of det Two identical rows (or columns) means A = 0 One row (column) a multiple of another then A = 0 Add a multiple of one row (column) to another leaves det unchanged 8

Det of triangular matrix is product of diagonal elements Multiply row (column) of A by λ multiplies det by λ Multiply every element of A by λ multiplies det by λ n (here n = ) AB = A B A + B = A + B (in general) Example Recall supply/demand equilibrium from eariler q s = a + Ap, q d = b + Bp, where a, b are constant vectors of length n and A, B are square n n matrices The equilibrium prices would then be a solution p to the equation (A B)p = b a, which yields p = (A B) (b a), provided A B is invertible Then Put in some numbers: suppose (as in earlier example) q s = 6 + 4p, q d = 50 p + p q s = 4 + p, q s = 0 + p p q s = ( ) 6 + 4 ( ) 4 0 p, q d = 0 Little of the geometry of by changing triangle area 9 Determinants and Inverses of Matrices ( ) ( ) 50 + p 0 Definition 9 The ijth minor M ij of a matrix A is the determinant of the matrix obtained by deleting the ith row and jth column of A The ijth cofactor of A is the C ij = ( ) i+j M ij (do cofactor pattern) Examples Definition 94 The determinant of A is any of the 6 sums: A = The other possible sums are always zero Examples + quick chat about volumes To find the inverse of a matrix with det = 0, Write down matrix of cofactors a ij C ij = a ij C ij i= j= Define adjoint adj(a) as transpose of matrix of cofactors A is divide by determinant Example 9

9 The Inverse of an n n matrix and its properties Can use cofactors to find the det and adjoint of higher order matrices just as before The cofactors are the determinants of the minors (up to sign), etc Becomes computationally time consuming however Easier (sometimes) method is to use two tricks: Diagonal method for determinant (A I) method for inverse Example A = 4 0 = 6 7 8 do both ways 6 9 Theorem 95 A, B non-singular AB non-singular and (AB) = B A (A ) = A (A ) T = (A T ) det(a ) = (det A) Proof Calculate B A (AB) and (AB)B A A A = I, etc Example Firm produces three outputs y, y, y with three inputs z, z, z The input-requirements matrix A = 0 4 relates the outputs to inputs according to z = Ay If firm wants to produce y =, y =, y = it needs z = (9, 7, 6) Suppose instead that firm has resources z = (7, 6, 9), then can produce y = A z = (5,, ) Here A = 6 0 4 Talk about complexity of say an administrator shutting down a company and trying to drain stock to zero while maximising profit 94 Cramer s rule Inverse matrices are useful for solving the equation Ax = b where A is square and invertible Cramer s rule is an alternative method that allows you to obtain as many or as few entries of x without having to compute all the others Simply use the formula: C C n x = A b A C n C nn b b n = b C + b C + + b n C n b C n + b C n + + b n C nn The upshot is that x i = A b C i + b C i + + b n C ni 0

for each i =,, n Otherwise said, x i is the determinant of the matrix A with its ith column replaced with b Example Calculate x if Ax = b where A = 0 and b = (,, 5) T Ans 6 Could have 4 obtained from A earlier Sim x = 9, x = 5 Original method is usually best when need or more values of x i Also has advantage that once you have A can change to any b and still solve Finding all solutions with Cramer is same number of calculations as by finding inverse Example Open Leontief Input-Output model Model economy as a basket of related goods, the outputs of each being required for the others For example suppose there are three sectors to consider; electronics, mining, automotive If these sectors are related by the following: Each $ of production in each of the industries in the left column requires the corresponding level of production in Electronics Mining/Drilling Automotive Electronics 04 0 0 the top row Let matrix of above Mining/Drilling 0 0 04 Automotive 0 0 0 coefficients be A: a ij denotes the monetary input required into industry j in order to produce one unit of industry i s output Let x be the vector of money outputs of each industry, then Ax is the vector of total production demand Let d be the demand of consumers from the industries Thus total demand on the industries is Ax + d In a closed economy this has to equal the supply demand, hence x = Ax + d The result is x = (I A) d For our example 06 0 0 (I A) = 0 08 04 0 0 09 954 057 00 = 086 55 0690 0747 045 64 If the public demands d = (5000, 000, 0000) then the total output of the industries must be x = (586, 758, 674) Of course could have used Cramer s rule if just one variable was required 0 Advanced topics in Linear Algebra 0 Vector Spaces Recall: a vector is a matrix with n rows and column Written v v T is a row vector Think of vectors geomterically as arrows to points Thus v = ( ) is the vector pointing from the origin in R to the point with co-ordinates (, ) Example of adding/subtracting vectors graphically + scaler multiplication Definition 0 The inner product of two vectors v, w is the product v w = v T w = v i w i The length of a vector is v = v + + v n = v v Ex of length calc

Definition 0 A list of n-dimensional vectors v,, v n are linearly dependent iff det(v,, v n ) = 0 Equivalently we have theat there exist scalers λ,, λ n not all zero such that λ v + + λ n v n = 0 Vectors are linearly independent otherwise Ie the only scalar solns to above eqn are all zero In -dim a pair of vectors is linearly independent if they point in separate directions (span parallelogram with non-zero area) In -dim need to span parallelepiped with non-zero volume ( vectors don t lie in a common plane) Examples (, ) T and (, ) T are lin ind (,, 4) T, (,, 0) T and (0,, ) T are lin ind (calc det), but replace v with (0, 0, ) T to get trio that is not lin ind Definition 0 A basis for a vector space V is a list of linearly independent vectors so that any vector in V can be generated as a linear combination of vectors in the basis Example See above Observe that in R n a basis contains exactly n linearly indep vectors How to solve? Suppose have a basis v,, v n of R n Given a vector u, how to we find the coefficients? Take dot products of u with each basis vector to get n eqns in n variables: Point is that it can be done! v v v v n v n v v n v n λ λ n = u v u v n Example Let u = (4, ), v = (, ) T and v = (, ) T We want to find λ, λ such that u = λ v + λ v Get, Easiest when a basis is orthonormal: ie v i v j = δ ij In such cases, u = λ v + + λ n v n = λ i = u v i Angles between vectors: The angle between two vectors ( π/ < θ π/) is related to the vectors by v w = v w cos θ Clear that vectors are orthogonal iff they have zero inner product Standard basis of R n is orthogonal - ex There are lots of vector spaces that are not R n Eg sets of functions, etc Bases/dimension defined similarly: dimension is simply the maximum number of linear indep vectors in the space Definition 04 The rank of a matrix is the maximum number of linearly independent columns It always equals the maximum number of linearly independent rows Example of a 4by5 of rank Theorem 05 An n n matrix is non-singular iff det A = 0 A exists rank(a) = n

0 The Eigenvalue Problem Eigenvalues and eigenvectors are convenient concepts to describe what a matrix does when it multiplies vectors You can say that a matrix stretches by a factor of in one direction and in another - this is exactly the eigendescription of a matrix Definition 06 Let A be n n A non-zero vector x R n is an eigenvector of A with eigenvalue λ if Ax = λx A stretches x by a ( factor ) of λ but ( doesn t change ( ) its direction Example: If A =, x 4 =, x ) =, then Ax = x Ax = 5x Characteristic equation: Write Ax = λx as (A λi)x = 0 We re looking for non-zero x solving this But such an x exists iff A λi is singular: iff det(a λi) = 0 p(λ) = det(a λi) = 0 is degree n polynomial If λ is a solution then A λ I is singular and so there exists a non-zero vector v with (A λ I)v = 0 Thus all solutions to p(λ) = 0 are eigenvalues and to each there corresponds at least one eigenvector Since p(λ) degree n, there are at most n eigenvalues - repeated roots mean that many matrices do not have this many Eigenvalues may also be complex Note that evectors are defined only up to non-zero scale Theorem 07 λ is an eigenvalue of A iff (A λi)x = 0 has anon-zero solution x A λi is singular det(a λi) = 0 Examples A = 7Id has λ = 7 v = anything ( ) ( ) ( ) A = has λ = 4,, v =, v = ( ) ( ) A = has λ =, v = Only one eigenvalue and one eigenvector 4 4 A = 0 has λ =,,, v = 0, 0 0 0, 0 5 A = 0 has λ =,, v =, 0 6 A = 0 has λ =, v = 0 0 0 0

Diagonalization Diagonalizing a square matrix is about finding a similar diagonal matrix Ie given A we want to find an invertible M and diagonal D such that A = MDM Theorem 08 Suppose that an n n matrix A has n independent eigenvectors v,, v n with eigenvalues λ,, λ n Then A = MDM where ( ) A = is diagonalizable ( ) 0 A = is diagonalizable A = 0 is diagonalizable 0 0 0 4 A = 0 is diagonalizable 0 M = (v,, v n ), D = diag(λ,, λ n ) ( ) is not diagonalizable, neither is 0 4 0 0 Definition 09 A square matrix is orthogonal if it satisfies Q T Q = I Theorem 00 If A is a symmetric matrix then it has n independent orthonormal eigenvectors Putting these as the columns of a matrix Q we have Q being orthogonal Moroever Q AQ is diagonal with e values down diagonal 0 Quadratic Forms Useful in calculus especially - will use later All need to know is how to make them out of matrices Definition 0 A quadratic form is a scalar expresiion q(x) = x T Ax where A is a square matrix Don t need A to be symmetric but can always choose it to be so Examples: Definition 0 A quadratic form q(x) = x T Ax is positive definite > 0 positive semi-definite 0 if q(x) negative definite < 0 negative semi-definite 0 for all x = 0 4

In the case that A is symmetric the above can be translated into statements about eigenvalues: A quadratic form q(x) = x T Ax is positive definite > 0 positive semi-definite 0 if the eigenvalues of A are all negative definite < 0 negative semi-definite 0 Theorem 0 Can write all quadratic forms using symmetric matrices A +/ (semi-) definiteness relates to signs of eigenvalues Definition 04 A principal minor of a square matrix A is the determinant of a matrix obtained by deleting some rows and columns from A - the rows and columns must be identical The leading principal minors of a square matrix A are the determinants of the matrices A = a, A = ( a a a a ), A =, A n = A - ie delete the last n k rows and columns Theorem 05 A symmetric matrix A is positive definite iff all leading principal minors are positive A symmetric matrix A is negative definite iff the leading principal minors alternate -/+/-/+ etc A symmetric matrix A is positive semi-definite iff some principal minors are zero and the rest are positive A symmetric matrix A is negative semi-definite iff A is positive semi-definite Examples q(x) = 7x xy + y is positive definite: leading principal minors 7, 6 q(x) = y + xy is positive semit-definite: principal minors 0,, 0 q(x) = x + y + z + 4xy + yz is positive semi-definite: principal minors,,, 0,,, 0 (by, by, by) Calculus for functions of n-variables Partial Differentiation Use (x,, x n ) to denote a point in R n and y = f (x) = f (x,, x n ) denotes a function of these variables Example y = 7x + sin(x x ) Continuity/differentiability is a little harder than in -dimension Give rough verbal descriptions of why Definition The partial derivative of the function f (x,, x n ) with respect to x i is the limit f f (x =,, x i + h,, x n ) f (x,, x n ) lim, x i h 0 h if the limit exists Often write f xi as a shorthand Idea: hold all variables constant and look at small changes in function related to small changes in one variable at a time Picture of surface and slopes Examples Let f (x, x, x ) = 7x + x x, then f x = 7, Let f (x, x ) = 7x + sin(x x f ), then x Notice how the chain-rule still applies f x = x and f x = x x = 7 + x cos(x x ) and f x = x cos(x x ) 5

Marginal-product functions Suppose that y = f (x,, x n ) represents the output of an company/industry/economy which depends on n inputs x,, x n The partial derivative f xi represents (approximately) the change in output that results from an increase of in the input x i, if all other inputs are unchanged This approximation will be better if the units are smaller f xi is thus the marginal product function with respect to the variable x i Examples Let y = 8x / x / be the production function (typically decreasing powers gives decreasing returns from more input - law of decreasing marginal productivity) Then the marginal production functions are f x = 6x / x /, f x = 9x x / General Cobb Douglas production function y = f (K, L) = AK α L β where A > 0, 0 < α, β < are constants dependent on the technology/industry Here K represents available capital and L available labor Then f K = αak α L β, f L = βak α L β Observe: as quantity of labor increases, the marginal benefit of the increase reduces; as capital increases the marginal benefit of it decreases (Ie doubling K does not result in twice the output) Conversely, increasing K leads to higher marginal product of L - labor benefits from extra capital Thus the inputs are complimentary The constant elasticity of substitution (CES) production function is the following y = A[δx r + ( δ)x r ] /r ], where A > 0, 0 < δ <, r > After some nasty calcs, observe that y x = δ A r ( y x ) r+, y x = δ A r ( y x ) r+ Time-dependence Suppose that the inputs x,, x n are dependent on time The chain rule then says that y (t) = f x x + + f xn x n Examples y = 7x + x where x = t, x = t Then calc Etc Second-order partial derivatives Can differentiate multiple times The partial derivative f xi x j is simply x j f xi differentials exist, one may keep differentiating forever Supposing that all Example Find derivs up to order of f (x, x ) = x x 4x x 6

Observe that the mixed partials satisfy Young s theorem: if the derivatives are continuous then the order one calculates derivatives is irrelevant; eg f x x = f x x Definition The gradient vector of a function of n variables f (x) is the vector (of length n) of partial derivatives f = Fact: gradient vector points orthogonally away from surface Example Definition The Hessian of a function f (x) is the n n matrix of second order partial derivatives f x x f x x n f = f xn x f xn x n Notice a function of n variables has n partial deriavtives of order and n of order (indeed n k of order k) Example Interpretation: f x + f x is approx of what f looks like in the x direction Ie f is the slope, and the second derivative indicates whether the function curves up or down Examples The Cobb-Douglas production function y = f (x, x ) = Ax αxβ with two inputs has Hessian ( ) α(α )Ax α f x β αβax α x β αβax α x β β(β )Ax α xβ Give 0 < α, β <, the signs are thus ( + + ) That f and f are negative says that as one increases one input, the marginal product of that variable decreases Conversely f > 0 says that as one increases on variable, the marginal product of the other increases First order total differential Used to writing dy dx = f (x) from which we get d f = f (x)dx This is useful for chain rule calculations, etc; d f = f (x)dx is the change of variable formula for integrals We think about this as d f measuring the size of a small change in f resulting from the small change dx in x An analogous thing can be done in multi-variable calculus: Definition 4 The first order total differential of a function y = f (x,, x n ) is dy = d f = f x dx + + f x n dx n = f dx + + f n dx n The total derivative allows us to calculate (approximately) the resulting small change in the value of the function f given small changes dx,, dx n in the inputs The approx is better if the changes are small compared to the size of f itself Example Let y = f (x, x ) = 60x / x / Then dy = 0x / x / dx + 0x / x / dx 7

Suppose that the inputs are currently x = 9, x = 8 and that we increase each by x by and x by (dx =, dx = ) Then the approximate change in y is dy = 0 + 0 = 0 + 0 = 50 4 Note f (9, 8) = 60 and f (0, 0) = 408775 60 + 50, so approximation is quite good Advantage is that the approx is easy to calculate once you have the partial derivative values, can plug in any changes in the inputs after this Implicit Functions Can use total derivatives to calculate slopes of implicitly defined functions F(x, y) = 0 F x dx + F y dy = 0 dy dx = F x F y Examples x + y = 5 at (, ) do as a function x sin y + y sin x = π at (π/, π/) can t repeat as a function Level curves Let y = f (x,, x n ) The set of points (x,, x n ) for which y is a constant is a level curve (n = ) or level set n Suppose y = f (x, x ) is constant Then 0 = d f = f dx + f dx The level curve thus has slope f / f Examples Find the slopes of the level curves of the function f (x, x ) = x x Here 0 = f dx + f dx = x x dx + x x dx dx dx = x x Isoquants If y = f (x, x ) is a production function for inputs x, x, then the level curves of y are described as isoquants curves of equal quantity Ie the level curve 00 = f (x, x ) is the curve whose points describe all the possible input combinations that will product 00 output It is common to want to substitute one input for another in order to produce the same outcome The rate at which one can do this is called the Marginal Rate of Technical Substitution and is defined as follows: MRTS = dx = f dx f Example Find the MRTS for the production function y = f (x, x ) = x / x / Here MRTS = x x If the inputs are currently x = x = 000, then MRTS = Notice (00 998) / = 99999 = (999 00) / Ie if delete one of x need to replace with two of x 8

Notion carries over to more variables Suppose have three variables y = f (x, x, x ) Then MRTS, = dx dx = f f is the marginal rate at which input can be substituted for input Example If y = f (x, x, x ) = x /4 x / x /, then the partial derivatives are f = x/ x / x /4, f = 4x/4 x / x / from which we get the marginal rates of technical susbstitution, f = 6x/4 x / x / MRTS, = f f = 4x x, MRTS, = f f = x x, MRTS, = f f = x x, MRTS, = f f = x 4x, MRTS, = f f = x x, MRTS, = f f = x x Ie if have equal quantities of all three inputs going into the system, then can replace one of x by 4/ of x, or of x with 4 of x, 9