Skip to content

Matrix

Definition: A matrix is a rectangular array of numbers, symbols, or expressions arranged in rows and columns.

General Form:

A=[a11a12a13a1na21a22a23a2na31a32a33a3nam1am2am3amn]A = \begin{bmatrix} a_{11} & a_{12} & a_{13} & \cdots & a_{1n} \\ a_{21} & a_{22} & a_{23} & \cdots & a_{2n} \\ a_{31} & a_{32} & a_{33} & \cdots & a_{3n} \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ a_{m1} & a_{m2} & a_{m3} & \cdots & a_{mn} \end{bmatrix}

Notation:

  • Matrix AA of order m×nm \times n (m rows, n columns)
  • Element in ithi^{th} row and jthj^{th} column: aija_{ij}
  • A=[aij]m×nA = [a_{ij}]_{m \times n}

Order/Dimension: m×nm \times n where mm = number of rows, nn = number of columns

Based on Order:

  • a) Row Matrix: Matrix with only one row (1×n1 \times n)

A=[2579]A = \begin{bmatrix} 2 & 5 & 7 & 9 \end{bmatrix}

  • b) Column Matrix: Matrix with only one column (m×1m \times 1)

B=[372]B = \begin{bmatrix} 3 \\ 7 \\ 2 \end{bmatrix}

  • c) Square Matrix: Matrix where m=nm = n (equal rows and columns)

C=[1234](2×2 square matrix)C = \begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix} \quad (2 \times 2 \text{ square matrix})

  • d) Rectangular Matrix: Matrix where mnm \neq n

Based on Elements:

  • e) Zero/Null Matrix: All elements are zero. Notation: OO or 00

O=[0000]O = \begin{bmatrix} 0 & 0 \\ 0 & 0 \end{bmatrix}

  • f) Diagonal Matrix: Square matrix where all non-diagonal elements are zero (aij=0 for ij)(a_{ij} = 0 \text{ for } i \neq j)

D=[500030007]D = \begin{bmatrix} 5 & 0 & 0 \\ 0 & 3 & 0 \\ 0 & 0 & 7 \end{bmatrix}

  • g) Scalar Matrix: Diagonal matrix where all diagonal elements are equal

S=[400040004]S = \begin{bmatrix} 4 & 0 & 0 \\ 0 & 4 & 0 \\ 0 & 0 & 4 \end{bmatrix}

  • h) Identity/Unit Matrix: Diagonal matrix where all diagonal elements are 1. Notation: II or InI_n

I3=[100010001]I_3 = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix}

  • i) Upper Triangular Matrix: Square matrix where all elements below the diagonal are zero (aij=0 for i>j)(a_{ij} = 0 \text{ for } i > j)

U=[257038004]U = \begin{bmatrix} 2 & 5 & 7 \\ 0 & 3 & 8 \\ 0 & 0 & 4 \end{bmatrix}

  • j) Lower Triangular Matrix: Square matrix where all elements above the diagonal are zero (aij=0 for i<j)(a_{ij} = 0 \text{ for } i < j)

L=[200530784]L = \begin{bmatrix} 2 & 0 & 0 \\ 5 & 3 & 0 \\ 7 & 8 & 4 \end{bmatrix}

A. Equality of Matrices

  • Two matrices AA and BB are equal if:
    1. They have the same order
    2. Corresponding elements are equal: aij=bija_{ij} = b_{ij} for all i,ji, j

B. Addition and Subtraction

  • Conditions: Matrices must have the same order

  • Addition: A+B=[aij+bij]A + B = [a_{ij} + b_{ij}]

  • Subtraction: AB=[aijbij]A - B = [a_{ij} - b_{ij}]

  • Properties:

    • Commutative: A+B=B+AA + B = B + A
    • Associative: (A+B)+C=A+(B+C)(A + B) + C = A + (B + C)
    • Additive Identity: A+O=AA + O = A
    • Additive Inverse: A+(A)=OA + (-A) = O
  • Example:

[1234]+[5678]=[681012]\begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix} + \begin{bmatrix} 5 & 6 \\ 7 & 8 \end{bmatrix} = \begin{bmatrix} 6 & 8 \\ 10 & 12 \end{bmatrix}

C. Scalar Multiplication

  • Definition: Multiplying every element by a scalar kk

kA=[kaij]kA = [k \cdot a_{ij}]

  • Properties:

    • k(A+B)=kA+kBk(A + B) = kA + kB
    • (k+m)A=kA+mA(k + m)A = kA + mA
    • (km)A=k(mA)(km)A = k(mA)
    • 1A=A1 \cdot A = A
    • 0A=O0 \cdot A = O
  • Example:

3[1234]=[36912]3\begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix} = \begin{bmatrix} 3 & 6 \\ 9 & 12 \end{bmatrix}

D. Matrix Multiplication

  • Condition: Number of columns in first matrix = number of rows in second matrix
  • If AA is m×nm \times n and BB is n×pn \times p, then ABAB is m×pm \times p
  • Formula:

(AB)ij=k=1naikbkj(AB)_{ij} = \sum_{k=1}^{n} a_{ik} \cdot b_{kj}

  • Properties:

    • Generally NOT commutative: ABBAAB \neq BA
    • Associative: (AB)C=A(BC)(AB)C = A(BC)
    • Distributive: A(B+C)=AB+ACA(B + C) = AB + AC
    • Identity: AI=IA=AAI = IA = A
    • (AB)T=BTAT(AB)^T = B^T A^T
  • Example:

[1234][5678]=[1×5+2×71×6+2×83×5+4×73×6+4×8]=[19224350]\begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix} \begin{bmatrix} 5 & 6 \\ 7 & 8 \end{bmatrix} = \begin{bmatrix} 1 \times 5 + 2 \times 7 & 1 \times 6 + 2 \times 8 \\ 3 \times 5 + 4 \times 7 & 3 \times 6 + 4 \times 8 \end{bmatrix} = \begin{bmatrix} 19 & 22 \\ 43 & 50 \end{bmatrix}

  • Definition: Matrix obtained by interchanging rows and columns

  • If A=[aij]m×nA = [a_{ij}]_{m \times n}, then AT=[aji]n×mA^T = [a_{ji}]_{n \times m}

  • Properties:

    • (AT)T=A(A^T)^T = A
    • (A+B)T=AT+BT(A + B)^T = A^T + B^T
    • (kA)T=kAT(kA)^T = kA^T
    • (AB)T=BTAT(AB)^T = B^T A^T
  • Example:

A=[123456]AT=[142536]A = \begin{bmatrix} 1 & 2 & 3 \\ 4 & 5 & 6 \end{bmatrix} \quad \Rightarrow \quad A^T = \begin{bmatrix} 1 & 4 \\ 2 & 5 \\ 3 & 6 \end{bmatrix}

Special Matrices Related to Transpose:

  • a) Symmetric Matrix: A square matrix where A=ATA = A^T (i.e., aij=ajia_{ij} = a_{ji})

S=[123256369]S = \begin{bmatrix} 1 & 2 & 3 \\ 2 & 5 & 6 \\ 3 & 6 & 9 \end{bmatrix}

  • b) Skew-Symmetric Matrix: A square matrix where A=ATA = -A^T (i.e., aij=ajia_{ij} = -a_{ji})

    • All diagonal elements must be zero.

K=[023204340]K = \begin{bmatrix} 0 & 2 & -3 \\ -2 & 0 & 4 \\ 3 & -4 & 0 \end{bmatrix}

  • Important Property: Any square matrix can be expressed as the sum of a symmetric and skew-symmetric matrix:

A=12(A+AT)+12(AAT)\boxed{A = \frac{1}{2}(A + A^T) + \frac{1}{2}(A - A^T)}

  • Symmetric part: 12(A+AT)\boxed{\frac{1}{2}(A + A^T)}
  • Skew-symmetric part: 12(AAT)\boxed{\frac{1}{2}(A - A^T)}

a) Minor of an element

  • Definition: The minor of an element aija_{ij} in a matrix AA, denoted by MijM_{ij}, is the determinant of the submatrix obtained by deleting the ii-th row and jj-th column from AA.

  • Notation: MijM_{ij} or Mij|M_{ij}|

  • For a 2×2 Matrix: A=[ab cd]A = \begin{bmatrix} a & b \ c & d \end{bmatrix} M11M_{11} = determinant after removing 1st row and 1st column = dd M12M_{12} = determinant after removing 1st row and 2nd column = cc M21M_{21} = determinant after removing 2nd row and 1st column = bb M22M_{22} = determinant after removing 2nd row and 2nd column = aa

  • For a 3×3 Matrix: A=[a11a12a13a21a22a23a31a32a33]A = \begin{bmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{bmatrix}

    M11=a22a23a32a33=a22a33a23a32M_{11} = \begin{vmatrix} a_{22} & a_{23} \\ a_{32} & a_{33} \end{vmatrix} = a_{22}a_{33} - a_{23}a_{32}

    M12=a21a23a31a33=a21a33a23a31M_{12} = \begin{vmatrix} a_{21} & a_{23} \\ a_{31} & a_{33} \end{vmatrix} = a_{21}a_{33} - a_{23}a_{31}

    M13=a21a22a31a32=a21a32a22a31M_{13} = \begin{vmatrix} a_{21} & a_{22} \\ a_{31} & a_{32} \end{vmatrix} = a_{21}a_{32} - a_{22}a_{31}

b) Cofactor of an Element

  • Definition: The cofactor of an element aija_{ij}, denoted by CijC_{ij} or AijA_{ij}, is given by:

Cij=(1)i+jMij\boxed{C_{ij} = (-1)^{i+j} \cdot M_{ij}}

  • where MijM_{ij} is the minor of element aija_{ij}.

  • Sign Pattern: The factor (1)i+j(-1)^{i+j} creates a checkerboard pattern of signs:

For any matrix: $$\begin{bmatrix}

  • & - & + & - & \cdots \
  • & + & - & + & \cdots \
  • & - & + & - & \cdots \
  • & + & - & + & \cdots \ \vdots & \vdots & \vdots & \vdots & \ddots \end{bmatrix}$$

For 2×2 Matrix: [++]\begin{bmatrix} + & - \\ - & + \end{bmatrix}

For 3×3 Matrix: [+++++]\begin{bmatrix} + & - & + \\ - & + & - \\ + & - & + \end{bmatrix} Properties of Minors and Cofactors

  • Property 1: The minor MijM_{ij} is always a determinant value (scalar).
  • Property 2: The cofactor CijC_{ij} differs from minor MijM_{ij} only by sign (1)i+j(-1)^{i+j}.
  • Property 3: When i+ji + j is even, Cij=MijC_{ij} = M_{ij} (positive sign)
  • Property 4: When i+ji + j is odd, Cij=MijC_{ij} = -M_{ij} (negative sign)
  • Property 5: For a diagonal matrix, all cofactors on the diagonal are products of the other diagonal elements.
  • Definition: Transpose of the cofactor matrix

  • Steps to Find Adjoint:

    1. Find the cofactor matrix C=[Cij]C = [C_{ij}]
    2. Take transpose: adj(A)=CT\text{adj}(A) = C^T
  • Properties:

    • Aadj(A)=adj(A)A=det(A)IA \cdot \text{adj}(A) = \text{adj}(A) \cdot A = \det(A) \cdot I
    • adj(AT)=(adj A)T\text{adj}(A^T) = (\text{adj }A)^T
    • adj(AB)=adj(B)adj(A)\text{adj}(AB) = \text{adj}(B) \cdot \text{adj}(A)
    • adj(kA)=kn1adj(A)\text{adj}(kA) = k^{n-1} \text{adj}(A) for n×nn \times n matrix
    • det(adj A)=[det(A)]n1\det(\text{adj }A) = [\det(A)]^{n-1} for n×nn \times n matrix
  • Example for 2×2 Matrix:

    If A=[abcd]\text{If } A = \begin{bmatrix} a & b \\ c & d \end{bmatrix}

    then adj(A)=[dbca]\boxed{\text{adj}(A) = \begin{bmatrix} d & -b \\ -c & a \end{bmatrix}}

  • Definition: A scalar value that can be computed from a square matrix

  • Notation: det(A)\det(A) or A|A|

  1. For 2×2 Matrix:

