Recently, it has been observed that finite impulse response controllers are an excellent basis for encrypted control, where privacy-preserving controller evaluations via special cryptosystems are the main focus. Beneficial properties of FIR filters are also well-known from digital signal processing, which makes them preferable over infinite impulse response filters in many applications. Their appeal extends to feedback control, offering design flexibility grounded solely on output measurements. However, designing FIR controllers is challenging, which motivates this work. To address the design challenge, we initially show that FIR controller designs for linear systems can equivalently be stated as static or dynamic output feedback problems. After focusing on the existence of stabilizing FIR controllers for a given plant, we tailor two common design approaches for output feedback to the case of FIR controllers. Unfortunately, it will turn out that the FIR characteristics add further restrictions to the LMI-based approaches. Hence, we finally turn to designs building on non-convex optimization, which provide satisfactory results for a selection of benchmark systems.

Finite impulse response (FIR) filters are well-established and widely used in signal processing [1]. In contrast, FIR feedback controllers (with orders larger than \( 0\) ) are relatively rare. However, it has recently been observed that FIR controllers are quite useful in the framework of encrypted control [2]. This young but emerging field of research deals with privacy-preserving networked control despite honest-but-curious platforms (e.g., clouds) or neighboring agents (see [3, 4] for an overview). More precisely, techniques such as homomorphic encryption are used to realize encrypted controller evaluations. While encrypted controllers offer many promising features, they are computationally demanding and subject to technical restrictions. For instance, an unlimited operating time is challenging to ensure for controllers with infinite impulse response (IIR). In contrast, it is straightforward and efficient to encrypt and operate FIR controllers using state-of-the-art homomorphic encryption schemes [2]. Theoretical insights on this fact can be found in [5].

Regarding the design of FIR controllers (and filters), one can distinguish direct and indirect approaches. The latter typically build on IIR controllers (or filters) and consider FIR approximations, e.g., via frequency response sampling or windowing [6]. However, FIR approximations tend to have high orders, and they are not suited to describe unstable dynamics, which are well-known disadvantages. Nevertheless, an indirect approach has also been utilized to realize the encrypted FIR controllers in [2].

Direct FIR controller designs are rare and often restrictive. For instance, an \( H_{\infty}\) -based design using heuristics is proposed in [7] and a convexified design for linear systems with uncertainties is considered in [8]. In the context of encrypted control, [9] derives a sufficient condition for the existence of FIR controllers for SISO systems, which can be verified using linear programming.

Given the limitations of existing design approaches, especially for the MIMO systems considered here, this paper investigates whether FIR controllers may benefit from tailoring established designs for output feedback [10]. The origin of this approach stems from our observation that FIR controller designs can be equivalently formulated as static or dynamic output feedback problems (as specified in Sect. 2). Hence, classical design approaches based on (convex) linear matrix inequalities (LMI) [11, 12] can, in principle, be used and may offer problem-specific solutions. Unfortunately, our analysis (see Sect. 4) shows that the FIR specifications lead to further restrictions, which we formalize as a helpful reference. Given the predominantly adverse theoretical outcomes, we conclude with a numerical case study, which confirms our results but also shows that heuristic non-convex designs can provide satisfactory designs albeit with higher computational efforts.

Roadmap. Section 2 deals with the problem specification and reformulations of the FIR controller design. The existence of stabilizing FIR controllers is analyzed in Section 3. Then, convexification strategies for FIR designs are investigated in Section 4. Finally, we provide a numerical case study and conclusions in Sections 5 and 6, respectively.

We consider linear discrete-time systems of the form

with the time step \( k\in\mathbb{N}\) , system state \( \mathbf{x}(k) \in \mathbb{R}^n\) , control input \( \mathbf{u}(k) \in \mathbb{R}^m\) , and output \( \mathbf{y}(k)\in \mathbb{R}^p\) . Throughout the paper, we make the following assumptions.

Now, our focus is on the design of a stabilizing FIR feedback controller of the form

where \( \mathbf{F}_i \) denotes a controller gain and \( \ell \in\mathbb{N}\) refers to the FIR’s order. For \( \ell=0\) , the controller (2) results in classical static output feedback \( \mathbf{F}_0\mathbf{y}(k)\) . Remarkably, also the cases with \( \ell>0\) can be interpreted in terms of static output feedback. To this end and inspired by [7], we introduce the augmented states and outputs

