Next Article in Journal
Inequalities in Triangular Norm-Based ∗-fuzzy ( L + ) p Spaces
Next Article in Special Issue
Complex Ginzburg–Landau Equation with Generalized Finite Differences
Previous Article in Journal
Optimal Replenishment Policy for Deteriorating Products in a Newsboy Problem with Multiple Just-in-Time Deliveries
Previous Article in Special Issue
A Meshless Method Based on the Laplace Transform for the 2D Multi-Term Time Fractional Partial Integro-Differential Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

State Feedback Regulation Problem to the Reaction-Diffusion Equation †

by
Francisco Jurado
* and
Andrés A. Ramírez
División de Estudios de Posgrado e Investigación, Tecnológico Nacional de México/I.T. La Laguna, Revolución Blvd. and Instituto Tecnológico de La Laguna Av., Torreón 27000, Mexico
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in the Proceedings of the 2019 16th International Conference on Electrical Engineering, Computing Science and Automatic Control (CCE), Mexico City, Mexico, 11–13 September 2019.
Mathematics 2020, 8(11), 1983; https://doi.org/10.3390/math8111983
Submission received: 13 September 2020 / Revised: 31 October 2020 / Accepted: 3 November 2020 / Published: 6 November 2020
(This article belongs to the Special Issue Applications of Partial Differential Equations in Engineering)

Abstract

:
In this work, we explore the state feedback regulator problem (SFRP) in order to achieve the goal for trajectory tracking with harmonic disturbance rejection to one-dimensional (1-D) reaction-diffusion (R-D) equation, namely, a partial differential equation of parabolic type, while taking into account bounded input, output, and disturbance operators, a finite-dimensional exosystem (exogenous system), and the state of the exosystem as the state to the feedback law. As is well-known, the SFRP can be solved only if the so-called Francis (regulator) equations have solution. In our work, we try with the solution of the Francis equations from the 1-D R-D equation following given criteria to the eigenvalues from the exosystem and transfer function of the system, but the state operator is here defined in terms of the Sturm–Liouville differential operator (SLDO). Within this framework, the SFRP is then solved for the 1-D R-D equation. The numerical simulation results validate the performance of the regulator.

1. Introduction

Some physical quantities from applications involving the diffusion and structural vibrations depend on both position and time. Systems whose dynamic evolves in an infinite-dimensional Hilbert space are modeled by partial differential equations (PDEs). Such systems are called infinite-dimensional systems. Because these systems reflect the spatial distribution of a physical quantity, these systems are called distributed parameter systems (DPSs). The aim when designing the controller for infinite-dimensional systems is that the control system is stable and robust in the presence of parametric uncertainties and external perturbations.
The design of classical controllers is commonly based on the transfer function of the system. Infinite-dimensional systems have transfer function. Unlike transfer functions from finite-dimensional systems, transfer functions of infinite-dimensional systems are not rational functions. If the transfer function from an infinite-dimensional system is provided, then the controller can be directly designed. A drawback from this last approach is that the controller will be infinite-dimensional and it may be approximated by a finite-dimensional system. For some practical applications, it could be that a transfer function is not available. Subsequently, a finite-dimensional model of the system must be drawn from which the controller design may be based. This latter approach is the most used methodology to design controllers of systems modeled by PDEs [1].
The main aim in the regulator problem is the design of a feedback law, such that the output of the system tracks a reference signal in presence of a external disturbance, where the reference signal and external disturbances are generated through an exosystem (also called exogenous system). The state feedback regulation problem (SFRP) and error feedback regulation problem (EFRP) are two approaches that may be used to solve the output regulation problem [2]. The regulator theory has a long history and has been used in both finite- and infinite-dimensional cases. In this work, we try with the SFRP under the assumption that the whole state from the exosystem is available to the controller.
In our work, the employed methodology follows the work for linear finite-dimensional systems [3,4,5,6]. In [5], the solvability for the multivariable linear regulator problem, which tries with the solvability for a pair of linear matrix equations, namely, the Francis (regulator) equations, was shown. Necessary and sufficient conditions for the solvability of the so-called Hautus equations were given in [7]. In [8] the Hautus conditions were stated, for finite-dimensional linear systems, claiming that any eigenvalue from the exosystem is not an invariant zero of the control system. In [3] the regulator problem was extended to finite-dimensional nonlinear systems where necessary and sufficient conditions, about the solvability of a pair of nonlinear regulator equations, are given to the solvability of the problem.
The regulator problem for linear finite-dimensional systems was extended to linear DPSs [9,10,11,12,13]. Control systems governed by a discrete spectral operator are presented in [11,12], where the so-called state operator comply with the property of spectrum decomposition [14,15] from which a controllability condition is determined by the spectrum, implying the stabilizability of the control system through a finite-dimensional controller. In [13], the geometric regulation approach was extended to linear DPSs with bounded input and output operators where both reference and disturbance signals considered in [11,12] are assumed as generated by a finite-dimensional exosystem. Additionally, simple criteria for the solvability of the Francis (regulator) equations were provided. Moreover, the regulator problem to set-point control with harmonic disturbance, such disturbance acting over a small spatial interval, for a one-dimensional (1-D) heat equation, the regulator problem for harmonic tracking of a damped wave equation, when considering full state feedback for both cases, as well as the regulator problem for harmonic tracking of a 1-D heat equation, but under the error feedback approach, these latter without considering disturbances, were solved. Neumann boundary conditions were considered to the cases mentioned above about the regulator problem for a 1-D heat equation. The case for regulation of linear systems with bounded input, output and disturbance operators, this last disturbance operator entering across the entire interval, was considered in [16] for a 1-D heat equation with Dirichlet condition at one boundary and Neumann condition at the other boundary. In [17] key results from [13], regarding the state feedback regulator problem, were extended to systems with unbounded control and observation operators. The regulator theory developed in such a work was applied to a one-dimensional (1-D) heat equation, a heat equation on a two-dimensional (2-D) unit rectangle, and on a damped Rayleigh beam.
Comprehensive reviews regarding contributions made to generalize the regulator problem to infinite-dimensional systems can be found in [17,18]. The output regulation problem for DPSs has been extensively studied for different classes of PDEs systems; for example, parabolic and hyperbolic PDEs [18,19,20], linearized Stokes flow, and steady-state incompressible Navier–Stokes fluid flow equations with spatially varying parameters [16], coupled wave equation [16], wave equations [21,22], beam and plate equations [16,17,23,24], and partial integro-differential equations [25]. In most of the works, the controller is done based on output feedback; these include system stabilization and availability of the state of the exosystem for feedback is not required.
Reaction-diffusion (R-D) equations have been playing an important role in applied mathematics. Dynamics from combustion processes, superconductivity in liquids, signal transmission across neurons, population growth, mass balance and mass transport, chemical kinetics, chemotaxis, metabolic processes in living cells, and the spread of diseases and viruses can be analyzed and simulated by R-D models [26,27,28,29,30,31,32,33,34]. In [35], boundary control to the R-D equation via the backstepping method was proposed. The main idea of using the backstepping method was about the use of a Volterra integral transformation along with feedback control that made that the R-D system behaves as the heat equation; this latter as an exponentially stable target system. To the best of our knowledge, there is no work regarding the solution of the SFRP for harmonic tracking with harmonic disturbance rejection to the 1-D R-D equation.
In this work, although it could be argued that the 1-D R-D equation is be among the simplest PDEs systems, our application is concerned with the SFRP to the design of a harmonic tracking regulator, but, in contrast with our previous work recently reported in [36], with harmonic disturbance rejection for this kind of system. It is well known that the SFRP is solved only if the Francis (regulator) equations have a solution. In our work, the main contribution is that the state operator to the Francis equations, corresponding to the 1-D R-D equation, here is defined in terms of the Sturm–Liouville differential operator (SLDO) to then solve them.
The manuscript is organized, as follows. In Section 1, the SFRP, modeling of DPS, control of PDEs, and applications of the R-D equation are summarized; the problem statement is formulated in Section 2; the regulator design is detailed at Section 3; the simulation results are included in Section 4; and, the conclusion is drawn at the end.

