Matrices are fundamental tools in biomedical engineering for organizing and manipulating data such as physiological measurements, imaging information, and system parameters. Applying matrices allows engineers to model complex biological systems, analyze signals, and design control systems for medical devices. This chapter introduces matrices, their types, notation, and basic operations, all contextualized for biomedical engineering applications in Kenya.

2.1 Introduction to matrices

Matrices are rectangular arrays of numbers arranged in rows and columns, used extensively in biomedical engineering for representing data sets and transformations. Understanding matrix types and operations enables efficient handling of multi-dimensional data from sources like ECG signals, MRI images, and patient monitoring systems.

2.1.1 Definition of a matrix

A matrix is a rectangular arrangement of numbers or elements in rows and columns enclosed by brackets. The size of a matrix is described by its number of rows \(m\) and columns \(n\), denoted as an \(m \times n\) matrix. Using matrices, biomedical engineers represent data such as sensor outputs or system coefficients compactly.

A matrix \(A\) with \(m\) rows and \(n\) columns is written as:

$$ 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} $$

where \(a_{ij}\) is the element in the \(i^{th}\) row and \(j^{th}\) column.

Worked Examples

Example 1: Represent a 2-channel ECG signal data for 3 time points as a matrix.

Given: ECG readings for 2 channels at 3 time points
Channel 1: 0.5 mV, 0.7 mV, 0.6 mV
Channel 2: 0.3 mV, 0.4 mV, 0.5 mV

Matrix \(E\):

$$ E = \begin{bmatrix} 0.5 & 0.7 & 0.6 \\ 0.3 & 0.4 & 0.5 \end{bmatrix} $$

Example 2: Write a \(3 \times 2\) matrix representing oxygen saturation levels (%) for 3 patients measured twice.

Given: Patient 1: 95, 96; Patient 2: 90, 92; Patient 3: 88, 89

Matrix \(O\):

$$ O = \begin{bmatrix} 95 & 96 \\ 90 & 92 \\ 88 & 89 \end{bmatrix} $$

Example 3: Identify the size of matrix

$$ M = \begin{bmatrix} 1 & 2 & 3 & 4 \\ 5 & 6 & 7 & 8 \end{bmatrix} $$

Matrix \(M\) has 2 rows and 4 columns, size \(2 \times 4\).

2.1.2 Different types of matrices: row, column, square, rectangular, diagonal, identity, zero

Matrices vary by their dimensions and element arrangements. Biomedical engineers must recognize matrix types to select appropriate operations for data processing and system modeling.

• Row matrix: Has only one row (\(1 \times n\)). Example: Sensor readings at one time instant.
• Column matrix: Has only one column (\(m \times 1\)). Example: Patient vital signs vector.
• Square matrix: Number of rows equals columns (\(n \times n\)). Used in system state representations.
• Rectangular matrix: Rows and columns differ (\(m \neq n\)). Common in data sets.
• Diagonal matrix: Square matrix with non-zero elements only on the main diagonal; zeros elsewhere. Represents scaling in systems.
• Identity matrix: Square matrix with ones on the main diagonal and zeros elsewhere. Acts as multiplicative identity.
• Zero matrix: All elements zero. Represents null data or transformations.

Worked Examples

Example 1: Identify matrix types

Given:

$$ A = \begin{bmatrix} 3 & 5 & 7 \end{bmatrix} $$

\(A\) is a row matrix (\(1 \times 3\)).

Example 2: Identify matrix types

Given:

$$ B = \begin{bmatrix} 2 \\ 4 \\ 6 \end{bmatrix} $$

\(B\) is a column matrix (\(3 \times 1\)).

Example 3: Check if matrix is diagonal or identity

Given:

$$ C = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix} $$

\(C\) is an identity matrix of size \(3 \times 3\).

Example 4: Determine if matrix is square or rectangular

Given:

$$ D = \begin{bmatrix} 7 & 8 \\ 9 & 10 \\ 11 & 12 \end{bmatrix} $$

Matrix \(D\) is rectangular \(3 \times 2\).

Example 5: Identify zero matrix

Given:

$$ Z = \begin{bmatrix} 0 & 0 \\ 0 & 0 \end{bmatrix} $$

Matrix \(Z\) is a zero matrix \(2 \times 2\).

