---
title: Robust PI-Lead Control Design
url: https://www.emergentmind.com/topics/robust-pi-lead-control
type: topic
---

# Robust PI-Lead Control Design

Searching arXiv for the cited papers to ground the article in current metadata.
{"query":"id:2511.21641 OR 2511.21641 robust PI-Lead control", "max_results": 5}
{"query":"2511.21641", "max_results": 5}
to=arxiv_search 招商总代ीjson
{"query":"2511.21641", "max_results": 5}
Robust PI-Lead control, in the literature summarized here, denotes a family of feedback designs that combine a proportional-integral compensator with a lead element in order to obtain adequate bandwidth, phase margin and robustness under model uncertainty or limited plant knowledge. In the linear model-free formulation, the target controller is \(C(s)=K_p\,\frac{T_i s+1}{T_i s}\) together with a standard lead element \(L(s)=\frac{\tau s+1}{\alpha \tau s+1}\), \(0<\alpha<1\), chosen from experimental observations only [2511.21641]. Related work treats the lead stage as a nonlinear reset element for precision mechatronics, derives discrete PI-Lead laws from incremental nonlinear dynamic inversion and time-delay control for flight-control systems, and embeds PI-Lead synthesis in robust PID parameter-space constructions [1805.09703].

## 1. Canonical form and operating assumptions

In the model-free formulation, the plant is assumed to be type–one integrating, open-loop stable, minimum-phase, and possibly subject to an unknown pure delay \(e^{-s\tau}\). The measured signals are a real-time measurable input \(u(t)\) and a noisy output \(x(t)\), and the plant is written as
\[
G(s)=\frac{x(s)}{u(s)}=\frac{\tilde G(s)}{s},
\]
with \(\tilde G(s)\) stable and minimum-phase [2511.21641]. This restriction is substantive: a common misconception is that model-free PI-Lead tuning is assumption-free, whereas the stated procedure is explicitly developed for this plant class.

The same source frames the controller design around the loop \(L(s)\,C(s)\,G(s)\), where the PI part provides the single-integrator structure and the lead part supplies additional phase margin. The design objective is not stated as exact pole placement but as a practical loop-shaping goal: adequate bandwidth, phase margin and robustness obtained experimentally, without internal model or identification. A plausible implication is that the method is intended for settings in which plant experiments are easier to obtain than a reliable parametric model, particularly motion-control systems.

A broader interpretation of robust PI-Lead control appears in related literature. In precision mechatronics, the lead stage can be replaced by a nonlinear reset element rather than a linear differentiator, while in flight control an incremental law can be mapped into a discrete PI-Lead increment. In robust PID theory, PI-Lead is treated as a constrained subfamily of PID controllers. This suggests that “robust PI-Lead control” is not a single synthesis doctrine but a structural class realized through several design paradigms.

## 2. Ultimate-sensitivity tuning without modeling

The model-free design in “Model-free practical PI-Lead control design by ultimate sensitivity principle” follows the ultimate sensitivity principle in Ziegler–Nichols style: close the loop with a pure gain \(K\), increase \(K\) until the closed-loop output exhibits sustained oscillations, denote the gain by \(K_u\), and denote the oscillation frequency by \(\omega_u\), or equivalently measure \(T_u=2\pi/\omega_u\) [2511.21641]. The proposed PI-Lead synthesis then proceeds in three steps.

The first step determines the integrator time constant \(T_i\). One sets \(K_p=1\), together with any pre-gain \(k\) needed to stabilize the initial loop, applies a step reference, and gradually decreases \(T_i\). At the first occurrence of sustained oscillations, \(T_i=\bar T_i\), one records
\[
\bar\omega^c_{pi}=\frac{1}{\bar T_i},
\qquad
\bar\omega_{gc}=\frac{2\pi}{\text{(measured oscillation period)}},
\]
and then chooses
\[
T_i=\frac{10}{\max\{\bar\omega_{gc},\bar\omega^c_{pi}\}}.
\tag{2.1}
\]
The stated interpretation is that the factor \(10\) shifts the PI corner one decade below the critical phase-drop region.