2. Problem Statement

2.1. Sturm–Liouville Boundary Value Problem

Let us consider a differential equation of the form
( p ( x ) y ) q ( x ) y + λ r ( x ) y = 0 , x ( 0 , 1 ) ,
with symmetric boundary conditions, also known as Sturm–Liouville boundary conditions,
α 1 y ( 0 ) + α 2 y ( 0 ) = 0 , β 1 y ( 1 ) + β 2 y ( 1 ) = 0 ,
at the endpoints. The Sturm–Liouville Equation (1) can be rewritten as
x i p ( x ) y x i q ( x ) y + λ r ( x ) y = 0
or, equivalently,
p ( x ) Δ y q ( x ) y + λ r ( x ) y = 0 ,
with Δ the Laplacian operator [26,37], assuming that p ( x ) is a C 1 function, q ( x ) and r ( x ) are continuous functions, p ( x ) > 0 , and r ( x ) > 0 [38]. In this case, the boundary value problem is said to be regular. It should be noticed that these latter conditions are satisfied in many problems in mathematical physics.
L = d d x p ( x ) d d x + q ( x )
is known as the SLDO [39].
Let us consider the linear homogeneous differential operator
L [ y ] = ( p ( x ) y ) + q ( x ) y .
Accordingly, the Equation (1) can be written as
L [ y ] = λ r ( x ) y .
Each one of the boundary conditions (2) involves only one of the boundary points, and that is why are said to be separated. In order to impose any restriction on y, at least one of the parameters α 1 and α 2 must be different from zero to the first boundary condition. Alike, at least one of the parameters β 1 and β 2 must differ from zero [40]. The boundary value problem consisting of (1) and (2) is called the Sturm-Liouville boundary value problem (SLBVP) [41].
Modeling that involves ordinary differential equations (ODEs) usually leads to initial value problems or boundary value problems, these latter written in the form of the SLBVPs, which are eigenvalue problems that involve a parameter λ related to frequencies, energies, or other physical quantities [42]. In fact, the next theorem summarizes the way about how to convert a linear second order differential equation A ( x ) y + B ( x ) y + C ( x ) y into one of the form ( 1 / r ( x ) ) ( ( p ( x ) y ) + q ( x ) y ) . If p ( x ) , q ( x ) , and r ( x ) are given, then A ( x ) = p ( x ) / r ( x ) , B ( x ) = p ( x ) / r ( x ) , and C ( x ) = q ( x ) / r ( x ) .
Theorem 1.
[43] Suppose that A ( x ) > 0 , B ( x ) and C ( x ) are analytic real-valued functions in the (finite or infinite) interval a < x < b , then exist functions p ( x ) > 0 , q ( x ) , and r ( x ) > 0 likewise analytic and real valued in a < x < b and
A ( x ) y + B ( x ) y + C ( x ) y = 1 r ( x ) ( p ( x ) y ) + q ( x ) y
identically in y.
Proof. □ 
See [43].
The expression on the right from this last equation is referred as the Sturmian form of the differential equation, also called self-adjoint form. Solutions of SLBVPs, called eigenfunctions, have many general properties in common, as the orthogonality property, useful in eigenfunctions expansions in terms of Fourier series, Legendre polinomials, Bessel functions, and other eigenfunctions [41,44].

