UNIT 4: CONTROL SYSTEM - EXAM-FOCUSED SHORT NOTES
1.0 SYSTEM MODELING & REPRESENTATION
1.1 Block Diagram Representation
-
Components:
-
Block: Represents a system/component with transfer function
G(s). -
Summing Point: Algebraic sum (
+or-) of input signals. -
Take-off Point: Point where signal is sampled for feedback or measurement.
-
Forward Path: Path from input
R(s)to outputC(s)containing blocks. -
Feedback Path: Path from output to summing point, transfer function
H(s).
-
-
Reduction Rules (Apply sequentially):
-
Series:
G₁G₂...Gₙ -
Parallel:
G₁ ± G₂ ± ... ± Gₙ -
Feedback:
T(s) = \frac{G(s)}{1 ± G(s)H(s)}(Negative feedback:1 + G(s)H(s)).
-
-
Open Loop vs. Closed Loop:
| Feature | Open Loop | Closed Loop | | :--- | :--- | :--- | | Feedback | Absent | Present | | Accuracy | Low, depends on calibration | High, corrects errors | | Stability | Inherently stable | Needs design for stability | | Complexity | Simple | Complex | | Sensitivity | High to parameter changes | Low to parameter changes |
[!TIP] Exam Focus: Block diagram reduction is foundational. Always simplify to a single forward path and single feedback loop before writing
T(s).
1.2 Signal Flow Graph (SFG) & Mason's Gain Formula
-
SFG: Graphical representation of linear algebraic equations.
-
Node: System variable (signal).
-
Branch: Directed edge with gain
G(s). -
Forward Path Gain: Product of branch gains along a path from input to output.
-
Loop Gain: Product of branch gains in a closed path.
-
Non-Touching Loops: Loops with no common nodes.
-
-
Mason's Gain Formula:
$$T(s) = \frac{C(s)}{R(s)} = \frac{\sum_{k=1}^{N} P_k \Delta_k}{\Delta}$$
Where:
* `P_k` = k-th forward path gain.
* `Δ` = `1 - (sum of all individual loop gains) + (sum of gain products of all possible two non-touching loops) - ...`
* `Δ_k` = Value of `Δ` for that part of the graph **not** touching the k-th forward path.
- Advantages over Block Diagram: Direct application, no need for successive reduction, handles complex interconnections easily.
[!TIP] Common Pitfall: While calculating
Δ_k, remove all nodes and branches of the k-th forward path, then computeΔfor the remaining subgraph.
1.3 Analogous Systems
-
Force-Voltage (F-V) / Direct Analogy:
| Mechanical (Translational) | Electrical (Voltage) | | :--- | :--- | | Force (F) | Voltage (V) | | Velocity (v) | Current (I) | | Displacement (x) | Charge (q) | | Mass (M) | Inductance (L) | | Damping (B) | Resistance (R) | | Spring (K) | Inverse Capacitance (1/C) |
-
Force-Current (F-I) / Inverse Analogy:
| Mechanical (Translational) | Electrical (Current) | | :--- | :--- | | Force (F) | Current (I) | | Velocity (v) | Voltage (V) | | Displacement (x) | Flux (φ) | | Mass (M) | Capacitance (C) | | Damping (B) | Conductance (1/R) | | Spring (K) | Inverse Inductance (1/L) |
[!TIP] Exam Tip: For rotational systems, replace Force
Fwith TorqueT, Velocityvwith Angular velocityω, Displacementxwith Angular displacementθ, MassMwith Moment of InertiaJ, DampingBwith Viscous friction coefficientB.
1.4 System Components & Actuators
-
AC Servomotor:
-
Principle: Two-phase induction motor. Two stator windings (90° apart): reference winding (voltage
V_r) and control winding (voltageV_c). Rotor is squirrel-cage. -
Assumptions:
-
Constant field flux (linear magnetization).
-
Small torque angle (linear torque-speed relation).
-
Negligible rotor inductance (
L_r ≈ 0).
-
-
Transfer Function Derivation:
Torque
T_e ∝ V_r V_c sin(θ) ≈ V_r V_c θ(for small θ).Mechanical equation:
J d²θ/dt² + B dθ/dt = T_e.After Laplace and linearization:
-
$$\frac{\Theta(s)}{V_c(s)} = \frac{K}{s(T_M s + 1)}$$
Where `K` = motor constant, `T_M = J/B` = mechanical time constant.
-
Stepper Motor:
-
Principle: Converts digital pulses into discrete angular rotations. Rotor moves in fixed step angles (e.g., 1.8°/step).
-
Types: Variable Reluctance (VR), Permanent Magnet (PM), Hybrid (HB).
-
Applications: Printers, plotters, CNC machines, robotics (open-loop position control).
-
-
Tacho-Generator:
-
Principle: DC generator where output voltage
V_t ∝ ω(angular speed). -
Use in Feedback: Provides a voltage proportional to shaft speed for speed control systems (velocity feedback). Transfer function:
V_t(s)/ω(s) = K_t.
-
2.0 TIME DOMAIN ANALYSIS & PERFORMANCE
2.1 Standard Test Input Signals
| Signal | Time Domain f(t) |
Laplace F(s) |
Use |
|---|---|---|---|
| Step | A u(t) |
A/s |
Test sudden change, steady-state error |
| Ramp | A t u(t) |
A/s² |
Test tracking of linearly increasing input |
| Parabolic | (A t²/2) u(t) |
A/s³ |
Test acceleration tracking |
| Impulse | δ(t) |
1 |
Test system's impulse response (inverse Laplace of G(s)) |
2.2 First-Order System Response
-
Transfer Function:
G(s) = \frac{K}{τs + 1} -
Unit Step Response (
R(s)=1/s):
$$c(t) = K \left(1 - e^{-t/τ}\right)$$
* **Time Constant `τ`**: Time to reach 63.2% of final value.
* **Rise Time `t_r`** (10% to 90%): `t_r ≈ 2.2τ`.
- Unit Ramp Response:
$$c(t) = K \left(t - τ + τ e^{-t/τ}\right)$$
* **Steady-State Error**: `e_ss = lim_{t→∞} (r(t) - c(t)) = τ`.
2.3 Second-Order System Response
- Standard Form:
$$G(s) = \frac{\omega_n^2}{s^2 + 2ζω_n s + \omega_n^2}$$
Where `ω_n` = natural frequency, `ζ` = damping ratio.
- Underdamped Response (
0 < ζ < 1):
$$c(t) = 1 - \frac{e^{-ζω_n t}}{\sqrt{1-ζ^2}} \sin(ω_d t + φ)$$
Where `ω_d = ω_n√(1-ζ²)` = damped frequency, `φ = arccos(ζ)`.
-
Key Performance Specifications:
- Maximum Overshoot
M_p:
- Maximum Overshoot
$$M_p = e^{\frac{-ζπ}{\sqrt{1-ζ^2}}} \times 100\%$$
\boxed{M_p = \exp\left(-\frac{ζπ}{\sqrt{1-ζ^2}}\right)}
* **Peak Time** `t_p`:
$$t_p = \frac{π}{ω_d} = \frac{π}{ω_n\sqrt{1-ζ^2}}$$
\boxed{t_p = \frac{π}{ω_n\sqrt{1-ζ^2}}}
* **Rise Time** `t_r` (0% to 100% for underdamped):
$$t_r ≈ \frac{π - φ}{ω_d}$$
* **Settling Time** `t_s` (2% criterion):
$$t_s ≈ \frac{4}{ζω_n}$$
* **Resonant Frequency** `ω_r` & **Resonance Peak** `M_r`:
$$ω_r = ω_n\sqrt{1-2ζ^2} \quad (for ζ < 1/\sqrt{2})$$
$$M_r = \frac{1}{2ζ\sqrt{1-ζ^2}}$$
-
Critically Damped (
ζ=1) & Overdamped (ζ>1):-
No overshoot.
-
Slower response than underdamped with same
ω_n. -
t_randt_slarger.
-
[!TIP] Derivation of
M_p: Setdc/dt = 0inc(t), solve fort = t_p, substitute back intoc(t_p).M_p = c(t_p) - 1.
2.4 Performance Specifications & ζ-ω_n Relationship
-
Trade-off: Higher
ζ→ LowerM_p, but longert_sandt_r. -
Effect of
ω_n: Higherω_n→ Faster response (smallert_p,t_r,t_s), but may increase sensitivity. -
Graph:
M_pvs.ζis exponentially decaying.M_p=16%corresponds toζ≈0.5;M_p=5%corresponds toζ≈0.7.
3.0 STABILITY ANALYSIS (TIME DOMAIN)
3.1 Routh-Hurwitz Criterion
-
Necessary & Sufficient Condition: All coefficients of characteristic polynomial
a_n s^n + ... + a_0 = 0must be positive AND no sign changes in the first column of Routh array. -
Routh Array Construction:
| sⁿ |
a_n|a_{n-2}|a_{n-4}| ... | | sⁿ⁻¹ |a_{n-1}|a_{n-3}|a_{n-5}| ... | | sⁿ⁻² |b_1 = (a_{n-1}a_{n-2} - a_n a_{n-3})/a_{n-1}|b_2| ... | | | sⁿ⁻³ |c_1 = (b_1 a_{n-3} - a_{n-1} b_2)/b_1| ... | | | | ... | ... | ... | | | -
Special Cases:
-
Row of Zeros: Indicates symmetrical roots on
jωaxis (marginal stability). Form Auxiliary Equation from row above, differentiate, replace row with coefficients. -
First Element Zero: Replace
ε → 0⁺in that element and continue.
-
-
Relative Stability: To check roots left of
s = -σline, substitutes = z - σand apply Routh to polynomial inz.
[!TIP] Exam Pattern: "Find range of K for stability" → Form characteristic equation → Build Routh array → Condition: All first column elements > 0 → Solve inequalities for
K.
3.2 Root Locus Technique
-
Definition: Plot of closed-loop pole locations as gain
Kvaries from0 → ∞. -
Construction Rules:
-
Branches: Equal to number of open-loop poles
n. -
Start/End Points: Start at OL poles (
K=0), end at OL zeros (K→∞) or ∞. -
Asymptotes (if
n > m):-
Number =
n - m -
Angles:
θ = \frac{(2q+1)180°}{n-m},q = 0,1,...,n-m-1 -
Centroid:
σ = \frac{\sum Re(poles) - \sum Re(zeros)}{n-m}
-
-
Angle of Departure/Arrival:
-
From pole
p_i:∠Departure = 180° - Σ(angles to other poles) + Σ(angles to zeros) -
To zero
z_i:∠Arrival = 180° - Σ(angles to other zeros) + Σ(angles to poles)
-
-
Breakaway/Break-in Points:
-
Solve
dK/ds = 0forK(s) = \frac{Product(s - zeros)}{Product(s - poles)}. -
Valid points: On real axis segment between poles/zeros.
-
-
Imaginary Axis Crossing: Use Routh criterion on characteristic equation
1 + KG(s)H(s) = 0or apply Angle Condition (∑∠G(jω)H(jω) = -180°).
-
-
Stability Judgment: Closed-loop system unstable if any root locus branch lies in Right-Half Plane (RHP) for
K > 0. -
Design Use: From root locus plot, read
Kfor desiredζ(draw line from origin at anglecos⁻¹(ζ)).
[!TIP] Sketching Order: 1) Mark OL poles/zeros. 2) Draw asymptotes. 3) Plot on real axis. 4) Find breakaway points. 5) Determine departure angles. 6) Sketch symmetric curves.
4.0 FREQUENCY DOMAIN ANALYSIS & STABILITY
4.1 Bode Plot (Asymptotic & Exact)
-
Magnitude Plot:
20 log₁₀ |G(jω)H(jω)|(dB) vslog ω. -
Phase Plot:
∠G(jω)H(jω)(degrees) vslog ω. -
Construction (Asymptotic):
-
Write
G(jω)H(jω)in standard form:K \frac{∏(1+jωT_z)}{∏(1+jωT_p)}. -
Magnitude:
-
Start at
20 log K(dB). -
For each pole at origin: slope
-20 dB/decstarting atω=1. -
For each zero at origin: slope
+20 dB/dec. -
For each simple pole at
ω = 1/T_p: slope changes by-20 dB/decatω = 1/T_p. -
For each simple zero at
ω = 1/T_z: slope changes by+20 dB/dec.
-
-
Phase:
-
Start:
0°(no zeros/poles at origin) or±90°n(n poles/zeros at origin). -
Each pole contributes
0° → -90°over 2 decades centered atω = 1/T_p. -
Each zero contributes
0° → +90°over 2 decades centered atω = 1/T_z.
-
-
-
Key Parameters:
-
Gain Crossover Frequency
ω_gc:|G(jω)H(jω)| = 1(0 dB). -
Phase Crossover Frequency
ω_pc:∠G(jω)H(jω) = -180°. -
Gain Margin (GM):
GM = 1/|G(jω_pc)H(jω_pc)|(linear) orGM = -20 log|G(jω_pc)H(jω_pc)|(dB). Positive GM → Stable. -
Phase Margin (PM):
PM = 180° + ∠G(jω_gc)H(jω_gc). Positive PM → Stable.
-
-
Stability Comment: For stable closed-loop, PM > 0 (typically > 30-50° for good transient response) and GM > 0 (typically > 6 dB).
[!TIP] Exact vs Asymptotic: Asymptotic is approximate. For accuracy, add correction factors at corner frequencies:
±3 dBfor simple pole/zero,±1 dBfor second-order factors.
4.2 Polar Plot (Nyquist Open-Loop)
-
Construction: Plot of
G(jω)H(jω)(real vs imaginary) asωvaries0 → ∞. -
Generalized Shapes:
| System Type | Polar Plot Shape (ω: 0→∞) | | :--- | :--- | | Type 0 | Starts on +real axis, ends at origin (0,0). | | Type 1 | Starts at ∞ on -90° line, ends at finite point on +real axis. | | Type 2 | Starts at ∞ on -180° line, ends at origin. |
-
Inverse Polar Plot: Plot of
1/[G(jω)H(jω)](used in Nyquist criterion). -
Comparison with Bode: Bode separates magnitude/phase; Polar combines them in one plot. Polar shows encirclements directly.
4.3 Nyquist Stability Criterion
-
Nyquist Contour: Semi-circle in RHP (
s = σ + jω,σfrom-∞to+∞,ωfrom+∞to-∞). -
Mapping Principle: Contour mapped by
F(s) = 1 + G(s)H(s). -
Encirclement
N: Number of clockwise encirclements of-1 + j0point by Nyquist plot ofG(s)H(s). -
Poles in RHP
P: Number of open-loop unstable poles. -
Closed-loop Unstable Poles
Z:\boxed{Z = P - N}
- Closed-loop stable iff
Z = 0(no poles in RHP).
- Closed-loop stable iff
-
Application:
-
Plot Nyquist of
G(s)H(s). -
Count
P(poles ofG(s)H(s)in RHP). -
Count
N(clockwise encirclements of -1 by Nyquist plot). -
Compute
Z. IfZ=0, closed-loop stable.
-
-
Gain Margin from Nyquist:
GM = 1/|G(jω_pc)H(jω_pc)|whereω_pcis frequency where plot crosses negative real axis. Distance from crossing point to -1.
[!TIP] Key:
Nis clockwise positive. ForP=0(stable open-loop),Nmust be0for stability (no encirclement of -1). ForP>0,Nmust equalPfor stability.
4.4 Gain Margin & Phase Margin
-
Gain Margin (GM): Factor by which gain can be multiplied before system becomes unstable.
GM = 1/|G(jω_pc)H(jω_pc)|. High GM → Less sensitive to gain changes. -
Phase Margin (PM): Additional phase lag required to reach -180° at
ω_gc.PM = 180° + ∠G(jω_gc)H(jω_gc). High PM → Less oscillatory, more damped response. -
Relationship: Both measure relative stability (how close to instability). Empirical:
PM ≈ ζ(for low ζ).PM > 30°usually desired.
5.0 STEADY-STATE ERROR ANALYSIS
5.1 Static Error Coefficients
- Position Error Constant
K_p:
$$K_p = \lim_{s→0} G(s)H(s)$$
- Velocity Error Constant
K_v:
$$K_v = \lim_{s→0} s G(s)H(s)$$
- Acceleration Error Constant
K_a:
$$K_a = \lim_{s→0} s^2 G(s)H(s)$$
- Derivation: From
E(s) = \frac{R(s)}{1 + G(s)H(s)}, apply Final Value Theorem (FVT) for standard inputs.
5.2 Steady-State Error e_ss
-
Type of System: Number of pure integrators (poles at origin) in
G(s)(open-loop).| System Type | Input |
e_ss(Unit Feedback) | | :--- | :--- | :--- | | Type 0 | StepA|\frac{A}{1+K_p}| | | RampAt|∞| | | ParabolicAt²/2|∞| | Type 1 | Step |0| | | Ramp |\frac{A}{K_v}| | | Parabolic |∞| | Type 2 | Step |0| | | Ramp |0| | | Parabolic |\frac{A}{K_a}| -
Generalized Error Coefficients:
$$e_{ss} = \frac{1}{1 + K_p} \text{ (step)}; \quad e_{ss} = \frac{1}{K_v} \text{ (ramp)}; \quad e_{ss} = \frac{1}{K_a} \text{ (parabolic)}$$
5.3 Limitations of Static Error Coefficients
-
Only valid for stable Type 0 & Type 1 systems.
-
Does not predict
e_ssfor unstable Type ≥1 systems (FVT not applicable). -
Inaccurate for inputs other than standard step/ramp/parabolic.
-
Does not give transient information.
[!TIP] Remember:
K_p,K_v,K_adepend only on open-loop TFG(s)H(s), not on closed-loop stability. Always check stability first (Routh/Root Locus).
6.0 CONTROLLERS & COMPENSATION DESIGN
6.1 Basic Controllers
| Controller | Transfer Function G_c(s) |
Effect on System |
|---|---|---|
| Proportional (P) | K_p |
Increases K_p, K_v, K_a proportionally. Reduces e_ss. No effect on system type. Can reduce stability if K_p too high. |
| Integral (I) | K_i/s |
Increases system type by 1. Eliminates e_ss for step (Type 0→1) or ramp (Type 1→2). Slows response, reduces PM. |
| Derivative (D) | K_d s |
Improves transient response (increases damping, reduces M_p, t_s). Increases high-frequency gain → noise amplification. |
| PID | K_p + K_i/s + K_d s = K_p(1 + 1/(T_i s) + T_d s) |
Combines benefits: reduces e_ss (via I), improves transient (via D). Tuning required. |
6.2 Need for Compensation
-
Improve transient response (increase PM,
ζ, reducet_s). -
Improve steady-state accuracy (increase
K_v,K_a). -
Reduce sensitivity to parameter variations.
-
Achieve controllability/observability in state-space.
6.3 Lead Compensation
- Phase Lead Network:
$$G_c(s) = \frac{1 + sT}{1 + αsT}, \quad 0 < α < 1$$
* Maximum phase lead: `φ_m = sin⁻¹\left(\frac{1-α}{1+α}\right)` at `ω_m = 1/(T√α)`.
* `α` determines `φ_m`; `T` sets `ω_m`.
-
Effect:
-
Increases Phase Margin (adds positive phase).
-
Increases bandwidth and speed (reduces
t_s). -
Increases
K_vslightly (if zero near origin).
-
-
Design Procedure (Bode):
-
From specs, determine required
PM_new = PM_required + safety (5-12°). -
Find
φ_max = PM_new - PM_current(from uncompensated Bode atω_gc). -
Determine
α = (1 - sin φ_max)/(1 + sin φ_max). -
Find
ω_m(where uncompensated phase =PM_new - φ_max - 90°? Actually:ω_gcshould be atω_mafter compensation. AdjustKso that|G_c(jω_m)G(jω_m)H(jω_m)| = 1). -
Calculate
T = 1/(ω_m √α). -
Final
Kchosen to meetK_vor gain crossover.
-
6.4 Lag Compensation
- Phase Lag Network:
$$G_c(s) = \frac{1 + sT}{1 + βsT}, \quad β > 1$$
* Maximum phase lag: `φ_m ≈ - (1 - α)/(1 + α)`? Actually for lag: `α = 1/β`, `φ_m = sin⁻¹\left(\frac{1-α}{1+α}\right)` but negative.
* `β` determines attenuation (and thus `K_v` improvement).
-
Effect:
-
Increases system type → Reduces
e_ss(improvesK_v). -
Slightly reduces PM (adds negative phase at
ω_gc). -
Lowers bandwidth.
-
-
Design Procedure (Bode):
-
Determine
KfromK_vrequirement:K_v = lim_{s→0} s G(s) = K(for Type 1). -
Draw Bode with this
K. CheckPM. IfPMinsufficient, design lead first, then add lag to adjustK_vwithout affectingPMmuch (place lag pole/zero near origin,ω_gc>>1/T).
-
6.5 Lag-Lead Compensation
-
Combined Network:
G_c(s) = G_lead(s) G_lag(s). -
When Used: Need simultaneous improvement in both transient (PM) and steady-state (
K_v). -
Design Steps:
-
Design lead to meet
PMspec (as in 6.3). -
Check
K_vfrom resultingK. If insufficient, design lag to increaseKwithout disturbingPM(place lag corner freq<< ω_gc). -
Alternatively, design simultaneously using Nichols or Bode.
-
6.6 Design Examples (Pattern)
-
Given:
G(s), specs:K_v = ?,PM ≥ ?°. -
Lead Design:
-
Find
KfromK_v(if Type 1:K_v = K). -
Plot Bode with
K. FindPM_current. -
φ_needed = PM_req + 5° - PM_current. -
α = (1 - sin φ_needed)/(1 + sin φ_needed). -
Find
ω_gc_desiredwhere uncompensated phase =-180° + PM_req + φ_needed. -
T = 1/(ω_gc_desired √α). -
Verify
PMat newω_gc.
-
-
Lag Design (if needed after lead):
-
β = K_v_required / K_v_current. -
Place
ω_z = ω_gc/10,ω_p = ω_z/β(so that|G_c| ≈ 1atω_gc).
-
7.0 STATE SPACE ANALYSIS
7.1 State Space Representation
-
State Variables: Minimum set
x₁(t), x₂(t), ..., xₙ(t)that completely describes system dynamics. -
State Vector:
x(t) = [x₁ x₂ ... xₙ]ᵀ. -
State Equations:
$$\dot{x}(t) = A x(t) + B u(t)$$
$$y(t) = C x(t) + D u(t)$$
Where `A` (n×n), `B` (n×r), `C` (m×n), `D` (m×r).
-
Modeling from Transfer Function:
-
Controllable Canonical Form: Companion matrix
A. -
Observable Canonical Form: Transpose of controllable form
Aᵀ.
-
-
Transfer Function Decomposition:
-
Direct Decomposition:
G(s) = C(sI - A)⁻¹B + D. -
Parallel Decomposition:
G(s) = Σ_{i=1}^n \frac{r_i}{s - p_i}(residuesr_ifrom partial fraction of(sI-A)⁻¹B).
-
7.2 Controllability & Observability
-
Controllability: Ability to transfer state from any initial to any final in finite time with suitable input.
-
Kalman's Rank Condition:
\boxed{\text{Rank } \mathcal{C} = n}
Where Controllability Matrix:
-
$$\mathcal{C} = [B \quad AB \quad A^2B \quad ... \quad A^{n-1}B]$$
-
Observability: Ability to determine initial state from output history.
-
Kalman's Rank Condition:
\boxed{\text{Rank } \mathcal{O} = n}
Where Observability Matrix:
-
$$\mathcal{O} = \begin{bmatrix} C \\ CA \\ CA^2 \\ \vdots \\ CA^{n-1} \end{bmatrix}$$
- Duality:
(A,B)controllable ↔(Aᵀ,Cᵀ)observable.
7.3 State Transition Matrix Φ(t)
-
Definition:
Φ(t) = e^{At}. -
Solution of State Equation (zero-input + zero-state):
$$x(t) = Φ(t)x(0) + \int_0^t Φ(t-τ) B u(τ) dτ$$
-
Properties:
-
Φ(0) = I -
Φ(t₁)Φ(t₂) = Φ(t₁ + t₂) -
Φ(-t) = [Φ(t)]⁻¹ -
dΦ(t)/dt = A Φ(t) = Φ(t) A -
Φ(t) = \mathcal{L}^{-1}[(sI - A)^{-1}]
-
-
Calculation:
-
Laplace Inverse:
Φ(t) = \mathcal{L}^{-1}[(sI - A)^{-1}]. -
Cayley-Hamilton: For 2×2/3×3, find
e^{At} = α₀(t)I + α₁(t)A + α₂(t)A²using eigenvalues.
-
7.4 Eigenvalues & Eigenvectors
-
Definitions: For matrix
A, scalarλand vectorvsatisfyAv = λv. -
Significance:
-
Eigenvalues
λ_i= closed-loop poles (forAmatrix of state-space). -
Determine stability (Re(λ) < 0 for all
i). -
Determine natural response modes (
e^{λ_i t}).
-
-
Calculation: Solve
det(A - λI) = 0forλ. For eachλ, solve(A - λI)v = 0forv. -
Modal Matrix
V: Matrix of eigenvectors[v₁ v₂ ... vₙ]. -
Diagonalization: If
Vinvertible,A = V Λ V⁻¹, whereΛ = diag(λ₁, λ₂, ..., λₙ). -
Relation to
Φ(t):
$$e^{At} = V e^{Λt} V^{-1} = V \begin{bmatrix} e^{λ₁t} & & \\ & \ddots & \\ & & e^{λₙt} \end{bmatrix} V^{-1}$$
[!TIP] Controllability/Observability Check: For
Ain controllable canonical form, always controllable. For observable canonical form, always observable.
8.0 ADVANCED TOPICS & SHORT NOTES
8.1 Linearization of Nonlinear Systems
- Concept: Approximate nonlinear system
ẋ = f(x,u)around an operating point(x₀, u₀)using Taylor series (ignore higher-order terms).
$$\delta\dot{x} = \frac{∂f}{∂x}\bigg|_{(x₀,u₀)} \delta x + \frac{∂f}{∂u}\bigg|_{(x₀,u₀)} \delta u$$
Where `δx = x - x₀`, `δu = u - u₀`.
- Effect of Feedback: Feedback linearization uses feedback to cancel nonlinearities, making overall system linear.
8.2 Effect of Pole Location on Stability & Response
-
Real Poles (
σ < 0):-
Left of origin → stable.
-
More negative (
σlarge magnitude) → faster response. -
Dominant pole (closest to origin) dictates response speed.
-
-
Complex Poles (
σ ± jω_d):-
σdetermines decay rate (stability margin). -
ω_ddetermines oscillation frequency. -
ζ = -σ/√(σ²+ω_d²)determines damping.
-
-
Pole-Zero Cancellation: Cancelling a pole with a zero can remove a mode but risks instability if cancelled pole is unstable or lightly damped.
8.3 Log Magnitude vs. Phase Plot
-
Plot of
20 log|G(jω)H(jω)|vs∠G(jω)H(jω)(both in dB/degrees). -
Relation to Nichols Chart: Nichols chart is a standardized version showing constant
M_p,ω_rcontours. -
Use: Directly read gain/phase margins from intersection with 0 dB and -180° lines. Assess stability margins visually.
8.4 Relative Stability
-
Beyond Absolute Stability: Measures how close system is to instability.
-
Quantification:
-
Gain Margin (GM): Factor of safety against gain increase.
-
Phase Margin (PM): Factor of safety against phase lag.
-
Delay Margin: Maximum time delay before instability.
-
-
Higher GM/PM → More Robustly Stable.
8.5 Comparative Analysis of Plots
| Plot | Advantages | Information |
|---|---|---|
| Bode | Handles magnitude/phase separately; easy to construct for complex systems; shows bandwidth. | ω_gc, PM, GM, slope changes. |
| Polar | Single plot; shows encirclements directly; good for gain/phase variations. | Overall shape, encirclements, GM (distance to -1). |
| Nyquist | Complete mapping of contour; handles unstable open-loop systems; rigorous stability test. | N, P, Z; GM from crossing point. |
[!TIP] Quick Recall: Bode = dB vs log ω (separate); Polar = Re vs Im (single curve); Nyquist = mapping of entire contour (includes infinite semicircle).
END OF UNIT 4 NOTES
Aligned with RGPV past papers (JUN 2025 - JUN 2022). Focus on derivations (overshoot, servomotor), rule applications (Routh, Root Locus, Mason), and design (Bode-based compensator).