key: cord-0054385-z5qryfop
authors: dos Santos Gomes, Daiana Caroline; de Oliveira Serra, Ginalber Luiz
title: Computational Approach For Real-Time Interval Type-2 Fuzzy Kalman Filtering and Forecasting via Unobservable Spectral Components of Experimental Data
date: 2021-01-03
journal: J Control Autom Electr Syst
DOI: 10.1007/s40313-020-00675-9
sha: 0c490ea5b7da082b6c22c94d08a7bb251c910c2a
doc_id: 54385
cord_uid: z5qryfop

In this paper, a methodology for design of Kalman filter, using interval type-2 fuzzy systems, in discrete time domain, via spectral decomposition of experimental data, is proposed. The adopted methodology consists of recursive parametric estimation of local state space linear submodels of interval type-2 fuzzy Kalman filter for tracking and forecasting of the dynamics inherited to experimental data, using an interval type-2 fuzzy version of Observer/Kalman Filter Identification (OKID) algorithm. The partitioning of the experimental data is performed by interval type-2 fuzzy Gustafson–Kessel clustering algorithm. The interval Kalman gains in the consequent proposition of interval type-2 fuzzy Kalman filter are updated according to unobservable components computed by recursive spectral decomposition of experimental data. Results illustrate the efficiency of proposed methodology, as compared to other approach widely cited in the literature, for filtering and tracking the state variables of Lorenz’s chaotic attractor in a noisy environment, and its applicability for adaptive and real-time forecasting the dynamic spread behavior of novel coronavirus 2019 (COVID-19) outbreak in state of Maranhão and Brazil.

In science and engineering is very common the solution of problems with stochastic nature such as prediction, separation and detection of signals in the presence of random noise (Mack and Habets 2020; Gomez-Garcia et al. 2020; Liu et al. 2019; Zhu et al. 2019 ). Kalman filter (KF) is the most well-known and used mathematical tool for stochastic estimation from noisy and uncertain measurements. It was proposed by Rudolph E. Kalman in 1960, who published his famous paper "A New Approach to Linear Filtering and Prediction Problem" (Kalman 1960 becoming a standard approach for optimal estimation. Since the time of its introduction, the Kalman filter has been the subject of extensive research and applications in the fields of orbit calculation, target tracking, integrated navigation, dynamic positioning, sensor data fusion, microeconomics, control, modeling, digital image processing, pattern recognition, image segmentation and image edge detection and others. This broad interest in KF is due to its optimality, convenient form for online real-time processing, easy formulation and implementation (Serra 2018) .

Fuzzy systems have been widely used for modeling (Wang and Wu 2019; Chin and Lin 2018; Kim et al. 2020a; Zhao and Lin 2019) and also for health area (Khodaei-Mehr et al. 2018; Pham and Berger 2011) . The successful applications of fuzzy systems are due to its structure based on rules, where the antecedent propositions of the rules define fuzzy operation regions and the consequent describes a corresponding physical behavior in those regions, and its capability of approximate functions as well as treat nonlinearities and uncertainties (Serra 2012) . Recently, type-2 fuzzy systems have been highlighted in several applications due to their better ability to deal with uncertain information (Evangelista and Serra 2019; Mendel 2019) . A special case within the study of type-2 fuzzy logic is the interval type-2 fuzzy sets that, by simplifying the membership functions assigned to the sets, define a footprint of uncertainty (FoU) in data processing, which is limited by upper and lower membership functions (Liang and Mendel 2000) .

With the emergence caused by COVID-19 and its rapid spread to several countries, authorities around the world have been implemented plans for fighting transmission of the virus based on guidelines provided by the World Health Organization (WHO) (Ryu and Chun 2020; WHO 2020a, b) . The impacts caused by COVID-19 pandemic have affected health, social, economic, political and cultural spheres, unprecedented in the recent worldwide history of epidemics (Chakraborty and Maity 2020) . In view to the problems faced to COVID-19 crisis, researchers from all scientific areas have proposed new studies and contributions as an effort to fight and understand the coronavirus disease Fishbane and Hirsch 2020; Sun et al. 2020; Brown et al. 2020; Huang et al. 2020; Kang et al. 2020) . In these contexts, computational modeling applied to analysis of epidemiological data has received increasing interest from the scientific community, aiming to characterize the dynamic evolution of infectious diseases outbreak and also helping in its control and prevention (Deeba et al. 2020; Price et al. 2019; Korcinska et al. 2020; Chen et al. 2020; van Gaalen et al. 2017; Sloan et al. 2020; van de Kassteele et al. 2019; Chowell et al. 2020) .

Associated to epidemiological data analysis, the proposal for adaptive and real-time modeling methodologies which take into account the processing of uncertainties inherited to these types of experimental data (underreportings, lack of information, incubation period of the virus, time to seek care and diagnosis) is still open. The proposed methodology, based on interval type-2 fuzzy Kalman filter, is useful for adaptive and real-time forecasting of uncertain experimental data. A new formulation of type-2 fuzzy version of Observer/Kalman Filter Identification (OKID) algorithm, is proposed, for updating, recursively, the consequent proposition of type-2 fuzzy Kalman filter using spectral components extracted from the experimental data. Interval fuzzy sets characterizing the antecedent of type-2 fuzzy Kalman filter inference system are estimated by, also proposed, a formulation of interval type-2 fuzzy version of Gustafson-Kessel clustering algorithm. Computational results of tracking the states of a nonlinear dynamic system with chaotic behavior in a noisy environment show the efficiency of proposed methodology as compared to other approach widely cited in the literature. Applicability of proposed methodology is illustrated through experimental results from real-time interval forecasting of the COVID-19's dynamic spread behavior in state of Maranhão and Brazil.