The second step determines \(K_p\) from transient overshoot and therefore from phase margin. With \(T_i\) fixed, one varies \(K_p\) until the closed loop exhibits a transient overshoot
\[
M=\frac{\max_t\,x(t)-x_{\rm ref}}{x_{\rm ref}}
\]
in the range \(30\%\text{–}40\%\). The overshoot is related to the damping ratio \(\zeta\) by
\[
M=\exp\!\Bigl(-\frac{\pi\zeta}{\sqrt{1-\zeta^2}}\Bigr),
\qquad \zeta\in(0,1),
\tag{2.3}
\]
and the corresponding loop phase margin is
\[
\varphi_m
=\arctan\!\Bigl(
\frac{2\zeta}{\sqrt{\sqrt{1+4\zeta^4}-2\zeta^2}}
\Bigr).
\tag{2.4}
\]
Ensuring \(M\approx30\%\text{–}40\%\) therefore guarantees \(\varphi_m\approx30^\circ\text{–}40^\circ\).

The third step introduces the lead element. The method fixes \(\alpha=0.1\), corresponding to approximately \(55^\circ\) maximum lead, and places the lead’s maximum-phase frequency at half to one decade above the PI corner, with the convenient choice
\[
\omega_{\max\varphi}=\frac{10^{1.5}}{T_i}.
\]
This yields
\[
\tau=\frac{T_i}{10^{1.5}},
\qquad
\alpha\tau=\frac{T_i}{10^{0.5}},
\]
and the unity-DC-gain lead element
\[
L(s)=
\frac{\bigl(10^{1.5}/T_i\bigr)s+1}
{\bigl(10^{0.5}/T_i\bigr)s+1}.
\tag{2.5}
\]

The illustrative experiment uses a noise-perturbed voice-coil actuator with gravity and friction, sampled at \(10\ \text{kHz}\). The pre-gain is \(k=300\), the step amplitude is \(x_{\rm ref}=9\ \text{mm}\), sustained oscillations appear at \(\bar T_i=0.031\ \text{s}\), and the measured values are \(\bar\omega^c_{pi}=32.26\ \text{rad/s}\) and \(\bar\omega_{gc}=20.73\ \text{rad/s}\), giving
\[
T_i=\frac{10}{32.26}\approx0.31\ \text{s}.
\]
Overshoot tuning gives \(M\approx35\%\) and \(K_p\approx450\), so that
\[
C(s)=450\,\frac{0.31\,s+1}{0.31\,s}
=\frac{139.5\,s+450}{0.31\,s},
\]
and the lead element becomes
\[
L(s)=\frac{0.031\,s+1}{0.0031\,s+1}.
\]

## 3. Loop-shaping criteria, robustness margins, and empirical behavior

The loop-shaping formulation defines the loop transfer \(L_{\rm loop}(s)=L(s)\,C(s)\,G(s)\) and the standard sensitivity functions
\[
S(j\omega)=\frac{1}{1+L_{\rm loop}(j\omega)},
\qquad
T(j\omega)=\frac{L_{\rm loop}(j\omega)}{1+L_{\rm loop}(j\omega)}.
\]
The stated robustness targets are phase margin \(\varphi_m\ge 30^\circ\), gain margin at least \(6\ \text{dB}\), and sensitivity peak
\[
M_s=\max_\omega |S(j\omega)|\lesssim 2,
\]
with the latter linked to disturbance rejection and robust stability [2511.21641].

The same work gives practical notes that clarify the intended experimental workflow. If necessary, a static pre-gain \(k\) may be used, for example to counteract gravity in a motion system. In Step 1, the reference is stepped and self-oscillations are observed at low frequencies, where they stand well above sensor noise. In Step 2, percent overshoot \(M\) is measured, and the paper states that noise has little effect on peak detection. In Step 3, the lead element is implemented directly from (2.5). No internal model or identification is used at any stage.

The experimental comparison reported for the voice-coil actuator states that PI-Lead, relative to plain PI, reduced overshoot by approximately \(50\%\) and improved settling under external perturbations. Against a classic Ziegler–Nichols PID tuned via \(K_p=0.6K_u\), \(T_i=0.5T_u\), and \(T_d=0.125T_u\), the PI-Lead gave a smoother, less aggressive action and superior robustness at low frequencies. A plausible implication is that the added lead is being used primarily to recover phase margin that would otherwise be traded away by aggressive integral action.

## 4. Reset-lead realizations and mitigation of the waterbed effect

A different line of work develops a nonlinear reset-lead element for precision mechatronics, starting from the observation that industrial PID consists of Lag, Lead, and Low Pass Filters, and that linear PID is inherently bounded by the waterbed effect, which creates a trade-off between precision tracking on one side and stability robustness on the other [1805.09703]. The reset controller is partitioned into a resetting part \(\Sigma_r\) and a non-resetting part \(\Sigma_{nr}\) in series. For the reset-lead, the resetting block is chosen as
\[
\Sigma_r(s)=\frac{1}{s/(\alpha\omega_d)+1},
\]
with state-space realization
\[
A_r=-\alpha\omega_d,\qquad B_r=\alpha\omega_d,\qquad C_r=1,\qquad D_r=0,
\]
and reset law
\[
x_r(t^+)=\gamma\,x_r(t)\qquad \text{whenever } e(t)=0,
\]
where \(\gamma\in[-1,1]\) tunes the degree of nonlinearity. The non-resetting block is
\[
\Sigma_{nr}(s)=K_p\cdot\Bigl(1+\frac{\omega_i}{s}\Bigr)\cdot\Bigl(1+\frac{s}{\omega_d}\Bigr)\cdot\frac{1}{1+s/\omega_l}.
\]
The overall open loop is \(L(s)=\Sigma_r(s)\,\Sigma_{nr}(s)\,G_p(s)\), and the first-harmonic describing function of the reset element is
\[
G_{DF}(j\omega)=C_r (j\omega I-A_r)^{-1} B_r (I+j\Theta_D(\omega)).
\]

The design procedure is explicit. One identifies the plant \(G_p(s)\); for the Lorentz stage,
\[
G_p(s)=\frac{0.5474}{0.5718\,s^2+0.95\,s+1}.
\]
Then one selects desired gain crossover \(\omega_c\) and phase margin \(\phi_m\). Two cases are given: Case A with \(\omega_c=150\ \text{Hz}\) and \(\phi_m\approx45^\circ\), and Case B with \(\omega_c=200\ \text{Hz}\) and \(\phi_m\approx45^\circ\). The remaining placements are \(\omega_i=\omega_c/10\), \(\omega_d=\omega_c/5\), and \(\omega_l=7\omega_c\). The reset-lead replaces the linear lead, \(\gamma\) is chosen to give the same \(\phi_m\) from describing-function analysis, with \(\gamma=0\) as an example, and \(\alpha\approx0.7\) compensates the reset-induced corner shift. Finally, \(K_p\) is scaled so that \(|L(j\omega_c)|\approx1\).

The robustness rationale is stated in contrast to the linear “tamed differentiator”
\[
D_{\rm lin}(s)=\frac{s/\omega_d+1}{s/\omega_t+1}.
\]
In the linear case, adding phase near \(\omega_c\) requires \(\omega_t\gg \omega_d\), which also raises high-frequency gain and hurts precision. By contrast, the reset-lead contributes \(+52^\circ\) of phase for \(\gamma=0\) without boosting gain at high \(\omega\), permits \(\omega_t=\omega_d\), and recovers low-frequency gain for tracking. The measured improvements are reported in two ways. For Case A, with the same \(\omega_c=150\ \text{Hz}\) and \(\phi_m=45^\circ\), RMS tracking error is reduced by \(34\%\), low-frequency \(-3\ \text{dB}\) sensitivity improves, and high-frequency \(|L|\) is approximately \(5\ \text{dB}\) lower than for linear PID. For Case B, with the same \(\phi_m=45^\circ\) and the same high-frequency \(|L|\) as PID, \(\omega_c\) increases from \(150\) to \(200\ \text{Hz}\), a \(33\%\) bandwidth increase.

The time-domain results use a 4th-order triangular reference with \(400\ \text{nm}\) amplitude and velocity, acceleration, and jerk limits. Case A reports PID RMS error \(41.5\ \text{nm}\) versus Reset-PID A \(23.1\ \text{nm}\), and Case B reports PID RMS \(41.5\ \text{nm}\) versus Reset-PID B \(15.7\ \text{nm}\). With \(50\ \text{nm}\) white noise injected, PID maximum error is approximately \(100\ \text{nm}\), while Reset-PID A is approximately \(60\ \text{nm}\). Closed-loop identification shows that \(S(j\omega)\) and \(T(j\omega)\) closely match describing-function predictions. A common misconception is that the lead element in PID must be implemented as a linear differentiator; the reset-lead results show a distinct nonlinear alternative.

## 5. Incremental PI-Lead synthesis via INDI and time-delay control

For flight-control systems, PI-Lead appears through the equivalence between incremental nonlinear dynamic inversion (INDI) and time-delay control (TDC). The plant is written as
\[
\dot\omega=f(\omega)+G(\omega)\delta,
\]
where \(\omega\in\mathbb{R}^3\) are body rates and \(\delta\in\mathbb{R}^3\) are control deflections [1701.08981]. The one-step nonlinear dynamic inversion law is
\[
\delta=G(\omega)^{-1}\bigl[\nu-f(\omega)\bigr],
\]
with virtual input \(\nu\) chosen from desired error dynamics. Taking a small time delay \(\lambda\) and performing a first-order Taylor increment around \((\omega_0,\delta_0)=(\omega(t-\lambda),\delta(t-\lambda))\) gives
\[
\dot\omega(t)\simeq \dot\omega_0+G_0\bigl[\delta(t)-\delta_0\bigr],
\qquad
\Delta\dot\omega\simeq G_0\Delta\delta,
\]
and therefore the standard INDI law
\[
\delta(t)=\delta(t-\lambda)+G_0^{-1}\bigl[\nu(t)-\dot\omega(t-\lambda)\bigr].
\]
In the alternative extended form \(\dot\omega=H(t)+\bar g\,\delta\), time-delay estimation of \(H(t)\) yields
\[
\delta(t)=\delta(t-\lambda)+\bar g^{-1}\bigl[\nu(t)-\dot\omega(t-\lambda)\bigr],
\]
so INDI and TDC deliver the same incremental law.

With sampled time \(\lambda=t_s\), the discrete law is
\[
\delta(k)=\delta(k-1)+\bar g^{-1}\bigl[\nu(k-1)-\dot\omega(k-1)\bigr].
\]
Choosing second-order desired error dynamics with \(e:=\omega_d-\omega\) and \(\nu=\dot\omega_d+k_P e+k_D\dot e\) leads to the standard discrete PI-Lead increment
\[
\delta(k)=\delta(k-1)+K\,t_s\Bigl[T_D\ddot e(k-1)+\dot e(k-1)+T_I^{-1}e(k-1)\Bigr].
\]
Term-by-term matching gives
\[
K=\frac{k_D}{\bar g\,t_s},
\qquad
T_I=\frac{k_D}{k_P},
\qquad
T_D=\frac{1}{k_D}.
\]
For closed-loop tuning in terms of natural frequency and damping ratio, one chooses
\[
s^2+2\zeta\omega_n s+\omega_n^2=0,
\]
so that
\[
k_P=\omega_n^2,\qquad k_D=2\zeta\omega_n,
\]
and therefore
\[
K=\frac{2\zeta\omega_n}{\bar g\,t_s},
\qquad
T_I=\frac{2\zeta}{\omega_n},
\qquad
T_D=\frac{1}{2\zeta\omega_n}.
\]
The damping ratio may be mapped to phase margin through
\[
PM \approx \arctan\!\Bigl[\frac{2\zeta}{\sqrt{\sqrt{1+4\zeta^4}-2\zeta^2}}\Bigr].
\]

The paper’s longitudinal pitch-rate example controls \(x_2=q\) to follow \(q_{\rm ref}\) and uses desired one-pole error dynamics \(\dot z_2+k_{P2}z_2=0\) with \(k_{P2}=50\ \text{rad/s}\). The incremental law becomes
\[
u(k)=u(k-1)+\bar g^{-1}\bigl[-k_{P2} z_2(k-1)-\dot q(k-1)+\dot q_{\rm ref}(k-1)\bigr].
\]
The equivalent PI parameters are \(K=(\bar g\,t_s)^{-1}\) and \(T_I=k_{P2}^{-1}=0.02\ \text{s}\). For \(t_s=1\ \text{ms}\) and \(\bar g=\hat g_2\), the mapping yields \(K\approx100/\hat g_2\); with \(M=2\) and \(\hat g_2\approx1.0\), the simulation uses \(K=100\), \(T_I=0.02\ \text{s}\), white-noise \(\sigma=10^{-3}\ \text{rad/s}\), and reports coincident INDI and PI-Lead closed-loop responses. The reported Bode plot shows a \(3\ \text{dB}\) bandwidth of approximately \(8\ \text{Hz}\) and \(PM\approx60^\circ\), consistent with \(\zeta\approx0.7\).

## 6. Robust stabilization regions and PI-Lead as a constrained PID family

A more formal robust-control treatment is provided by robust PID theory for finite sets of SISO LTI plants \(G_\nu(s)=B_\nu(s)/A_\nu(s)\). The central result is that the set of all continuous-time PID controllers
\[
K(s)=k_P+\frac{k_I}{s}+k_D s
\]
that simultaneously stabilize all plants is exactly the union of convex polygonal slices in \((k_P,k_I,k_D)\)-space, each slice taken at fixed \(k_P\) [1303.0425]. If
\[
\mathcal S(k_P)=\{(k_I,k_D)\mid K(s)\text{ with gains }(k_P,k_I,k_D)\text{ is Hurwitz for every plant}\},
\]
then
\[
\mathcal S=\bigcup_{k_P}\{k_P\}\times \mathcal S(k_P),
\]
and each \(\mathcal S(k_P)\) is a convex polygon whose edges arise from singular frequencies on the imaginary axis.

The continuous-time construction uses
\[
p_\nu(s)=A_\nu(s)(k_I+k_P s+k_D s^2)+B_\nu(s).
\]
On \(s=j\omega\), with \(A_\nu(j\omega)=R_A(\omega)+jI_A(\omega)\) and \(B_\nu(j\omega)=R_B(\omega)+jI_B(\omega)\), the condition \(p_\nu(j\omega)=0\) decouples into
\[
H(\omega;k_P,k_I,k_D)
:=R_A(\omega)k_I-\omega^2R_A(\omega)k_D-\omega I_A(\omega)k_P+R_B(\omega)=0,
\]
\[
G(\omega;k_P,k_I,k_D)
:=I_A(\omega)k_I-\omega^2I_A(\omega)k_D+\omega R_A(\omega)k_P+I_B(\omega)=0.
\]
At a singular frequency \(\omega'\), one obtains the straight-line relation
\[
k_I-\omega'^2 k_D=-\frac{R_B(\omega')}{R_A(\omega')}
\]
in the \((k_I,k_D)\)-plane and the \(k_P\)-plot relation
\[
k_P=-\frac{I_B(\omega')}{\omega' R_A(\omega')}.
\]
The algorithm computes common stabilizing \(k_P\)-intervals, solves \(k_P(\omega)=k_{P_i}\) for singular frequencies, forms boundary lines, computes transition signs, intersects the associated half-spaces, and returns the robustly stabilizing region as a union of polygonal slices.

Within this framework, a PI-Lead controller is written as
\[
K(s)=\frac{k_I+(k_P+k_I T_L)s+k_P T_L s^2}{s(1+\alpha T_L s)},
\]
and embedded in standard PID form \(Q(s)=k_I+k_P s+k_D s^2\) through
\[
k_D=k_P T_L,\qquad k_P^*=k_P+k_I T_L,\qquad k_I^*=k_I.
\]
The PI-Lead family therefore satisfies the linear constraints
\[
k_D-T_L k_P=0,
\qquad
k_P^*-k_P-T_L k_I=0.
\]
The implementation proposed in the source is to compute the robust PID region \(\mathcal S\), impose these linear equations to reduce dimensionality to a one-parameter curve, and intersect that curve with \(\mathcal S\) to obtain admissible PI-Lead gains. A common misconception is that robust PI-Lead tuning is necessarily heuristic; this parameter-space construction shows a complete picture of all PI-Lead controllers that robustly stabilize the given plant family.

The illustrative example uses two second-order plants,
\[
G_1(s)=\frac{1}{s^2+2\zeta_1\omega_{n1}s+\omega_{n1}^2},
\qquad
G_2(s)=\frac{1}{s^2+2\zeta_2\omega_{n2}s+\omega_{n2}^2},
\]
with \((\zeta_1,\omega_{n1})=(0.5,1)\) and \((\zeta_2,\omega_{n2})=(0.7,2)\). The common stabilizing interval is reported as approximately \([0.5,3.0]\). Choosing \(T_L=0.5\) and \(k_P=1.5\) gives \(k_D=0.75\), and a permissible interval such as \(k_I\in[0.8,1.4]\) is obtained by enforcing the polygon inequalities. The Robsin toolbox is reported to implement singular-frequency detection, polygon construction, \(k_P\)-interval gridding, and robust region visualization.

Robust PI-Lead control therefore spans purely experimental tuning for type-one integrating plants, nonlinear reset realization for precision mechatronics, incremental discrete synthesis for flight control, and exact robust-stabilization region computation for plant families. This suggests a common structural theme—PI action combined with lead phase advance—while the notion of robustness varies across the literature, ranging from empirical phase-margin targets and low-frequency disturbance rejection to Lyapunov-based verification and simultaneous stabilization of multiple LTI plants.

Source: https://www.emergentmind.com/topics/robust-pi-lead-control