2.2. Abstract Control System Model

Let us consider an abstract distributed parameters control system in the Hilbert space Z in the form
z t = A z ( t ) + B i n u ( t ) + B d d ( t ) ,
z ( 0 ) = z 0 , z 0 Z ,
y ( t ) = C z ( t ) ,
where the state operator A refers to an unbounded densely defined operator, z ( t ) stands for the state of the system, z 0 is for the initial condition, z t is the derivative of z ( t ) with respect to time, u ( t ) U is the control input, and y ( t ) Y is the measured output. U , Y are Hilbert spaces of either finite or infinite dimension. B i n is an input operator, B d L ( U , Z ) is a disturbance operator, d ( t ) refers to the disturbance, and C is an output operator.
The operator A is given in terms of a linear elliptic partial differential operator L in the Hilbert space Z = L 2 ( Ω ) of infinite dimension with bounded domain Ω R n with piecewise smooth boundary. To the second order case,
L = i , j = 1 n x j a i j ( x ) x i a ( x )
where a ( x ) 0 for all x Ω , a i j C ( Ω ¯ ) , a i j ( x ) = a j i ( x ) , and Ω ¯ is the closure of Ω . If there exist constants 0 < c 1 < c 2 < , by uniform ellipticity it means that
c 1 | ξ | 2 i , j = 1 n a i j ( x ) ξ i ξ j c 2 | ξ | 2
for all ξ R n and x Ω ¯ , where | · | represents the Euclidean norm in R n .
In our work, from (5) and (11), the state operator is defined here as given in the form of the SLDO, namely,
A = A + λ ,
with
A = d 2 / d x 2 .
The output operator C L ( Z , Y ) is a set of bounded operators C i given by
y i ( t ) = C i z = 1 | Ω i | Ω i z ( x , t ) d x ,
for some Ω j of the domain Ω , with Lebesgue measure
| Ω i | = Ω i d x > 0 .
In the same setting,
y i ( t ) = C i z = z , Ψ i = Ω z ( x , t ) Ψ i ( x ) d x
with
Ψ i ( x ) = 1 | Ω i | 1 Ω i ( x ) L 2 ( Ω )
and indicator function
1 Ω i ( x ) = 1 , x Ω i 0 , x Ω i .
Accordingly,
y = C z = [ C 1 ( z ) C 2 ( z ) C n c ( z ) ] T .

2.3. Exogenous System (Exosystem)

Let us consider a finite-dimensional exosystem
d w ( t ) d t = S w ( t ) ,
y r ( t ) = Q w ( t ) ,
d ( t ) = P w ( t ) ,
w ( 0 ) = w 0 ,
which generates both reference output y r ( t ) and disturbance d ( t ) , W is the state space of the exosystem, S L ( W ) , Q L ( W , Y ) , and P L ( W , Z ) .
Let us consider the system given by
z t ( x , t ) = A z ( x , t ) ,
z ( x , 0 ) = z 0 ( x ) .
If there exist positive constants M and α , such that
z ( · , t ) = e A t z 0 M e α t z 0 t 0 ,
with z 0 Z , then the system (20) and (21) is exponentially stable. From the above, it is said that the state operator A generates an exponentially stable C 0 semigroup in Z [45], i.e., A is stable. Accordingly, here it is assumed that the uncontrolled system (20) and (21), i.e., to the case when u = 0 and d = 0 , is exponentially stable. Additionally, it is assumed that the exogenous system (16)–(19) is neutrally stable. To the linear case, in the Lyapunov sense, this is equivalent to the origin being stable implying that σ ( S ) i R and S has no nontrivial Jordan blocks. ρ ( T ) refers to the resolvent set of an operator T and σ ( T ) refers to the spectrum of T.
Let us consider the error signal
e ( t ) = y ( t ) y r ( t ) ,
= C z ( t ) Q w ( t ) .
The main task for the regulator is to force the output of the system to track a reference signal in presence of a disturbance d ( t ) , i.e., e ( t ) 0 as t . Accordingly, the problem is stated as follows.
Problem 1.
State Feedback Regulator Problem (SFRP).
The SFRP consists in to find a control law
u ( t ) = Γ w ( t )
in function of the state of the exosystem with Γ L ( W , U ) , such that for the system
d z ( t ) d t = A z ( t ) + ( B d P + B i n Γ ) w ( t ) ,
d w ( t ) d t = S w ( t ) ,
which corresponds to the interconnection of (8)–(10) with (16)–(19), the error norm
e ( t ) = y ( t ) y r ( t ) ,
= C z ( t ) Q w ( t ) ,
satisfies
e ( t ) 0 as t ,
for any z 0 Z and w 0 W .
Being e ( t ) a finite dimension vector, then all l p norms in (30) are equivalent. Because exponential stability for the system (20) and (21) has been assumed, then a state feedback control law is not involved. In this work, our main focus is with tracking and harmonic disturbance rejection, so, the stabilization problem is out of scope of our proposal. Thus, the solvability of the SFRP is stated in the next theorem.
Theorem 2.
If there exist mappings Π L ( W , Z ) and Γ L ( W , U ) , with r a n k ( Π ) D ( A ) , satisfying the Francis (regulator) equations
Π S = A Π + B i n Γ + B d P ,
C Π = Q ,
then the linear SFRP has solution. The feedback control law that solves the SFRP is given by
u ( t ) = Γ w ( t ) .
Proof. 
The proof can be carried out along the same lines as in [16].
From this last theorem, from the assumptions that e A t is an exponentially stable semigroup and that the exosystem is neutrally stable, if (31) holds then e ( t ) 0 as t for all z 0 Z , w 0 W , if and only if [ C Π Q ] = 0 . □

