Session 12 - Engineering With Matrices

Skill Block 1

Skill Block 2

Worked Example

Environmental Engineering – Fluid Flow in Reservoirs

Context: This problem utilizes a $3 \times 3$ matrix inverse to solve an environmental tracking scenario.
Problem: The concentration levels of a tracking dye in three cascading retention ponds ( $C_1, C_2, C_3$ ) satisfy the following steady-state system:

$\begin{cases} 4C_1 - C_2 = 8 \\ -C_1 + 3C_2 - C_3 = 12 \\ -C_2 + 2C_3 = 4 \end{cases}$

Using a matrix calculator, find the inverse of the coefficient matrix ( $A^{-1}$ ) and solve for the concentration vector $x = A^{-1}b$ .

Worked Solution:

Step 1: Write the System as a Matrix Equation ( $Ax = b$ )

To construct the matrix form, align the concentration variables ( $C_1$ , $C_2$ , and $C_3$ ) in order for all three equations. If a variable is missing from a specific pond's equation, insert a 0 coefficient to preserve the alignment:

Pond 1: $4C_1 - 1C_2 + 0C_3 = 8$
Pond 2: $-1C_1 + 3C_2 - 1C_3 = 12$
Pond 3: $0C_1 - 1C_2 + 2C_3 = 4$

Now we separate this system into the standard $Ax = b$ matrix components:

Coefficient Matrix ( $A$ ):

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

Variable Column Vector ( $x$ ):

x = \begin{bmatrix} C_1 \\ C_2 \\ C_3 \end{bmatrix}

Constant Column Vector ( $b$ ):

b = \begin{bmatrix} 8 \\ 12 \\ 4 \end{bmatrix}

Combined matrix representation:

\begin{bmatrix} 4 & -1 & 0 \\ -1 & 3 & -1 \\ 0 & -1 & 2 \end{bmatrix} \begin{bmatrix} C_1 \\ C_2 \\ C_3 \end{bmatrix} = \begin{bmatrix} 8 \\ 12 \\ 4 \end{bmatrix}

Step 2: Compute the Inverse Matrix ( $A^{-1}$ )

Following the problem instructions to use a matrix calculator to find $A^{-1}$ , entering matrix $A$ provides the following exact fraction outputs and corresponding decimal approximations:

A^{-1} = \begin{bmatrix} \frac{5}{18} & \frac{1}{9} & \frac{1}{18} \\ \frac{1}{9} & \frac{4}{9} & \frac{2}{9} \\ \frac{1}{18} & \frac{2}{9} & \frac{11}{18} \end{bmatrix} \approx \begin{bmatrix} 0.2778 & 0.1111 & 0.0556 \\ 0.1111 & 0.4444 & 0.2222 \\ 0.0556 & 0.2222 & 0.6111 \end{bmatrix}

Step 3: Solve for the Concentration Vector ( $x = A^{-1}b$ )

Now, multiply the inverse matrix $A^{-1}$ by the constants vector $b$ to isolate the unknown concentrations:

x = A^{-1}b

x = \begin{bmatrix} \frac{5}{18} & \frac{1}{9} & \frac{1}{18} \\ \frac{1}{9} & \frac{4}{9} & \frac{2}{9} \\ \frac{1}{18} & \frac{2}{9} & \frac{11}{18} \end{bmatrix} \begin{bmatrix} 8 \\ 12 \\ 4 \end{bmatrix}

Perform row-by-column matrix multiplication:

Top Row ( $C_1$ ):

\left(\frac{5}{18} \cdot 8\right) + \left(\frac{1}{9} \cdot 12\right) + \left(\frac{1}{18} \cdot 4\right) = \frac{40}{18} + \frac{24}{18} + \frac{4}{18} = \frac{68}{18} = \frac{34}{9}

Middle Row ( $C_2$ ):

\left(\frac{1}{9} \cdot 8\right) + \left(\frac{4}{9} \cdot 12\right) + \left(\frac{2}{9} \cdot 4\right) = \frac{8}{9} + \frac{48}{9} + \frac{8}{9} = \frac{64}{9}

Bottom Row ( $C_3$ ):

\left(\frac{1}{18} \cdot 8\right) + \left(\frac{2}{9} \cdot 12\right) + \left(\frac{11}{18} \cdot 4\right) = \frac{8}{18} + \frac{48}{18} + \frac{44}{18} = \frac{100}{18} = \frac{50}{9}

This gives the final solution vector:

x = \begin{bmatrix} \frac{34}{9} \\ \frac{64}{9} \\ \frac{50}{9} \end{bmatrix}

