Nonlinear observability algorithms with known and unknown inputs: analysis and implementation

TL;DR

This paper compares two symbolic algorithms, ORC-DF and FISPO, for analyzing the observability of nonlinear systems with known and unknown inputs, highlighting their differences, strengths, and limitations, and implementing them in a MATLAB toolbox.

Contribution

It provides an analytical comparison of ORC-DF and FISPO algorithms, and integrates both into the STRIKE-GOLDD toolbox for easier application.

Findings

FISPO is more generally applicable to nonlinear models.

ORC-DF can be more efficient for affine models with known inputs.

The choice of algorithm depends on model characteristics.

Abstract

The observability of a dynamical system is affected by the presence of external inputs, either known (such as control actions) or unknown (disturbances). Inputs of unknown magnitude are especially detrimental for observability, and they also complicate its analysis. Hence the availability of computational tools capable of analysing the observability of nonlinear systems with unknown inputs has been limited until lately. Two symbolic algorithms based on differential geometry, ORC-DF and FISPO, have been recently proposed for this task, but their critical analysis and comparison is still lacking. Here we perform an analytical comparison of both algorithms and evaluate their performance on a set of problems, discussing their strengths and limitations. Additionally, we use these analyses to provide insights about certain aspects of the relationship between inputs and observability. We find that, while ORC-DF and FISPO follow a similar approach, they differ in key aspects that can have a substantial influence on their applicability and computational cost. The FISPO algorithm is more generally applicable, since it can analyse any nonlinear ODE model. The ORC-DF algorithm analyses models that are affine in the inputs, and if those models have known inputs it is sometimes more efficient. Thus, the optimal choice of a method depends on the characteristics of the problem under consideration. To facilitate the use of both algorithms we implement the ORC-DF algorithm in a new version of STRIKE-GOLDD, a MATLAB toolbox for structural identifiability and observability analysis. Since this software tool already had an implementation of the FISPO algorithm, the new release allows modellers and model users the convenience of choosing between different algorithms in a single tool, without changing the coding of their model.