Abstract

This paper considers the parameter identification of Wiener systems with colored noise. The difficulty in the identification is that the model is nonlinear and the intermediate variable cannot be measured. Particle swarm optimization is an artificial intelligence evolutionary method and is effective in solving nonlinear optimization problem. In this paper, we obtain the identification model of the Wiener system and then transfer the parameter identification problem into an optimization problem. Then, we derive a particle swarm optimization iterative (PSOI) identification algorithm to identify the unknown parameter of the Wiener system. Furthermore, a gradient iterative identification algorithm is proposed to compare with the particle swarm optimization iterative algorithm. Numerical simulation is carried out to evaluate the performance of the PSOI algorithm and the gradient iterative algorithm. The simulation results indicate that the proposed algorithms are effective and the PSOI algorithm can achieve better performance over the gradient iterative algorithm.

1. Introduction

Almost all practical systems are nonlinear [13]. Many identification methods have been developed for linear systems [4, 5], bilinear systems [68], and nonlinear systems [9]. The Wiener models are a typical class of nonlinear systems and are widely used in industrial production process [10, 11]. The Wiener nonlinear system consists of a dynamic linear subsystem and a static nonlinear subsystem and has the characteristics of complex structure between subsystems [12, 13]. One of the difficulties in identifying Wiener nonlinear model parameters is that the intermediate variable (the output of the linear subsystem) cannot be measured, and the identification issues for Wiener systems have attracted great attention [14].

The iterative identification method is generally used to identify the system with unknown item in the model information vector [1517]. The basic idea of iterative identification is to estimate the unknown items in the information vector by using the iterative parameter estimation of the previous step [18, 19]. The iterative identification method is an important branch of system identification, which can be realized by using gradient search, least squares principle, and Newton optimization [2022].

The particle swarm optimization algorithm is an evolutionary computing technique which is based on the simulation of birds’ flock [23, 24]. The basic idea of particle swarm optimization algorithm is to find the optimal solution through collaboration and information sharing among individuals in the group [25]. This algorithm has attracted the attention of academia with the advantages of easy implementation, high precision, and fast convergence [26]. Compared with the conventional optimization methods, it has excellent optimized performances and characteristics [27]. The particle swarm optimization algorithm has been widely used in function optimization, system identification, and fuzzy control [2830]. Recently, Chen and Wang proposed a stochastic gradient algorithm and a particle swarm optimization algorithm to estimate all the unknown parameters of the Hammerstein system [31]. In this paper, we use the particle swarm optimization algorithm and the gradient iterative algorithm to identify the unknown parameters of the Wiener systems with colored noise.

The rest of the paper is organized as follows. Section 2 gives the system description for the Wiener model. Section 3 gives the particle swarm optimization algorithm for Wiener nonlinear systems. Section 4 derives a gradient iterative algorithm for the discussed system. Section 5 provides an example for illustrating the results in this paper. Finally, some conclusions are given in Section 6.

2. System Description

Consider the Wiener system shown in Figure 1 with the following expressions: where , , and are polynomials in the shift operator with

Assume that the degrees and are known and and for

Define the linear subsystem output as and the noise model output as

The static nonlinear block is a nonlinear function where the basis are known nonlinear functions of , the unknown parameters are the coefficients of the nonlinear functions and assume that the degree is known. Without loss of generality, let the first coefficient of nonlinear block be unity and rewrite the as

In the above equations, and are the system input and output, respectively, and is a Gaussian distributed white noise with zero mean and variance . From (3), we have where

From (4), we can obtain where

Thus, the Wiener nonlinear system model can be written as follows: where

3. The Particle Swarm Optimization Algorithm

With the development of optimization theory, some new intelligent algorithms have been proposed to solve the problem of traditional system identification, such as the genetic algorithm [32], the ant colony algorithm [33], and the particle swarm algorithm [34, 35], these algorithms enrich the system identification technology. Particle swarm optimization algorithm is a nature-inspired evolutionary algorithm, and it has been successful in solving a wide range of real-value optimization problems [36]. In the following, the particle swarm optimization algorithm is used to identify the unknown parameters of the Wiener nonlinear systems.

Suppose that the search space is -dimensional and a particle swarm consists of particles.

Define the information vector as

Let represent the data length. Define the stacked output vector and the stacked information matrix as

Define the independent position of each particle and the independent velocity as follows:

Let denote the estimates of at iteration Define as the best position of each particle at iteration

Let and denote the estimates of and at iteration

Then, the estimates can be obtained as follows:

According to the basic principle of the particle swarm algorithm, the best position of each particle satisfies the following cost function:

Let denote the global best position of all the particles where satisfies

Define According to (7), we can obtain the estimation of

According to the principle of particle swarm optimization, each particle goes to a new position and a new velocity at iteration as follows:

Replacing and with and in (11), we can obtain the estimate Thus, we can obtain the particle swarm optimization iterative (PSOI) identification algorithm as follows:

The steps of the PSOI algorithm are listed as follows: (1)Let set the initial values as , , and Set the initial factor and give a small positive number . Set , , and .(2)Collect the input and output data and , , form by (27) and by (26). Construct and by (31) and (32), respectively.(3)Update the velocity of each particle , according to (25).(4)Update the position of each particle by (24).(5)Compute the best position of each particle by (33).(6)Determine the best position of all the particles by (34).(7)Compute by (36). Form by (35), compute by (28).(8)Compare and : if then terminate the procedure and obtain the estimate ; otherwise, increase by 1 and go to Step 2.

The flowchart of PSOI algorithm is shown in Figure 2.

Remark 1. The major factors that influence the performance of the particle swarm optimization include and . and are positive constants between 0 and 2. is the step size that adjusts the particle to its own best position. is the step size that regulates the particle to the global best position. is called the inertia factor and is an important adjusting parameter of the PSOI algorithm. A larger can facilitate global optimization; otherwise, a smaller one can facilitate local optimization. It can be chosen as a constant between 0.1 and 0.9 generally. and are two independent random numbers uniformly distributed in the range of [0, 1].

4. Gradient Iterative Algorithm

The gradient search is a very basic and ancient search method [37, 38]. It is widely used in parameter identification of nonlinear systems [3941]. In the following, based on the gradient search principle, a gradient iterative identification algorithm for Wiener nonlinear model is derived.

Consider the latest group data from to and define the stacked output vector the stacked information matrix , and the stacked noise vector as follows:

From (11), we have

Define the criterion function

Let as an iterative variable and is the th iterative estimation of parameter vector at time For the optimization problem (39), the gradient iterative algorithm is obtained by using the negative gradient search where is the iterative step-size. However, in the upper formula of (40), the gradient iterative estimate impossible to calculate because the stacked information vector contains unknown intermediate variables and The solution is to replace the unknown variables and by and , respectively. Let , , and denote the estimates of , , and at iteration , respectively

Let denote the estimates of at iteration and let denote the estimates of at iteration

Thus, can be calculated by the following:

Define as the estimates of at iteration and then the estimates of can be obtained by the following:

Replacing by in (40) and rewriting give

Thus, we can obtain the gradient iterative (GI) estimation algorithm for Wiener models

The steps of the GI algorithm are listed as follows: (1)Let , , and and give a small positive number .(2)Collect the input and output data and and form by (49) and by (51).(3)Form by (56), compute by (55), and form by (52).(4)Compute by (54) and form by (53).(5)Form by (50) and form by (48).(6)Choose a suitable by (57) and compute by (47).(7)Compare and if then terminate the procedure and obtain ; otherwise, increase by 1 and go to Step 2.

5. Examples

Consider the following Wiener nonlinear systems:

In simulation, the input is taken as an uncorrelated stochastic signal sequence with zero mean and unit variance and as a Gaussian white noise sequence with zero mean and variance Applying the GI algorithm and the PSOI algorithm to estimate the parameters of this system, the parameter estimates and their errors are shown in Tables 1 and 2 and Figures 3 and 4. In the PSOI algorithm simulation, each swarm generation contains 10 position particles. The coefficients and in (25) are set to 0.8, 1.2, and 1.8, respectively.

From the simulation results in Tables 1 and 2 and Figures 3 and 4, we can draw the following conclusions: (i)As increases, the parameter estimation errors given by the GI algorithm and PSOI algorithm gradually become smaller (see Tables 1 and 2).(ii)The PSOI algorithm has a faster convergence rate than the GI algorithm (see Figures 3 and 4).(iii)The PSOI algorithm has a higher estimation accuracy than the GI algorithm, which can be seen from Tables 1 and 2.

6. Conclusions

In this paper, we derived the particle swarm optimization iterative algorithm and the gradient iterative algorithm for Wiener nonlinear systems. Compared with the gradient iterative algorithm, the particle swarm optimization algorithm has a higher estimation accuracy and has a faster convergence rate. The proposed approaches in the paper can be combined with other mathematical tools [4247] to study the performances of some parameter estimation algorithms and can be applied to other multivariable systems with different structures and disturbance noises and other literature [4852].

Data Availability

The relevant data can be obtained by email to the author.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China (61403217) and the Jiangsu Province Postdoctoral Research Funding Plan (1601129B).