In the last years, studies involving the integration of fuzzy systems and Kalman filters have been proposed in the literature (Pires and Serra 2019; Eyoh et al. 2018; Gil et al. 2019; Bouhentala et al. 2019; Hwang et al. 2019) . In , fuzzy sets are combined with a novel trace optimization method based on extended Kalman filter with nested probabilistic-numerical linguistic information applied in tracking a maneuvering target. Based on limited and uncertain information from different sensors, the proposed methodology is able to merge this information and is efficient when applied in the problem of trace optimization of an unknown maneuvering target in Sichuan province in China. In (Asl et al. 2020) , an optimization methodology of an adaptive unscented Kalman filter (UKF) is presented by means of an evolutionary fuzzy algorithm named fuzzy adaptive grasshopper optimization algorithm. The efficiency of the method is verified by application to different benchmark functions, such as a robotic manipulator and a servo-hydraulic system, performing better when compared to previous versions of UKF. In (Matía et al. 2019) , a reformulation in the uncertainty representation in fuzzy Kalman filters (FKF) is proposed, to solve the problem of uncertainty propagation that occurred in other approaches present in the literature for the FKF. Trapezoidal possibility distributions is used to represent fuzzy variables, defining regions of possibility and impossibility in asymmetric sets, which contributes to model sensor's uncertainty with more accuracy. In (Yang et al. 2019 ), a strategy based on fuzzy logic control to treat pathological symptoms of movement disorder with higher performance is presented. A slowly varying hidden variable in a neural network is estimated using an unscented Kalman filter and is used as a feedback variable to enhance control performance. The presented design, with enhanced control performance and higher hardware efficiency, has significant potential for clinical treatment of movement disorders and opens a new perspective for the applications in the fields of neural control engineering and brain-machine interfaces. In (Benhamida et al. 2019 ), an improved model predictive direct torque control (MPDTC) based on an adaptive fuzzy logic modulator and an extended Kalman filter is introduced. The methodology consists of an adaptive fuzzy logic-based duty cycle vector modulation, and two voltage vectors (VV) instead of a single VV are applied during the whole control cycle to promote the torque control performance and reduce its undulations. Moreover, a design of an extended Kalman filter estimator for drive systems is proposed. Despite the extensive literature in this context, there are still many fields to be explored regarding the association of Kalman filters and fuzzy logic.

Several mathematical, computational and statistical methods have already been proposed and widely applied in the prediction of infectious diseases worldwide (Rajaei et al. 2019; Watkins et al. 2020; Martins et al. 2018; He et al. 2018) . In (Weng et al. 2020 ), a seasonal autoregressive integrated moving average (SARIMA) model was used to predict the incidence of Chlamydia Trachomatis (CT) infection in Shenzhen city, China. The proposed model breaks down the time series data corresponding to monthly cases of CT infection into behavioral patterns such as trend, seasonal and random, in order to understand the epidemiological behavior of the disease and obtain more accurate predictions. In (Koolhof et al. 2020) , negative binomial regression models were developed to forecast human infections by Ross River virus (RRV), which is Australia's most epidemiologically important mosquito-borne disease. The model uses data from climatic, environmental and oceanographic variables in order to understand the factors associated with RRV transmission in epidemiologically important regions in the state of Victoria and to establish an early warning forecasting system. In (Hranac et al. 2019) , ensemble niche models to predict spatiotemporally varying bat birthing to examine how birthing cycles of African fruit bats, molossid bats, and non-molossid microbats inform the spatiotemporal occurrence of Ebola Virus Disease (EVD) spillover. The model uses pool sparse data and predicted the risk of EVD spillover at locations of the two 2018 EVD outbreaks in the Democratic Republic of the Congo. In (Chowell et al. 2020 ), an ensemble model was introduced for sequential forecasting that weights a set of plausible models and use a frequentist computational bootstrap approach to evaluate its uncertainty. The feasibility of the approach was demonstrated using simple dynamic differential equation models and the trajectory of outbreak scenarios of the Ebola Forecasting Challenge, providing forecasting horizon of 1-4 weeks. In (Stocks et al. 2020 ), a novel multilayered dynamic transmission model is presented for hepatitis C virus HCV transmission within a people who inject drugs (PWID) population. The approach based on this model is able to estimate the number of undiagnosed PWIDs, the true incidence, the average time until diagnosis, the reproduction numbers and associated uncertainties, using routine surveillance data.

Recently, with the beginning of the COVID-19 epidemic outbreak, several researchers have proposed model-based data analysis approaches applied to novel Coronavirus 2019 (Zhong et al. 2020; Lin et al. 2020; Mohd and Sulayman 2020; Duan and Zhang 2020; Chintalapudi et al. 2020; Fredj and Chrif 2020) . The objective of these studies is to characterize the evolution of pandemic in certain regions and, thus, to contribute for the requirements adopted to contain the contamination by virus and allocation of resources. In this sense, a mathematical model for forecasting the transmission dynamics of COVID-19 in Korea is proposed in (Kim et al. 2020b ). The mathematical model is based on SEIR model (Susceptible-Exposed-Infectious-Recovered) and takes account the behavioral changes that were implemented in the population since local transmission began. The study is able to predict the final size and the timing of the end of the epidemic as well as the maximum number of isolated individuals using daily confirmed cases comparing epidemiological parameters between the national level and the Daegu/Gyeongbuk area. The study proposed in (Park et al. 2020 ) addresses the role of asymptomatic carriers in transmission poses challenges for control of the COVID-19 pandemic. A mathematical framework is used to evaluate expected effects of asymptomatic transmission on the basic reproduction number R 0 (i.e., the expected number of secondary cases generated by an average primary cases in a fully susceptible population) and the fraction of new secondary cases attributable to asymptomatic individuals. The proposed methodology models the viral spread using a renewal equation framework, which allows to model the current incidence of infected individuals as a function of previous incidence and how infectiousness of an infected individual varies over the course of their infection. This analysis shows that understanding the temporal course of asymptomatic transmission can be important for assessing the importance of this route of transmission and for disease dynamics. In (Kanagarathinam and Sekar 2020), a SEIR model for COVID-19 was developed to show the importance of estimation of reproduction number (R 0 ). The work is focused on predicting the 2019 Novel Coronavirus outbreak at the early stage in India based on estimation of reproduction number (R 0 ) and help to take active measures prior to the spread of 2019-nCoV disease. The data of daily newly infective cases in India have been used as samples to obtain the reproduction number. The maximum likelihood approach has been carried out to the distribution of R 0 . In (Rustam et al. 2020) , several supervised machine learning models are compared to predict the number of patients affected by COVID-19. The four standard forecast models were used: linear regression (LR), least absolute shrinkage and selection operator (LASSO), support vector machine (SVM) and exponential smoothing (ES). Three types of predictions are made by each of the models, such as the number of newly infected cases, the number of deaths and the number of recoveries in next 10 days.

Differently from aforementioned approaches and others ones found from the literature, the scope of this paper outlines the integration of Kalman filter and interval type-2 fuzzy systems for tracking and forecasting the COVID-19 dynamic spread behavior. The design of interval type-2 fuzzy Kalman filter, according to proposed methodology, is based on spectral unobservable components and uncertainty regions extracted from experimental data.

The originality of proposed methodology is outlined by the following main contributions:

-Applicability for adaptive and real-time forecasting of uncertain experimental data; -A new formulation of interval type-2 fuzzy version of Gustafson-Kessel clustering algorithm for estimating and shaping the interval fuzzy sets characterizing the antecedent proposition of type-2 fuzzy Kalman filter inference system; -A new formulation of type-2 fuzzy version of Observer/ Kalman Filter Identification (OKID) algorithm, for recursive updating the consequent proposition of type-2 fuzzy Kalman filter based on unobservable spectral components extracted from the experimental data.