abcd=adbc\begin{vmatrix} a & b \\ c & d \end{vmatrix} = ad - bc

  1. For 3×3 Matrix:

abcdefghi=a(eifh)b(difg)+c(dheg)\begin{vmatrix} a & b & c \\ d & e & f \\ g & h & i \end{vmatrix} = a(ei - fh) - b(di - fg) + c(dh - eg)

  • Expansion by Cofactors:

det(A)=j=1naijCij(along any row or column)\det(A) = \sum_{j=1}^{n} a_{ij} C_{ij} \quad \text{(along any row or column)}

  • where Cij=(1)i+jMijC_{ij} = (-1)^{i+j} M_{ij} (cofactor)

  • and MijM_{ij} = minor (determinant of submatrix after removing ithi^{th} row and jthj^{th} column)

  • Properties of Determinants:

    • det(I)=1\det(I) = 1
    • det(O)=0\det(O) = 0
    • det(AT)=det(A)\det(A^T) = \det(A)
    • det(kA)=kndet(A)\det(kA) = k^n \det(A) for n×nn \times n matrix
    • det(AB)=det(A)det(B)\det(AB) = \det(A) \cdot \det(B)
    • If two rows/columns are identical, det(A)=0\det(A) = 0
    • If we interchange two rows/columns, determinant changes sign
    • Multiplying a row/column by kk multiplies determinant by kk
    • Adding a multiple of one row to another doesn’t change determinant
  • Singular Matrix: det(A)=0\det(A) = 0

  • Non-singular Matrix: det(A)0\det(A) \neq 0

  • Definition: A1A^{-1} is the inverse of AA if AA1=A1A=IA \cdot A^{-1} = A^{-1} \cdot A = I

  • Formula:

A1=1det(A)adj(A)\boxed{A^{-1} = \frac{1}{\det(A)} \cdot \text{adj}(A)}

  • Conditions for Inverse to Exist:

    • Matrix must be square
    • det(A)0\det(A) \neq 0 (non-singular)
  • Properties:

    • (A1)1=A(A^{-1})^{-1} = A
    • (AB)1=B1A1(AB)^{-1} = B^{-1} A^{-1}
    • (AT)1=(A1)T(A^T)^{-1} = (A^{-1})^T
    • det(A1)=1det(A)\det(A^{-1}) = \frac{1}{\det(A)}
    • (kA)1=1kA1(kA)^{-1} = \frac{1}{k} A^{-1}
  • Example:

A=[2314]A = \begin{bmatrix} 2 & 3 \\ 1 & 4 \end{bmatrix}

det(A)=2×43×1=83=5\det(A) = 2 \times 4 - 3 \times 1 = 8 - 3 = 5

adj(A)=[43 12]\text{adj}(A) = \begin{bmatrix} 4 & -3 \ -1 & 2 \end{bmatrix}

A1=15[4312]=[4/53/51/52/5]A^{-1} = \frac{1}{5}\begin{bmatrix} 4 & -3 \\ -1 & 2 \end{bmatrix} = \begin{bmatrix} 4/5 & -3/5 \\ -1/5 & 2/5 \end{bmatrix}

Used to solve systems of equations and find inverse:

  1. Row Interchange: RiRjR_i \leftrightarrow R_j (swap two rows)
  2. Row Scaling: RikRiR_i \to kR_i (multiply a row by non-zero scalar)
  3. Row Addition: RiRi+kRjR_i \to R_i + kR_j (add multiple of one row to another)

Elementary Matrices: Matrices obtained by performing one elementary operation on the identity matrix

  • Definition: The rank of a matrix is the maximum number of linearly independent rows (or columns)

  • Methods to Find Rank:

    1. Convert to row echelon form using elementary row operations
    2. Count the number of non-zero rows==
  • Properties:

    • rank(A)min(m,n)\text{rank}(A) \leq \min(m, n) for m×nm \times n matrix
    • rank(A)=rank(AT)\text{rank}(A) = \text{rank}(A^T)
    • rank(AB)min(rank(A),rank(B))\text{rank}(AB) \leq \min(\text{rank}(A), \text{rank}(B))
    • If AA is invertible n×nn \times n matrix, rank(A)=n\text{rank}(A) = n (full rank)
  • Row Echelon Form: ⭐

    • All non-zero rows are above zero rows
    • Leading entry of each non-zero row is to the right of the leading entry of the row above
  • Reduced Row Echelon Form:

    • Row echelon form
    • Leading entry in each non-zero row is 1
    • All entries above and below each leading 1 are zero
10. System of Linear Equations (❌ Not all important for GATE)
Section titled “10. System of Linear Equations (❌ Not all important for GATE)”
  • General Form:

