Linear Algebra Notes

Matrices and Systems of Equations

Introduction

  • In 1949, economist Wassily Leontief created a model of equations to understand the interdependence between different sectors of an economy.
  • His original problem involved 500 equations with 500 unknowns, which was computationally intensive for the time, so he reduced it to 42 equations and 42 unknowns.
  • Solving the reduced problem required 56 hours of computation.
  • Leontief's model was a system of linear equations, where unknown quantities are related through linear combinations.
  • The goal was to solve the equations simultaneously to find the unknown values that satisfy the system.
  • Linear models have become vital in various fields like physics, mathematics, engineering, biology, and social sciences due to their simplicity and increased computational capacity.
  • Systems of linear equations represent mathematical models where several unknown quantities are related linearly through various equations.
  • These equations are solved simultaneously to find the values that satisfy all the equations at the same time.
  • Examples of applications include:
    • Water supply networks: Determining pressures and flow rates to plan expansion and ensure proper water distribution.
    • Structural calculations: Modeling and analyzing the resistance and behavior of structures under different loads and forces.
    • Chemical reactions: Determining the amounts of substances present before and after a reaction.
    • Electrical networks: Analyzing the current and voltage at different points.
    • Logistics: Planning and optimizing the production and transportation of goods.
    • Transportation: Planning and optimizing routes, schedules, and vehicle usage.
    • Classification problems: Grouping elements based on their characteristics.
    • Image processing: Reconstructing images or creating graphical animations.
    • Finance: Optimizing portfolios and valuing financial options.
    • Modeling physical phenomena: Simulating the movement of bodies in space and fluid flow in pipes.
    • Game theory: Modeling interactions between economic agents.
    • Modeling biological systems: Simulating population dynamics and interactions between species.
    • Design of mechanical, electrical, and electronic systems; signal processing.