In this section, the proposed methodology for designing the interval type-2 fuzzy Kalman filter computational model from experimental data is presented. Formulations for pre-processing the experimental data by spectral analysis, parametric estimation of interval type-2 fuzzy Kalman filter antecedent proposition, parametric estimation of interval type-2 fuzzy Kalman filter consequent proposition and its recursive updating mechanism are addressed.

The singular spectral analysis technique is a mathematical tool for analyzing and decomposing complex time series into simpler components within the original data. Such unobservable components have relevant characteristics about the corresponding time series behavior (Elsner 2002 ).

Let the initial experimental data set referring to p outputs of the dynamic system under analysis, with N b samples, given by:

where y k ∈ R p , with k = 1, . . . , N b , is the output vector of the dynamic system at instant of time k. From this initial data set, a trajectory matrix H is defined, for each of the dimensions of Y, considering a set of ρ delayed vectors with dimension δ, which is an integer number defined by user with 2 ≤ δ ≤ N b − 1 and ρ = N b − δ + 1, given by:

y 1 y 2 y 3 · · · y ρ y 2 y 3 y 4 · · · y ρ+1 . . . . . . . . . . . . . . .

and the covariance matrix S is obtained as follows:

Applying the singular value decomposition (SVD) procedure to matrix S, is obtained a set of eigenvalues in decreasing order such that σ 1 ≥ σ 2 ≥ · · · ≥ σ δ ≥ 0 with their respective eigenvectors φ 1 , φ 2 , . . . , φ δ . Considering d = max {ς, such that σ ς > 0}, and V ς = H T φ ς / √ σ ς with ς = 1, . . . , d, the singular value decomposition of the trajectory matrix H, can be rewritten as:

where the matrix H ς | ς=1,...,d is elementary (it has rank equal to 1) and is given by:

The regrouping of H ς | ς=1,...,d into ξ linearly independent matrices terms I j | j=1,...,ξ , such that ξ ≤ d, results in H = I 1 + I 2 + · · · + I ξ

The unobservable components α j | j=1,...,ξ extracted from experimental data, resulted from the matrices I j | j=1,...,ξ , are given by:

where δ * = min(δ, ρ), ρ * = max(δ, ρ) and N b = δ + ρ − 1.

After the initialization of spectral analysis algorithm in the training step, the next steps will be repeated for each time instant k = N b + 1, N b + 2, . . . , as formulated in sequel. The value of ρ is increased by:

The covariance matrix is updated, recursively, as follows:

where ϒ k = ψ k ψ T k ∈ R δ×δ with ψ k = y ρ , y ρ+1 , . . . , y k T ∈ R δ×1 . Applying SVD procedure to covariance matrix S k , the set of eigenvalues and their respective eigenvectors are updated such that the dynamic system output y k can be rewritten by: 

. . , represents the samples of extracted unobservable components at instant k. The recursive singular spectral analysis, according to proposed methodology, is implemented as described in Algorithm 1.

The spectral components extracted by singular spectral analysis in section 2.1 that presents smaller eigenvalues are considered as residuals and they are used as data set for recursive estimation of measurement noise covariance R, inherited to experimental data. The initial covariance of measurement noise v, is given by:

where N b is the length of experimental data set to be used in training step for initial parameterization of interval type-2 fuzzy Kalman filter and v is the spectral component α j considered as residual. The recursive updating of covariance R, at instants of time k = N b + 1, N b + 2, . . . , is given by:

Once the covariance R of the residual v is computed, an adaptive factor X k , with k = N b +1, N b +2, . . . , is proposed,

..,ξ %Training step;

Step 1: Compute ρ = N b − δ + 1;

Step 2: Construct the trajectory matrix H -Eq. (2);

Step 3: Compute the covariance matrix S -Eq. (3);

Step 4: Apply the SVD method in the covariance matrix S and obtain the set of eigenvalues σ 1 ≥ σ 2 ≥ · · · ≥ σ δ ≥ 0 with their respective eigenvectors φ 1 , φ 2 , . . . , φ δ ;

Step 5: Rewrite the SVD of the matrix S in the form of Eq. (4);

Step 6: Regroup the matrices H ς | ς =1,...,d in ξ linearly independent matrix terms I j | j=1,...,ξ -Eq. (6);

Step 7: Compute the unobservable components α j k | j=1,...,ξ -Eq. (7); %Recursive step;

Step 1: Update ρ -Eq. (8);

Step 2: Update the covariance matrix S k -Eq. (9);

Step 3: Update the set of eigenvalues

Step 4: Rewrite the sample y k in the form of Eq. (10);

Step 5: (11); end to be used as weighting factor for implementing the recursive updating mechanism of interval type-2 fuzzy Kalman filter, as follows:

The algorithm for recursive estimation of noise covariance R and weighting factor X , according to proposed methodology, is implemented as described in Algorithm 2.

The adopted structure of interval type-2 fuzzy Kalman filter presents the i| [i=1,2,...,c] -th fuzzy rule, given by: Compute the initial measurement noise covariance R -Eq. (12);

Compute the initial adaptive factor X k -Eq. (14); (14); end with n-th order, m inputs, p outputs, where Z k is the linguistic variable of the antecedent; W i is the interval type-2 fuzzy set; x i k ∈ R n is the estimated interval states vector of the nonlinear dynamic system; ŷ i k ∈ R p is the estimated interval output vector and u k ∈ R m is the input signal. The matrices

are, respectively, state matrix, input matrix, output matrix, direct transmission matrix and Kalman gain matrix, which are uncertain parameters that describe the dynamics of the nonlinear system within a region of uncertainty. The residual

where y k ∈ R p is the real output of the dynamic system and ŷ i k is the interval estimated output by i-th linear submodel, and the X k is the weighting factor, as formulated in section 2.1.3, which corresponding to the level of filtering to be applied by interval type-2 fuzzy Kalman filter on the experimental data.

The interval type-2 fuzzy Kalman filter approximates the dynamic behavior inherited to experimental data through the weighted sum of Kalman filters defined in the consequent proposition of type-2 fuzzy Kalman filter rules, according to normalized interval activation degrees, of each i-th rule, as follows:

corresponds to lower and upper activation degrees in i-th cluster, respectively, and c is the number of rules of interval type-2 fuzzy Kalman filter, such that:

The partitioning of experimental data implies to computing of operating regions and, necessarily, the number of rules of the interval type-2 fuzzy Kalman filter. An interval type-2 fuzzy version of Gustafson-Kessel clustering algorithm was proposed, which is formulated in the sequel. Given the experimental data set Z, previously collected, of size n × N b , choose the number of clusters 1 < c < N b ; the initial interval partition matrix U (0) ∈ R c×N b , the termination tolerance E > 0 and the interval weighting exponent m = m, m , where m and m correspond to, respectively, weighting exponent of upper and lower membership functions. Repeat for l = 1, 2, . . .

Step 1 -Compute the centers of the clusters:

Step 2 -Compute the covariance matrices of the clusters:

Step 3 -Compute the distances:

Step 4 -Update the partition matrix:

where

The interval type-2 fuzzy Gustafson-Kessel clustering algorithm, according to proposed methodology, is implemented as described in Algorithm 3.

Algorithm 3: Interval Type-2 Fuzzy Gustafson-Kessel Clustering Algorithm

Step 1: Compute the centers of the clusters v i (l) -Eq. (19);

Step 2: Compute the covariance matrices of the clusters F i -Eq. (20);

Step 3 -Compute the distances D i k F i -Eq. (21); Step 4: Update the partition matrix U (l) -Eq. (22)-(24); end until U (l) − U (l−1) < E;

Once the interval type-2 fuzzy Kalman filter is obtained by black box modeling approach, it is necessary to estimate, from experimental data, the matrices that compose the interval type-2 fuzzy Kalman filter inference system consequent proposition, as formulated in Eq. (15). In this sense, an interval type-2 fuzzy OKID (Observer/Kalman Filter Identification) algorithm, is proposed. The interval type-2 membership values from partitions defined on experimental data, which were estimated by interval type-2 Gustafson-Kessel clustering algorithm, have been considered as weighting criteria for computing the consequent proposition of interval type-2 fuzzy Kalman filter inference system. The interval type-2 fuzzy OKID algorithm is formulated in the sequel.

Let the experimental data set Z, such that

where α * k corresponds to spectral components extracted from the experimental data. Choose an appropriate number of Markov parameters q, through the following steps:

Step 1 -Compute the matrix of regressors , given by:

Step 2 -Compute the interval observer Markov parameters

is the diagonal weighting matrix of the i-th fuzzy rule obtained from the interval type-2 GK fuzzy clustering algorithm and is the interval observer Markov parameters of i-th rule such that Manipulating Eq. (26):

Assuming

rewriting as:

Equation (32) is solved by QR factorization method, which is numerically robust since it avoids matrix inverse operations. 

where Q i is an orthogonal matrix, such that Q i −1 = Q i T and R i is an upper triangular matrix. Because the matrix R i is upper triangular, Eq. (33) The temporal behavior of a measurement noise covariance R k associated to state variable x 1 and b the weighting factor X k used for implementing the recursive updating mechanism of interval type-2 fuzzy Kalman filter

Step 3 -Compute the observer gain and system Markov parameters:

Thus, the system Markov parameters are obtained as follows: 

Step 4 -Construct the Hankel matrix H i ( j − 1) ∈ R γ p×βm :

where γ and β are sufficiently large arbitrary integers defined by user.

Step 5 -For j = 1, decompose H i (0) using Singular Value Decomposition:

where i ∈ R α p×α p and i ∈ R βm×βm are orthogonal matrices and i ∈ R α p×βm is the diagonal matrix of singular values defined as:

such that n is the number of significant singular values and determines the minimum order of the type-2 fuzzy Kalman filter. Thus, the size of matrices in Eq. (45) is reduced to the minimum order, as follows:

where i n ∈ R α p×n , i n ∈ R βm×n , i n ∈ R n×n are the resulting matrices after the reduction to minimum order.

Step 6 -Compute the observability and controllability matrices: 

is the observability matrix and

is the controllability matrix.

Step 7 -Compute the matrices that make up the consequent proposition of interval type-2 fuzzy Kalman filter:

Step 8 -Compute the interval Kalman gain matrix:

where Y i o j is the observer gain Markov parameters, P i β is the observability matrix and K i k is the interval Kalman gain matrix. Manipulating Eq. (56):

Assuming

(57) is rewritten as follows:

Equation (58) is solved by QR factorization method being applied to A i and obtaining the interval Kalman gain matrix K i k in the same way as done in

Step 2 for determining interval Markov parameters.

Inference System After the initial estimation of interval type-2 fuzzy Kalman filter, the Kalman filters, into consequent proposition of interval type-2 fuzzy Kalman filter inference system, are updated recursively at instants of time k = N b + 1, k = N b + 2, . . . , to each new sample from experimental data set. Considering the regressors vector, at instant k, given by

the interval observer Markov parameters Y i are obtained by recursive updating of Eq. (32), as follows:

Once U i k and ℵ i k have been updated, and applying the QR factorization in U i k , the i-th Markov parameters are updated. The consequent proposition of the type-2 fuzzy Kalman filter is updated recursively by repeating the Step 3 to Step 7. Similarly, the interval type-2 fuzzy Kalman gain matrix K i k is obtained by recursive updating of Eq. (58), as follows:

Once A i k and N i k have been updated, and applying the QR factorization method in A i k , the interval type-2 fuzzy Kalman gain matrix is updated.

The interval type-2 fuzzy Observer/Kalman Filter Identification algorithm, according to proposed methodology, is implemented as described in Algorithm 4.

In this section, computational results from filtering and tracking the states of a nonlinear dynamic system with chaotic behavior in a noisy environment, are shown to illustrate the efficiency of proposed methodology in comparison to approach in (Páramo-Carranza et al. 2017) . Also, experimental results from real-time interval tracking and forecasting the COVID-19's dynamic spread behavior in state of Maranhão and Brazil are shown to demonstrate the applicability of proposed methodology. 

In this section, computational results for filtering and tracking the states of a nonlinear dynamic system, with chaotic for i = 1 to c do

Step 1: Compute the interval observer Markov parameters Y i -Eq. (31)-(33);

Step 2: Compute the interval system Markov parameters Y i -Eq. (38)-(40);

Step 3: Compute the interval observer gain Markov parameters Y i o -Eq. (41)-(43);

Step 4: Construct the Hankel matrices H i (0) and H i (1) -Eq. (44);

Step 5: Decompose H i (0) using SVD method -Eq. (45) and determine the minimum order of realization n for the application -Eq. (46);

Step 6: Compute the observability matrix P Step 1: Construct the regressors vector λ k -Eq. (59); for i = 1 to c do

Step 2: Update the interval Markov parameters through Eq. (60)-(61);

Step 3: Repeat Step 3 to Step 7 described in training step;

Step 4: Update the interval Kalman gain matrix through Eq. (62)-(63) end end tor, described by (Huang et al. 2015 ):

where the parameters = 10, = 8/3 and = 28 provide a chaotic behavior.