2.1.3 Notation and elements of a matrix

Matrix elements are denoted by \(a_{ij}\), where \(i\) indicates the row and \(j\) the column. This notation is critical for referencing entries in biomedical data matrices such as pixel intensities in medical images or coefficients in system models.

The general matrix \(A\) is expressed as:

$$ A = [a_{ij}]_{m \times n} $$

where \(i = 1, 2, \ldots, m\) and \(j = 1, 2, \ldots, n\).

Worked Examples

Example 1: Write the element in the 2nd row, 3rd column of

$$ A = \begin{bmatrix} 1 & 2 & 3 \\ 4 & 5 & 6 \end{bmatrix} $$

Element \(a_{23} = 6\).

Example 2: Express the element \(a_{31}\) of matrix

$$ B = \begin{bmatrix} 7 & 8 \\ 9 & 10 \\ 11 & 12 \end{bmatrix} $$

Element \(a_{31} = 11\).

Example 3: Given matrix

$$ C = \begin{bmatrix} 2 & 4 & 6 \\ 8 & 10 & 12 \\ 14 & 16 & 18 \end{bmatrix} $$

Find the element \(a_{22}\):

\(a_{22} = 10\).

Example 4: Given matrix \(D = [d_{ij}]_{3 \times 3}\), with

$$ D = \begin{bmatrix} d_{11} & d_{12} & d_{13} \\ d_{21} & d_{22} & d_{23} \\ d_{31} & d_{32} & d_{33} \end{bmatrix} $$

If \(d_{ij} = i + j\), find \(d_{23}\):

$$ d_{23} = 2 + 3 = 5 $$

Example 5: For matrix \(E = [e_{ij}]_{2 \times 4}\), where \(e_{ij} = 3i, j\), find \(e_{14}\):

$$ e_{14} = 3(1) - 4 = 3 - 4 = -1 $$

2.1.4 Basic operations: addition, subtraction, scalar multiplication

Matrix operations allow biomedical engineers to combine and scale data sets, such as overlaying sensor readings or adjusting signal amplitudes.

• Addition: Two matrices of the same size are added by summing corresponding elements.
• Subtraction: Similar to addition, subtract corresponding elements.
• Scalar multiplication: Multiply every element in a matrix by a scalar (real number).

Conditions: For addition and subtraction, matrices must be the same size. Scalar multiplication applies to any matrix.

Formulas:

• Addition: \((A + B)_{ij} = a_{ij} + b_{ij}\)
• Subtraction: \((A - B)_{ij} = a_{ij} - b_{ij}\)
• Scalar multiplication: \((kA)_{ij} = k \times a_{ij}\), where \(k\) is a scalar.

Worked Examples

Example 1: Add matrices

$$ A = \begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix}, \quad B = \begin{bmatrix} 5 & 6 \\ 7 & 8 \end{bmatrix} $$

$$ A + B = \begin{bmatrix} 1+5 & 2+6 \\ 3+7 & 4+8 \end{bmatrix} = \begin{bmatrix} 6 & 8 \\ 10 & 12 \end{bmatrix} $$

Example 2: Subtract matrices

$$ C = \begin{bmatrix} 9 & 8 \\ 7 & 6 \end{bmatrix}, \quad D = \begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix} $$

$$ C - D = \begin{bmatrix} 9-1 & 8-2 \\ 7-3 & 6-4 \end{bmatrix} = \begin{bmatrix} 8 & 6 \\ 4 & 2 \end{bmatrix} $$

Example 3: Multiply matrix by scalar

$$ k = 3, \quad E = \begin{bmatrix} 1 & 0 \\ -1 & 2 \end{bmatrix} $$

$$ kE = \begin{bmatrix} 3 \times 1 & 3 \times 0 \\ 3 \times (-1) & 3 \times 2 \end{bmatrix} = \begin{bmatrix} 3 & 0 \\ -3 & 6 \end{bmatrix} $$

Example 4: Add two \(3 \times 2\) matrices

$$ F = \begin{bmatrix} 2 & 4 \\ 6 & 8 \\ 10 & 12 \end{bmatrix}, \quad G = \begin{bmatrix} 1 & 3 \\ 5 & 7 \\ 9 & 11 \end{bmatrix} $$