{a11x1+a12x2++a1nxn=b1a21x1+a22x2++a2nxn=b2am1x1+am2x2++amnxn=bm\begin{cases} a_{11}x_1 + a_{12}x_2 + \cdots + a_{1n}x_n = b_1 \\ a_{21}x_1 + a_{22}x_2 + \cdots + a_{2n}x_n = b_2 \\ \vdots \\ a_{m1}x_1 + a_{m2}x_2 + \cdots + a_{mn}x_n = b_m \end{cases}

  • Matrix Form: AX=B\boxed{AX = B}
  • where AA is coefficient matrix, XX is variable matrix, BB is constant matrix

Methods of Solution

  • a) Matrix Inversion Method (when AA is square and non-singular): ⭐

X=A1B\boxed{X = A^{-1}B}

  • b) Cramer’s Rule (for square systems with det(A)0\det(A) \neq 0):

xi=det(Ai)det(A)x_i = \frac{\det(A_i)}{\det(A)} - where AiA_i is AA with ii-th column replaced by BB

  • c) Gaussian Elimination:

    • Convert augmented matrix [AB][A|B] to row echelon form
    • Use back substitution
  • d) Gauss-Jordan Elimination:

    • Convert to reduced row echelon form
    • Read off solutions directly

Consistency of System: ( ⭐GATE 2025 )

  • Augmented Matrix: [AB][A|B]
  • Using rank:
    • Consistent (has solution): rank(A)=rank([AB])\text{rank}(A) = \text{rank}([A|B])
      1. Unique solution: rank(A)=rank([AB])=n\text{rank}(A) = \text{rank}([A|B]) = n (number of variables)
      2. Infinite solutions: rank(A)=rank([AB])<n\text{rank}(A) = \text{rank}([A|B]) < n
    • Inconsistent (no solution): rank(A)<rank([AB])\text{rank}(A) < \text{rank}([A|B])
11. Special Products and Matrices ( ❌Not important for GATE)
Section titled “11. Special Products and Matrices ( ❌Not important for GATE)”

a) Idempotent Matrix: A2=AA^2 = A

  • Example: Identity matrix

b) Nilpotent Matrix: Ak=OA^k = O for some positive integer kk

  • Smallest such kk is called index of nilpotency

c) Involutory Matrix: A2=IA^2 = I

  • AA is its own inverse: A1=AA^{-1} = A

d) Orthogonal Matrix: AAT=ATA=IAA^T = A^T A = I

  • A1=ATA^{-1} = A^T
  • det(A)=±1\det(A) = \pm 1
  • Columns (and rows) form orthonormal sets

e) Hermitian Matrix (for complex matrices): A=AˉTA = \bar{A}^T (conjugate transpose)

f) Unitary Matrix: AAˉT=IA\bar{A}^T = I

Definition: Sum of diagonal elements of a square matrix

tr(A)=i=1naii\boxed{\text{tr}(A) = \sum_{i=1}^{n} a_{ii}}

Trace–Eigenvalue Relation: (⭐GATE 2023)

  • Summation of the eigenvalue is equal to the trace of the matrix tr(A)=i=1nλi\boxed{\text{tr}(A)=\sum_{i=1}^{n}\lambda_i}

  • Holds for any square matrix (diagonalizable or not)

  • Follows from characteristic polynomial

  • Works over real or complex field

  • Frequently used to find unknown eigenvalues using trace shortcut

Properties:

  • tr(A+B)=tr(A)+tr(B)\text{tr}(A + B) = \text{tr}(A) + \text{tr}(B)
  • tr(kA)=ktr(A)\text{tr}(kA) = k \cdot \text{tr}(A)
  • tr(AB)=tr(BA)\text{tr}(AB) = \text{tr}(BA)
  • tr(AT)=tr(A)\text{tr}(A^T) = \text{tr}(A)
  • tr(ABC)=tr(CAB)=tr(BCA)\text{tr}(ABC) = \text{tr}(CAB) = \text{tr}(BCA) (cyclic property)
13. Applications of Matrices ( ❌Not important for GATE)
Section titled “13. Applications of Matrices ( ❌Not important for GATE)”
  1. Solving Linear Equations: Systems of equations in science, engineering, economics
  2. Computer Graphics: Transformations (rotation, scaling, translation)
  3. Cryptography: Encoding and decoding messages
  4. Economics: Input-output models, Leontief models
  5. Network Analysis: Graph theory, connectivity
  6. Quantum Mechanics: State vectors, operators
  7. Statistics: Correlation matrices, covariance matrices
  8. Machine Learning: Data representation, neural networks

Determinant (2×2):

A=abcd=adbc\boxed{ |A| = \begin{vmatrix} a & b \\ c & d \end{vmatrix} = ad - bc}

Inverse:

A1=1Aadj(A)\boxed{A^{-1} = \frac{1}{|A|} \cdot \text{adj}(A)}

Adjoint (2×2):

adj[abcd]=[dbca]\boxed{\text{adj}\begin{bmatrix} a & b \\ c & d \end{bmatrix} = \begin{bmatrix} d & -b \\ -c & a \end{bmatrix}}

Product Rules:

(AB)T=BTAT\boxed{(AB)^T = B^T A^T}

(AB)1=B1A1\boxed{(AB)^{-1} = B^{-1} A^{-1}}

Rank Property:

rank(A)=rank(AT)\boxed{\text{rank}(A) = \text{rank}(A^T)}

Solution of AX=BAX = B:

X=A1B(when A is invertible)\boxed{X = A^{-1}B \quad \text{(when } A \text{ is invertible)}}


Eigenvector: A non-zero vector v\vec{v} is called an eigenvector of a square matrix AA if there exists a scalar λ\lambda such that:

Av=λv\boxed{A\vec{v} = \lambda\vec{v}}

Eigenvalue: The scalar λ\lambda is called the eigenvalue corresponding to the eigenvector v\vec{v}.

  • v\vec{v}: Eigen Vector
  • λ\lambda: Eigen Value

