The Root Locus Interactive Calculator analyzes the path of closed-loop system poles as gain varies, enabling control engineers to design stable feedback systems with desired transient response characteristics. Root locus analysis is fundamental to classical control theory, providing visual insight into how system stability, damping, and natural frequency change with controller gain. This calculator supports both continuous-time and discrete-time systems, computing pole locations, stability margins, and performance metrics across the entire gain range.
📐 Browse all free engineering calculators
System Diagram
Root Locus Calculator
Root Locus Equations
Characteristic Equation
1 + KG(s)H(s) = 0
Where:
- K = proportional gain parameter (variable)
- G(s) = forward path transfer function
- H(s) = feedback path transfer function
- s = complex frequency variable (σ + jω)
Angle Condition
∠G(s)H(s) = (2k + 1)×180°, k = 0, ±1, ±2, ...
Or equivalently:
Σ(angles from zeros) - Σ(angles from poles) = (2k + 1)×180°
This condition determines which points in the s-plane lie on the root locus.
Breakaway/Break-in Points
dK/ds = 0
Equivalently:
N'(s)D(s) - N(s)D'(s) = 0
Where:
- N(s) = numerator polynomial of G(s)H(s)
- D(s) = denominator polynomial of G(s)H(s)
- Prime (′) indicates derivative with respect to s
Angle of Departure from Complex Pole
θd = 180° + Σφzi - Σφpi
Where:
- θd = angle of departure (degrees)
- φzi = angle from i-th zero to the pole in question
- φpi = angle from i-th other pole to the pole in question
Second-Order Dominant Poles
s = -ζωn ± jωn√(1 - ζ²)
Performance Relationships:
Ts = 4/(ζωn), Mp = 100×exp(-πζ/√(1-ζ²))
Where:
- ζ = damping ratio (dimensionless, 0 to 1 for underdamped)
- ωn = undamped natural frequency (rad/s)
- Ts = 2% settling time (seconds)
- Mp = percent overshoot (%)
Gain Margin
GMdB = 20 log10(Kcrit/Kcurrent)
Where:
- Kcrit = critical gain at stability boundary
- Kcurrent = current operating gain
- GMdB = gain margin in decibels
Theory & Engineering Applications
Root locus analysis, developed by Walter R. Evans in 1948, provides a graphical method for examining how the poles of a closed-loop transfer function move in the complex plane as a system parameter (typically gain K) varies from zero to infinity. Unlike frequency-domain methods such as Bode plots, root locus directly reveals transient response characteristics including natural frequency, damping ratio, settling time, and overshoot by mapping pole locations to time-domain performance metrics. This technique has become indispensable in classical control system design, particularly for single-input single-output (SISO) systems where proportional feedback dominates.
Fundamental Root Locus Properties
The root locus satisfies both magnitude and angle conditions derived from the characteristic equation 1 + KG(s)H(s) = 0. The angle condition—that the phase of G(s)H(s) must equal an odd multiple of 180 degrees—determines which points in the s-plane lie on the locus, independent of gain value. The magnitude condition then determines the specific gain value at any point on the locus. A non-obvious but crucial insight is that root locus segments on the real axis exist only to the left of an odd number of real-axis poles and zeros. This property, often overlooked in introductory treatments, directly follows from the angle contribution of real-axis singularities: each real pole or zero to the right of a test point contributes 0° while those to the left contribute ±180°.
The number of locus branches equals the number of open-loop poles, with each branch originating at an open-loop pole (when K = 0) and terminating either at an open-loop zero or at infinity along an asymptote. For systems with more poles than zeros (the typical case), the asymptotes radiate from the centroid σa = (Σpoles - Σzeros)/(n - m) at angles θk = (2k+1)×180°/(n-m) where n is the number of poles, m is the number of zeros, and k = 0, 1, 2, ..., (n-m-1). These asymptotic directions are precisely the angles at which the locus approaches infinity, providing critical information about high-gain behavior without detailed calculation.
Breakaway and Break-in Point Analysis
Breakaway points, where the root locus leaves the real axis, and break-in points, where branches return to the real axis, occur where dK/ds = 0. While textbooks present this condition formulaically, the physical interpretation reveals design insight: at breakaway, the system transitions from having all real poles (overdamped response) to having complex conjugate poles (underdamped oscillatory response). This transition point often represents the maximum real-part magnitude achievable, which directly determines the fastest possible settling time for a given system structure. For a system with transfer function K(s+z)/[(s+p₁)(s+p₂)], the breakaway point can be found by solving (p₁+p₂-2z) ± √[(p₁-p₂)² + 4z(z-p₁-p₂)] / 2 = 0, but designers should recognize that pole-zero configurations with breakaway points far into the left half-plane enable aggressive control without oscillation.
A practical limitation often ignored in academic treatments: breakaway point calculations become numerically unstable for higher-order systems (fifth order and above) because the derivative condition produces high-degree polynomials with closely-spaced or multiple roots. In industrial practice, control engineers typically use computational tools to plot the locus directly rather than solving for breakaway points analytically. Furthermore, not all solutions to dK/ds = 0 represent actual breakaway points—they must lie on a valid locus segment, satisfying the angle condition independently.
Worked Example: Proportional Control of a DC Motor Position System
Consider a permanent magnet DC motor with armature-controlled position feedback. The open-loop transfer function relating shaft angle to armature voltage is:
G(s) = 25 / [s(s + 4)(s + 6)]
We wish to design a proportional controller with gain K such that the closed-loop system achieves 16.3% overshoot (corresponding to ζ = 0.5) with the fastest possible settling time. The characteristic equation is:
1 + K × 25 / [s(s+4)(s+6)] = 0
s³ + 10s² + 24s + 25K = 0
Step 1: Determine Root Locus Structure
The system has three poles at s = 0, s = -4, and s = -6, with no finite zeros. The locus begins at these three poles. Since n - m = 3 - 0 = 3, three branches approach infinity along asymptotes at angles (2k+1)×180°/3 = 60°, 180°, and 300° (equivalently, -60°). The asymptote centroid is located at σa = (0 - 4 - 6 - 0)/(3 - 0) = -3.33.
Step 2: Find Breakaway Point
On the real axis between s = 0 and s = -4, the locus breaks away toward complex values. Setting dK/ds = 0:
From K = -[s(s+4)(s+6)]/25, taking the derivative:
dK/ds = -[(3s² + 20s + 24)]/25 = 0
Solving: 3s² + 20s + 24 = 0 yields s = [-20 ± √(400-288)]/6 = [-20 ± 10.58]/6
This gives s = -1.57 and s = -5.10. The breakaway point at s = -1.57 lies between the poles at 0 and -4, confirming it's valid. At s = -1.57, the gain is K = -[(-1.57)(-1.57+4)(-1.57+6)]/25 = 0.534.
Step 3: Locate Poles for ζ = 0.5
For a damping ratio of 0.5, the dominant poles must lie on rays from the origin at angle θ = arccos(0.5) = 60° from the negative real axis. The dominant pole locations are s = -σ ± jω where ω/σ = tan(60°) = √3, so ω = σ√3.
To find the intersection of the root locus with the ζ = 0.5 line, we check the angle condition. Testing s = -2.5 + j4.33 (magnitude 5.0 rad/s):
- Angle from pole at s = 0: arctan(4.33/-2.5) = 120°
- Angle from pole at s = -4: arctan[4.33/(-2.5+4)] = arctan(4.33/1.5) = 70.9°
- Angle from pole at s = -6: arctan[4.33/(-2.5+6)] = arctan(4.33/3.5) = 51.0°
- Total angle: 120° + 70.9° + 51.0° = 241.9° ≈ 240° (close to 180° + 60° = 240°)
This confirms s ≈ -2.5 ± j4.33 lies on the locus. The corresponding gain K is found from the magnitude condition:
K = |s(s+4)(s+6)|/25 = |(-2.5+j4.33)(1.5+j4.33)(3.5+j4.33)|/25
Computing magnitudes: |s| = 5.0, |s+4| = √(1.5² + 4.33²) = 4.58, |s+6| = √(3.5² + 4.33²) = 5.56
K = (5.0 × 4.58 × 5.56)/25 = 5.09
Step 4: Verify Performance Metrics
With ωn = 5.0 rad/s and ζ = 0.5:
- Damped natural frequency: ωd = ωn√(1-ζ²) = 5.0 × √0.75 = 4.33 rad/s ✓
- Settling time (2%): Ts = 4/(ζωn) = 4/(0.5 × 5.0) = 1.6 seconds
- Peak time: Tp = π/ωd = π/4.33 = 0.726 seconds
- Percent overshoot: Mp = 100 × exp(-πζ/√(1-ζ²)) = 100 × exp(-π×0.5/0.866) = 16.3% ✓
The third pole (non-dominant) at K = 5.09 is found from the characteristic equation s³ + 10s² + 24s + 127.25 = 0. Using the known roots s = -2.5 ± j4.33, the sum of all roots is -10, so the third pole is at s = -10 - (-2.5) - (-2.5) = -5.0. This third pole is fast compared to the dominant pair (settling in 0.8 seconds versus 1.6 seconds), validating the dominant pole approximation.
Design Conclusion: A proportional gain K = 5.09 achieves the target damping ratio with settling time of 1.6 seconds. Increasing gain further would move poles higher into the left half-plane, reducing settling time but creating excessive overshoot. This represents the optimal tradeoff for proportional control; improved performance requires adding lead compensation.
Discrete-Time Root Locus
For digital control systems, root locus analysis extends to the z-plane where the unit circle (|z| = 1) replaces the imaginary axis as the stability boundary. The discrete-time characteristic equation 1 + KG(z)H(z) = 0 produces a root locus with the same angle condition, but pole locations map differently to time-domain behavior. A continuous-time pole at s = σ + jω maps to z = e^(σT_s) × e^(jωT_s) where T_s is the sampling period. Constant damping ratio lines in the s-plane become logarithmic spirals in the z-plane, while constant natural frequency lines become radial lines from the origin. A practical complication: the bilinear transformation w = (z-1)/(z+1) maps the unit circle to the imaginary axis, allowing use of continuous-time root locus techniques, but pole/zero locations and gain values transform nonlinearly, often confusing less experienced designers.
Industrial Applications and Multi-Loop Systems
While root locus excels for SISO systems, modern control applications often involve multiple feedback loops where the technique extends through sequential loop closure. In cascaded control architectures (common in process industries), the inner loop is closed first, its characteristic equation becoming the "plant" for the outer loop's root locus analysis. For a two-loop system with inner loop transfer function K₁G₁(s)/(1+K₁G₁(s)), the outer loop root locus uses this closed inner-loop expression, revealing how inner-loop bandwidth affects outer-loop stability margins. Robotics applications frequently employ this structure: a fast inner torque loop stabilizes actuator dynamics while a slower outer position loop achieves trajectory tracking. A key limitation emerges in multi-input multi-output (MIMO) systems where multiple gains vary simultaneously—classical root locus cannot capture cross-coupling effects, necessitating state-space techniques or multi-variable Nyquist criteria.
Aerospace control systems leverage root locus for flutter suppression and autopilot design, where gain scheduling accommodates parameter variations across the flight envelope. The fundamental insight: plotting multiple root loci for different operating points (Mach number, altitude, angle of attack) reveals whether a single fixed-gain controller maintains adequate stability margins. In automotive electronic stability control, root locus guides selection of yaw rate feedback gain that balances responsiveness against noise sensitivity, with higher gains improving disturbance rejection but moving closed-loop poles closer to the imaginary axis, reducing phase margin. For more on related control applications, explore our collection of engineering calculators.
Practical Applications
Scenario: Tuning a Robotic Arm Controller
Marcus, a robotics engineer at an industrial automation company, is designing a position controller for a six-axis robotic arm used in precision assembly. The arm's joint actuators have been modeled with a third-order transfer function, and Marcus needs to select a proportional gain that achieves less than 5% overshoot while minimizing settling time. Using this root locus calculator in "Second-Order System Analysis" mode, he inputs the identified natural frequency (ω_n = 12.3 rad/s) and tests various damping ratios. By targeting ζ = 0.7 (corresponding to 4.6% overshoot), the calculator reveals the closed-loop pole locations at -8.61 ± j8.77 rad/s and predicts a settling time of 0.465 seconds. Marcus then verifies this gain value (K = 47.2) experimentally, finding that the arm reaches its target position in 0.48 seconds with 4.3% overshoot—close enough to the prediction to proceed with integration testing. The root locus analysis saved Marcus three days of trial-and-error tuning that would have been required using empirical methods alone.
Scenario: Analyzing Flight Control Stability Margins
Dr. Sarah Chen, a flight controls engineer at an aerospace company, is conducting stability margin analysis for a new unmanned aerial vehicle's pitch attitude controller. The baseline controller has been flight-tested at moderate gains, but program requirements demand verification of gain margin to ensure safe operation under sensor noise and aerodynamic uncertainties. Using the "Gain Margin & Stability" calculator mode, Sarah inputs the current gain (K = 8.5), crossover frequency (ω_pc = 4.7 rad/s), measured phase margin (52°), and the critical gain where the system becomes marginally stable (K_crit = 34.2, determined from wind tunnel frequency response tests). The calculator computes a gain margin of 12.05 dB, indicating the system can tolerate a four-fold gain increase before instability—well above the 6 dB minimum required by military specifications. This analysis convinces the certification authority that the controller has adequate robustness, allowing the program to proceed to autonomous flight testing without requiring expensive redesign of the control laws.
Scenario: Process Control Loop Optimization
Jennifer, a process control engineer at a chemical manufacturing plant, is tasked with improving the performance of a temperature control loop on a large reactor vessel. The existing PID controller exhibits slow response (settling time of 18 minutes) and occasional oscillations when disturbances occur. After performing a step test to identify the process dynamics, Jennifer obtains dominant pole estimates at -0.083 ± j0.042 rad/min. Using the "Dominant Pole Location" calculator mode, she enters these values along with the measured settling time and peak time. The calculator reveals that the damping ratio is only 0.41, explaining the oscillatory behavior, and the natural frequency is just 0.093 rad/min. To achieve a settling time under 12 minutes with acceptable damping (ζ = 0.6), Jennifer uses the calculator to determine that the dominant poles need to move to -0.167 ± j0.148 rad/min. This corresponds to increasing the controller gain by a factor of 2.3, which she implements carefully through a series of small adjustments while monitoring reactor stability. The optimized loop now settles in 10 minutes with no overshoot, increasing production throughput by 7% because batch cycles complete faster.
Frequently Asked Questions
Free Engineering Calculators
Explore our complete library of free engineering and physics calculators.
Browse All Calculators →🔗 Explore More Free Engineering Calculators
About the Author
Robbie Dickson — Chief Engineer & Founder, FIRGELLI Automations
Robbie Dickson brings over two decades of engineering expertise to FIRGELLI Automations. With a distinguished career at Rolls-Royce, BMW, and Ford, he has deep expertise in mechanical systems, actuator technology, and precision engineering.