Once that the problem of interest, in this paper, is based on the time series related to state variables x 1 , x 2 and x 3 from Lorenz's chaotic attractor, the variable u k , in Eq. 15 of proposed methodology, is considered as white noise signal with low amplitude. A data set from Lorenz's chaotic attractor, with total length of 10000 samples by a sampling period of T = 1 ms, was generated. The first 4000 samples were used in training step for initial parameterization of interval type-2 fuzzy Kalman filter. The unobservable components associated to state variable x 1 of Lorenz's chaotic attractor were extracted by singular spectral analysis approach, according to section 2.1, for pre-processing of data set. The noisy state variable x 1 was decomposed into 2 spectral components, where the component with residual behavior was used for estimation of measurement noise covariance R k and weighting factor X , and the other component, assuming correlated to nominal dynamic of state variable x 1 , was used for parameterizing the interval type-2 fuzzy Kalman filter. The partitions of computational data from noisy state variable x 1 were defined by interval type-2 fuzzy Gustafson-Kessel clustering algorithm, so the antecedent proposition, the rules number and consequent proposition, of the interval type-2 fuzzy Kalman filter, could be estimated successfully. For implementing the interval type-2 fuzzy clustering algorithm, the following parameters were adopted: number of clusters c = 2, interval weighting exponent m = [1.5, 2.0] and termination tolerance E = 10 −5 . According to computational data of Lorenz's chaotic attractor state variables x 2 and x 3 , the unobservable components from noisy state variable x 1 , the interval type-2 fuzzy normalized membership values and the parameterization of interval type-2 fuzzy Kalman filter by training and recursive steps were performed. The parametric estimation of consequent proposition of interval type-2 fuzzy Kalman filter inference system, in Eq. 15, took into account the partitions on noisy state variable x 1 as weighting criterion, and the parameters values: q = 1, γ = 10 e β = 10. The signal-to-noise ratio (SNR) in the state variable x 1 was varied, so that the ultimate goal was for filtering and tracking, from noisy computational data set, the nominal state variables x 1 , x 2 and x 3 as accurate as possible. The mean square error (MSE) values, used as validation criterion, for different levels of SNR, considering the effect of 100 realizations in the filtering and tracking the nominal state variables x 1 , x 2 and x 3 , are shown in Figs. 1, 2, 3, respectively. It is observed the superior performance of interval type-2 fuzzy Kalman filter, based on proposed methodology, for filtering and tracking the Lorenz's chaotic attractor state variables, considering the lower values of MSE obtained by upper and lower limits of the interval type-2 fuzzy Kalman filter as compared to values of MSE obtained by approach in (Páramo-Carranza et al. 2017) .

Considering the particular case of S N R = 0.65, the dynamic behavior for noisy and nominal cases of Lorenz's chaotic attractor state variable x 1 is shown in Fig. 4 . According to experimental data of state variable x 1 shown in Fig. 4 , the pre-processed unobservable components in Fig. 5 and the interval type-2 fuzzy normalized membership values in Fig. 6 , the parametric estimation of the type-2 fuzzy Kalman filter was computed by training and recursive steps. The com- Fig. 19 The temporal behavior of spectral unobservable components α j | j=1,...,7 , which were extracted from experimental data of daily deaths, in the state of Maranhão ponent with residual behavior was used for estimation of measurement noise covariance R k and weighting factor X , and the other component, assuming correlated to nominal dynamic of state variable x 1 , was used for parameterizing the interval type-2 fuzzy Kalman filter. As it can be seen in Fig. 5 , the unobservable component α 1 presenting a residual behavior was used for estimation of measurement noise covariance R k and weighting factor X , and the unobservable component α 2 assuming correlated to nominal dynamic of state variable x 1 was used for parameterizing the interval type-2 fuzzy Kalman filter. Based on estimation of measurement noise covariance R k associated to state variable x 1 , the weighting factor X k to be used as weighting factor for implementing the recursive updating mechanism of interval type-2 fuzzy Kalman filter is given by:

The temporal behavior of the measurement noise covariance R k and the weighting factor X k , which are associated with filtering and tracking the state variable x 1 from Lorenz chaotic attractor, by interval type-2 fuzzy Kalman filter, is shown in Fig. 7 . The confidence region, as shown in Fig. 8 , created by interval type-2 fuzzy Kalman filter taking into account uncertainties, estimated by interval type-2 membership functions shown in Fig. 6 , inherited to noisy data from Lorenz's chaotic attractor, illustrates its efficiency for filtering and tracking the state variables x 1 , x 2 and x 3 . The estimation of interval type-2 fuzzy Kalman gain matrices K i | i=1,...,2 , during recursive updating of interval type-2 fuzzy Kalman filter for filtering and tracking the state variables x 1 , x 2 and x 3 of Lorenz's chaotic attractor, is shown in Fig. 9 . The recursive estimation of the interval type-2 fuzzy matrices A i , B i , C i and D i , with i = 1, . . . , 2, in the consequent proposition of the interval type-2 fuzzy Kalman filter inference system, during recursive updating of interval type-2 fuzzy Kalman filter for filtering and tracking the state variables x 1 , x 2 and x 3 of Lorenz's chaotic attractor, is shown in Figs. 10, 11 , 12, 13. The upper and lower instantaneous activation degrees related to interval type-2 fuzzy Kalman filter inference system, for filtering and tracking the state variables 

In this section, a more detailed discussion on the results shown in Figs. 1, 2, 3, 8, 16 and Table 1 Interval type-2 fuzzy matrix C i , during recursive updating of interval type-2 fuzzy Kalman filter for tracking and forecasting the COVID-19 dynamic spread experimental data related to daily deaths reports within period ranging from 19 of May 2020 to 25 of September, in state of Maranhão: a Rule 1, b Rule 2, c Rule 3 from discretization of the Lorenz chaotic nonlinear model in Eq. 64, the linearization at specific operation points of state variables x 1 , x 2 and x 3 , and the definition of membership functions with pre-established format for partitioning the dynamic behavior of states x 1 , x 2 and x 3 . The adopted discretization for fuzzy Kalman filter design is based on the current sample (at instant k) of the dynamic behavior Fig. 28 Instantaneous normalized fuzzy activation degrees of type-2 fuzzy Kalman inference system, during its recursive updating for tracking and forecasting the COVID-19 dynamic spread experimental data related to daily deaths reports within period ranging from 19 of May Fig. 29 Real-time efficiency, based on VAF validation criterion, of interval type-2 fuzzy Kalman filter, during its recursive updating for tracking and forecasting the COVID-19 dynamic spread experimental data related to daily deaths reports within period ranging from 19 of May 2020 to 29 of September, in state of Maranhão Fig. 30 The experimental data of daily deaths reports within period from 29 of February 2020 to 18 of May 2020, in Brazil formance for high errors (Franklin et al. 1997; Serra 2018) . The adopted linearization considers the application of partial derivatives to Lorenz chaotic attractor in specific operation points of the state variables x 1 , x 2 and x 3 , which implies a filter design dynamically limited to specified operating points (Luenberger 1979; Chen 1999; Antsaklis and Liu 2003) . Thus, as the dynamic behavior of Lorenz chaotic attractor, in a noisy environment, gets out the operating points of linearized state variables x 1 , x 2 and x 3 , the approximation errors increase in prediction and updating steps of discrete fuzzy Kalman filter and polarize the estimation of state variables x 1 , x 2 and x 3 . The adopted partitions of the discourse universe Fig. 31 The efficiency of unobservable components, according to VAF criterion, for representing the experimental data of the daily deaths reports, in Brazil of state variable x 1 consider operation regions with preestablished membership functions from expert knowledge, discarding any similarity degree from dynamic behavior of state variable x 1 , which implies the increasing of filtering errors (Wang 1997; Babuska 1998; Sato-Ilic and Jain 2006) .