Systems of Linear Equations

  • Definition 1.1 (System of Linear Equations): A linear equation in the variables x<em>1,,x</em>nx<em>1, …, x</em>n can be expressed as a linear combination equal to a constant, bb:
    a<em>1x</em>1++a<em>nx</em>n=b.a<em>1x</em>1 + … + a<em>nx</em>n = b.
  • The variables are multiplied by constants (coefficients a<em>1,,a</em>na<em>1, …, a</em>n) and summed.
  • The independent term and coefficients can be real or complex.
  • A system of linear equations is a collection of one or more equations involving the same variables:
    {a<em>11x</em>1+a<em>12x</em>2++a<em>1nx</em>n=b<em>1 a</em>21x<em>1+a</em>22x<em>2++a</em>2nx<em>n=b</em>2  a<em>m1x</em>1+a<em>m2x</em>2++a<em>mnx</em>n=bm \begin{cases} a<em>{11}x</em>1 + a<em>{12}x</em>2 + … + a<em>{1n}x</em>n = b<em>1 \ a</em>{21}x<em>1 + a</em>{22}x<em>2 + … + a</em>{2n}x<em>n = b</em>2 \ … \ a<em>{m1}x</em>1 + a<em>{m2}x</em>2 + … + a<em>{mn}x</em>n = b_m \ \end{cases}
  • If the independent terms are all zero (b<em>1=b</em>2==bn=0b<em>1 = b</em>2 = … = b_n = 0), the system is homogeneous.
  • Example 1.1 (Water Supply Network): In a water supply network with 4 nodes, the flow through each pipe is determined by establishing equilibrium at each node (flow in = flow out). Example is shown to derive the equations.

Solution of a System of Equations

  • Definition 1.2 (Solution of a System of Linear Equations): A list of values (x<em>1,,x</em>n)T(x<em>1, …, x</em>n)^T that satisfies all equations in the system.
  • A system can have:
    • No solution (incompatible).
    • Exactly one solution (compatible determined).
    • Infinitely many solutions (compatible indeterminate).
  • Definition 1.3 (Equivalent Systems of Linear Equations): Two systems are equivalent if they have the same set of solutions.
  • Theorem 1.1: Basic operations that produce equivalent systems include:
    • Multiplying an equation by a non-zero constant.
    • Adding or subtracting equations.
    • Replacing an equation with a linear combination of the other equations.
2x2 Systems
  • Linear equations with two variables represent lines in the plane.
  • Systems of two equations with two unknowns can represent:
    • Parallel lines (incompatible, no solution).
    • Intersecting lines (compatible determined, one solution).
    • Overlapping lines (compatible indeterminate, infinite solutions).
  • Examples with equations and graphs are given for each case.

Matrix Representation of a System

  • A linear system of mm equations with nn variables can be represented by a matrix AA of coefficients (size m<br/>nm <br />\newline n):
    {a<em>11x</em>1+a<em>12x</em>2++a<em>1nx</em>n=b<em>1 a</em>21x<em>1+a</em>22x<em>2++a</em>2nx<em>n=b</em>2  a<em>m1x</em>1+a<em>m2x</em>2++a<em>mnx</em>n=b<em>m A=(a</em>11amp;a<em>12a</em>1n a<em>21a</em>22amp;amp;a<em>2n  a</em>m1amp;a<em>m2a</em>mn)\begin{cases} a<em>{11}x</em>1 + a<em>{12}x</em>2 + … + a<em>{1n}x</em>n = b<em>1 \ a</em>{21}x<em>1 + a</em>{22}x<em>2 + … + a</em>{2n}x<em>n = b</em>2 \ … \ a<em>{m1}x</em>1 + a<em>{m2}x</em>2 + … + a<em>{mn}x</em>n = b<em>m \ \end{cases} \Rightarrow A = \begin{pmatrix} a</em>{11} &amp; a<em>{12} & … & a</em>{1n} \ a<em>{21} & a</em>{22} &amp; … &amp; a<em>{2n} \ … & … & … & … \ a</em>{m1} &amp; a<em>{m2} & … & a</em>{mn} \end{pmatrix}
  • Using the vector of independent terms b=(b<em>1,b</em>2,,b<em>m)Tb = (b<em>1, b</em>2, …, b<em>m)^T and the vector of variables x=(x</em>1,x<em>2,,x</em>n)Tx = (x</em>1, x<em>2, …, x</em>n)^T, the system can be written as Ax=bAx = b, where A=(a<em>1,a</em>2,,a<em>n)A = (a<em>1, a</em>2, …, a<em>n) and a</em>ia</em>i are the columns of matrix AA.
  • Definition 1.4 (Product of a Matrix by a Vector): Ax=A(x<em>1 x</em>2  x<em>n)=(a</em>11amp;a<em>12a</em>1n a<em>21a</em>22amp;amp;a<em>2n  a</em>m1amp;a<em>m2a</em>mn)(x<em>1 x</em>2  x<em>n)=x</em>1a<em>1+x</em>2a<em>2++x</em>nan.Ax = A \begin{pmatrix} x<em>1 \ x</em>2 \ … \ x<em>n \end{pmatrix} = \begin{pmatrix} a</em>{11} &amp; a<em>{12} & … & a</em>{1n} \ a<em>{21} & a</em>{22} &amp; … &amp; a<em>{2n} \ … & … & … & … \ a</em>{m1} &amp; a<em>{m2} & … & a</em>{mn} \end{pmatrix} \begin{pmatrix} x<em>1 \ x</em>2 \ … \ x<em>n \end{pmatrix} = x</em>1a<em>1 + x</em>2a<em>2 + … + x</em>na_n.
  • The augmented matrix combines the coefficient matrix and the vector of independent terms:
    (Ab)=(a<em>11a</em>12amp;amp;a<em>1nb</em>1 a<em>21a</em>22amp;amp;a<em>2nb</em>2 amp;amp;amp;amp; a<em>m1a</em>m2amp;amp;a<em>mnb</em>m)(A|b) = \begin{pmatrix} a<em>{11} & a</em>{12} &amp; … &amp; a<em>{1n} & b</em>1 \ a<em>{21} & a</em>{22} &amp; … &amp; a<em>{2n} & b</em>2 \ … &amp; … &amp; … &amp; … &amp; … \ a<em>{m1} & a</em>{m2} &amp; … &amp; a<em>{mn} & b</em>m \end{pmatrix}
  • Definition 1.5 (Rank of a Matrix): The rank of a matrix is the maximum number of linearly independent rows or columns.
  • Definition 1.6: The number of degrees of freedom is given by nRank(A)n - Rank(A).
Compatibility of a System
  • The system of equations can be expressed as a linear combination of the columns of the coefficient matrix:
    x<em>1(a</em>11 a<em>21  a</em>m1)+x<em>2(a</em>12 a<em>22  a</em>m2)++x<em>n(a</em>1n a<em>2n  a</em>mn)=(b<em>1 b</em>2  b<em>m)x</em>1a<em>1+x</em>2a<em>2++x</em>nan=b.x<em>1 \begin{pmatrix} a</em>{11} \ a<em>{21} \ … \ a</em>{m1} \end{pmatrix} + x<em>2 \begin{pmatrix} a</em>{12} \ a<em>{22} \ … \ a</em>{m2} \end{pmatrix} + … + x<em>n \begin{pmatrix} a</em>{1n} \ a<em>{2n} \ … \ a</em>{mn} \end{pmatrix} = \begin{pmatrix} b<em>1 \ b</em>2 \ … \ b<em>m \end{pmatrix} \Rightarrow x</em>1a<em>1 + x</em>2a<em>2 + … + x</em>na_n = b.
  • The vector of independent terms, bb, is a linear combination of the columns of AA.
  • The solution consists of the coefficients x<em>1,,x</em>nx<em>1, …, x</em>n that allow obtaining bb as a linear combination of the columns of AA.
  • The system is compatible if bb is a linear combination of the columns of AA.
  • Theorem 1.2 (Rouche-Frobenius): A system of linear equations Ax=bAx = b, where AA is an m<br/>nm <br />\newline n matrix, is compatible if the rank of AA equals the rank of the augmented matrix. If they coincide and the rank is nn (the number of unknowns), the system is compatible determined.

Solving Linear Systems

Gaussian Elimination
  • Gaussian elimination transforms the original system into an equivalent one with an upper triangular coefficient matrix (coefficients below the diagonal are zero).
  • Starting from the last rows, the variables are easily solved to obtain the solution.
  • Permissible operations include:
    • Interchanging equations.
    • Multiplying an equation by a non-zero constant.
    • Replacing an equation with a linear combination of itself and another equation.
  • The operations are performed on the augmented matrix.
  • Example 1.2 (Compatible Determined System): Example of solving a 3x3 system using gaussian elimination, showing steps on the augmented matrix and the resulting unique solution x<em>1=3,x</em>2=12,x3=2x<em>1 = 3, x</em>2 = 12, x_3 = 2.
  • Example 1.3 (Compatible Indeterminate System): Example of solving a 3x3 system with infinite solutions, showing steps on the augmented matrix and the resulting solution set depending on a parameter ρ\rho: x<em>1=5ρ,x</em>2=10+ρ,x3=ρx<em>1 = 5 - \rho, x</em>2 = 10 + \rho, x_3 = \rho.
  • Example 1.4 (Incompatible System): Example of solving a 3x3 system with no solution, showing steps on the augmented matrix and the resulting contradiction 0=80 = -8.
Gauss-Jordan Elimination
  • Gauss-Jordan elimination makes zeros above and below the diagonal until a reduced echelon matrix is reached.
  • Definition 1.7 (Echelon Form): A rectangular matrix is in echelon form if:
    • All non-zero rows are above rows of all zeros.
    • The leading entry of a row is to the right of the leading entry of the row above it.
    • Entries below a leading entry are zero.
      (amp;amp;amp;amp; 0amp;amp;amp;amp; 0amp;0amp;amp;amp; 0amp;0amp;0amp;0amp;0 0amp;0amp;0amp;0amp;0 )\begin{pmatrix} * &amp; * &amp; * &amp; * &amp; * \ 0 &amp; * &amp; * &amp; * &amp; * \ 0 &amp; 0 &amp; * &amp; * &amp; * \ 0 &amp; 0 &amp; 0 &amp; 0 &amp; 0 \ 0 &amp; 0 &amp; 0 &amp; 0 &amp; 0 \ \end{pmatrix}
  • If a matrix in echelon form also satisfies:
    • The leading entry of each non-zero row is 1.
    • Each leading 1 is the only non-zero entry in its column.
      Then it is in reduced echelon form.
      (1amp;0amp;0amp;amp; 0amp;1amp;0amp;amp; 0amp;0amp;1amp;amp; 0amp;0amp;0amp;0amp;0 0amp;0amp;0amp;0amp;0 )\begin{pmatrix} 1 &amp; 0 &amp; 0 &amp; * &amp; * \ 0 &amp; 1 &amp; 0 &amp; * &amp; * \ 0 &amp; 0 &amp; 1 &amp; * &amp; * \ 0 &amp; 0 &amp; 0 &amp; 0 &amp; 0 \ 0 &amp; 0 &amp; 0 &amp; 0 &amp; 0 \ \end{pmatrix}
  • The leading entries form a staircase pattern.
  • Every matrix is row equivalent to a unique reduced echelon matrix.
  • Definition 1.8 (Pivot): A pivot is a non-zero number on the main diagonal used to create zeros below it.
  • Example 1.5 (Water Supply Network Solution): Solving the system from Example 1.1 using Gauss-Jordan to derive the general solution depending on parameter ρ\rho:
    (x<em>1 x</em>2 x<em>3 x</em>4)=(70 30 20 0)+ρ(1 1 1 1)\begin{pmatrix} x<em>1 \ x</em>2 \ x<em>3 \ x</em>4 \end{pmatrix} = \begin{pmatrix} 70 \ 30 \ 20 \ 0 \end{pmatrix} + \rho \begin{pmatrix} 1 \ -1 \ 1 \ 1 \end{pmatrix}

Matrices

  • Matrices are fundamental in linear algebra and used in science, engineering, computing, and other fields.
  • A matrix is an ordered collection of numbers arranged in rows and columns, used to represent systems of linear equations, linear transformations, data, and other mathematical objects.
  • Definition 1.9 (Matrix): A matrix stores data in a rectangular form using rows and columns. If AA is an mm row and nn column matrix, we denote it by A=(a<em>ij)A = (a<em>{ij}), where a</em>ija</em>{ij} is the coefficient in row ii, column jj.

Matrix Operations

Sum
  • Definition 1.10 (Sum): If A=(a<em>ij)A = (a<em>{ij}) and B=(b</em>ij)B = (b</em>{ij}) are matrices of the same dimensions, the sum is A+B=(a<em>ij+b</em>ij)A + B = (a<em>{ij} + b</em>{ij}).
  • Example 1.6 (Matrix Sum): Calculating the total sales of a company by adding the sales from two warehouses, represented by matrices AA and BB.
Scalar Product
  • Definition 1.11 (Scalar Product): If A=(a<em>ij)A = (a<em>{ij}) is a matrix and ρ\rho is a scalar (R or C), the product of a matrix by a scalar is ρA=(ρa</em>ij)\rho A = (\rho a</em>{ij}).
  • Example 1.7 (Scalar Product): Increasing salaries of employees by 10% using scalar product on the matrix of salaries AA.
Properties of Matrix Operations
  • Theorem 1.3: If AA, BB, and CC are matrices of the same dimension, with scalars ρ<em>1\rho<em>1 and ρ</em>2\rho</em>2, then:
    • Commutative: A+B=B+AA + B = B + A.
    • Associative: A+(B+C)=(A+B)+CA + (B + C) = (A + B) + C.
    • Additive identity: A+0=AA + 0 = A, where 00 is the zero matrix.
    • Distributive: (ρ<em>1+ρ</em>2)A=ρ<em>1A+ρ</em>2A(\rho<em>1 + \rho</em>2)A = \rho<em>1A + \rho</em>2A.
    • Distributive: ρ<em>1(A+B)=ρ</em>1A+ρ1B\rho<em>1(A + B) = \rho</em>1A + \rho_1B.
    • Associative: ρ<em>1(ρ</em>2A)=(ρ<em>1ρ</em>2)A\rho<em>1(\rho</em>2A) = (\rho<em>1\rho</em>2)A.
Matrix Product
  • Definition 1.12 (Matrix Product): Given matrices AA and BB such that the number of columns of AA equals the number of rows of BB, the product is defined as AB=(Ab<em>1,Ab</em>2,,Ab<em>n)AB = (Ab<em>1, Ab</em>2, …, Ab<em>n), where the columns are the product of matrix AA with the columns of BB (as in Definition 1.5). The element (i,j)(i, j) of the product is (AB)</em>ij=a<em>i1b</em>1j+a<em>i2b</em>2j++a<em>inb</em>mj(AB)</em>{ij} = a<em>{i1}b</em>{1j} + a<em>{i2}b</em>{2j} + … + a<em>{in}b</em>{mj}.
  • Example 1.8 (Transformations in the Plane): The effect of some matrices on vectors in the plane:
    • The matrix B=(1amp;0 0amp;1)B = \begin{pmatrix} 1 &amp; 0 \ 0 &amp; -1 \end{pmatrix} calculates the symmetric vector with respect to the x-axis.
      (1amp;0 0amp;1)(x<em>1 x</em>2)=(x<em>1 x</em>2)\begin{pmatrix} 1 &amp; 0 \ 0 &amp; -1 \end{pmatrix} \begin{pmatrix} x<em>1 \ x</em>2 \end{pmatrix} = \begin{pmatrix} x<em>1 \ -x</em>2 \end{pmatrix}
    • The matrix A=(1amp;0 0amp;1)A = \begin{pmatrix} -1 &amp; 0 \ 0 &amp; 1 \end{pmatrix} calculates the symmetric vector with respect to the y-axis.
      (1amp;0 0amp;1)(x<em>1 x</em>2)=(x<em>1 x</em>2)\begin{pmatrix} -1 &amp; 0 \ 0 &amp; 1 \end{pmatrix} \begin{pmatrix} x<em>1 \ x</em>2 \end{pmatrix} = \begin{pmatrix} -x<em>1 \ x</em>2 \end{pmatrix}
    • To obtain the symmetric point with respect to the origin, multiply:
      (1amp;0 0amp;1)(x<em>1 x</em>2)=(x<em>1 x</em>2)\begin{pmatrix} -1 &amp; 0 \ 0 &amp; -1 \end{pmatrix} \begin{pmatrix} x<em>1 \ x</em>2 \end{pmatrix} = \begin{pmatrix} -x<em>1 \ -x</em>2 \end{pmatrix}
Properties of the Matrix Product
  • Theorem 1.4: If AA, BB, and CC are matrices (with suitable dimensions for the operations) and ρ\rho is a scalar, then:
    • A(BC)=(AB)CA(BC) = (AB)C.
    • A(B+C)=AB+ACA(B + C) = AB + AC.
    • (A+B)C=AC+BC(A + B)C = AC + BC.
    • ρ(AB)=ρ(A)B=A(ρB)\rho(AB) = \rho(A)B = A(\rho B).
  • The matrix product is not commutative (ABBAAB \neq BA).
Special Matrices
  • Definition 1.13:
    • Identity: A square matrix with ones on the diagonal and zeros elsewhere.
      I=(1amp;0amp;amp;0amp;0 0amp;1amp;amp;0amp;0 amp;amp;amp;amp; 0amp;0amp;amp;1amp;0 0amp;0amp;amp;0amp;1 )I = \begin{pmatrix} 1 &amp; 0 &amp; … &amp; 0 &amp; 0 \ 0 &amp; 1 &amp; … &amp; 0 &amp; 0 \ … &amp; … &amp; … &amp; … &amp; … \ 0 &amp; 0 &amp; … &amp; 1 &amp; 0 \ 0 &amp; 0 &amp; … &amp; 0 &amp; 1 \ \end{pmatrix}
    • Zero: A matrix where all coefficients are zero.
      0=(0amp;0amp;amp;0amp;0 0amp;0amp;amp;0amp;0 amp;amp;amp;amp; 0amp;0amp;amp;0amp;0 0amp;0amp;amp;0amp;0 )0 = \begin{pmatrix} 0 &amp; 0 &amp; … &amp; 0 &amp; 0 \ 0 &amp; 0 &amp; … &amp; 0 &amp; 0 \ … &amp; … &amp; … &amp; … &amp; … \ 0 &amp; 0 &amp; … &amp; 0 &amp; 0 \ 0 &amp; 0 &amp; … &amp; 0 &amp; 0 \ \end{pmatrix}
    • Transpose: Given an m<br/>nm <br />\newline n matrix AA, its transpose ATA^T is an n<br/>mn <br />\newline m matrix whose columns are the rows of AA.
    • Orthogonal: A matrix AA is orthogonal if AAT=ATA=IAA^T = A^TA = I.
    • Symmetric: A matrix AA such that A=ATA = A^T (coefficients a<em>ij=a</em>jia<em>{ij} = a</em>{ji} coincide).
  • Theorem 1.5: If AA is an mnm \newline n matrix, then:
    • AI<em>n=AAI<em>n = A, where I</em>nI</em>n is the n-sized identity matrix.
    • I<em>mA=AI<em>mA = A, where I</em>mI</em>m is the m-sized identity matrix.

Elementary Matrices

  • Operations in Gaussian elimination can be performed by multiplying by elementary matrices, obtained from the identity matrix by:
    • EijE_{ij}: Interchanging rows ii and jj.
    • Ei(λ)E_i(\lambda): Multiplying row ii by λ\lambda.
    • Eij(λ)E_{ij}(\lambda): Adding λ\lambda times row ii to row jj.
  • Example 1.9: Effects from premultiplying matrix A by elementary matrices:
    • Swapping two rows.
    • Multiplying a row by a scalar.
    • Adding a multiple of one row to another row.
  • Elementary matrices that premultiply a matrix perform the same row operations used to construct the elementary matrix. Post-multiplication performs column operations.
  • Gaussian elimination becomes a sequence of elementary matrices premultiplying the augmented matrix.
  • Definition 1.14 (Euclidean Scalar Product): Given vectors a=(a<em>1,,a</em>n)Ta = (a<em>1, …, a</em>n)^T and b=(b<em>1,,b</em>n)Tb = (b<em>1, …, b</em>n)^T, the Euclidean scalar product is ab=aTb=(a<em>1,,a</em>n)(b<em>1  b</em>n)=a<em>1b</em>1++a<em>nb</em>na • b = a^T b = (a<em>1, …, a</em>n) \begin{pmatrix} b<em>1 \ … \ b</em>n \end{pmatrix} = a<em>1b</em>1 + … + a<em>nb</em>n.
  • The element (i,j)(i, j) of the matrix product ABAB is the scalar product of the ii-th row of AA and the jj-th column of BB.
  • Definition 1.15 (Power of a Matrix): If AA is a square matrix and kk a positive integer, then AkA^k means multiplying k copies of AA. A0=IA^0 = I.
  • Example 1.10 (Bird Population): Shows example of modeling a bird population with two age groups by expressing current and next-year populations with matrix operations
  • Theorem 1.6: Given suitably sized matrices A,BA, B, and CC compatible for multiplication:
    • (A+B)T=AT+BT(A + B)^T = A^T + B^T.
    • (AB)T=BTAT(AB)^T = B^TA^T.
    • (AT)T=A(A^T)^T = A.

Matrix Inverse

  • Definition 1.16 (Matrix Inverse): Given a square matrix AA, its inverse, denoted A1A^{-1}, is a matrix such that AA1=A1A=IAA^{-1} = A^{-1}A = I. If AA is orthogonal (ATA=IA^TA = I), then AT=A1A^T = A^{-1}.
  • Definition 1.17 (Invertible Matrix): An n×nn \times n matrix AA is invertible if its inverse exists. If it exists, it is unique. A non-invertible matrix is singular, otherwise, it is non-singular.
  • Theorem 1.7: If AA is an invertible matrix, the following statements are equivalent:
    • The system Ax=bAx = b has a unique solution, x=A1bx = A^{-1}b.
    • The homogeneous system Ax=0Ax = 0 has the unique solution x=0x = 0.
    • The rows and columns of AA are linearly independent.
  • Theorem 1.8:
    • If AA is invertible, so is A1A^{-1}, and (A1)1=A(A^{-1})^{-1} = A.
    • If AA and BB are invertible and have compatible dimensions, then (AB)1=B1A1(AB)^{-1} = B^{-1}A^{-1}.
    • If AA is invertible, so is ATA^T, and (AT)1=(A1)T(A^T)^{-1} = (A^{-1})^T.
    • If AA is invertible, (Ak)1=(A1)k(A^k)^{-1} = (A^{-1})^k.
  • Example 1.11 (Electrical Circuit): Example of equations derived from Ohm and Kirchhoff laws and then used for expressing the solution matrix notation
Finding Inverse using Gaussian Elimination
  • For an invertible coefficient matrix from Example 1.11, the inverse is calculated by applying Gaussian elimination until reaching the identity matrix, performing the same operations on the identity matrix to obtain the inverse.
Alternative Formulas to Calculate the Inverse
  • For a 2x2 matrix A=(aamp;b camp;d)A = \begin{pmatrix} a &amp; b \ c &amp; d \end{pmatrix}, if adbc0ad - bc \neq 0, then A1=1adbc(damp;b camp;a)A^{-1} = \frac{1}{ad - bc} \begin{pmatrix} d &amp; -b \ -c &amp; a \end{pmatrix}.
  • The term adbcad - bc is the determinant of the matrix.

Determinants

Introduction

  • In linear algebra history, determinants appeared for resolution of system of equations even before the matrix concept.
  • Defined as a number associated with a square matrix, allowing us to check for invertibility or not. Vital in the resolution and interpretation of linear system of equations.

Determinant of a Matrix

  • From solving a 2-equation 2-unknown system. If a<em>110a<em>{11} \neq 0, a</em>11x<em>1+a</em>12x<em>1=b</em>1;a<em>21x</em>1+a<em>22x</em>1=b2a</em>{11}x<em>1+ a</em>{12}x<em>1 = b</em>1 ; a<em>{21}x</em>1+ a<em>{22}x</em>1 = b_2
  • After Gauss elimination, we can see now the solution is: x<em>1=b</em>1a<em>22b</em>2a<em>11a</em>22a<em>11a</em>12a<em>21;x</em>2=b<em>2a</em>11b<em>1a</em>21a<em>22a</em>11a<em>12a</em>21x<em>1 = \frac{b</em>1a<em>{22} - b</em>2a<em>{11}}{a</em>{22}a<em>{11} - a</em>{12}a<em>{21}} ; x</em>2 = \frac{b<em>2a</em>{11} - b<em>1a</em>{21}}{a<em>{22}a</em>{11} - a<em>{12}a</em>{21}}
  • Note that there is single solution if a<em>22a</em>11a<em>12a</em>210a<em>{22}a</em>{11} - a<em>{12}a</em>{21} \neq 0.
  • Definition 2.1: Given a matrix A, The determinant is. Det(A)=A=a<em>22a</em>11a<em>12a</em>21Det(A)=|A| =a<em>{22}a</em>{11} - a<em>{12}a</em>{21}
  • Definition 2.2: For any matrix A, define it's minor associated with a<em>ija<em>{ij}, denoted by A</em>ij|A</em>{ij}|, which is the determinant of the submatrix AijA_{ij}, obtained by removing the i row and j column.
  • and C<em>ij=(1)i+jA</em>ijC<em>{ij} = (-1)^{i+j}|A</em>{ij}| is called Cofactor
  • The definition allows defining it recursively of ANY square matrix:
  • Definition 2.3: A determinant of A is the sum of entries of any row of A multiplied by the cofactors, det(A)=A=<em>j=1na</em>ijC<em>ij=a</em>i1C<em>i1+a</em>i2C<em>i2++a</em>inC<em>indet(A) = |A| = \sum<em>{j=1}^n a</em>{ij}C<em>{ij} = a</em>{i1}C<em>{i1} + a</em>{i2}C<em>{i2} + … + a</em>{in}C<em>{in}, or the sum of any columndet(A)=A=</em>i=1na<em>ijC</em>ij=a<em>1jC</em>1j+a<em>2jC</em>2j++a<em>njC</em>njdet(A) = |A| = \sum</em>{i=1}^n a<em>{ij}C</em>{ij} = a<em>{1j}C</em>{1j} + a<em>{2j}C</em>{2j} + … + a<em>{nj}C</em>{nj}
  • Example 2.1 In the previous example of electrical circuit, can verify it is invertible because determinant is not zero: |A| = 8 !=0 (development with first row cofactors)
  • To compute determinants of 3 by 3 can use Sarrus rule. for any A=(a<em>ij)A =(a<em>{ij}) then det(A)=;a</em>11a<em>22a</em>33+a<em>12a</em>23a<em>31+a</em>13a<em>21a</em>32a<em>31a</em>22a<em>13a</em>32a<em>23a</em>11a<em>33a</em>21a12det(A) = ; a</em>{11}a<em>{22}a</em>{33} + a<em>{12}a</em>{23}a<em>{31} + a</em>{13}a<em>{21}a</em>{32} - a<em>{31}a</em>{22}a<em>{13} - a</em>{32}a<em>{23}a</em>{11} - a<em>{33}a</em>{21}a_{12}

Additional Properties

  • If A is Diagonal, then the Determinant is just the diagonal elements multiplied.
  • If A is a squared matrix, then:
    • If a multiple is added to a row the determinant does not change.
    • If two rows are swapped, then determinant times -1.
    • If any row is multiplied by a constant, then the determinant is multiplied by the constant.
    • Also det(A)=det(AT)det(A) = det(A^T)
    • If another B exists $det(AB) = det(A)det(B)
    • But det(A+B)det(A)+det(B)det(A+B) \neq det(A) + det(B)
    • Given an scalar rho, det(ρA)=ρndet(A)det(\rho A ) = \rho^n det(A)

More Theorems

  • Squared Matrix A is invertible if and only if the Determinant is Not Zero
  • If A is a Squared matrix and if one of those conditions is valid, then the Determinant equals 0 :
    • The row is fully zero.
    • Two Rows equal.
    • One row as linear combination of the rest of rows.
  • Example 2.2 A message is encoded numerically. Use Matrix multiplication to hide the message.

Rank Of a Matrix

  • Definition 2.4: The range is the biggest size minor after trying the non-null minors submatrix
  • Theorem 2.5: Range does not change if:
    • Rows or Columns are interchanged.
    • Rows or Columns are transposed.
    • Row/Column times scalar number.
    • Row/Column sum of the Row/Column times scalar number..

Computing the Range

  • Orlado's Method.
    • Gauss triangulation.

Use for Inverses With Determinants

  • Definition 2.5: Given A, its adjoint matrix adj(A)<em>ij=(1)i+jA</em>ijadj(A)<em>{ij} = (-1)^{i+j}|A</em>{ij}|.
  • Theorem 2.6: Given A, then the Invert is A1=1Aadj(A)TA^{-1} = \frac{1}{|A|}adj(A)^T

Vectorial Spaces

Introduction

  • Vector operations can be represented as directions with Length and Direction.
    Operations such as sums and scaling functions (multiplying the vector with a scalar number are defined
  • Definition 3.1 Vectorial Space : Given a set V operations internal + and external *. Such that an element in K(R or C) times V results again in V, then it says V is a vectorial space or vector space, where:
    • Respecting Internal Operation (+) : u+v element of V ; u + v = v + u ; u+(v+w) = (u+v)+w
      V has a Zero such that u+0=u for all Elements in V, V has a negative such that u + (-u) = 0
    • Respecting External Operation (*) : a * u is in Va *(u + v) = a * u + a * v ; (a + B)* u = a * u + B * u
      a(B * u) = (aB)*u
      The Neutral elements results again in u, 1 * u = u , The opposed of elements.
  • A Vectorial space always has One 0 and just One opposed element

Examples and Non-Examples

  • Example 3.1: Real Numbers N-tuple with sum and Scalar Product, Complex numbers (same), Also the set fo Matrix mXn has vector behavior, Polynomials of Grade minor than or Equal to 2, or Continuos functions on a point a/b .
  • Example 3.2: Matrix 2X2 is not valid, Polynomias Grade EQUAL to 2 NOT valid.
    The zero is no valid.

Sub-Spaces

  • Definition 3.2: A not empty V subset of, where the Internal Sum and External escalar product are defined is a sub-space if they have vectorial Space with the SAME operations.
  • Checking Vectorial Space W!=0, W is subspace of V if, If u, v inside W, then u + V in W, If inside W scalar a*u inside W.
    Zero Siempre in subspace, If a vector and its opposite exist in the subspace then by add they give zero.
  • Example 3.3: Test if W is a subspace in R3. A test example. Comply those two theoream valid then W is Subspace.
    Theorem3.2 The generated combination for ANY number of vectors inside a vectorial Space ,is still a vectorial Sub-space.

Nul Space

Example 3.7 A flux across streets needs solving:
System description - Question- What King of Solutions has this system? What grades of freedom ,And what if you let variables free?.
-Its solution results in the next vectorial subSpace.
Definition 3.3 The nul spacio of a matrix of dimensions m x n, N ul(A), are all the set of a Homogeneus system Ax =0
Theorem 3.3 Te nul spacio of a matrix has dimension m x n, N ul(A), is a vectorial subSpace f Kn

Linear Dependency

  • **Definition 3.4 Combination lineal de vectores **. A Vector inside V is combination lineal. v = r1u1 + r2 u2...
  • Definition 3.5 vectores Lineal Independientes They are lineally independent if, in any combination lineal that equals the zero vector ,the sole possible results are Zero Coefficient.
  • If not this way, those are define s Lineal Dependency.

Base and Dimensions

  • Definition 3.6 Generator System. Are those that form a Base, when a vector is combination Lineal.
  • Definition 3.7 Base. Set fo Vectors inside vectorial space called System, when
    Are independent. Generate
    Then any generator vector to find a BAse is needed find the Higher number of independent Vectors
    Theorem 3.4 Any basis inside some vectorial space, has a Same Number of Vectors ,all this
    Definition 3.8 The canonical Base of vectorial Space ,such base is:
    Definition 3.9 Coordinates. For vectorial V, If B vectors are Base ,and V can be described v=p1u1+pnunv = p1u1…+ pnun Then p its coordinates given B
    If you have and, you need a Base, it requires finding what maximal independents vectors exist. and calculate de Range with new vetors, adding.
    , If any, that is not, add others vectos find the Lineal independents for generate the B.

Intersections and adds Inside vectoriale Space.

  • Add : vector that is described sum in subspaces. for those Subdpace if any, It will give its generating set,
  • Sum:. Find independent base there inside the sum set, has high independents.

Applications For linear spaces transformations