3. Regulator Design

Let us consider the 1-D R-D equation that is given by
z t ( x , t ) = D z x x ( x , t ) + λ z ( x , t ) + B i n u + B d d ,
z ( 0 , t ) = 0 , ( Dirichlet BC )
z x ( 1 , t ) = 0 , ( Neumann BC )
z ( x , 0 ) = ϕ ( x ) ,
y ( t ) = C z ( t ) ,
where D z x x ( x , t ) is referred as the diffusion term, with diffusion coefficient (constant) D > 0 , and λ z ( x , t ) is the reaction (source) term, with λ an arbitrary constant. Refer to [37] to have an idea about typical values that D may take when trying with heat conduction or molecular diffusion.
The system (34)–(38) is formulated in the form (8)–(10) in the Hilbert state space Z = L 2 ( 0 , 1 ) . z x refers to the partial derivative with respect to space and z x x refers to the second partial derivative with respect to space also. The maximal elliptic operator is given by L = d 2 / d x 2 in (11) with domain D ( L ) = H 2 ( 0 , 1 ) , the Sobolev space of functions ζ Z with d ζ / d x both continuous on ( 0 , 1 ) and d 2 ζ / d x 2 Z . The state operator (12) is a self-adjoint operator in Z , i.e.,
A [ ζ ] = ζ + λ ζ ,
where
D ( A ) = { ζ H 2 ( 0 , 1 ) : ζ ( 0 ) = 0 , ζ ( 1 ) = 0 } Z .
The spectrum of A denoted by
σ ( A ) = { λ k } k = 0 ,
where λ k = μ k 2 with μ k = k 1 2 π , is purely discrete with a corresponding set of orthonormal eigenvectors
ζ k ( x ) = 2 sin ( μ k x ) , k = 1 , 2 , .
The operator A is assumed as an infinitesimal generator of an exponentially stable C 0 semigroup in Z in terms of the eigenfunction expansion
e A t ζ = j = 0 e λ j t ζ , ζ j ζ j .
It should be noticed that the system under consideration is a single-input/single-output system with scalar input and output B i n and C, respectively, and B d is for a bounded disturbance. Additionally, the input, output, and disturbance operators are bounded operators acting in the interior of the domain.
The input to the system is spatially uniform over a small interval about a fixed point x i n = x 0 ( 0 , 1 ) , where B i n u = b ( x ) u with
b ( x ) = 1 2 ν 0 1 [ x 0 ν 0 , x 0 + ν 0 ] ( x )
and
1 [ a , b ] ( x ) = 1 , x [ a , b ] 0 , x [ a , b ] .
The input operators B i n and B d are given as
B i n u ( t ) = j = 1 n i n B i n j u j ( t ) , B d d ( t ) = j = 1 n d B d j d j ( t ) ,
where u j ( t ) and d j ( t ) are scalar control inputs and disturbances, respectively. B i n j and B d j ( x ) are characteristic functions of a bounded subset of Ω , namely
B i n j ( x ) = 1 | Ω j | 1 Ω j ( x ) .
To guarantee that B i n j Z , here it is assumed that | Ω j | > 0 .
The output is the average transport reaction over a small interval about a point x o u t = x 1 ( 0 , 1 ) , i.e.,
C ϕ = 0 1 c ( x ) ϕ ( x ) d x
with
c ( x ) = 1 2 ν 1 1 [ x 1 ν 1 , x 1 + ν 1 ] ( x ) .
Because C ζ = ζ , c , C is a bounded linear observation functional on Z .
In our work, d ( t ) = A d sin ( β t ) R entering across the whole interval is considered, so B d = 1 .
In our proposal, the Francis Equations (31) and (32) take the form
Π S w = A Π w + B i n Γ w + B d P w ,
C Π w = Q w = w 1 ,
with
Π = [ Π 1 Π 2 Π 3 Π 4 ] , w ( t ) = w 1 w 2 w 3 w 4 = A r sin ( α t ) A r cos ( α t ) A d sin ( β t ) A d cos ( β t ) ,
S = 0 α 0 0 α 0 0 0 0 0 0 β 0 0 β 0 , Q = [ 1 0 0 0 ] , P = [ 0 0 1 0 ] .
Because of the block diagonal structure of S, the Francis equations can be decoupled into two separate parts, the first one to work with the harmonic tracking task and the last one to work with the rejection of a harmonic disturbance.

3.1. Harmonic Tracking