On the other hand, the methodology for interval type-2 fuzzy Kalman filter design, proposed in this paper, considers the pre-processing of the noisy dynamic data of Lorenz chaotic attractor state variables, the computation of an interval adaptive similarity measure related to dynamic behavior of state variable x 1 resulting in membership functions with different formats and orientations for estimating regional models useful to characterize the dynamic behavior of Lorenz chaotic attractor. The adopted pre-processing method for fuzzy Kalman filter design is based on spectral decomposition of the noisy dynamic data of state variable x 1 , resulting in unobservable spectral components with interpretable data patterns which allows to discard the spectral components associated with noise characteristic for better mathematical characterization of Lorenz chaotic attractor in the operation regions, recursively updated at each sampling time k, in order to reduce the filtering errors (Golyandina and Zhigljavsky 2013; Abdollahzade et al. 2015; Hossein and Saeid 2015; Juang 1994; Qi et al. 2019) . For partitioning the dynamic data of state variable x 1 an interval adaptive similarity mechanism is adopted to define interval operation regions represented by interval membership functions, with different formats and orientations, adapted to the topological structure associated to the variability of dynamic data, in order to minimize the filtering errors (Babuska 1998; Höppner et al. 1999 ; Chaomurilige and Yang 2017).

The temporal behavior of unobservable spectral components α j | j=1,...,10 , which were extracted from experimental data of daily deaths, in Brazil 

In this section, experimental results for forecasting analysis of dynamic spread behavior taking into account the experimental data of daily deaths reports caused by COVID-19 in state of Maranhão and Brazil are presented.

The experimental data of daily deaths reports in the period ranging from 18 of March 2020 to 18 of May 2020 in state of Maranhão are shown in Fig. 17 , which were obtained from the database provided by the State Health Organization of Maranhão 1 . Once that the problem of interest, in this paper, Fig. 34 The confidence region generated by interval type-2 fuzzy Kalman filter for tracking the experimental data of daily deaths reports, from 29 of February 2020 to 18 of May 2020, in Brazil is based on the time series related to daily deaths reports in state of Maranhão, the variable u k , in Eq. 15 of proposed methodology, is considered as white noise signal with low amplitude.

The pre-processing of experimental data by singular spectral analysis was able to extract the unobservable components associated to daily deaths reports. The variance accounted for (VAF) was considered as criterion for evaluating the appropriate number of these components, within a range from 2 to 15 ones, for best representation of the experimental data, as shown in Fig. 18 . As it can be seen, considering the cost-benefit balance for computational practical application of the proposed methodology, the appropriated number of unobservable components was ξ = 7, with VAF value of 99.92% in efficiency to represent as accurately as possible the experimental data and, at same time, reducing the computational load of interval type-2 fuzzy Kalman filter algorithm. These spectral unobservable components, which were extracted from daily deaths reports in state of Maranhão, are shown in Fig. 19 .

The partitions of experimental data related to daily deaths reports were defined by interval type-2 fuzzy Gustafson-Kessel clustering algorithm, as shown in Fig. 20 , so the antecedent proposition, the rules number and consequent proposition, of the type-2 fuzzy Kalman filter, could be estimated successfully. For implementing the proposed type-2 fuzzy clustering algorithm, the following parameters were adopted: number of clusters c = 3, interval weighting exponent m = [2.0, 2.3] and termination tolerance E = 10 −5 .

The implementation of interval type-2 fuzzy OKID algorithm, for parametric estimation of consequent proposition in the type-2 fuzzy Kalman filter inference system, in Eq. 15, took into account the partitions on daily deaths reports, in Fig. 20 , as weighting criterion, and the parameters values: q = 1, γ = 5 e β = 5. According to experimental data of daily deaths reports in state of Maranhão shown in Fig. 17 , the pre-processed unobservable components shown in Fig. 19 and the interval type-2 fuzzy normalized membership values shown in Fig. 20 , the initial parametric estimation of the type-2 fuzzy Kalman filter was computed by training step. The confidence region, as shown in Fig. 21 , created by initial estimation of interval type-2 fuzzy Kalman filter taking into account uncertainties estimated by interval type-2 membership functions shown in Fig. 20 , inherited to experimental data ranging from 18 of March 2020 to 18 of May 2020, illustrates its efficiency for tracking the experimental data of daily deaths reports in the state of Maranhão. From this confidence region, shown in Fig. 21 , interval normal distribution projections were estimated, delimiting upper and lower limits for forecasting the future daily deaths reports in the state of Maranhão. The efficiency of interval type-2 fuzzy Kalman filter based on its initial estimation by training step from experimental data of daily deaths reports ranging from 18 of March 2020 to 18 of May 2020, for forecasting the future (validation) experimental data of daily deaths reports within the period ranging from 19 of May 2020 to 27 of May 2020, is shown in Fig. 22a . From 05 of May 2020 to 17 of May 2020, it was established a lockdown in metropolitan city of São Luís, and other four municipalities as well, in state of Maranhão, with positive reflexes on the numbers of confirmed cases and daily deaths reports. Due to this social context changing established by government in state of Maranhão, it was necessary a recursive updating of interval type-2 fuzzy Kalman filter on 24 of June 2020, for proposals of new tracking and forecasting future daily deaths reports within the period ranging from 25 of June 2020 to 03 of July 2020, as shown in Fig. 22b . The effect of lockdown can be observed within the period ranging from 28 of May 2020 to 01 of June 2020. The new recursive updating of the interval type-2 fuzzy Kalman filter was on 28 of July 2020, for proposals of new tracking and forecasting future daily deaths reports within the period ranging from 29 of July 2020 to 25 of August 2020, as shown in Fig. 22c . The new recursive updating of the interval type-2 fuzzy Kalman filter was on 25 of August 2020, for proposals of new tracking and forecasting future daily deaths reports within the period ranging from 26 of August 2020 to 25 of September 2020, as shown in Fig. 22d . The new recursive updating of the interval type-2 fuzzy Kalman filter was on 25 of September 2020, for proposals of new tracking and forecasting future daily deaths reports within the period ranging from 26 of September 2020 to ahead, as shown in Fig. 22e . It can be seen an efficiency in the adaptability of interval normal distribution projections created in real time by interval type-2 fuzzy Kalman filter, which illustrates its applicability for tracking and forecasting the COVID-19 dynamic spread experimental data related to 

The experimental data corresponding to daily deaths reports within the period ranging from 29 of February 2020 to 18 of May 2020, in Brazil, are shown in Fig. 30 , which were extracted from official report by Ministry of Health of Brazil 2 . Once that the problem of interest, in this paper, is based on the time series related to daily deaths reports in Brazil, the variable u k , in Eq. 15 of proposed methodology, is considered as white noise signal with low amplitude. The pre-processing of experimental data by singular spectral analysis was able to extract the unobservable components associated to daily deaths reports. The variance accounted for (VAF) was considered as criterion for evaluating the appropriate number of these components, within a range from 2 to 15 ones, for best representation of the experimental data, as shown in Fig. 31 . As it can be seen, considering the cost-benefit balance for computational practical application of the proposed methodology, the appropriated number of unobservable components was ξ = 10, with VAF value of 99.98% in efficiency to represent as accurately as possible the experimental data and, at same time, reducing the computational load of interval type-2 fuzzy Kalman filter algorithm. These spectral unobservable components, which were extracted from daily deaths reports in Brazil, are shown in Fig. 32 .

The partitions of experimental data related to daily deaths reports were defined by interval type-2 fuzzy Gustafson-Kessel clustering algorithm, as shown in Fig. 33 , so the antecedent proposition, the rules number and consequent proposition, of the type-2 fuzzy Kalman filter, could be estimated successfully. For implementing the proposed type-2 fuzzy clustering algorithm, the following parameters were adopted: number of clusters c = 3, interval weighting exponent m = [1.5, 2.3] and termination tolerance E = 10 −5 .

The implementation of interval type-2 fuzzy OKID algorithm, for parametric estimation of consequent proposition in the type-2 fuzzy Kalman filter inference system, in Eq. 15, took into account the partitions on daily deaths reports, in Fig 33, as weighting criterion, and the parameters values: q = 1, γ = 15 e β = 15. According to experimental data of daily deaths reports in Brazil shown in Fig. 30 , the pre-processed unobservable components shown in Fig. 32 and the interval type-2 fuzzy normalized membership values shown in Fig. 33 , the initial parametric estimation of the type-2 fuzzy Kalman filter was computed by training step. The confidence region, as shown in Fig. 34 , created by initial estimation of interval type-2 fuzzy Kalman filter taking into account uncertainties estimated by interval type-2 member-ship functions shown in Fig. 33 , inherited to experimental data ranging from 29 of February 2020 to 18 of May 2020, illustrates its efficiency for tracking the experimental data of daily deaths reports in Brazil.

From this confidence region, shown in Fig. 34 , interval normal distribution projections were estimated, delimiting upper and lower limits for forecasting the future daily deaths reports in Brazil. The efficiency of interval type-2 fuzzy Kalman filter based on its initial estimation by training step from experimental data of daily deaths reports ranging from 29 of February 2020 to 18 of May 2020, for forecasting the future (validation) experimental data of daily deaths reports within the period ranging from 19 of May 2020 to 27 of May 2020, is shown in Fig. 35a . The new recursive updating of interval type-2 fuzzy Kalman filter on 24 of June 2020, for proposals of new tracking and forecasting future daily deaths reports within the period ranging from 25 of June 2020 to 30 of June 2020, is shown in Fig. 35b . The new recursive updating of the interval type-2 fuzzy Kalman filter was on 23 of July 2020, for proposals of new tracking and forecasting future daily deaths reports within the period ranging from 24 of July 2020 to 13 of August 2020, as shown in Fig. 35c . The new recursive updating of the interval type-2 fuzzy Kalman filter was on 28 of August 2020, for proposals of new tracking and forecasting future daily deaths reports within the period ranging from 29 of August 2020 to 25 of September 2020, as shown in Fig. 35d . The new recursive updating of the interval type-2 fuzzy Kalman filter was on 25 of September 2020, for proposals of new tracking and forecasting future daily deaths reports within the period ranging from 26 of September 2020 to ahead, as shown in Fig. 35e . It can be seen an efficiency in the adaptability of interval normal distribution projections created in real time by interval type-2 fuzzy Kalman filter, which illustrates its applicability for tracking and forecasting the COVID-19 dynamic spread experimental data related to daily deaths reports in Brazil. The estimation of interval type-2 fuzzy Kalman gain matrices K i | i=1,...,3 , during recursive updating of interval type-2 fuzzy Kalman filter, within period ranging from 19 of May 2020 to 25 of September, in Brazil, is shown in Fig. 36 . The recursive estimations of the interval type-2 fuzzy matrices A i , B i , C i and D i , with i = 1, . . . , 3, in the consequent proposition of the interval type-2 fuzzy Kalman filter inference system, during recursive updating of interval type-2 fuzzy Kalman filter, within period ranging from 19 of May 2020 to 25 of September, are shown in Figs. 37, 38, 39, 40 . The upper and lower instantaneous activation degrees related to interval type-2 fuzzy Kalman filter inference system, during its training and recursive step within period ranging from 29 of February 2020 to 25 of September, in Brazil, are shown in Fig. 41 . The efficiency of interval type-2 fuzzy Kalman filter, during its recursive updating for tracking and forecasting the COVID-19 dynamic spread experimental data related to daily deaths reports within period ranging from 19 of May 2020 to 29 of September, in Brazil, was validated through Variance Accounted For (VAF) criterion, as shown in Fig. 42 .

According to results from real-time interval tracking and forecasting the COVID-19's dynamic spread behavior, based on experimental data of daily deaths reports, in Brazil, the following auxiliary information, and helpful for Ministry of Health of Brazil, it can be inferred that:

-The proposed methodology can be helpful to help the Ministry of Health of Brazil as auxiliary information in the diagnosis regarding the flexibility of social activities, in Brazil. Figure 35 depicts a time-based displacement of the interval projections as well as their amplitudes, as the daily deaths reports are processed by interval type-2 fuzzy Kalman filter. From the recursive updating on 25 of September, the forecasting interval projections indicate the month of January 2021 as more adequate for reevaluating the requirements on flexibility of social activities, in Brazil.

In this paper, an approach for experimental data-based design of interval type-2 fuzzy Kalman filters was proposed. The computational results showed efficiency of designed interval type-2 fuzzy Kalman filter, due to its spectral pre-processing analysis mechanism, for filtering and tracking data set of Lorenz's chaotic attractor in a noisy environment, as compared to another approach widely cited in the literature. The experimental results showed the applicability of designed interval type-2 fuzzy Kalman filter, due to its recursive updating mechanism, for adaptive and real-time forecasting the COVID-19 spread dynamics related to daily deaths reports in state of Maranhão and Brazil. From the recursive updating on 25 of September 2020, it could be inferred by forecasting interval projections the month of January 2021 as more adequate for reassessment the requirements on flexibility of social activities in Brazil (containing the 27 states and federal district) and the month of November as more adequate for reassessment the requirements on flexibility of social activities in state of Maranhão.

For further works, the formulation and applicability of proposed methodology in the context of evolving interval type-2 fuzzy systems are of particular interest.

(2020b) Infection prevention and control for the safe management of a dead body in the context of covid-19: interim guidance

A new hybrid enhanced local linear neuro-fuzzy model based on the optimized singular spectrum analysis and its application for nonlinear and chaotic time series forecasting

Stability and control of dynamical systems with applications

Fuzzybased parameter optimization of adaptive unscented kalman filter: Methodology and experimental validation

Fuzzy modeling for control

Torque ripple minimization in predictive torque control method of PMSM drive using adaptive fuzzy logic modulator and EKF estimator

Interval-valued membership function estimation for fuzzy modeling. Fuzzy Sets and Systems

Hydroxychloroquine vs. azithromycin for hospitalized patients with suspected or confirmed COVID-19 (HAHPS): Protocol for a pragmatic, open label, active comparator trial

COVID-19 outbreak: Migration, effects on society, global environment and prevention. Science of The Total Environment

Deterministic annealing gustafson-kessel fuzzy clustering algorithm

Linear system theory and design

Epidemiological features and time-series analysis of influenza incidence in urban and rural areas of shenyang, china

Robust genetic algorithm and fuzzy inference mechanism embedded in a sliding-mode controller for an uncertain underwater robot

COVID-19 virus outbreak forecasting of registered and recovered cases after sixty day lockdown in italy: A data driven model approach

Real-time forecasting of epidemic trajectories using computational dynamic ensembles. Epidemics, 30, 100379

Global transmission and evolutionary dynamics of the chikungunya virus

ARIMA modelling and forecasting of irregularly patterned COVID-19 outbreaks using japanese and south korean data

Analysis of time series structure: SSA and related techniques

Multivariable statespace recursive identification algorithm based on evolving type-2 neural-fuzzy inference system

Hybrid learning for interval type-2 intuitionistic fuzzy logic systems as applied to identification and prediction problems

COVID-19 with different severities: A multicenter study of clinical features

Erythropoiesis-stimulating agent treatment in patients with COVID-19

Digital Control of Dynamic Systems

Novel corona virus disease infection in tunisia: Mathematical model and the impact of the quarantine strategy

Adaptive neuro-fuzzy control for discrete-time nonaffine nonlinear systems

Singular spectrum analysis for time series

Lossy signal-interference filters and applications

Modeling the 2016-2017 yemen cholera outbreak with the impact of limited medical resources

Singular spectrum analysis of biomedical signals

Predicting ebola virus disease risk and the role of african bat birthing

Novel grid multiwing butterfly chaotic attractors and their circuit design

SARS-CoV-2 viral load in clinical samples from critically ill patients

On-line obstacle detection, avoidance, and mapping of an outdoor quadrotor using EKF-based fuzzy tracking incremental control

Fuzzy cluster analysis: methods for classification, data analysis and image recognition

Applied system identification

A new approach to linear filtering and prediction problems

Estimation of reproduction number (ro) and early prediction of 2019 novel coronavirus disease (COVID-19) outbreak in india using statistical computing approach

Nowcasting the number of new symptomatic cases during infectious disease outbreaks using constrained p-spline smoothing

Optimal neuro-fuzzy control of hepatitis c virus integrated by genetic algorithm

Decentralized h ∞ sampled-data fuzzy filter for nonlinear interconnected oscillating systems with uncertain interconnections

Prediction of COVID-19 transmission dynamics using a mathematical model considering behavior changes. Epidemiology and Health, e2020026

The forecasting of dynamical ROSS river virus outbreaks: Victoria, Australia. Epidemics, 30, 100377

Detection of norovirus infections in denmark

Interval type-2 fuzzy logic systems: theory and design

A conceptual model for the coronavirus disease 2019 (COVID-19) outbreak in wuhan, china with individual reaction and governmental action

A robust localization method for unmanned surface vehicle (USV) navigation using fuzzy adaptive kalman filtering

Introduction to dynamic systems: theory, models, and applications

Deep filtering: Signal extraction and reconstruction using complex time-frequency filters

Computational models for trapping ebola virus using engineered bacteria

The fuzzy kalman filter: Improving its implementation by reformulating uncertainty representation. Fuzzy Sets and Systems

Comparing the performance potentials of interval and general type-2 rule-based fuzzy systems in terms of sculpting the state space

Unravelling the myths of r0 in controlling the dynamics of COVID-19 outbreak: A modelling perspective

Discrete-time kalman filter for takagi-sugeno fuzzy models

The time scale of asymptomatic transmission affects estimates of epidemic potential in the COVID-19 outbreak

Automated detection of white matter changes in elderly people using fuzzy, geostatistical, and information combining models

Methodology for evolving fuzzy kalman filter identification

Using routine testing data to understand circulation patterns of influenza a, respiratory syncytial virus and other respiratory viruses in victoria, australia

Fuzzy system identification and adaptive control

Control of malaria outbreak using a non-linear robust strategy with adaptive gains

COVID-19 future forecasting using supervised machine learning models

An interim review of the epidemiological characteristics of 2019 novel coronavirus

Innovations in fuzzy clustering: theory and applications

Frontiers in advanced control systems

Kalman filters: Theory for advanced applications

Spatial and temporal clustering of patients hospitalized with laboratory-confirmed influenza in the united states

Dynamic modeling of hepatitis c transmission among people who inject drugs. Epidemics, 30, 100378

A qualitative study on the psychological experience of caregivers of COVID-19 patients

A course in fuzzy systems and control

Tracking a maneuvering target by multiple sensors using extended kalman filter with nested probabilistic-numerical linguistic information

Robust guaranteed cost sampled-data fuzzy control for uncertain nonlinear time-delay systems

Robust economic model predictive control of continuous-time epidemic processes

Time series analysis and forecasting of chlamydia trachomatis incidence using surveillance data from 2008 to 2019 in shenzhen, china

Design of hidden-property-based variable universe fuzzy control for movement disorders and its efficient reconfigurable implementation

Wavelet-TSK-type fuzzy cerebellar model neural network for uncertain nonlinear systems

Early prediction of the 2019 novel coronavirus outbreak in the mainland china based on simple mathematical model

Signal detection in generalized gaussian distribution noise with nakagami fading channel

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

Acknowledgements This work was developed in the Laboratory of Computational Intelligence Applied to Technology at Federal Institute of Education, Science and Technology of Maranhão. The authors are grateful to Coordination for the Improvement of Higher Education Personnel (CAPES) for financial support and to the Master and Doctorate Program in Electrical Engineering at Federal University of Maranhão (PPGEE-UFMA) for their support in the development of this research.