which lead to the state space representation

where the augmented matrices are given by

Clearly, for the augmented system (3), the FIR controller (2) results in static output feedback of the form

Now, instead of augmenting the system, one can also reformulate (2) in terms of the linear dynamic controller

where \( \hat{\mathbf{x}}\) is the controller state and \( \mathbf{H}\) , \( \mathbf{G}\) , \( \mathbf{E}\) , and \( \mathbf{D}\) denote the controller parameters. In fact, specifying

renders (6) equivalent to (2). Both reformulations allow specifying our focus on stabilizing FIR controllers with a given order \( \ell\) . In fact, based on (3), we are aiming for \( \mathbf{\mathcal{F}}_\ell\) leading to Schur stable

Analogously, building on (6), we are looking for \( \mathbf{E}\) and \( \mathbf{D}\) resulting in a Schur stable matrix

Remarkably, the closed-loop matrices (8) and (9) are equal. We formalize this observation in the following lemma in order to simplify the discussion below.

Before we dive into possible design strategies for (2), we derive a necessary condition for the existence of a stabilizing FIR controller. This condition builds on the trivial observation that all eigenvalues of \( \mathbf{H}\) in (7) are zero. As a consequence, \( \mathbf{H}\) is Schur stable and the dynamic controller in (6) is (asymptotically) stable. Since this controller is an equivalent representation of (2), we conclude that the FIR controller design can only be successful if the system is stabilizable by a stable controller. Such systems are called strongly stabilizable and they possess the parity interlacing property (PIP) [13]. Now, for single-input single-output (SISO) systems, we find the following condition.

Variants of this theorem also exist for multiple-input multiple-output (MIMO) systems but, to the best of the authors’ knowledge, only for the continuous-time case [13, p. 84]. While a technical extension to the discrete-time MIMO case is beyond the scope of this paper, it is still clear by definition that strong stabilizability is a necessary condition for the existence of a stabilizing FIR controller. In order to highlight this point and for later validation, we consider the two SISO systems

which we specify via transfer functions for convenience. Taking into account that \( \infty\) counts as an unstable zero for strictly proper systems, it is easy to see that \( G_1\) possesses the PIP whereas \( G_2\) does not. Now, while only \( G_1\) is strongly stabilizable, classical stabilizability is given for both systems. For instance, the controllers

stabilize \( G_1\) and \( G_2\) , respectively. Here, \( C_1\) is a (stable) FIR controller, whereas \( C_2\) is neither FIR nor stable. In fact, due to Theorem 1, it is impossible to find a stabilizing (and by construction stable) FIR controller for \( G_2\) (as later confirmed by our numerical studies in Section 5).

Next, we motivate that strong stabilizability is not only necessary but also sufficient for the existence of stabilizing FIR controllers. In fact, strong stabilizability guarantees the existence of a stable dynamic controller of the form (6) but without the structural restrictions in (7). Nonetheless, it is well-known that such dynamics with a Schur stable \( \mathbf{H}\) can be approximated with arbitrary precision by FIR dynamics. To see this, note that the output of (6) for an initial controller state \( \hat{\mathbf{x}}(0)=\mathbf{0}\) is given by

Computing FIR approximations using a rectangular window method leads to the specifications \( \mathbf{F}_0 :=\mathbf{D}\) and \( \mathbf{F}_i := \mathbf{E} \mathbf{H}^{i-1} \mathbf{G}\) for every \( i \in \{1,\dots,\ell\}\) (see [6, Chapter 5.4]). Then, the difference between (12) and (2) is \( \sum_{i=\ell}^{k-1} \mathbf{E} \mathbf{H}^i \mathbf{G}\mathbf{y}(k-1-i)\) . Clearly, for sufficiently large FIR orders \( \ell\) , the difference tends to \( 0\) (and the same observation holds for \( \hat{\mathbf{x}}(0) \neq \mathbf{0}\) if also \( k\) is large enough). As a consequence, stable dynamic controllers suggest the existence of stabilizing FIR controllers. Nevertheless, in this paper, we aim for a direct design of FIR controllers and not for the approximation of stable IIR controllers, since this often requires unnecessarily large orders \( \ell\) . Still, large \( \ell\) can be beneficial for a direct design, as indicated by the following lemma. In fact, the feasibility of a controller design may only be obtained but never lost when increasing \( \ell\) .

In summary, designing stabilizing FIR controllers is only meaningful for strongly stabilizable systems and for these, it is reasonable to search for the smallest feasible FIR order \( \ell\) .

A stabilizing FIR controller is characterized by a Schur stable matrix \( \mathbf{\Phi}_\ell\) , as in (10). Clearly, based on Lyapunov’s theory, Schur stability holds if and only if there exists a positive definite \( \mathbf{\mathcal{P}}_\ell\) satisfying

Unfortunately, finding suitable \( \mathbf{\mathcal{F}}_\ell\) and \( \mathbf{\mathcal{P}}_\ell\) or even deciding whether they exist or not is significantly harder than solving the related problem

associated with static state feedback (i.e., \( \mathbf{u}(k)=\mathbf{\mathcal{K}}_\ell \tilde{\mathbf{x}}(k)\) ). In fact, the nonlinear coupling of \( \mathbf{\mathcal{K}}_\ell\) and \( \mathbf{\mathcal{P}}_\ell\) in (14) can be resolved by carrying out a congruence transformation with \( \mathbf{\mathcal{W}}_\ell:=\mathbf{\mathcal{P}}_\ell^{-1}\) , applying the Schur complement, and using the substitution \( \mathbf{\mathcal{L}}_\ell=\mathbf{\mathcal{K}}_\ell \mathbf{\mathcal{P}}_\ell^{-1}\) . Then, we obtain the LMI

which solves (14) by means of \( \mathbf{\mathcal{K}}_\ell=\mathbf{\mathcal{L}}_\ell\mathbf{\mathcal{W}}_\ell^{-1}\) . Clearly, if the resulting \( \mathbf{\mathcal{K}}_\ell\) is such that an \( \mathbf{\mathcal{F}}_\ell\) exists which satisfies

then \( \mathbf{\mathcal{P}}_\ell\) and \( \mathbf{\mathcal{F}}_\ell\) also solve (13). Remarkably, a decomposable solution for (14) as in (16) exists whenever (13) is solvable. This trivially follows from the fact that a solution to (13) in terms of \( \mathbf{\mathcal{P}}_\ell\) and \( \mathbf{\mathcal{F}}_\ell\) implies a solution to (14) in terms of \( \mathbf{\mathcal{P}}_\ell\) and \( \mathbf{\mathcal{K}}_\ell=\mathbf{\mathcal{F}}_\ell\mathbf{\mathcal{C}}_\ell\) , where \( \mathbf{\mathcal{K}}_\ell\) is decomposable by construction.

Now, it is interesting to note that (14) and consequently (15) are always feasible for our setup. In fact, the feasibility of (14) holds once the pair \( (\mathbf{\mathcal{A}}_\ell,\mathbf{\mathcal{B}}_\ell)\) is stabilizable, which is given according to the following lemma in combination with Assumption 1.

Given the guaranteed feasibility and convexity of (15), one could aim for solving (15) and trying to derive a solution to (13) via (16). This is expedient in special cases such as, e.g., \( \mathbf{C}=\mathbf{I}_n\) . However, in general, the ability to find a suitable \( \mathbf{\mathcal{F}}_\ell\) requires explicitly considering the constraint (16) while solving (14) or, equivalently, to solve (13) directly. Unfortunately, both problems are typically hard [15]. Hence, classical solution strategies build on different convexifications of (13) that enable an efficient numerical solution. Two popular ones are discussed next.

4.1 Convexification for static output feedback

4.2 Convexification for dynamic output feedback

The analysis in the previous sections showed that, on the one hand, the existence of a stabilizing FIR controller can easily be checked by testing whether the system is strongly stabilizable (see Section 3). On the other hand, since classical design approaches based on convex reformulations typically fail (see Section 4 and the examples studied below), directly designing a FIR controller appears to be hard. Yet, there exist various numerical tools tailored for the design of static output feedback which is compatible with a FIR controller via (5). This allows us to validate and discuss our results based on selected examples. For instance, the HiSyn toolbox [18] provides a useful toolkit. Moreover, we will consider genetic algorithms (GA), which also have been applied successfully to design output feedback [19] and even classical FIR filters [20, 21]. Specifically, we are using the GA instance from MATLAB’s global optimization toolbox. Recalling that we aim for stabilizing FIR controllers, we simply consider

as the goal for the numerical design, where \( \rho\left(\mathbf{\Phi}_\ell(\mathbf{\mathcal{F}}_{\ell})\right)\) denotes the spectral radius of \( \mathbf{\Phi}_\ell\) , depending on the choice of the controller parameters \( \mathbf{\mathcal{F}}_\ell\) . Note, in this context, that the HiSyn toolbox directly supports the design criterion (23). Still, minimizing the spectral radius is non-trivial. In fact, it has been pointed out in [22] that the non-smooth function \( \rho\) may even have a locally unbounded gradient with respect to \( \mathbf{\mathcal{F}}_\ell\) . Both characteristics can be handled by the HiSyn toolbox and GA. In fact, the HiSyn toolbox, while using gradient information, treats discontinuities carefully, and GA are gradient-free optimization algorithms by design.

We illustrate our findings and the numerical approaches with four examples. First, we pick up the two examples in (11) (from now on denoted as System 1 and 2, respectively). Then, we consider two examples from the literature. Specifically, the linearized batch reactor (System 3) in [23, p. 62] with a sampling period of \( 0.1\) is given by

which serves as a benchmark and is also studied in [2]). Lastly, we consider the system investigated for the FIR design in [7, p. 5] (System 4), with

The selection captures different cases relevant for our problem of interest. In fact, all systems but System 2 (i.e., \( G_2\) ) are strongly stabilizable. Systems 1 and 4 can be stabilized by classical static output feedback (i.e., a FIR controller with order \( \ell=0\) ), whereas this is impossible for System 2, and it seems to be not possible for System 3.

Now, for all examples and independent of the previous characterization, we repeat the following two steps: i) We investigate whether a convex FIR design is possible based on (18), where we note that considering \( \ell=0\) is sufficient according to Lemma 4; ii) We apply the HiSyn toolbox as well as GA to (approximately) solve (23) for increasing \( \ell \in \{0,1,\dots,5\}\) . Regarding step ii), accounting for the non-convexity of the problem, we perform \( 10\) optimization runs with randomly selected initial guesses for each \( \ell\) and each solver. In this context, if a feasible solution has been found for \( \ell-1\) , we use the warm start \( \mathbf{\mathcal{F}}_{\ell}:=\begin{pmatrix} \mathbf{\mathcal{F}}_{\ell-1}^\ast & \mathbf{0} \end{pmatrix}\) suggested by Lemma 2 as a seed for the random initialization.

The central results of our numerical analysis are illustrated in Figure 1 in terms of the optimized spectral radii resulting for the different systems, FIR orders, and solvers. We initially note that the results confirm that a stabilizing FIR controller cannot be found for System 2, as predicted by Theorem 1. Moreover, it is confirmed that Systems 1 and 4 can be stabilized by classical static output feedback since stabilizing FIR controllers are found even for \( \ell=0\) . Remarkably, a convex design of such controllers via (18) is only possible for System 4 but not for System 1, which highlights the restrictiveness underlying this convexification. At this point, it is interesting to note that the convex approach in [9], restricted to SISO systems, is likewise infeasible for System 1 (and 2). As expected, the LMI (18) is further infeasible for Systems 2 and 3 and each \( \ell \in \mathbb{N}\) . However, as apparent from Figure 1, a stabilizing FIR controller can be found for System 3 and \( \ell\geq 1\) by both the HiSyn toolbox and GA. Finally, one can make two observations. First, the spectral radius is typically decreasing with \( \ell\) approximately until \( \ell=n\) despite the additional design freedom (see Lemma 2). Second, the HiSyn toolbox provides better or equally good performance, while the running time is about \( 15\) times faster than GA due to the aforementioned use of gradients in the optimization procedure.

Motivated by applications in encrypted control, we considered the design of stabilizing FIR controllers in this paper. We found that the existence of such controllers is closely related to the notion of strong stabilizability (see Section 3). Regarding the design of FIR controllers, one can, in principle, make use of reformulations in terms of augmented static output feedback or structurally restricted dynamic output feedback (Section 2). However, classical convexifications based on these reformulations suffer from the FIR’s special structure, as shown in Section 4. Therefore, these design approaches are impractical (or even impossible) for a reliable controller synthesis. Nonetheless, designs can, e.g., be achieved by non-convex optimization algorithms as shown in Section 5.

Future research directions are twofold. First, we will investigate further classical design approaches in the light of FIR controllers. Second, applying FIR controllers in privacy-preserving (encrypted) control as low-level primitives is of interest.