Accordignly, along the same lines from [16], for the first part on which it is desirable to track a harmonic signal y r ( t ) = A r sin ( α t ) consider an exosystem given as
w t = S w , w ( 0 ) = 0 A r , S = 0 α α 0 ,
with solution
w ( t ) = A r sin ( α t ) A r cos ( α t ) .
Taking P = [ 0 0 ] and Q = [ 1 0 ] then
y r ( t ) = Q w = A r sin ( α t ) .
In this case W = R 2 , so Π = [ Π 1 Π 2 ] and Γ = [ Γ 1 Γ 2 ] R 2 with Π j Z .
The Francis equations applied to the vector w = [ w 1 w 2 ] T W give the following system
Π S w = A Π w + B i n Γ w ,
C Π w = Q w .
From these last equations, the regulator equation on the left can be written as
α Π 1 w 2 α Π 2 w 1 = A Π 1 w 1 + A Π 2 w 2 + B i n Γ 1 w 1 + B i n Γ 2 w 2 .
Because (46) must be fulfilled for all w, let us consider first the case w 1 = 1 and w 2 = 0 and then w 1 = 0 and w 2 = 1 yielding
α Π 2 A Π 1 = B i n Γ 1 ,
α Π 1 A Π 2 = B i n Γ 2 .
Remarking that the exogenous system (16)–(19) is neutrally stable, multiplying (48) by i = 1 to then adding the result to (47) it yields
( i α I A ) Π 1 + i ( i α I A ) Π 2 = B i n i Γ 2 + B i n Γ 1 .
Because i α ρ ( A ) , multiplying from the left by ( i α I A ) 1 both sides of (49) results
Π 1 + i Π 2 = ( i α I A ) 1 B i n ( Γ 1 + i Γ 2 ) .
Applying C to both sides of (50) and recalling that C Π w = Q w , with C Π 1 = 1 and C Π 2 = 0 , then Q = [ 1 0 ] . Consequently,
1 = C ( i α I A ) 1 B i n ( i Γ 2 + Γ 1 ) = G ( i α ) ( Γ 1 + i Γ 2 ) .
Rewriting (51) in terms of G ( i α ) = Re ( G ( i α ) ) + i Im ( G ( i α ) ) , it yields
1 = ( Re ( G ( i α ) ) + i Im ( G ( i α ) ) ) ( Γ 1 + i Γ 2 ) .
From the above, matching real and imaginary parts,
1 = ( Re ( G ( i α ) ) Γ 1 Im ( G ( i α ) ) Γ 2 , 0 = ( Im ( G ( i α ) ) Γ 1 + Re ( G ( i α ) ) Γ 2 .
Hence,
Γ 1 = Re ( G ( i α ) ) | G ( i α ) | 2 , Γ 2 = Im ( G ( i α ) ) | G ( i α ) | 2 .
Accordingly, from the fact that
G ( i α ) 1 = 1 G ( i α ) = G ( i α ) ¯ | G ( i α ) | 2 = Re ( G ( i α ) ) i Im ( G ( i α ) ) | G ( i α ) | 2 ,
the desired control is given by
Γ = [ Γ 1 Γ 2 ] = [ Re ( G ( i α ) 1 ) Im ( G ( i α ) 1 ) ] .
Here, it has been assumed that the system is real, i.e., G ( s ¯ ) = G ( s ) ¯ for all s ρ ( A ) . It should be noticed that G ( i α ) 0 for solvability.

3.2. Harmonic Disturbance Rejection

Now, for the last part where is desirable to reject a sinusoidal disturbance such as d ( t ) = A d sin ( β t ) , let us chose an exosystem that is governed by the harmonic oscillator
w t = S w , w ( 0 ) = 0 A d , S = 0 β β 0 ,
with solution
w ( t ) = A d sin ( β t ) A d cos ( β t ) .
Taking Q = [ 0 0 ] and P = [ 1 0 ] , then
y r ( t ) = Q w = 0 .
Again, to this case W = R 2 so, we look for Π = [ Π 3 Π 4 ] , where Π j Z , and Γ = [ Γ 3 Γ 4 ] R 2 .
The Francis equations that are applied to the vector w = [ w 3 w 4 ] T W result in the system
Π S w = A Π w + B d P w + B i n Γ w ,
C Π w = 0 .
Expanding the above matrix multiplications, the regulator equation from the left becomes
β Π 3 w 4 β Π 4 w 3 = A Π 3 w 3 + A Π 4 w 4 + B d w 3 + B i n Γ 3 w 3 + B i n Γ 4 w 4 .
Because (55) has to hold for w, now consider the case for w 3 = 1 and w 4 = 0 to then consider that for w 3 = 0 and w 4 = 1 resulting
β Π 4 A Π 3 = B d + B i n Γ 3 ,
β Π 3 A Π 4 = B i n Γ 4 .
By multiplying (57) by i = 1 and adding the result to (56), we get
( i β I A ) Π 3 + i ( i β I A ) Π 4 = B d + B i n i Γ 4 + B i n Γ 3 .
Noting that i β ρ ( A ) , applying ( i β I A ) 1 to both sides of (58), we get
Π 3 + i Π 4 = ( i β I A ) 1 B d + ( i β I A ) 1 B i n ( i Γ 4 + Γ 3 ) .
By multiplying both sides of (59) by C, from (54) implies
C Π 3 = 0 , C Π 4 = 0 ,
thus, we get
0 = C ( i β I A ) 1 B d + C ( i β I A ) 1 B i n ( i Γ 4 + Γ 3 ) = G B d ( i β ) + G ( i β ) ( i Γ 4 + Γ 3 )
where we have used the notation
G B d ( s ) = C ( s I A ) 1 B d .
Finally, solving for Γ we have
Γ 3 = Re ( G ( i β ) 1 ) Re ( G B d ( i β ) ) + Im ( G ( i β ) 1 ) Im ( G B d ( i β ) ) ,
Γ 4 = Re ( G ( i β ) 1 ) Im ( G B d ( i β ) ) Im ( G ( i β ) 1 ) Re ( G B d ( i β ) ) ,
where
Γ = [ Γ 3 Γ 4 ] .
Consequently, by combining (52) and (60) we get
Γ = [ Γ 1 Γ 2 Γ 3 Γ 4 ] .
Although it could be argued that the 1-D R-D equation was so easy to find an explicit formula for the transfer function G ( s ) = C ( s I A ) 1 B i n , then there is no reason for which the state feedback regulator problem to this type of PDE had not been reported in the literature in spite of its importance to a myriad of applications.

4. Simulation Results

In our numerical simulation, we have set A d = 2 , A r = 2 , α = 4 , β = 4 , x 0 = 0.75 , x 1 = 0.25 and ν 0 = ν 1 = 0.25 . From Figure 1, Figure 2 and Figure 3, it can be seen that the controlled output y ( t ) tracks the reference signal y r ( t ) from an initial condition φ ( x ) = 4 cos ( π x ) . Additionally, it can be seen that the error signal e ( t ) tends to zero as time tends to infinity. The performance of the regulator is shown for different λ values. The convergence of y ( t ) to y r ( t ) is faster for small λ values. From Figure 4, Figure 5 and Figure 6, the solution surface is shown for every λ value. In mass balance the source term λ z can be interpreted as the rate of production (destruction) of species per unit volume. From (34), the sign from the diffusion term agrees with the observation that mass flows from high to low mass concentration. In general, pressure gradients, temperature, and external forces affect the mass flux but their effects can be neglected and take the diffusion coefficient D > 0 to be constant. In fact, the units may be chosen in a convenient way in order to make D = 1 . In our work, the diffusion coefficient is set to the unit value.

5. Conclusions

In [16], from an abstract control system model for the 1-D heat equation, the SFRP was solved in terms of the transfer function of the system, but with feedback control law in function of the state of the exosystem. The SFRP to the 1-D heat equation was solved for harmonic tracking and rejection of a constant disturbance. The SLDO (5) is present for most of the R-D equation types [34]. In our work, we propose the state operator for an abstract control system model to the 1-D R-D equation, in contrast with that in [16], namely, (13), as given in the form of the SLDO (12). Accordingly, the 1-D R-D equation is characterized, along the same lines as in [16], in terms of the Francis (regulator) equations, but with state operator (12) and then these are solved. The SFRP to the 1-D R-D equation is solved for harmonic tracking with harmonic disturbance rejection. The simulation results validate the performance of the regulator, i.e., the error tends to zero as time tends to infinity. From all of the above, we conclude that our proposal performs well, i.e., the SFRP for the 1-D R-D equation has solution. As future work, our proposal may be extended to those approaches in [16], and so on, on which the abstract control system model to the R-D equation could falls.

Author Contributions

Conceptualization, F.J.; methodology, F.J. and A.A.R.; software, A.A.R.; validation, F.J. and A.A.R.; formal analysis, F.J. and A.A.R.; investigation, F.J.; data curation, A.A.R.; writing—original draft preparation, F.J.; writing—review and editing, F.J.; visualization, A.A.R.; resources, F.J.; supervision, F.J.; funding acquisition, F.J.; project administration, F.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Tecnológico Nacional de México (TecNM) under project 8007.20-P and, partially, under grant number 30844 from EDD 2019 program.

Acknowledgments

This work was supported by CONACYT-SENER, México, through grant number 862135.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DPSsDistributed Parameters Systems
EFRPError Feedback Regulation Problem
ODEsOrdinary Differential Equations
PDEsPartial Differential Equations
R-DReaction-Diffusion
SFRPState Feedback Regulation Problem
SLBVPSturm-Liouville Boundary Value Problem
SLDOSturm-Liouville Differential Operator

References

  1. Morris, K. Control of Systems Governed by Partial Differential Equations in The Control Handbook, 2nd ed.; Levine, W.S., Ed.; CRC Press: Boca Raton, FL, USA, 2011. [Google Scholar]
  2. Isidori, A. Nonlinear Control Systems; Springer: London, UK, 1995. [Google Scholar]
  3. Isidori, A.; Byrnes, C.I. Output Regulation of Nonlinear Systems. IEEE Trans. Autom. Control 1990, 35, 131–140. [Google Scholar] [CrossRef]
  4. Francis, B.A.; Wonham, W.M. The Internal Model Principle of Control Theory. Automatica 1976, 12, 457–465. [Google Scholar] [CrossRef]
  5. Francis, B.A. The Linear Multivariable Regulator Problem. SIAM J. Control Optim. 1977, 15, 486–505. [Google Scholar] [CrossRef]
  6. Wonham, W.M. Linear Multivariable Control: A Geometric Approach; Springer: New York, NY, USA, 1985. [Google Scholar]
  7. Hautus, M.L.J. Linear Matrix Equations with Applications to the Regulator Problem. In Outils et modèles Mathématiques Pour L’automatique: L’analyse de Systèmes et le Traitement du Signal; Editions du Centre National de la Recherche Scientifique: Paris, France, 1982. [Google Scholar]
  8. Knobloch, H.W.; Isidori, A.; Flockerzi, D. Topics in Control Theory; Birkhäuser Verlag: Basel, Switzerland, 1993. [Google Scholar]
  9. Pohjolainen, S.A. On The Asymptotic Regulation Problem for Distributed Parameter Systems. Control Distrib. Parameter Syst. 1983, 197–201. [Google Scholar] [CrossRef]
  10. Pohjolainen, S.A. Robust Multivariable PI–Controller for Infinite Dimensional Systems. IEEE Trans. Autom. Control 1982, 27, 17–30. [Google Scholar] [CrossRef]
  11. Schumacher, J.M. Dynamic Feedback in Finite- and Infinite-Dimensional Linear Systems; Mathematisch Centrum Tracts Centrum Voor Wiskunde en Informatica: Amsterdam, The Netherlands, 1981; Volume 143. [Google Scholar]
  12. Schumacher, J.M. Finite-Dimensional Regulators for a Class of Infinite–Dimensional Systems. Syst. Control Lett. 1983, 3, 7–12. [Google Scholar] [CrossRef] [Green Version]
  13. Byrnes, C.I.; Laukó, I.G.; Gilliam, D.S.; Shubov, V.I. Output Regulation for Linear Distributed Parameter Systems. IEEE Trans. Autom. Control 2000, 45, 2236–2252. [Google Scholar]
  14. Müller, P.H.; Kato, T. Perturbation theory for linear operators (Die Grundlehren der mathematischen Wissenschaften in Einzeldarstellungen mit besonderer Berücksichtigung der Anwendungsgebiete, Band 132) XX+ 592 S. m. 3 Fig. Springer–Verlag. Preis geb. DM 79, 20. Z. Angew. Math. Und Mech. 1967, 47. [Google Scholar] [CrossRef]
  15. Curtain, R.F.; Zwart, H. An Introduction to Infinite-Dimensional Linear Systems Theory; Springer: New York, NY, USA, 1995. [Google Scholar]
  16. Aulisa, E.; Gilliam, D. A Practical Guide to Geometric Regulation for Distributed Parameter Systems; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  17. Natarajan, V.; Gilliam, D.S.; Weiss, G. The State Feedback Regulator Problem for Regular Linear Systems. IEEE Trans. Autom. Control 2014, 59, 2708–2723. [Google Scholar] [CrossRef]
  18. Xu, X.; Dubljevic, S. Output and Error Feedback Regulator Designs for Linear Infinite-Dimensional Systems. Automatica 2017, 83, 170–178. [Google Scholar] [CrossRef]
  19. Deutscher, J. A Backstepping Approach to the Output Regulation of Boundary Controlled Parabolic PDEs. Automatica 2015, 57, 56–64. [Google Scholar] [CrossRef]
  20. Deutscher, J. Finite-Time Output Regulation for Linear 2 × 2 Hyperbolic Systems using Backstepping. Automatica 2017, 75, 54–62. [Google Scholar] [CrossRef]
  21. Zhou, H.C.; Guo, B.Z. Performance Output Tracking for One-Dimensional Wave Equation Subject to Unmatched General Disturbance and Non-Collocated Control. Eur. J. Control 2018, 39, 39–52. [Google Scholar] [CrossRef]
  22. Humaloja, J.; Kurula, M.; Paunonen, L. Approximate Robust Output Regulation of Boundary Control Systems. IEEE Trans. Autom. Control 2019, 64, 2210–2223. [Google Scholar] [CrossRef] [Green Version]
  23. Rebarber, R.; Weiss, G. Internal Model based Tracking and Disturbance Rejection for Stable Well-Posed Systems. Automatica 2003, 39, 1555–1569. [Google Scholar] [CrossRef]
  24. Jin, F.F.; Guo, B.Z. Boundary Output Tracking for an Euler–Bernoulli Beam Equation with Unmatched Perturbations from a Known Exosystem. Automatica 2019, 109, 108507. [Google Scholar] [CrossRef]
  25. Xu, X.; Dubljevic, S. Output Regulation for a Class of Linear Boundary Controlled First–Order Hyperbolic PIDE Systems. Automatica 2017, 85, 43–52. [Google Scholar] [CrossRef]
  26. Logan, J.D. An Introduction to Nonlinear Partial Differential Equations; John Wiley & Sons: Hoboken, NJ, USA, 2008. [Google Scholar]
  27. Smoller, J. Shock Waves and Reaction-Diffusion Equations; Springer: New York, NY, USA, 1983. [Google Scholar]
  28. Murray, J.D. Mathematical Biology I: An Introduction; Springer: Berlin/Heidelberg, Germany, 2002. [Google Scholar]
  29. Murray, J.D. Mathematical Biology II: Spatial Models and Biomedical Applications; Springer: Berlin/Heidelberg, Germany, 2003. [Google Scholar]
  30. Tveito, A.; Winther, R. Introduction to Partial Differential Equations A Computational Approach; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  31. Banks, H.T.; Tran, H.T. Mathematical and Experimental Modeling of Physical and Biological Processes; CRC Press: Boca Raton, FL, USA, 2009. [Google Scholar]
  32. Gahn, M.; Neuss-Radu, M.; Knabner, P. Homogenization of Reaction-Diffusion Processes in a Two-Component Porous Medium with Nonlinear Flux Conditions at the Interface. SIAM J. Appl. Math. 2016, 76, 1819–1843. [Google Scholar] [CrossRef]
  33. Cardone, G.; Perugia, C.; Timofte, C. Homogenization Results for a Coupled System of Reaction-Diffusion Equations. Nonlinear Anal. 2019, 188, 236–264. [Google Scholar] [CrossRef] [Green Version]
  34. Karafyllis, I.; Krstic, M. Input-to-State Stability for PDEs; Springer: Cham, Switzerland, 2019. [Google Scholar]
  35. Krstic, M.; Smyshlyaev, A. Boundary Control of PDEs A Course on Backstepping Designs; SIAM: Philadelphia, PA, USA, 2008. [Google Scholar]
  36. Ramírez, A.A.; Jurado, F. Tracking Regulator Design with Disturbance Rejection to the Reaction–Diffusion Equation. In Proceedings of the 2019 16th International Conference on Electrical Engineering, Computing Science and Automatic Control (CCE), Mexico City, Mexico, 11–13 September 2019. [Google Scholar]
  37. Seidman, T.I. Control of the Heat Equation in The Control Handbook, 2nd ed.; Levine, W.S., Ed.; CRC Press: Boca Raton, FL, USA, 2011. [Google Scholar]
  38. Strauss, W.A. Partial Differential Equations an Introduction; John Wiley & Sons: Hoboken, NJ, USA, 2008. [Google Scholar]
  39. Mathews, J.; Walker, R.L. Mathematical Methods of Physics; The Benjamin/Cummings Publishing Co.: Menlo Park, CA, USA, 1970. [Google Scholar]
  40. Boyce, W.E.; DiPrima, R.C.; Meade, D.B. Elementary Differential Equations and Boundary–Value Problems; John Wiley & Sons: Hoboken, NJ, USA, 2017. [Google Scholar]
  41. Kreyszig, E. Advanced Engineering Mathematics; John Wiley & Sons: Hoboken, NJ, USA, 2006. [Google Scholar]
  42. Farlow, S.J. Partial Differential Equations for Scientists and Engineers; Dover Publications: New York, NY, USA, 1993. [Google Scholar]
  43. Holland, S.S., Jr. Applied Analysis by the Hilbert Space Method, an Introduction with Applications to the Wave, Heat, and Schrödinger Equations; Dover Publications: Mineola, NY, USA, 1990. [Google Scholar]
  44. Kaplan, K. Advanced Calculus; Pearson Education: Boston, MA, USA, 2003. [Google Scholar]
  45. Pazy, A. Semigroups of Linear Operators and Applications to Partial Differential Equations; Springer: New York, NY, USA, 1983. [Google Scholar]
Figure 1. Performance of the regulator for λ = 2.5 with Γ 1 = 1.8424 , Γ 2 = 1 , Γ 3 = 0.5 , and Γ 4 = 0 .
Figure 1. Performance of the regulator for λ = 2.5 with Γ 1 = 1.8424 , Γ 2 = 1 , Γ 3 = 0.5 , and Γ 4 = 0 .
Mathematics 08 01983 g001
Figure 2. Performance of the regulator for λ = 0 , i.e., to the case for which the R-D equation is reduced to the Fick’s second law of diffusion (also called simply as diffusion equation), for whose case Γ 1 = 2.4674 , Γ 2 = 1 , Γ 3 = 0.5 , and Γ 4 = 0 .
Figure 2. Performance of the regulator for λ = 0 , i.e., to the case for which the R-D equation is reduced to the Fick’s second law of diffusion (also called simply as diffusion equation), for whose case Γ 1 = 2.4674 , Γ 2 = 1 , Γ 3 = 0.5 , and Γ 4 = 0 .
Mathematics 08 01983 g002
Figure 3. Performance of the regulator for λ = 2.5 , for whose case Γ 1 = 3.0924 , Γ 2 = 1 , Γ 3 = 0.5 , and Γ 4 = 0 .
Figure 3. Performance of the regulator for λ = 2.5 , for whose case Γ 1 = 3.0924 , Γ 2 = 1 , Γ 3 = 0.5 , and Γ 4 = 0 .
Mathematics 08 01983 g003
Figure 4. Plot of the solution surface for λ = 2.5 .
Figure 4. Plot of the solution surface for λ = 2.5 .
Mathematics 08 01983 g004
Figure 5. Solution surface for the case in which the reaction term is omitted, i.e., to the case for which λ = 0 .
Figure 5. Solution surface for the case in which the reaction term is omitted, i.e., to the case for which λ = 0 .
Mathematics 08 01983 g005
Figure 6. Plot of the solution surface for λ = 2.5 .
Figure 6. Plot of the solution surface for λ = 2.5 .
Mathematics 08 01983 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jurado, F.; Ramírez, A.A. State Feedback Regulation Problem to the Reaction-Diffusion Equation. Mathematics 2020, 8, 1983. https://doi.org/10.3390/math8111983

AMA Style

Jurado F, Ramírez AA. State Feedback Regulation Problem to the Reaction-Diffusion Equation. Mathematics. 2020; 8(11):1983. https://doi.org/10.3390/math8111983

Chicago/Turabian Style

Jurado, Francisco, and Andrés A. Ramírez. 2020. "State Feedback Regulation Problem to the Reaction-Diffusion Equation" Mathematics 8, no. 11: 1983. https://doi.org/10.3390/math8111983

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop