CELE Engineering Mathematics — Advanced Engineering MathematicsStudy Notes

A complex number z represents a two-dimensional quantity in the form z = a + bi, where a is the real part, b is the imaginary part, and i² = −1. In engineering, complex numbers model AC quantities (voltage, current, impedance) and oscillatory systems. Understanding multiple representations of the same number is crucial for efficient computation. **Rectangular Form:** z = a + bi - a = Re(z) (real component) - b = Im(z) (imaginary component) - Useful for addition and subtraction **Polar Form:** z = r∠θ = r(cos θ + i sin θ) - r = |z| = magnitude (modulus) - θ = argument (angle from positive real axis, measured counterclockwise) - Useful for multiplication, division, and powers **Exponential Form:** z = re^(iθ) (Euler's formula) - Compact representation, simplifies exponential operations - Connection: e^(iθ) = cos θ + i sin θ **Conversion Formulas:** Magnitude: |z| = √(a² + b²) Argument: θ = arctan(b/a), with careful attention to quadrant - If a > 0, b > 0: First quadrant, θ = arctan(b/a) - If a < 0, b > 0: Second quadrant, θ = 180° − arctan|b/a| - If a < 0, b < 0: Third quadrant, θ = 180° + arctan|b/a| - If a > 0, b < 0: Fourth quadrant, θ = 360° − arctan|b/a| or −arctan|b/a| From polar to rectangular: - a = r cos θ - b = r sin θ **Engineering Context:** In AC circuit analysis per IEEE and PRC standards, impedance Z = R + jX (using j for the imaginary unit in engineering) combines resistance R and reactance X. The magnitude |Z| gives the total impedance, while the angle determines the phase shift between voltage and current.

Operations on complex numbers follow rules that depend on the form used. In general, rectangular form is preferred for addition and subtraction (combine real and imaginary parts separately), while polar form is preferred for multiplication and division (multiply/divide magnitudes, add/subtract angles). **Addition and Subtraction (Rectangular):** z₁ ± z₂ = (a₁ ± a₂) + i(b₁ ± b₂) Example: (3 + 4i) + (1 − 2i) = 4 + 2i Example: (5 − 3i) − (2 + i) = 3 − 4i **Multiplication:** Rectangular method: (a + bi)(c + di) = ac − bd + i(ad + bc) This follows FOIL, remembering that i² = −1. Polar method (superior for most engineering calculations): z₁ · z₂ = r₁r₂∠(θ₁ + θ₂) Multiply magnitudes, add angles. **Division:** Rectangular method (requires rationalization by conjugate): (a + bi)/(c + di) = [(a + bi)(c − di)]/[(c + di)(c − di)] = [(ac + bd) + i(bc − ad)]/(c² + d²) Polar method (superior for engineering): z₁/z₂ = (r₁/r₂)∠(θ₁ − θ₂) Divide magnitudes, subtract angles. **De Moivre's Theorem (Powers and Roots):** For integer power n: zⁿ = [r∠θ]ⁿ = rⁿ∠(nθ) For the n-th root (n distinct values): z^(1/n) = r^(1/n)∠[(θ + 360°k)/n], where k = 0, 1, 2, ..., n−1 This is fundamental for solving equations like "find all fifth roots" or analyzing vibration modes with complex eigenvalues. **Conjugate Operations:** Complex conjugate of z = a + bi is z* = a − bi Properties: |z*| = |z|, arg(z*) = −arg(z), z · z* = |z|² Used in division and in phasor analysis for AC circuits.

A matrix is a rectangular array of numbers organized in rows and columns. For civil engineering applications, matrices organize data in structural analysis (stiffness matrices), finite element methods, transformation matrices for coordinate changes, and solution of simultaneous equations. **Matrix Notation and Dimensions:** A matrix A with m rows and n columns is denoted A(m×n) or Aₘₙ. Element in row i, column j is denoted aᵢⱼ. Example 3×2 matrix: A = [1 2] [3 4] [5 6] **Types of Matrices:** - Square matrix: m = n (important for determinants and eigenvalues) - Identity matrix I: diagonal matrix with 1's on diagonal, 0's elsewhere I₃ = [1 0 0] [0 1 0] [0 0 1] - Zero matrix: all elements are 0 - Diagonal matrix: non-zero only on main diagonal - Symmetric matrix: A = Aᵀ (transpose equals original; common in structural stiffness) - Skew-symmetric matrix: A = −Aᵀ **Matrix Addition and Subtraction:** Only matrices of same dimensions can be added/subtracted. Add/subtract element-wise. A + B = [aᵢⱼ + bᵢⱼ] Example: [1 2] + [4 0] = [5 2] [3 4] [1 2] [4 6] **Scalar Multiplication:** cA = [c·aᵢⱼ] (multiply each element by scalar c) Example: 2·[1 2] = [2 4] [3 4] [6 8] **Matrix Multiplication:** Product AB is defined only if number of columns in A equals number of rows in B. If A is (m×n) and B is (n×p), then AB is (m×p). Element (AB)ᵢⱼ = Σₖ aᵢₖ·bₖⱼ (dot product of row i of A with column j of B) CRITICAL: Matrix multiplication is NOT commutative. AB ≠ BA in general. Example: A(2×3) · B(3×2) = C(2×2) A = [1 2 3], B = [4 5] [4 5 6] [6 7] [8 9] C₁₁ = (1)(4) + (2)(6) + (3)(8) = 4 + 12 + 24 = 40 C₁₂ = (1)(5) + (2)(7) + (3)(9) = 5 + 14 + 27 = 46 C₂₁ = (4)(4) + (5)(6) + (6)(8) = 16 + 30 + 48 = 94 C₂₂ = (4)(5) + (5)(7) + (6)(9) = 20 + 35 + 54 = 109 C = [40 46] [94 109] **Transpose of a Matrix:** Aᵀ is obtained by swapping rows and columns. If A is (m×n), then Aᵀ is (n×m). In structural analysis, transpose is used when converting between different coordinate systems. **Determinant (for Square Matrices Only):** The determinant det(A) or |A| is a scalar that provides critical information: - If det(A) ≠ 0: Matrix is invertible (has unique solution in systems of equations) - If det(A) = 0: Matrix is singular (no unique solution) - Magnitude relates to volume scaling in transformations **2×2 Determinant:** det([a b]) = ad − bc [c d] **3×3 Determinant (Cofactor Expansion along first row):** det(A) = a₁₁·det([a₂₂ a₂₃]) − a₁₂·det([a₂₁ a₂₃]) + a₁₃·det([a₂₁ a₂₂]) [a₃₂ a₃₃] [a₃₁ a₃₃] [a₃₁ a₃₂] **Properties of Determinants:** - det(AB) = det(A)·det(B) - det(Aᵀ) = det(A) - Swapping two rows changes sign of determinant - Multiplying a row by scalar c multiplies determinant by c - If two rows are identical, determinant is 0 - det(I) = 1 (identity matrix) **Matrix Inverse:** For square matrix A with det(A) ≠ 0: A·A⁻¹ = A⁻¹·A = I For 2×2 matrix: A⁻¹ = (1/det A)·[d -b] [-c a] where A = [a b] [c d] For larger matrices, use cofactor method, row reduction, or numerical software. **Engineering Significance:** In finite element analysis and structural stiffness methods, the global stiffness matrix K is square and (if properly constrained) has non-zero determinant, allowing unique solution K·D = F (stiffness × displacement = force).

Cramer's Rule provides an elegant determinant-based method for solving systems of linear equations Ax = b, where A is the coefficient matrix, x is the unknown vector, and b is the right-hand side vector. This method is particularly useful for 2×2 and 3×3 systems common in structural analysis and circuit problems. **General Form of Linear System:** a₁₁x₁ + a₁₂x₂ + ... + a₁ₙxₙ = b₁ a₂₁x₁ + a₂₂x₂ + ... + a₂ₙxₙ = b₂ ... aₙ₁x₁ + aₙ₂x₂ + ... + aₙₙxₙ = bₙ In matrix form: Ax = b **Cramer's Rule:** If det(A) ≠ 0 (matrix is non-singular), then: xᵢ = det(Aᵢ)/det(A) where Aᵢ is the matrix obtained by replacing the i-th column of A with the vector b. **Conditions for Use:** 1. System must have exactly as many equations as unknowns (square system) 2. det(A) ≠ 0 (ensures unique solution) 3. If det(A) = 0 and det(Aᵢ) ≠ 0 for some i: no solution exists 4. If det(A) = 0 and det(Aᵢ) = 0 for all i: infinitely many solutions (system is dependent) **Advantages:** - Conceptually elegant; shows each unknown as ratio of determinants - Useful for hand calculation of small systems - Reveals when unique solution exists (det(A) ≠ 0) **Disadvantages:** - Computationally expensive for large systems (many determinants to calculate) - Numerically unstable with ill-conditioned matrices - For large systems, Gaussian elimination is preferred **Engineering Applications:** - Finding unknown forces in indeterminate trusses (method of joints) - Analyzing multi-loop circuits using node voltage method - Solving for unknown reactions and internal forces in beams - Finite element analysis with small element matrices

Vectors are mathematical quantities with both magnitude and direction. In civil engineering, vectors represent forces, displacements, velocities, and stress resultants. Two primary vector operations—dot product and cross product—are essential for force equilibrium, moment calculations, and three-dimensional geometry. **Vector Representation:** A vector in 3D space is written as: v = (vₓ, vᵧ, vᵢ) or v = vₓi + vᵧj + vᵢk where i, j, k are unit vectors along x, y, z axes respectively. Magnitude (length): |v| = √(vₓ² + vᵧ² + vᵢ²) **Dot Product (Scalar Product):** The dot product of vectors a and b produces a scalar (single number): a · b = aₓbₓ + aᵧbᵧ + aᵢbᵢ (component form) Alternatively: a · b = |a||b|cos(θ) where θ is the angle between the vectors. **Properties of Dot Product:** - Commutative: a · b = b · a - Distributive: a · (b + c) = a · b + a · c - If a · b = 0: vectors are perpendicular (θ = 90°) - If a · b > 0: acute angle (θ < 90°) - If a · b < 0: obtuse angle (θ > 90°) - a · a = |a|² (magnitude squared) **Engineering Applications of Dot Product:** - Finding component of force in a given direction: F_dir = (F · u)u, where u is unit vector - Determining if forces/displacements are perpendicular - Calculating work: W = F · d = |F||d|cos(θ) - Finding angle between two members in a structure **Cross Product (Vector Product):** The cross product of vectors a and b produces a vector perpendicular to both: a × b = |i j k | |aₓ aᵧ aᵢ | |bₓ bᵧ bᵢ | Expanding: a × b = (aᵧbᵢ − aᵢbᵧ)i − (aₓbᵢ − aᵢbₓ)j + (aₓbᵧ − aᵧbₓ)k Alternatively: |a × b| = |a||b|sin(θ) where θ is the angle between vectors. Direction: Use right-hand rule. If fingers curl from a to b, thumb points in direction of a × b. **Properties of Cross Product:** - Non-commutative: a × b = −(b × a) - Distributive: a × (b + c) = a × b + a × c - If a × b = 0: vectors are parallel or antiparallel (θ = 0° or 180°) - |a × b| equals area of parallelogram formed by a and b - Perpendicular to both a and b - i × j = k, j × k = i, k × i = j **Engineering Applications of Cross Product:** - Calculating moment (torque): M = r × F (position × force) - Finding area of irregular surfaces - Determining normal to a plane defined by two vectors - Analyzing three-dimensional stress and strain states - Determining whether forces cause rotation about a point **Angle Between Two Vectors:** From dot product definition: cos(θ) = (a · b)/(|a||b|) θ = arccos[(a · b)/(|a||b|)] This is useful for finding angles between members in trusses or angle of incidence in wind load analysis.

Understanding vector magnitude and direction is fundamental for translating physical quantities into mathematical form. Direction cosines provide a systematic way to express vector direction in 3D space, essential for analyzing loads on inclined surfaces and member orientations in 3D trusses. **Vector Magnitude:** For vector v = (vₓ, vᵧ, vᵢ): |v| = √(vₓ² + vᵧ² + vᵢ²) Magnitude is always non-negative and represents the "length" of the vector in physical units (forces in kN, displacements in mm, etc.). **Unit Vector:** A unit vector has magnitude 1 and indicates direction only. Any vector can be converted to a unit vector by dividing by its magnitude: û = v/|v| = (vₓ/|v|, vᵧ/|v|, vᵢ/|v|) Unit vector is dimensionless and purely directional. Unit vectors along coordinate axes: i = (1, 0, 0), j = (0, 1, 0), k = (0, 0, 1) **Direction Cosines:** For a vector v making angles α, β, γ with the positive x, y, z axes respectively: cos(α) = vₓ/|v| = component of unit vector in x-direction cos(β) = vᵧ/|v| = component of unit vector in y-direction cos(γ) = vᵢ/|v| = component of unit vector in z-direction These are called direction cosines. They satisfy the fundamental relationship: cos²(α) + cos²(β) + cos²(γ) = 1 This identity is crucial for checking if direction is properly specified; if the sum of squares doesn't equal 1, direction specification is incomplete or inconsistent. **Converting Between Forms:** 1. From vector components to direction cosines: Calculate |v|, then divide each component by magnitude. 2. From direction cosines to vector: If direction cosines are known and magnitude |v| is known: vₓ = |v|cos(α) vᵧ = |v|cos(β) vᵢ = |v|cos(γ) 3. From two points to direction: If vector extends from point A(x₁, y₁, z₁) to B(x₂, y₂, z₂): v = (x₂−x₁, y₂−y₁, z₂−z₁) Then calculate magnitude and direction cosines as above. **Engineering Applications:** - Specifying cable or member orientation in 3D structures - Determining components of inclined loads (e.g., wind on sloped roof) - Converting between local and global coordinate systems - Analyzing strut forces in space frames - Finite element orientation matrices