Novel Direct and Accurate Identification of Kalman Filter for General Systems Described by a Box-Jenkins Model

A novel robust Kalman filter (KF)-based controller is proposed for a multivari-able system to accurately track a specified trajectory under unknown stochastic disturbance and measurement noise. The output is a sum of uncorrelated signal, disturbance and measurement noise. The system model is observable but not controllable while the signal one is controllable and observable. An emulator-based two-stage identification is employed to obtain a robust model needed to design the robust controller. The system and KF are identified and the signal and output error estimated. From the identified models, minimal realizations of the signal and KF, the disturbance model and whitening filter are obtained using balanced model reduction techniques. It is shown that the signal model is a transfer matrix relating the system output and the KF residual, and the residual is the whitened output error. The disturbance model is identified by inverse filtering. A feedback-feedforward controller is designed and implemented using an internal model of the reference driven by the error between the reference and the signal estimate, the feedforward of reference and output error. The successful evaluation of the proposed scheme on a simulated autonomously-guided drone gives ample encouragement to test it later, on a real one.


Introduction
In conventional Kalman filter applications, the system involved is typically linearized and then identified. Based on the identified system, the Kalman filter is then identified. In this chapter, we propose a novel approach in that (a) the system is represented by a more general model, termed multi-input multi-output Box-Jenkins (MIMO BJ model) which subsumes all previous classical models, such as ARMA models and their derivatives, and (b) the associated Kalman filter identification is carried out directly, i.e. it does not necessitate the prior identification of the system involved. The various tools involved in our proposed approach are all explained below. consistent, unbiased, and efficient. The fundamental requirement of identification is that the leftover signal from identification, namely the residual is a zero-mean white noise process that contains no information. To meet this requirement, both the input and output of the system are filtered. Among the class of all linear whitening filters, the Kalman filter is the best. The system model is indirectly identified by minimizing the residual generated by the Kalman filter instead of the equation error. The Subspace Method (SM) uses the structure of the state-space model of the Kalman filter, whereas the prediction error method (PEM), which is the gold standard in system identification, is developed from the residual model [1,[9][10][11].

Emulator-based two-stage identification
The static and dynamic behavior of a physical system change as a result of variations in the parameters of some of its subsystems such as sensors, actuators, plant, disturbance models, and controllers. As the parameters of these subsystems are not generally accessible to generate data, instead, emulators, which are hardware or software devices, are connected in cascade to the output, input or both, of the subsystems. An emulator is a transfer function block which mimics the variations in the associated subsystems including the disturbance model. An emulator takes the form of a static gain or an all-pass filter to induce gain or phase variations in the subsystem it is connected to. Emulator parameters are perturbed to mimic various normal and abnormal, or faulty, operating scenarios resulting from variations in these subsystems. The emulator-generated data is employed in (a) the identification of robust systems and signal models and their associated Kalman filters using the two-stage identification scheme [2,3,[6][7][8].
A two-stage identification is used in various applications including the nonparametric identification of impulse response, estimation of Markov parameters in the SM, in model predictive control, identification of a signal model and in system identification. The use of the two-stage identification is inspired by the seminal paper by [12] for an accurate estimation of the parameters of an impulse response from measurements in an additive white noise. It is shown via simulation that the variance of the parameter estimation error approaches the Cramer-Rao lower bound [13]. Further, it is shown analytically that using a high-order model (with an order several times larger than the true order) improves significantly the accuracy of the parameter estimates. The two-stage scheme has not received much attention in system identification although it has been mentioned as an alternative scheme to the PEM [1,14], and has been successfully employed in identification in [15][16][17].
It should be emphasized that the prediction error method (PEM), viewed as a gold standard for system identification, is not geared for the estimation of the signal buried in the output, i.e. it is developed for the ARMA model and not for the Box-Jenkins one. A two-stage identification of the Box-Jenkins model is proposed as the system model is observable but not controllable while the signal model is both controllable and observable: • In the first stage, the robust system model and the associated Kalman filter are identified using the emulator-generated data using PEM, and the signal and the output error are both estimated. Further, the whitening filter that relates the output error and the residual is obtained.
The high-order for the first stage and the reduced-order for the second are both selected using the Akaike information criterion (AIC) and are cross-checked by verifying the whiteness of the associated residual. The two-stage identification has also been successfully employed in the identification model.
The question arises as to how to obtain the system model and the signal model from the identified high-order model Kalman filters. A key property of the Kalman filter is established here, namely that the transfer matrix of the signal and the system is the matrix fraction description model derived from the Kalman filter residual model of the system. This property is exploited to derive the signal and the system transfer matrices. The state-space models of the signal and the system models are derived from the identified state-space models of the Kalman filters. Thus, the proposed scheme identifies (a) the Kalman filter for the system, (b) the Kalman filter for the signal first. Then, the system model and the signal model are separately obtained.
The proposed scheme is further extended to identify the signal model to complement the PEM. In the first stage, a very high-order model is identified using PEM. In the second stage, the signal model is identified using a balanced model reduction of the high-order identified model obtained in the first stage. The PEM and state-space method (SM) are both tailored to identify the signal model and estimate the signal by employing the proposed version of the two-stage identification scheme. The results of the comparison of the performance of these methods in identifying the system and signal models are presented.

Highlights of the contributions
• The Auto-Regressive (AR), The Moving Average (MA), and the Auto-Regressive and Moving Average (ARMA) models are all special cases of the proposed Box-Jenkins model. As this model is more general and hence has wider applications, including robust controller design; estimation of latent variables; monitoring of the status of the system, fault diagnosis, development of condition-based maintenance programs and design of faulttolerant systems; filtering of signals, speech enhancement, noise and echo cancelation in communication; 2-D image filtering and tracking of moving objects.
• The state-space models of the system and the signal models are derived from the identified Kalman filters, by invoking the (causal) invertibility of the output error and the residual [5].
• An efficient scheme to monitor the status of the system may be implemented from the proposed scheme. First the status of the system is monitored by analyzing the residual of the Kalman filter of the system model. If there is a variation, then the residual of the Kalman filter of the signal model is analyzed to ascertain whether a fault has occurred.
• In practice, disturbances are inevitable, and can negatively affect the system performance. When the system is in an abnormal state, it is not in general easy to determine whether the abnormal operation is the result of variations in the disturbance or the occurrence of a fault. The proposed scheme provides a simple solution by analyzing the residuals of both Kalman filters as the residual of the Kalman filter for the system captures the variations in both the system and disturbance models, while that of the Kalman filter for the signal, and captures only the variations in the signal model. This is crucial for reducing false alarm, and all its concomitant risks and costs, resulting from variations in the disturbance and not in the signal model [19].
• The PEM (SM) may be tailored to identify the signal model and estimate the signal itself by using the proposed two-stage identification scheme.

Applications
Applications include monitoring the status of the system and the signal models, distinguishing between the variations in the disturbance model and those in the signal model to help diagnose a fault in the system and ensure a low false alarm probability, estimating the latent variable, namely the signal, developing a framework for applications including robust controller design; fault diagnosis; speech and biological signal processing; tracking of moving objects, design of soft sensors to replace maintenance-prone hardware sensors, evaluate and monitor product quality, meeting the ever-increasing need for fault-tolerant systems for mission-critical systems found in aerospace, the nuclear power systems, and autonomous vehicles.

Problem formulation
The output y k ð Þ ∈ R q is an additive sum of the signal s k ð Þ ∈ R q , disturbance, d k ð Þ ∈ R q and the measurement noise v k ð Þ ∈ R q where R is real scalar field.
Where the signal and the disturbance models are: Where u k ð Þ ∈ R q is the input; w k ð Þ ∈ R p is zero-mean white noise process that generates the disturbance d k ð Þ ∈ R p , and is uncorrelated with the measurement are qxp transfer matrix of order n s and n w respectively; ϑ k ð Þ¼d k ð Þ þ v k ð Þ is the output error. The signal model G s z ð Þ is formed of cascade and parallel combinations of the subsystems such as actuators, plant and the sensors. Let the state space model of the signal and the disturbance models be respectively Figure 1 shows the input-output model relating the input, the signal model, the signal, the disturbance model, the disturbance, the measurement noise and the output.
Linear regression model: System: signal, the disturbance and the measurement noise.
Using, (1)-(3), the expression for the linear regression becomes: are respectively the denominator and numerator polynomials. The model is termed Box-Jenkins model. Note that the model that generates the equation error υ z ð Þ is a Moving Average (MA) model, whereas the one that generates the output error ϑ k ð Þ is an Auto-Regressive Moving-Average (ARMA) model.
Augmented state-space model: The augmented state-space representation of the multi-input and multi-output ( Þ representing a p-input, q-output system, is given by: Ã T ∈ R q are respectively the state, the input and output; n ¼ n s þ n w is the order; w k ð Þ ∈ R n and v k ð Þ ∈ R q are respectively the disturbances and measurement noise; The augmented transfer matrix is not a minimal realization of the system output model as there is (stable) pole-zero cancelation since the polynomial D w z ð Þ, which is common to both the numerator N sw z ð Þ and the denominator D sw z ð Þ. In other words, N sw z ð Þ and D sw z ð Þ are not coprime. The signal model . Assumptions: It is assumed that (a) the disturbance w k ð Þ and the measurement noisev k ð Þ are independent zero-mean Gaussian white noise processes with unknown but finite covariance, , respectively, and are inaccessible, (b) A; C ð Þis observable, (c) the signal and disturbance models are both minimal, A s ; B s ; C s ; D s ð Þand A w ; B w ; C w ; D w ð Þ are both controllable and observable, (d) The initial conditions x 0 ð Þ, w k ð Þ and v k ð Þ are mutually uncorrelated. However, the signal s z ð Þ and the disturbance w z ð Þ may have spectral overlap, (e) the output error is bounded.

Kalman filter
Predictor form: A robust Kalman filter of the identified system A 0 ; B 0 ; C 0 ; D 0 À Á relating the system input u z ð Þ and system output y k ð Þ to the estimated output y k ð Þ is:x the residual or the innovation sequence; the Kalman gain K 0 ∈ R nxq ensures the asymptotic stability of the Kalman filter, i.e. (A 0 À K 0 C 0 ) is strictly Hurwitz having all its eigenvalues strictly inside the unit circle.
Innovation form: There is duality between the predictor, and the innovation forms of the Kalman filter [5]. The output y k ð Þ and the residual e k ð Þ are (causally) invertible. In other words, e k ð Þ can be generated from the output y k ð Þ and r k ð Þ using the (causal) predictor form, and y k ð Þ can be generated from e k ð Þ and r k ð Þ using the innovation form. The Kalman filter given by (7) is termed the predictor form and can be expressed in an alternative form, termed the innovation form, given by: Figure 2 shows the system and the Kalman filter which embodies the system model A; B; C ð Þ . The inputs to the Kalman filter are the input r k ð Þ and the output y k ð Þwhich is corrupted by the noise v k ð Þ and affected by the disturbance w k ð Þ.

Residual model
The frequency-domain expression relating the input u z ð Þ ∈ R p and the output y z ð Þ ∈ R q to the residual e z ð Þ ∈ R q is given by the following model termed the residual model: where D z ð Þ and N z ð Þ are matrix polynomials, F z ð Þ is the scalar characteristic polynomial termed Kalman polynomial, F z ð Þ ¼ zI À A 0 þK 0 C 0 ; (10) ; a ijℓ and b ijℓ are the coefficients of the polynomials D ij z ð Þ and N ij z ð Þ, respectively. The rational polynomials F À1 z ð ÞD z ð Þ and F À1 z ð ÞN z ð Þ associated with the system output y z ð Þ and the input u z ð Þ are termed as an output IIR filter, and an input IIR filter, respectively. The estimate of the Kalman filterŷ k ð Þ is: The residual model of the Kalman filter forms the backbone of the proposed identification scheme.

The key properties of the Kalman filter
The map relating the signal and its model, and the output IIR filter and an input IIR filter of residual model is developed next.
The following lemmas are developed by invoking the key property namely that the residual is a zero-mean white noise process if and only if there is no mismatch between the actual model of the system and its identified model embodied in the Kalman filter [4], that is, the identified model embodied in the Kalman filter is identical to that of the actual model:

Derivation of the signal and the signal model
The following Lemma 1 shows that (a) the estimate of the signal model is the matrix fraction description relating the transfer matrices relating the residual of the Kalman filter to the input, and the output of the system; (b) the estimate of the signal is its output generated by the system input; and the Kalman filter whitens the output error.
Lemma 1 (a) The left-matrix description of the MIMO signal model derived from the state- ð Þ and the left-matrix description of the Kalman filter derived from the residual Proof: (a) Consider the residual model (9). Substituting for y z ð Þ yields: Since the residual is a zero-mean, white noise process and is uncorrelated with u z ð Þ andυ z ð Þ, correlating both sides with the input u z À1 ð Þ yields: Assuming that the input correlation is not identically equal to zero, i.e. E u z ð Þu z À1 ð Þ ½ 6 ¼ 0 yields: where and

Derivation of the output error and its model
The following Lemma 2 shows that the output error is the difference between the output and the estimate of the signal; the estimate of the output error model is the matrix fraction description of the transfer matrices relating the residual of the Kalman filter to the input; the estimate of the output error is obtained as its output when its input is the residual. It is assumed that the transfer matrix of the Kalman filter F À1 z ð ÞD z ð Þ ¼ I À C 0 zI À A 0 þK 0 C 0 À Á À1 K 0 relating the output y z ð Þ and the residual e z ð Þ is minimum-phase, that is both the numerator and the denominator polynomials are asymptotically stable.

Lemma 2.
If the matrix D z ð Þis invertible, the output error ϑ z is termed the disturbance model estimate which generates the output error when excited by the residual.
As the input w k ð Þ driving the disturbance model is not accessible, then by substituting the actual input w k ð Þ by the residual e k ð Þ, although both are zero-mean white noise processes, only the denominator polynomial of the disturbance model can be identified. Hence the term "disturbance model estimate".
Minimum realization of the output error model is obtained using balanced model reduction method by treating e z ð Þ as the input andθ k ð Þ as the output of a model [7].

Minimal realization of the signal model
There are two approaches to identifying the signal model and the signal. One approach is by deriving them from the residual model of the Kalman filter as shown in Lemma 1 given by (12) and the other approach is to invoke the duality between the predictor form (7) and the innovation form (8) of the Kalman filter. The latter approach may be more convenient.
In view of (12), the system model G z ð Þ and the signal model G s z ð Þ is derived from the identified Kalman filter (7) by simply replacing the transition matrix of the Kalman filter A 0 À K 0 C 0 by the system transition matrixA 0 .

Proof
There is duality between the predictor form (7) and the innovation form (8) of the Kalman filter [5]. The output y k ð Þand the residual e k ð Þ are (causally) invertible. In other words, e k ð Þ can be generated from the output y k ð Þ and r k ð Þ using the (causally) predictor form, and y k ð Þ can be generated from e k ð Þ and r k ð Þ using the innovation form [3]. Moreover, s z ð ÞN s z ð Þ respectively, as shown in (4), implying that The minimum realization of the system A 0 s ; B 0 s ; C 0 s ; D 0 s À Á is obtained from the balanced model reduction method by treating u k ð Þ as the input andŝ k ð Þ as the output of a model [7].

Emulator-based two-stage identification
An identified model at each operating point characterizes the behavior of the system in the neighborhood of that point. In practice, however, the system model may be perturbed because of variations in the parameters of that system. To overcome this problem, the system model is identified by performing a number of emulator parameter-perturbed experiments proposed in [7][8][9]. Each experiment consists of perturbing one or more emulator parameters. A robust model is identified as the best fit to the input-output data from the set of emulated perturbations. The robust model thus obtained characterizes the behavior of the system over wider operating regions (in the neighborhood of the operating point) whereas the conventional model characterizes the behavior merely at the nominal operating point (that is, the conventional approach assumes that the model of the system remains unperturbed at every operating point). In [7][8][9], it is theoretically shown that the identification errors resulting from the variations in the emulator parameters are significantly lower compared to those of the conventional ones based on performing a single experiment (that is, without using emulators). The emulator-based identification scheme is inspired from the model-free artificial neural network approach which captures the static and dynamic behaviors by presenting the neural network with data covering likely operating scenarios. The PEM identifies the robust model of the plant, and the Kalman filter associated with the plant is then derived from the identified model without any a-priori knowledge of the statistics, such as covariance of the disturbance and measurement noise affecting the input-output data.
An accurate emulator-based model identification scheme is proposed and employed here. An emulator, which is modeled as a product of first-order all-pass filters and which induces phase and gain changes, is connected in cascade to the input, output or both, of the signal model to emulate a set of likely operating regimes around the nominal operating point. The identified model is obtained as the best fit over all emulated operating regions, thereby ensuring both accuracy and robustness of the identified model.

Two-stage identification
• In the first stage, a robust model of the system A 0 ; B 0 ; C 0 ; D 0 À Á and its associated Kalman filter are identified using PEM from the set of the emulator-generated input-output data. Then the estimate s 0 k ð Þ of the signal s k ð Þ and the estimateθ k ð Þ of the output error ϑ k ð Þ are derived.
• In the second stage, using the key properties established in Lemmas 1-3, the robust signal model A 0 s ; B 0 s ; C 0 s ; D 0 s À Á and its associated Kalman filter are obtained using balanced model reduction method and the PEM.
Akaike Information Criterion: To select an appropriate order for the identified system model in the first stage, and for the signal model in the second stage, the widely popular Akaike Information Criterion (AIC) is used, which weights both the parameter estimation error and the complexity of the model so as to arrive at an optimal order [1].

Signal model and the Kalman filter
Similar to the Kalman filter for the system (7), the Kalman filter for the signal is: Wherex s k ð Þ ∈ R n s ;ŝ k ð Þ ∈ R n s ; the residual e s k ð Þ ¼ e s1 k ð Þ e s2 k ð Þ e s3 k ð Þ … e sq k ð Þ Â Ã T ∈ R q is the residual; and K s ∈ R n s xq is the Kalman gain.
Status monitoring: The residuals e k ð Þ and e s k ð Þ of the Kalman filters (7) and (20) are employed to monitor the status of the overall system and to detect and isolate faults in the signal and disturbance models and the sensors. The proposed scheme provides a sound framework for developing fault-tolerant systems and conditionbased maintenance systems as well.

Evaluation on the illustrative example
The proposed two-stage identification scheme and the key properties of the Kalman filter established in the lemmas in Sections 2.3.1-2.3.3 are verified using the illustrative example given in Section 3.1. The results of this illustration are shown below in Figure 3a and b.
Subfigures A and B, of Figure 3a compare the true step response of the signal and its Kalman filter estimate; subfigures C and D show the output error ϑ k ð Þ and its estimate.
Remarks: These subfigures confirm the accuracy of the estimates of the signal and the output error (18) established in Lemmas 1 and 2. Subfigures A and B, of Figure 3b show the autocorrelation of the equation error whereas subfigures C and D show the autocorrelations of the residual.
Moreover, these subfigures clearly confirm that the equation error is a colored noise that is whitened by the KF, thus confirming (17) of Lemma 1 and making the KF residual a zero-man white noise process. Table 1 compares the true and estimated poles of the signal and disturbances models. The estimated poles are obtained from the model reduction techniques employed in the second stage of the two-stage identification scheme.
Remarks: The estimated poles are close to the true ones, especially those of the signal.

Evaluation of the proposed scheme
The management of leakage faults in fluid systems is becoming increasingly important in recent years from the point of view of economy, potential hazards, pollution, and conservation of scarce resources. Leakage in pipes and storage tanks occurs due to faulty joints, aging, excessive loads, holes caused by corrosion and accidents and the like. The process control system is a MIMO system that exhibits turbulence and is modeled as a combination of a signal, (which includes an ideal noise-free height, flow rate, and control input), a disturbance that includes effects of turbulence, and a measurement noise. The augmented model of the signal, and the disturbance, whose output is a sum of the signal, the disturbance, and the measurement noise described by Box-Jenkins model. The transfer matrices of the It is assumed that the disturbance is a Gaussian stochastic process and the measurement noise is a zero-mean Gaussian white noise process. The measurement output is, in general, an additive combination of the signal, disturbance and measurement noise. The output error, which is a sum of the disturbance and measurement noise, is assumed to be bounded. The signal and the disturbance are both modeled as outputs of linear time-invariant systems driven by some known input, and a Gaussian zero-mean white noise process, respectively. It is assumed that the signal, disturbance and measurement noise are mutually uncorrelated with each other.

Physical two-tank fluid system
3 the actuator G 1 ¼ G 0 1 γ a where G 0 1 is the fault-free transfer function, and subsystem 4 the leakage fault sensor gain γ ℓ . The fault-free cases correspond toγ si ¼ 1 : i ¼ 0, 1, 2, γ a ¼ 1 and γ ℓ ¼ 1. The various subsystems and sensor blocks are all shown in Figure 4b. The first two blocks G 0 and G 1 ¼ G 0 1 γ a , represent the controller and the actuator sub-systems, respectively. The leakage is modeled by the gain γ ℓ which is used to quantify the amount of flow lost from the first tank. Thus, the net outflow from tank 1 is quantified by the gain (1 À γ ℓ ). Since the two blocks G 0 2 and (1 À γ ℓ ) cannot be dissociated from each other, they are fused into a single block labeled G 2 ¼ G 0 2 1 À γ ℓ ð Þ. The physical two-tank system is controlled using LABVIEW which acquires the flow rate, and the height sensor outputs. The controller is implemented in LABVIEW and the controller output drives the actuator, namely the DC motor and pump combination. A fault in the sensor is introduced by including the emulator block, γ si : i ¼ 0, 1, 2 in the control input, flow rate, the height sensors, respectively in LABVIEW software. Similarly, an actuator fault is introduced by including an emulator γ a between the controller output and the input to the DC motor. The leakage fault is simulated by opening the drainage valve of the first tank. The amount by which the valve is opened is modeled by the emulator γ ℓ .
The height, flow rate, and control input profiles under various types of faults, are shown in Figure 5.    The height, the flow rate and the control input profiles under various types of faults are shown in Figure 5a for the nonlinear model, namely the dead-band effect of the actuator on the flow measurements. The measurement outputs are corrupted by the disturbance and measurement noise. Figure 5b show the outputs of the twostage identification of the linearized signal model. Subfigures A, B and C on the top show height profiles, and subfigures D, E and F in the middle show the flow rate profiles, and G, H, and I at the bottom show the control input profiles under leakage, actuator and sensor faults, respectively. The faults are induced by varying the appropriate emulator parameters to 0.25, 0.5 and 0.75 times the nominal values to represent 'small', 'medium' and 'large' faults. However, by its control design objective, the closed-loop PI controller will hide any fault that may occur in the system and hence will make it difficult to detect it. Also, the physical system exhibits a highly-nonlinear behavior. The flow rate saturates at 4.5 ml/s. The deadband effect in the actuator exhibits itself as a delay in the output response and saturation of the flow.
Remarks: The two-stage identification is employed to estimate the height, flow rate and the control input; their estimates are shown in Figure 5b. Comparing subfigures D, E, F confirms the superior performance of the identified estimates, thanks to the use of emulators. Figure 6a shows the residuals and their test statistics, and Figure 6b shows the autocorrelations of the residuals when the system is subject to leakage, actuator, and sensor faults of various degrees such a small, medium and large fault sizes. Subfigures A, B, and C; D, E, and F; and G, H, and I of Figure 6a show the residuals and their statistics when there is a leakage, actuator and sensor faults, respectively. Subfigures A, B, and C; D, E, and F; and G, H, and I of Figure 6b show the corresponding auto-correlations for different fault types.
Remarks: The Bayes decision strategy was employed to assert the fault type, i.e., to decide whether it is either a leakage or an actuator or sensor fault, respectively, using the fault isolation scheme proposed in [8]. The variance of the residual, which is the maximum value of the autocorrelation function evaluated at the origin (i.e. at zero delay), indicates the fault size.
The proposed Kalman-filter-based scheme can detect and isolate small and nascent faults and estimate the fault size. Thanks to the emulator-generated data, it can also provide an accurate prognosis of the status of the system.

Conclusions
Emulator-based identification of a wider class of multiple-input and multipleoutput system governed by Box-Jenkins model and the associated Kalman filter directly from the input-output data without a-priori knowledge of the disturbance and measurement noise statistics, and the establishment of key properties of estimation of the signal, the output error and their models are developed. The applications include monitoring the status of the system including faults, distinguishing between the variations in the disturbance model and those in the signal model to help diagnose a fault in the system and ensure a low false alarm probability, developing a framework for controlling autonomous vehicles, and meeting the ever-increasing need for fault-tolerant systems. The proposed emulator-based two-stage identification and estimation of the signal and its model were evaluated physical laboratory-scale process control system so as to estimate the signal corrupted by disturbance such as turbulence. Thanks to emulator-based identification, the estimates of the signal were accurate, the detection and isolation of leakage faults were promising and, as such, provide sufficient encouragement and impetus to try the proposed scheme on real-life processes in our future work.