Final Engineering Interpretation

Interpreting each computed value back into the environmental tracking context:

Dye concentration in Pond 1 ( $C_1$ ): $\frac{34}{9} \approx \mathbf{3.78 \text{ units}}$
Dye concentration in Pond 2 ( $C_2$ ): $\frac{64}{9} \approx \mathbf{7.11 \text{ units}}$
Dye concentration in Pond 3 ( $C_3$ ): $\frac{50}{9} \approx \mathbf{5.56 \text{ units}}$

Active Practice Problem

Civil Engineering – Indeterminate Foundation Supports

Context: This problem highlights what happens when a system lacks a unique solution due to redundant structural constraints.
Problem: A concrete foundation beam rests on three vertical supports. Due to static indeterminacy, the vertical equilibrium and moment balances only yield two equations for the three unknown reaction forces ( $R_1, R_2, R_3$ ):

$\begin{cases} R_1 + R_2 + R_3 = 80 \\ 2R_1 + 4R_2 + 6R_3 = 240 \end{cases}$

Explain mathematically why this system cannot yield a single, unique solution vector. What does this condition mean for a civil engineer trying to calculate the exact load on each individual support?

Workspace

Solution

Step 1: Write the System in Matrix Form ( $Ax = b$ )

To see how the numbers interact, let's group the system of equations into the standard matrix structure:

\begin{cases} R_1 + R_2 + R_3 = 80 \\ 2R_1 + 4R_2 + 6R_3 = 240 \end{cases}

Coefficient Matrix ( $A$ ): Pulling the weights from the reactions gives a rectangular layout with 2 rows and 3 columns.

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

Variable Vector ( $x$ ): The 3 unknown foundation forces.

x = \begin{bmatrix} R_1 \\ R_2 \\ R_3 \end{bmatrix}

Constant Vector ( $b$ ): The loading conditions.

b = \begin{bmatrix} 80 \\ 240 \end{bmatrix}

Step 2: The Mathematical Explanation (Why it lacks a unique solution)

There are two primary reasons why this matrix configuration cannot isolate a single unique answer:

1. The System is Underdetermined (Non-Square Matrix)

A unique solution requires a perfect one-to-one balance between constraints and unknowns. Because matrix $A$ is a rectangular $2 \times 3$ shape, it represents fewer equations (2 constraints) than variables (3 unknowns).

2. The Matrix Inversion Method Fails

To compute a unique solution vector via $x = A^{-1}b$ , the coefficient matrix must be square ( $n \times n$ ) to calculate a determinant and find an inverse ( $A^{-1}$ ). Because $A$ is not square, $\det(A)$ is undefined and $A^{-1}$ does not exist.

3. The Equations "Overlap" (Consistent, Dependent System)

If you reduce the rows to see how they interact algebraically:

Multiply the first equation by 2: $2R_1 + 2R_2 + 2R_3 = 160$
Subtract this from the second equation: $(2R_1 + 4R_2 + 6R_3) - (2R_1 + 2R_2 + 2R_3) = 240 - 160$
This simplifies to: $2R_2 + 4R_3 = 80 \implies R_2 + 2R_3 = 40$

We can choose any value for $R_3$ (making it a "free variable"), and calculate valid matching numbers for $R_1$ and $R_2$ . Because the equations are mathematically compatible but leave a variable free, the system falls under the Consistent, Dependent category, meaning it contains infinitely many solutions.

Step 3: Physical Engineering Interpretation

For a civil engineer, this mathematical condition represents static indeterminacy.

When a foundation beam rests on three rigid supports, basic rigid-body physics equations (balancing upward forces and balancing rotational moments) only give you 2 equations of equilibrium. This is physically not enough information to pinpoint how the load splits among the columns.

What it means for the designer:

Statics alone is insufficient: You cannot solve the true load on each individual support using simple force balances.
Material properties matter: In reality, the exact distribution of those 80 units of force depends heavily on physical variables missing from basic linear setups —such as how much the concrete beam flexes, the elasticity of the support columns, and whether the ground soil undergoes minor settling underneath one of the locations.
Next Steps: To get a single exact target answer, the engineer must upgrade their structural model to incorporate material deformation equations (compatibility equations) alongside these base matrix constraints.

Extra Practice & Extensions

Context: This problem requires evaluating the determinant to catch conflicting physical models.
Problem: An engineer attempts to calibrate a dual-mass dampening system and sets up the following matrix equation:

$\begin{bmatrix} 3 & -6 \\ -1 & 2 \end{bmatrix} \begin{bmatrix} x_1 \\ x_2 \end{bmatrix} = \begin{bmatrix} 9 \\ 4 \end{bmatrix}$
1. Compute the determinant of the coefficient matrix, $\det(A)$ .
2. Using the criteria for types of systems, classify this system as Consistent Independent, Consistent Dependent, or Inconsistent.
3. Interpret what this mathematical result indicates about the physical parameters chosen for the dampening system.

Skill Block 3

Worked Example

Industrial Engineering – Production Line Optimization

Context: A multi-step optimization problem requiring both representation and numerical inversion.
Problem: A manufacturing plant uses two primary automated assembly machines ( $M_1$ and $M_2$ ) to produce two types of medical sensors (Sensor A and Sensor B).
- Each unit of Sensor A requires 3 hours on $M_1$ and 2 hours on $M_2$ .
- Each unit of Sensor B requires 4 hours on $M_1$ and 5 hours on $M_2$ .
- The factory has exactly 70 hours of available time on $M_1$ and 70 hours of available time on $M_2$ this week.
1. Set up a system of linear equations and write it in the matrix form $Ax = b$ .
2. Compute $A^{-1}$ and solve the matrix equation to find exactly how many units of Sensor A and Sensor B should be produced to completely utilize the available machine hours.

Worked Solution:

Here is the complete step-by-step solution to represent and solve the production line optimization problem using the matrix inverse method.

Part 1: Set up the System and Matrix Form ( $Ax = b$ )

First, define the decision variables for the manufacturing line:

$x_1 =$ number of units of Sensor A to produce
$x_2 =$ number of units of Sensor B to produce

Using the time allocations specified for each machine, we establish the following linear system:

\begin{cases} 3x_1 + 4x_2 = 70 \quad \text{(Machine 1 hours)} \\ 2x_1 + 5x_2 = 70 \quad \text{(Machine 2 hours)} \end{cases}

To express this production system in the standard matrix form $Ax = b$ :

Coefficient Matrix ( $A$ ): Captures the hourly processing rates for both sensors on each machine.

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

Variable Column Vector ( $x$ ): The column containing our production unknowns.

x = \begin{bmatrix} x_1 \\ x_2 \end{bmatrix}

Constant Column Vector ( $b$ ): The maximum available operating capacity for the week.

b = \begin{bmatrix} 70 \\ 70 \end{bmatrix}

Combining these terms yields the full matrix setup:

\begin{bmatrix} 3 & 4 \\ 2 & 5 \end{bmatrix} \begin{bmatrix} x_1 \\ x_2 \end{bmatrix} = \begin{bmatrix} 70 \\ 70 \end{bmatrix}

Part 2: Compute the Inverse Matrix ( $A^{-1}$ )

To isolate our production target variables, we need to find the inverse of the $2 \times 2$ coefficient matrix.

Calculate the determinant ( $\det(A)$ ):

\det(A) = (3)(5) - (4)(2)

\det(A) = 15 - 8 = 7

Because $\det(A) \neq 0$ , the matrix is non-singular, confirming a single unique production schedule exists.

Form the Adjugate Matrix:

Swap the elements on the main diagonal (3 and 5 swap places) and switch the signs of the off-diagonal elements (4 and 2 become $-4$ and $-2$ ):

\text{Adjugate}(A) = \begin{bmatrix} 5 & -4 \\ -2 & 3 \end{bmatrix}

Scale by $1/\det(A)$ :

A^{-1} = \frac{1}{7} \begin{bmatrix} 5 & -4 \\ -2 & 3 \end{bmatrix} = \begin{bmatrix} \frac{5}{7} & -\frac{4}{7} \\ -\frac{2}{7} & \frac{3}{7} \end{bmatrix}

Part 3: Solve the Matrix Equation ( $x = A^{-1}b$ )

Now, multiply the calculated inverse matrix by our target hour vector $b$ . Keep in mind that the order of operations matters here because matrix multiplication is non-commutative ( $A^{-1}b \neq bA^{-1}$ ):

x = A^{-1}b

x = \frac{1}{7} \begin{bmatrix} 5 & -4 \\ -2 & 3 \end{bmatrix} \begin{bmatrix} 70 \\ 70 \end{bmatrix}

Perform row-by-column multiplication:

Top Row Calculation ( $x_1$ ): $(5)(70) + (-4)(70) = 350 - 280 = 70$
Bottom Row Calculation ( $x_2$ ): $(-2)(70) + (3)(70) = -140 + 210 = 70$

Divide the calculated rows by our remaining scaling factor of 7:

x = \frac{1}{7} \begin{bmatrix} 70 \\ 70 \end{bmatrix} = \begin{bmatrix} 10 \\ 10 \end{bmatrix}

Final Engineering Interpretation

Interpreting each specific variable's value back into our plant optimization context:

Sensor A Target ( $x_1$ ): 10 units
Sensor B Target ( $x_2$ ): 10 units

By scheduling exactly 10 units of Sensor A and 10 units of Sensor B for production, the plant will completely utilize the 70 hours of operational capacity on both machine $M_1$ and machine $M_2$ without leaving any idle machine time.

Active Practice Problem

Aerospace Engineering – Flight Trim Configurations

Context: This problem tests the critical step of checking the determinant before executing an automated simulation solution.
Problem: An automated flight control system calculates three elevator trim parameters ( $p_1, p_2, p_3$ ) using the following matrix array:

$A = \begin{bmatrix} 1 & 2 & 3 \\ 2 & 4 & 6 \\ 3 & -1 & 4 \end{bmatrix}, \quad b = \begin{bmatrix} 10 \\ 20 \\ 5 \end{bmatrix}$
1. Before attempting to compute $A^{-1}$ , calculate $\det(A)$ .
2. Based on your determinant calculations, will the flight control computer be able to find a unique baseline trim setting using the matrix inverse method ( $x = A^{-1}b$ )? Why or why not?
3. Identify if this specific system has infinitely many solutions or no solution by checking how the rows interact.

Solution Key:

Part 1: Calculate the Determinant, $\det(A)$

Given the system's coefficient matrix:

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

To calculate the determinant of a $3 \times 3$ matrix, we can use cofactor expansion along the first row:

\det(A) = 1 \cdot \det\begin{bmatrix} 4 & 6 \\ -1 & 4 \end{bmatrix} - 2 \cdot \det\begin{bmatrix} 2 & 6 \\ 3 & 4 \end{bmatrix} + 3 \cdot \det\begin{bmatrix} 2 & 4 \\ 3 & -1 \end{bmatrix}

Now, compute the $2 \times 2$ determinants using the $ad - bc$ cross-multiplication method:

First component: $\det\begin{bmatrix} 4 & 6 \\ -1 & 4 \end{bmatrix} = (4)(4) - (6)(-1) = 16 + 6 = 22$
Second component: $\det\begin{bmatrix} 2 & 6 \\ 3 & 4 \end{bmatrix} = (2)(4) - (6)(3) = 8 - 18 = -10$
Third component: $\det\begin{bmatrix} 2 & 4 \\ 3 & -1 \end{bmatrix} = (2)(-1) - (4)(3) = -2 - 12 = -14$

Substitute these values back into the main expansion equation:

\det(A) = 1(22) - 2(-10) + 3(-14)

\det(A) = 22 + 20 - 42

\det(A) = 42 - 42 = 0

Part 2: Determine if a Unique Solution Exists via the Inverse Method

Based on the calculated determinant, no, the flight control computer will not be able to find a unique baseline trim setting using the matrix inverse method ( $x = A^{-1}b$ ).

Why? A system only possesses a single unique solution if its coefficient matrix is invertible. Matrix inversion requires dividing by the determinant. Because $\det(A) = 0$ , the matrix $A$ is classified as singular (non-invertible). Attempting to evaluate $x = A^{-1}b$ would force an impossible mathematical division-by-zero error inside the simulation script.

Part 3: Identify the System Type by Analyzing Row Interactions

Since $\det(A) = 0$ , the system must be classified as either Consistent Dependent (infinitely many solutions) or Inconsistent (no solution). To identify which one it is, let's look closely at how the rows of the equations interact:

Equation 1 (Row 1): $1p_1 + 2p_2 + 3p_3 = 10$
Equation 2 (Row 2): $2p_1 + 4p_2 + 6p_3 = 20$
Equation 3 (Row 3): $3p_1 - 1p_2 + 4p_3 = 5$

Notice the exact proportional relationship between Row 1 and Row 2:

\text{Row 2} = 2 \times \text{Row 1}

2(1p_1 + 2p_2 + 3p_3) = 2(10) \implies 2p_1 + 4p_2 + 6p_3 = 20

Because the scaling applies perfectly across both the coefficients and the target constant on the right side, the equations do not contradict each other. Instead, Equation 2 is completely redundant—it provides zero new physical constraints to the system.

According to our classification framework:

Condition: $\det(A) = 0$ and the constraints perfectly overlap rather than contradict.
System Type: Consistent, Dependent.
Meaning: The system has infinitely many solutions.

Flight System Diagnostic Interpretation

For an aerospace engineer calibrating automated control surfaces, this diagnostic output reveals a critical design trait:

Under-constrained Configuration: Because one equation is redundant, you effectively only have 2 unique physical constraints governing 3 independent trim parameters ( $p_1, p_2, p_3$ ).
Physical Impact: The flight control computer cannot pick out one perfect baseline configuration on its own because there is a whole mathematical line of infinite combinations that will technically satisfy this balance. To fix this "free variable" problem, the software must either introduce an additional flight constraint (like prioritizing the setting that consumes the lowest battery power) or a human operator must manually freeze one of the three trim parameters to lock down a single baseline behavior.

Show your steps and calculations here.

MATH-182 (Calculus II) Curriculum

Session 12 - Engineering With Matrices

Micro-Lecture

Engineering Context: Engineering Calculus II Application

Skill Block 1

Regroup 1

Skill Block 2

Worked Example

Step 1: Write the System as a Matrix Equation ( $Ax = b$ )

Step 2: Compute the Inverse Matrix ( $A^{-1}$ )

Step 3: Solve for the Concentration Vector ( $x = A^{-1}b$ )

Final Engineering Interpretation

Active Practice Problem

Workspace

Solution

Step 1: Write the System in Matrix Form ( $Ax = b$ )

Step 2: The Mathematical Explanation (Why it lacks a unique solution)

1. The System is Underdetermined (Non-Square Matrix)

2. The Matrix Inversion Method Fails

3. The Equations "Overlap" (Consistent, Dependent System)

Step 3: Physical Engineering Interpretation

What it means for the designer:

Extra Practice & Extensions

Regroup 2

Skill Block 3

Worked Example

Part 1: Set up the System and Matrix Form ( $Ax = b$ )

Part 2: Compute the Inverse Matrix ( $A^{-1}$ )

Part 3: Solve the Matrix Equation ( $x = A^{-1}b$ )

Final Engineering Interpretation

Active Practice Problem

Part 1: Calculate the Determinant, $\det(A)$

Part 2: Determine if a Unique Solution Exists via the Inverse Method

Part 3: Identify the System Type by Analyzing Row Interactions

Flight System Diagnostic Interpretation

Regroup 3

Synthesis Wrap-up

MATH-182 (Calculus II) Curriculum

Session 12 - Engineering With Matrices

Micro-Lecture

Engineering Context: Engineering Calculus II Application

Skill Block 1

Regroup 1

Skill Block 2

Worked Example

Step 1: Write the System as a Matrix Equation (Ax=bAx = bAx=b)

Step 2: Compute the Inverse Matrix (A−1A^{-1}A−1)

Step 3: Solve for the Concentration Vector (x=A−1bx = A^{-1}bx=A−1b)

Final Engineering Interpretation

Active Practice Problem

Workspace

Solution

Step 1: Write the System in Matrix Form (Ax=bAx = bAx=b)

Step 2: The Mathematical Explanation (Why it lacks a unique solution)

1. The System is Underdetermined (Non-Square Matrix)

2. The Matrix Inversion Method Fails

3. The Equations "Overlap" (Consistent, Dependent System)

Step 3: Physical Engineering Interpretation

What it means for the designer:

Extra Practice & Extensions

Regroup 2

Skill Block 3

Worked Example

Part 1: Set up the System and Matrix Form (Ax=bAx = bAx=b)

Part 2: Compute the Inverse Matrix (A−1A^{-1}A−1)

Part 3: Solve the Matrix Equation (x=A−1bx = A^{-1}bx=A−1b)

Final Engineering Interpretation

Active Practice Problem

Part 1: Calculate the Determinant, det⁡(A)\det(A)det(A)

Part 2: Determine if a Unique Solution Exists via the Inverse Method

Part 3: Identify the System Type by Analyzing Row Interactions

Flight System Diagnostic Interpretation

Regroup 3

Synthesis Wrap-up

Step 1: Write the System as a Matrix Equation ( $Ax = b$ )

Step 2: Compute the Inverse Matrix ( $A^{-1}$ )

Step 3: Solve for the Concentration Vector ( $x = A^{-1}b$ )

Step 1: Write the System in Matrix Form ( $Ax = b$ )

Part 1: Set up the System and Matrix Form ( $Ax = b$ )

Part 2: Compute the Inverse Matrix ( $A^{-1}$ )

Part 3: Solve the Matrix Equation ( $x = A^{-1}b$ )

Part 1: Calculate the Determinant, $\det(A)$