$$ F + G = \begin{bmatrix} 2+1 & 4+3 \\ 6+5 & 8+7 \\ 10+9 & 12+11 \end{bmatrix} = \begin{bmatrix} 3 & 7 \\ 11 & 15 \\ 19 & 23 \end{bmatrix} $$

Example 5: Multiply matrix by scalar \(k = -2\)

$$ H = \begin{bmatrix} 3 & -1 \\ 0 & 4 \end{bmatrix} $$

$$ kH = \begin{bmatrix} -2 \times 3 & -2 \times (-1) \\ -2 \times 0 & -2 \times 4 \end{bmatrix} = \begin{bmatrix} -6 & 2 \\ 0 & -8 \end{bmatrix} $$

2.1.5 Special Matrices

Certain matrices have properties that simplify calculations in biomedical engineering tasks such as image processing and system control.

• Symmetric matrix: Equal to its transpose, \(A = A^T\). Used in covariance matrices in signal analysis.
• Skew-symmetric matrix: \(A = -A^T\). Applied in rotational transformations.
• Orthogonal matrix: \(A^T A = I\). Represents rotations preserving length.
• Singular matrix: Non-invertible matrix with zero determinant. Important in system solvability.
• Diagonalizable matrix: Can be represented as \(A = PDP^{-1}\), where \(D\) is diagonal. Used in system eigenvalue analysis.

Worked Examples

Example 1: Check if matrix is symmetric

$$ A = \begin{bmatrix} 2 & 3 \\ 3 & 4 \end{bmatrix} $$

Transpose \(A^T = \begin{bmatrix} 2 & 3 \\ 3 & 4 \end{bmatrix}\), which equals \(A\), so \(A\) is symmetric.

Example 2: Check if matrix is skew-symmetric

$$ B = \begin{bmatrix} 0 & 5 \\ -5 & 0 \end{bmatrix} $$

Transpose \(B^T = \begin{bmatrix} 0 & -5 \\ 5 & 0 \end{bmatrix} = -B\), so \(B\) is skew-symmetric.

Example 3: Verify orthogonality

$$ C = \begin{bmatrix} 0 & 1 \\ -1 & 0 \end{bmatrix} $$

Calculate \(C^T C\):

$$ C^T = \begin{bmatrix} 0 & -1 \\ 1 & 0 \end{bmatrix} $$

$$ C^T C = \begin{bmatrix} 0 & -1 \\ 1 & 0 \end{bmatrix} \begin{bmatrix} 0 & 1 \\ -1 & 0 \end{bmatrix} = \begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix} = I $$

Example 4: Determine if matrix is singular by finding determinant

$$ D = \begin{bmatrix} 4 & 2 \\ 6 & 3 \end{bmatrix} $$

Calculate determinant:

$$ |D| = (4)(3) - (2)(6) = 12 - 12 = 0 $$

Since determinant is zero, \(D\) is singular.

Example 5: Diagonalize matrix

Given:

$$ A = \begin{bmatrix} 4 & 1 \\ 0 & 3 \end{bmatrix} $$

Eigenvalues: \(\lambda_1 = 4\), \(\lambda_2 = 3\), so \(A\) is diagonalizable as:

$$ A = PDP^{-1}, \quad D = \begin{bmatrix} 4 & 0 \\ 0 & 3 \end{bmatrix} $$

\(P\) is formed by eigenvectors of \(A\).

Chapter Summary

This chapter introduced the concept of matrices, defining them and exploring various types such as row, column, square, rectangular, diagonal, identity, and zero matrices. The notation and elements of matrices were explained, followed by basic operations including addition, subtraction, and scalar multiplication. Special matrices like identity, diagonal, symmetric, and skew-symmetric matrices were also covered. The chapter then detailed matrix multiplication, outlining the rules and properties including associative and distributive laws, and discussed the transpose of a matrix with its characteristics. Determinants for 2x2 and 3x3 matrices were calculated, and methods for finding inverses, such as the adjoint method and Gauss-Jordan elimination, were examined along with the conditions necessary for an inverse to exist. The use of matrices to represent systems of linear equations was introduced, followed by the application of Cramer's rule and the inverse matrix method to solve these systems. Finally, the inverse determinant method was presented as another approach for solving linear systems efficiently.