Geometric Interpretation: When matrix AA acts on eigenvector v\vec{v}, the result is a vector in the same direction (or opposite direction) as v\vec{v}, scaled by factor λ\lambda.

  • To find eigenvalues, we solve:

Av=λvA\vec{v} = \lambda\vec{v} Avλv=0A\vec{v} - \lambda\vec{v} = \vec{0} (AλI)v=0(A - \lambda I)\vec{v} = \vec{0}

  • For non-trivial solution (where v0\vec{v} \neq \vec{0}):

det(AλI)=0\boxed{\det(A - \lambda I) = 0}

  • This is called the Characteristic Equation of matrix AA.

  • Characteristic Polynomial: p(λ)=det(AλI)\boxed{p(\lambda) = \det(A - \lambda I)}

3. Finding Eigenvalues and Eigenvectors ⭐⭐⭐
Section titled “3. Finding Eigenvalues and Eigenvectors ⭐⭐⭐”

Step-by-Step Process:

  1. Step 1: Find eigenvalues by solving det(AλI)=0\det(A - \lambda I) = 0
  2. Step 2: For each eigenvalue λi\lambda_i, find eigenvectors by solving (AλiI)v=0(A - \lambda_i I)\vec{v} = \vec{0}

Example: Find eigenvalues and eigenvectors of A=[4213]A = \begin{bmatrix} 4 & 2 \\ 1 & 3 \end{bmatrix}

Solution:

  1. Step 1: Find eigenvalues

    AλI=[4λ213λ]A - \lambda I = \begin{bmatrix} 4-\lambda & 2 \\ 1 & 3-\lambda \end{bmatrix}

    det(AλI)=(4λ)(3λ)2=λ27λ+10=0\det(A - \lambda I) = (4-\lambda)(3-\lambda) - 2 = \lambda^2 - 7\lambda + 10 = 0

    (λ5)(λ2)=0(\lambda - 5)(\lambda - 2) = 0 \quad, Eigenvalues: λ1=5,λ2=2\lambda_1 = 5, \lambda_2 = 2

  2. Step 2: Find eigenvectors

    • For λ1=5\lambda_1 = 5: (A5I)v=[1212][v1v2]=[00](A - 5I)\vec{v} = \begin{bmatrix} -1 & 2 \\ 1 & -2 \end{bmatrix} \begin{bmatrix} v_1 \\ v_2 \end{bmatrix} = \begin{bmatrix} 0 \\ 0 \end{bmatrix} \quad, Solution: v1=[21]\vec{v}_1 = \begin{bmatrix} 2 \\ 1 \end{bmatrix} (or any scalar multiple)

    • For λ2=2\lambda_2 = 2: (A2I)v=[22 11][v1 v2]=[0 0](A - 2I)\vec{v} = \begin{bmatrix} 2 & 2 \ 1 & 1 \end{bmatrix}\begin{bmatrix} v_1 \ v_2 \end{bmatrix} = \begin{bmatrix} 0 \ 0 \end{bmatrix} \quad, Solution: v2=[1 1]\vec{v}_2 = \begin{bmatrix} 1 \ -1 \end{bmatrix} (or any scalar multiple)

  1. Property 1: Sum of eigenvalues = Trace of matrix

i=1nλi=tr(A)=i=1naii\sum_{i=1}^{n} \lambda_i = \text{tr}(A) = \sum_{i=1}^{n} a_{ii}

  1. Property 2: Product of eigenvalues = Determinant of matrix

i=1nλi=det(A)\prod_{i=1}^{n} \lambda_i = \det(A)

  1. Property 3: If λ\lambda is an eigenvalue of AA, then:
  • λk\lambda^k is an eigenvalue of AkA^k (⭐GATE 2023)
  • λ+c\lambda + c is an eigenvalue of A+cIA + cI
  • kλk\lambda is an eigenvalue of kAkA
  • 1λ\frac{1}{\lambda} is an eigenvalue of A1A^{-1} (if AA is invertible and λ0\lambda \neq 0)
  1. Property 4: Eigenvalues of AA and ATA^T are the same

  2. Property 5: For triangular (upper or lower) and diagonal matrices, eigenvalues are the diagonal elements

  3. Property 6: If AA is singular, then λ=0\lambda = 0 is an eigenvalue

  4. Property 7: Eigenvalues of a real symmetric matrix are always real

  5. Property 8: Eigenvalues of an orthogonal matrix have absolute value 1

  6. Property 9: Eigenvalues of a Hermitian matrix are real

  7. Property 10: Eigenvalues of a skew-Hermitian matrix are purely imaginary or zero

  1. Property 1: Eigenvectors corresponding to distinct eigenvalues are linearly independent

  2. Property 2: If v\vec{v} is an eigenvector, then kvk\vec{v} (for any non-zero scalar kk) is also an eigenvector with the same eigenvalue

  3. Property 3: For a real symmetric matrix, eigenvectors corresponding to distinct eigenvalues are orthogonal

  4. Property 4: The eigenvectors corresponding to an eigenvalue with algebraic multiplicity mm span a subspace of dimension m\leq m

6. Algebraic and Geometric Multiplicity ( ❌Not important for GATE)
Section titled “6. Algebraic and Geometric Multiplicity ( ❌Not important for GATE)”

Algebraic Multiplicity: The number of times an eigenvalue appears as a root of the characteristic polynomial

Geometric Multiplicity: The dimension of the eigenspace (number of linearly independent eigenvectors) corresponding to an eigenvalue

Geometric MultiplicityAlgebraic Multiplicity\text{Geometric Multiplicity} \leq \text{Algebraic Multiplicity}

Example:

A=[500050003]A = \begin{bmatrix} 5 & 0 & 0 \\ 0 & 5 & 0 \\ 0 & 0 & 3 \end{bmatrix}

  • Characteristic equation: (5λ)2(3λ)=0(5-\lambda)^2(3-\lambda) = 0
  • Eigenvalues: λ1=5\lambda_1 = 5 (algebraic multiplicity 2), λ2=3\lambda_2 = 3 (algebraic multiplicity 1)
  • For λ1=5\lambda_1 = 5: geometric multiplicity = 2 (two independent eigenvectors)

A matrix AA is diagonalizable if there exists an invertible matrix PP such that:

P1AP=D\boxed{P^{-1}AP = D}

where DD is a diagonal matrix containing eigenvalues of AA, and columns of PP are the corresponding eigenvectors.

D=[λ1000λ2000λn]D = \begin{bmatrix} \lambda_1 & 0 & \cdots & 0 \\ 0 & \lambda_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & \lambda_n \end{bmatrix}

Condition for Diagonalization: An n×nn \times n matrix is diagonalizable if and only if it has nn linearly independent eigenvectors.

Guaranteed Diagonalizable:

  • Matrices with nn distinct eigenvalues
  • Real symmetric matrices
  • Hermitian matrices
  • Orthogonal matrices
  • Unitary matrices

Applications of Diagonalization:

  1. Computing Powers: Ak=PDkP1A^k = PD^kP^{-1}
  2. Matrix Functions: f(A)=Pf(D)P1f(A) = Pf(D)P^{-1}
  3. Solving Differential Equations
8. Similarity of Matrices ( ❌Not important for GATE)
Section titled “8. Similarity of Matrices ( ❌Not important for GATE)”

Two matrices AA and BB are similar if there exists an invertible matrix PP such that:

B=P1AP\boxed{B = P^{-1}AP}

Properties of Similar Matrices:

  • Same eigenvalues
  • Same determinant
  • Same trace
  • Same rank
  • Same characteristic polynomial
9. Cayley-Hamilton Theorem ( ❌Not important for GATE)
Section titled “9. Cayley-Hamilton Theorem ( ❌Not important for GATE)”

Theorem: Every square matrix satisfies its own characteristic equation.

If p(λ)=det(AλI)p(\lambda) = \det(A - \lambda I) is the characteristic polynomial, then:

p(A)=Op(A) = O

Example: If A=[1234]A = \begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix}

  • Characteristic equation: λ25λ2=0\lambda^2 - 5\lambda - 2 = 0

  • Then: A25A2I=OA^2 - 5A - 2I = O

  • Applications:

    • Finding matrix inverse
    • Computing higher powers of matrices
    • Simplifying matrix polynomials

Finding Inverse using Cayley-Hamilton:

  • If characteristic equation is λn+cn1λn1++c1λ+c0=0\lambda^n + c_{n-1}\lambda^{n-1} + \cdots + c_1\lambda + c_0 = 0

  • Then: An+cn1An1++c1A+c0I=OA^n + c_{n-1}A^{n-1} + \cdots + c_1A + c_0I = O

  • If c00c_0 \neq 0 (i.e., det(A)0\det(A) \neq 0):

A1=1c0(An1+cn1An2++c1I)A^{-1} = -\frac{1}{c_0}(A^{n-1} + c_{n-1}A^{n-2} + \cdots + c_1I)

10. Special Eigenvalue Problems ( ❌Not important for GATE)
Section titled “10. Special Eigenvalue Problems ( ❌Not important for GATE)”

Symmetric Matrices:

  • All eigenvalues are real
  • Eigenvectors corresponding to distinct eigenvalues are orthogonal
  • Always diagonalizable

Orthogonal Matrices (ATA=IA^TA = I):

  • All eigenvalues have absolute value 1
  • λ=1|\lambda| = 1 or λ=eiθ\lambda = e^{i\theta}

Positive Definite Matrices:

  • All eigenvalues are positive
  • All leading principal minors are positive

Idempotent Matrices (A2=AA^2 = A):

  • Eigenvalues are either 0 or 1

Nilpotent Matrices (Ak=OA^k = O):

  • All eigenvalues are 0

Involutory Matrices (A2=IA^2 = I):

  • Eigenvalues are either 1 or -1
11. Power Method (Finding Dominant Eigenvalue) ( ❌Not important for GATE)
Section titled “11. Power Method (Finding Dominant Eigenvalue) ( ❌Not important for GATE)”

The power method is an iterative technique to find the largest (dominant) eigenvalue and its corresponding eigenvector.

  • Algorithm:

    1. Start with initial guess x0\vec{x}_0 (non-zero)
    2. Iterate: xk+1=AxkAxk\vec{x}_{k+1} = \frac{A\vec{x}_k}{|A\vec{x}_k|}
    3. Eigenvalue: λxkTAxkxkTxk\lambda \approx \frac{\vec{x}_k^T A \vec{x}_k}{\vec{x}_k^T \vec{x}_k}
  • Converges to: Dominant eigenvalue (eigenvalue with largest absolute value)


System of Linear Equations (Advanced Concepts)

Section titled “System of Linear Equations (Advanced Concepts)”

System of equations:

{a11x1+a12x2++a1nxn=b1a21x1+a22x2++a2nxn=b2am1x1+am2x2++amnxn=bm\begin{cases} a_{11}x_1 + a_{12}x_2 + \cdots + a_{1n}x_n = b_1 \\ a_{21}x_1 + a_{22}x_2 + \cdots + a_{2n}x_n = b_2 \\ \vdots \\ a_{m1}x_1 + a_{m2}x_2 + \cdots + a_{mn}x_n = b_m \end{cases}

Matrix form: AX=BAX = B

where:

A=[a11a12a1na21a22a2nam1am2amn],X=[x1x2xn],B=[b1b2bm]A = \begin{bmatrix} a_{11} & a_{12} & \cdots & a_{1n} \\ a_{21} & a_{22} & \cdots & a_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ a_{m1} & a_{m2} & \cdots & a_{mn} \end{bmatrix}, \quad X = \begin{bmatrix} x_1 \\ x_2 \\ \vdots \\ x_n \end{bmatrix}, \quad B = \begin{bmatrix} b_1 \\ b_2 \\ \vdots \\ b_m \end{bmatrix}

Augmented Matrix: [AB][A | B] ⭐⭐

  1. Homogeneous System: AX=0AX = 0 (where B=0B = 0) ⭐

    • Always has at least the trivial solution X=0X = 0
    • Non-trivial solutions exist if and only if det(A)=0\det(A) = 0 (for square matrices) ⭐
  2. Non-homogeneous System: AX=BAX = B (where B0B \neq 0)

For Square Systems (m=nm = n):

Determinant ConditionDeterminantSolution Type
det(A)0\det(A) \neq 0Non-zeroUnique solution: X=A1BX = A^{-1}B
det(A)=0\det(A) = 0 and rank(A)=rank([AB])\text{rank}(A) = \text{rank}([A\|B])ZeroInfinite solutions
det(A)=0\det(A) = 0 and rank(A)<rank([AB])\text{rank}(A) < \text{rank}([A\|B])ZeroNo solution (inconsistent)

For General Systems (m×nm \times n):

Let r=rank(A)r = \text{rank}(A) and r=rank([AB])r' = \text{rank}([A|B])

Rank ConditionSolution Type
r<rr < r'Inconsistent (no solution)
r=r=nr = r' = nUnique solution
r=r<nr = r' < nInfinite solutions (with nrn - r free parameters)
4. Rouché-Capelli Theorem ( ❌Not important for GATE)
Section titled “4. Rouché-Capelli Theorem ( ❌Not important for GATE)”

A system AX=BAX = B is consistent if and only if:

rank(A)=rank([AB])\boxed{\text{rank}(A) = \text{rank}([A|B])}

Number of free variables = nrank(A)n - \text{rank}(A)

System: AX=0AX = 0

Properties:

  • Always has trivial solution X=0X = 0
  • Has non-trivial solution if and only if det(A)=0\det(A) = 0 (for square AA)
  • For m×nm \times n system with m<nm < n: always has infinite solutions
  • Solution space forms a vector space (null space or kernel of AA)

Dimension of solution space: nrank(A)n - \text{rank}(A) (called nullity of AA)

Rank-Nullity Theorem: For m×nm \times n matrix AA:

rank(A)+nullity(A)=n\text{rank}(A) + \text{nullity}(A) = n

6. Solution Methods ( ❌ Not all important for GATE)
Section titled “6. Solution Methods ( ❌ Not all important for GATE)”

A. Cramer’s Rule

  • For square system with det(A)0\det(A) \neq 0:

xi=det(Ai)det(A)x_i = \frac{\det(A_i)}{\det(A)}

  • where AiA_i is the matrix obtained by replacing the ii-th column of AA with BB.

  • Example: Solve {2x+3y=8 x+4y=9\begin{cases} 2x + 3y = 8 \ x + 4y = 9 \end{cases}

    A=[23 14],B=[8 9]A = \begin{bmatrix} 2 & 3 \ 1 & 4 \end{bmatrix}, \quad B = \begin{bmatrix} 8 \ 9 \end{bmatrix}

    det(A)=83=5\det(A) = 8 - 3 = 5

    A1=[83 94],det(A1)=3227=5A_1 = \begin{bmatrix} 8 & 3 \ 9 & 4 \end{bmatrix}, \quad \det(A_1) = 32 - 27 = 5

    x=55=1x = \frac{5}{5} = 1

    A2=[28 19],det(A2)=188=10A_2 = \begin{bmatrix} 2 & 8 \ 1 & 9 \end{bmatrix}, \quad \det(A_2) = 18 - 8 = 10

    y=105=2y = \frac{10}{5} = 2

  • Limitation: Computationally expensive for large systems

B. Matrix Inversion Method

  • For square system with det(A)0\det(A) \neq 0:

X=A1B \boxed{X = A^{-1}B}

  • Steps:
    1. Find A1A^{-1} using A1=1det(A)adj(A)A^{-1} = \frac{1}{\det(A)} \text{adj}(A)
    2. Compute X=A1BX = A^{-1}B

C. Gaussian Elimination

  • Process:

    1. Form augmented matrix [AB][A|B]
    2. Use elementary row operations to convert to row echelon form
    3. Use back substitution to find solutions
  • Elementary Row Operations:

    • RiRjR_i \leftrightarrow R_j (interchange rows)
    • RikRiR_i \to kR_i (multiply row by non-zero constant)
    • RiRi+kRjR_i \to R_i + kR_j (add multiple of one row to another)
  • Row Echelon Form:

    • All zero rows are at the bottom
    • Leading coefficient (pivot) of each non-zero row is to the right of the pivot above it
  • Example: {x+2y+z=9 2xy+3z=8 3x+y+2z=13\begin{cases} x + 2y + z = 9 \ 2x - y + 3z = 8 \ 3x + y + 2z = 13 \end{cases}

    Augmented matrix: [1219213831213]\left[\begin{array}{ccc|c} 1 & 2 & 1 & 9 \\ 2 & -1 & 3 & 8 \\ 3 & 1 & 2 & 13 \end{array}\right]

    After row operations: [1219051100024]\left[\begin{array}{ccc|c} 1 & 2 & 1 & 9 \\ 0 & -5 & 1 & -10 \\ 0 & 0 & -2 & -4 \end{array}\right]

    Back substitution: z=2,y=2,x=1z = 2, y = 2, x = 1

D. Gauss-Jordan Elimination

  • Process:

    1. Form augmented matrix [AB][A|B]
    2. Convert to reduced row echelon form (RREF)
    3. Read solutions directly
  • Reduced Row Echelon Form:

    • Row echelon form
    • Leading coefficient (pivot) in each row is 1
    • Each pivot is the only non-zero entry in its column
  • Example: Same system as above

    After Gauss-Jordan: [100101020012]\left[\begin{array}{ccc|c} 1 & 0 & 0 & 1 \\ 0 & 1 & 0 & 2 \\ 0 & 0 & 1 & 2 \end{array}\right]

    Solution: x=1,y=2,z=2x = 1, y = 2, z = 2

E. LU Decomposition

  • Decompose AA into product of lower triangular LL and upper triangular UU:

A=LUA = LU

  • Process:

    1. Find LL and UU such that A=LUA = LU
    2. Solve LY=BLY = B (forward substitution)
    3. Solve UX=YUX = Y (backward substitution)
  • Advantage: Efficient for solving multiple systems with same AA but different BB

7. Particular and General Solutions ( ❌ Not important for GATE)
Section titled “7. Particular and General Solutions ( ❌ Not important for GATE)”

For non-homogeneous system AX=BAX = B:

General Solution = Particular Solution + Homogeneous Solution

X=Xp+XhX = X_p + X_h

where:

  • XpX_p is any particular solution of AX=BAX = B
  • XhX_h is the general solution of AX=0AX = 0

Example: If AX=BAX = B has particular solution Xp=[120]X_p = \begin{bmatrix} 1 \\ 2 \\ 0 \end{bmatrix}

and homogeneous solution Xh=t[111]X_h = t\begin{bmatrix} 1 \\ -1 \\ 1 \end{bmatrix} (where tt is any real number)

Then general solution: X=[120]+t[111]=[1+t2tt]X = \begin{bmatrix} 1 \\ 2 \\ 0 \end{bmatrix} + t\begin{bmatrix} 1 \\ -1 \\ 1 \end{bmatrix} = \begin{bmatrix} 1+t \\ 2-t \\ t \end{bmatrix}

8. Vector Space Interpretation ( ❌ Not important for GATE)
Section titled “8. Vector Space Interpretation ( ❌ Not important for GATE)”

Column Space (Range/Image of AA): C(A)=Ax:xRnC(A) = {A\vec{x} : \vec{x} \in \mathbb{R}^n}

  • Dimension = rank(A)\text{rank}(A)

Null Space (Kernel of AA): N(A)=x:Ax=0N(A) = {\vec{x} : A\vec{x} = \vec{0}}

  • Dimension = nullity(A)=nrank(A)\text{nullity}(A) = n - \text{rank}(A)

Row Space: Space spanned by rows of AA

  • Dimension = rank(A)\text{rank}(A)

Left Null Space: Null space of ATA^T

Fundamental Theorem:

rank(A)+nullity(A)=n\text{rank}(A) + \text{nullity}(A) = n

dim(C(A))+dim(N(AT))=m\dim(C(A)) + \dim(N(A^T)) = m

9. Least Squares Solution ( ❌ Not important for GATE)
Section titled “9. Least Squares Solution ( ❌ Not important for GATE)”

For overdetermined systems (m>nm > n, more equations than unknowns) with no exact solution:

Normal Equation:

ATAX=ATBA^TAX = A^TB

Least Squares Solution:

X=(ATA)1ATBX = (A^TA)^{-1}A^TB

This minimizes AXB2|AX - B|^2

Application: Best fit line, curve fitting, data regression

10. Condition Number and Stability ( ❌ Not important for GATE)
Section titled “10. Condition Number and Stability ( ❌ Not important for GATE)”

Condition Number: Measures sensitivity of solution to changes in AA or BB

κ(A)=AA1\kappa(A) = |A| \cdot |A^{-1}|

For eigenvalue-based definition:

κ(A)=λmaxλmin\kappa(A) = \frac{\lambda_{\max}}{\lambda_{\min}}

Well-conditioned: κ(A)\kappa(A) is small (close to 1)

Ill-conditioned: κ(A)\kappa(A) is large

Singular matrix: κ(A)=\kappa(A) = \infty

11. Important GATE Concepts Summary ( ❌ Not important for GATE)
Section titled “11. Important GATE Concepts Summary ( ❌ Not important for GATE)”

For Homogeneous Systems (AX=0AX = 0):

  • Always consistent (at least trivial solution)
  • Non-trivial solution exists     \iff det(A)=0\det(A) = 0 (for square AA)
  • Number of free variables = nrank(A)n - \text{rank}(A)

For Non-homogeneous Systems (AX=BAX = B):

  • Consistent     \iff rank(A)=rank([AB])\text{rank}(A) = \text{rank}([A|B])
  • Unique solution     \iff rank(A)=rank([AB])=n\text{rank}(A) = \text{rank}([A|B]) = n
  • Infinite solutions     \iff rank(A)=rank([AB])<n\text{rank}(A) = \text{rank}([A|B]) < n

Key Formulas:

  • Cramer’s Rule: xi=det(Ai)det(A)x_i = \frac{\det(A_i)}{\det(A)}
  • Matrix Inversion: X=A1BX = A^{-1}B
  • Rank-Nullity: rank(A)+nullity(A)=n\text{rank}(A) + \text{nullity}(A) = n

Solution Methods Priority (for GATE):

  1. Cramer’s Rule (for 2×2 or 3×3 with unique solution)
  2. Matrix Inversion (for small square non-singular systems)
  3. Gaussian Elimination (for general systems)
  4. Rank method (for consistency check)