The problem of signal tracking in discrete linear time invariant systems, in the presence of a disturbance signal in the plant, is solved using a new zero-variation methodology. A discrete-time dynamic output feedback controller is designed in order to minimize the ℋ∞ norm between the exogen input and the output signal of the system, such that the effect of the disturbance is attenuated. Then, the zeros modification is used to minimize the ℋ∞ norm from the reference input signal to the error signal. The error is taken as the difference between the reference and the output signal. The proposed design is formulated in linear matrix inequalities (LMIs) framework, such that the optimal solution of the stated problem is obtained. The method can be applied to plants with delay. The control of a delayed system illustrates the effectiveness of the proposed method.
1. Introduction
In a control systems theory, the design of controller using pole placement of closed loop discrete-time systems can be easily done. In [1] a controller using pole placement is used to obtain an exact plot of complementary root locus, of biproper open-loop transfer functions, using only well-known root locus rules. However, the problem of zero placement is not very much studied by the control researchers. In [2] a discrete-time pole placement is obtained by a control design technique that uses simple and multirate sample. The methodology proposed in [3] preserves the ℋ2 state feedback controller optimality by pole placement in a Z plain region specified in design. In the field of discrete-time systems pole placement we find [4], where discrete adaptative controllers are designed considering arbitrary zero location. Also, in [5] a class of nonmodeled dynamics is controlled using a zero placement.
In [6] a methodology is proposed using zero and pole placement for discrete-time systems, to obtain the signal tracking and disturbance rejection, respectively. However, for the signal tracking problem, when a state feedback estimator is proposed, a modification occurs in ℋ∞-norm value obtained with the initial controller that provides the disturbance rejection. The methodology proposed in this paper has the advantage of maintaining the ℋ∞-norm value obtained with the initial controller for the signal tracking problem.
The problem of signal tracking, in the presence of disturbance signal for continuous-time plant, was solved in [7], using a zero variation methodology. A methodology with a simpler mathematic formulation is proposed in [8]. The signal tracking problem with disturbance rejection in discrete-time systems is solved by an analytic method in [9], however, the mathematic formulation is complex and a frequency selective tracking is not presented as proposed in this manuscript. In [10] the use of linear matrix inequalities (LMIs) is considered in design of controllers, filters and stability study. Also, in [11], LMIs are used for the design of a dynamic output feedback controller in order to guarantee the asymptotic stability of a continuous-time system and minimize the upper bound of a given quadratic cost function. Furthermore, the LMI formulation has been used in several engineering problems (see, e.g., [12–20]).
This manuscript proposes a formulation of a signal tracking with disturbance rejection optimization problem for discrete-time systems in the linear matrix inequalities framework, such that the optimal solution of the stated control problem is obtained. The proposed method is simpler than the other tracking techniques, and the main result is that when the problem is feasible the optimal solution is obtained with small computation effort, as the LMIs can be solved using linear programming algorithms, with polynomial convergence. The software MATLAB [21] is used to find the LMI solutions, when the problem is feasible. The control of a delayed system illustrates the effectiveness of the proposed method.
2. Statement of the Problems
Consider a controllable and observable linear time-invariant multi-input multi-output (MIMO) discrete-time system,x(k+1)=Ax(k)+Buu(k)+Bww(k),y(k)=C1x(k),z(k)=C2x(k),x(0)=0,k∈[0;∞),
where A∈ℜn×n, Bu∈ℜn×p, Bw∈ℜn×q, C1∈ℜm×n, C2∈ℜm×n, x(k) is the state vector, y(k) is the output vector, u(k) is the control input and w(k) is the disturbance input (exogenous input).
Problem 1.
The disturbance rejection problem for discrete-time systems, using the dynamic output feedback of the system described in (2.1), is the following: minimize the upper bound of ℋ∞ norm from the exogenous input w(k) to the output z(k). In this context, the objective is to design a controller ℋ∞, Kc(z), that attenuates the effect of disturbance signal in the output of the system. And in the tracking process, it is needful to design the controllers M and N that minimize the ℋ∞ norm between the reference input r(k) and the tracking error r(k)-z(k).
Remark 2.1.
The problems of weighted disturbance reduction and weighted reference tracking are closely related to the above one and will be addressed in Section 5.
Remark 2.2.
The block diagram of the control process used in this manuscript to solve Problem 1 is given in Figure 1, where Kc(z)=Cc(zI-Ac)-1Bc is a ℋ∞ dynamic output feedback controller, the controllers M and N were used to solve the zero variation problem in order to obtain a tracking system, r(k) is the reference input signal.
Discrete-time optimal control system with pole placement and zero modification.
The state space equation of the control system shown in Figure 1 can be written as:[x(k+1)xc(k+1)]=[ABuCcBcC1Ac][x(k)xc(k)]+[BuNM]r(k)+[Bw0]w(k),y(k)=[C10][x(k)xc(k)],z(k)=[C20][x(k)xc(k)],e(k)=r(k)-z(k)=r(k)-[C20][x(k)xc(k)].
Rewriting the system (2.2) in a compact form, it follows that:x¯(k+1)=Amx¯m(k)+Bmr(k)+Bnw(k),e(k)=-Cmx¯(k)+Dmr(k),z(k)=Cmx¯(k),
wherex¯m(k)=[x(k)xc(k)],Am=[ABuCcBcC1Ac],Dm=1,Bm=[BuNM],Bn=[Bw0],Cm=[C20].
Using the 𝒵-transform in order to solve the system (2.3), consider the initial conditions equal to zero. One obtains the transfer function between input signals (reference input and exogenous input) and the measured output of the system as showed in the following equation:Z(z)=Cm(zI-Am)-1BmR(z)+Cm(zI-Am)-1BnW(z).
For the transfer function from w(Z) to z(Z) in (2.6), the minimization of ℋ∞ norm is obtained with the initial design of ℋ∞ controller, that implies in the minimization of the perturbation effect in to the system output.
Figure 1 shows the addition of the term Mr(k) in the structure of the Kc(z) controller. The purpose of the controller M is only to change the zeros of the transfer function from r(k) to u(k) and it does not change the poles obtained in the initial design of Kc(z). The transfer function from W(z) to Z(z) is not changed by N or M, according to (2.4) and (2.5). In this way the performance of the ℋ∞ norm controller is not affected.
For the optimal tracking design, the relation between error signal and reference signal described in (2.7) is considered, making the perturbation signal W(z) equal to zero in (2.6),Hm(z)=E(z)R(z)=-Cm(zI-Am)-1Bm+Dm.
In this case, using the zero modification one can design a tracking system that minimizes the ℋ∞ norm between the reference input r(k) and the tracking error r(k)-z(k). In Section 4, motivated by the work in [22], we show that M and N modify the zeros from r(k) to u(k). The process of the zeros modification does not interfere in the disturbance rejection. Therefore in agreement with (2.6), Bm has no influence on the transfer function from W(z) to Z(z). In (2.6) one uses the zeros location, by the specifications of the N and M in Bm, in the process of minimization of the ℋ∞ norm of the transfer function between the reference signal and the tracking error.
3. ℋ∞ Dynamic Output Feedback Controller Design
The following theorem leads to a new method to design the Kc(z) in a LMI framework, and the goal is to attenuate the effects of exogenous signal in the output of discrete-time systems. By using [23], a pole placement constraint region with radius r and center in (-q,0) is required and used in this work to provide the designer with an expedite way to keep the controller gains within appropriate bounds.
Theorem 3.1.
Consider the system (2.1) with dynamic output feedback by the ℋ∞ upper bound controller, Kc(z). Then the optimal solution of the ℋ∞ norm between the input w(k) and the output z(k), with pole placement in a region of radius r and center in (-q,0) can be obtained from the solution of the following LMI optimization problem,
‖H‖∞2=minμs.t.[RI0BwAR+BuCjAIS0SBwAjSA+BjC200IDnC2RC2Bw′Bw′SDn′μI00RA′+Cj′Bu′Aj′RC2′0RIA′A′S+C2′Bj′C2′0IS]>0,[-Rr-IrAR+BuCj+RqA+Iq-Ir-SrAj+IqSA+BjC1+SqRA′+Cj′Bu′+RqAj′+Iq-Rr-IrIq+A′A′S+C1′Bj′+Sq-Ir-Sr]<0[RIIS]>0,
where, R=R′∈ℜn×n, S=S′∈ℜn×n, Aj, Bj, Cj and Dj in (3.1), (3.2), and (3.3) are the set of LMIs optimization variables. The radius r and center in (-q,0) are the pole placement constraints illustred in Figure 2
For solution of (3.4): ψE′=I-RS, where ψ and E′ can be obtained by L-U decomposition of I-RS [24].
The ℋ∞ compensator dinamic matrices Kc(z)=Cc(zI-Ac)-1Bc+Dc is obtained with the solution of equation (3.4).
In (3.1), (3.2) and (3.3), one has the following change of variables:Cj=Ccψ′,Bj=EBc,Aj=(SA+EBcC1)R+(SBuCc+EAc)ψ′.
The pole placement region in the 𝒵 plane.
Proof.
A realization of dynamic output feedback system Rwz is as followings:
Rwz=Ccl+(zI-Acl)-1Bcl,
where
Acl=[ABuCcBcC1Ac],Bcl=[Bw0],Ccl=[C20].
The optimization problem below described in LMI framework [25], is used to design the ℋ∞ compensator with pole placement constraints
‖H‖∞2=minμs.t.[Q̌0BclAclQ̌0IDCclQ̌Bcl′D′μI0Q̌Acl′Q̌Ccl′0Q̌]>0,[-rpQ̌AclQ̌+qQ̌Q̌q+Q̌Acl′-rpQ̌]<0,Q̌=Q̌′>0,μ>0.
However, a high computational effort is needed to solve this problem, because the optimization problems (3.7) and (3.8) are described as a solution of BMIs. Then, using a linear transformation, the problem can be easily solved, based on LMI framework. First, the matrix Q̌ and its inverse are considered as follows
Q̌=[Rψψ′J],Q̌-1=[SEE′S],
where R=R′∈ℜn×n, S=S′∈ℜn×n, and
Q̌Γ2=Γ1,withΓ1=[RIψ′0],Γ2=[IS0E′].
The condition (3.3) is obtained by considering the Lyapunov matrix Q̌>0 and premultiplying and postmultiplying Q̌ by Γ2′ and Γ2, respectively. After, the condition (3.1) is obtained premultiplying and postmultiplying the inequation (3.7) by (3.11) and (3.12), respectively. Then, pre and postmultiplying (3.8) by (3.13) and (3.14), respectively, results in condition (3.2)
[Γ2′0000I0000I0000Γ2′],[Γ20000I0000I0000Γ2],[Γ2′00Γ2′],[Γ200Γ2].
The zero modification is showed in the next section.
Remark 3.2.
In this paper, the methodology adopted to solve the pole location problem affords the designer an expedite way to keep the controller gains within appropriate bounds, a key requisite for implementation purposes.
4. Zeros Variation
Motivated by the work in [22] the present paper uses the zero variation in order to obtain the global optimum of ℋ∞ norm to solve the tracking problem.
A controller that modifies the zeros from r(k) to u(k) is designed considering the system (A,Bu,C1) as shown in Figure 3. The zeros of closed-loop system are allocated at arbitrary places according to the M and N values, where M∈ℜn×1 and N∈ℜ1×1.
Zeros variation of the controlled system.
The plant is described byx(k+1)=Ax(k)+Bu(k),y(k)=C1x(k).
The 𝒵 transform of (4.1) with zero initial condition, is[zI-A]X(z)=BuU(z),Y(z)=C1X(z).
A zero of the system is a value of z such that the system output is zero even with a nonzero state-and-input combination. Thus if we are able to find a nontrivial solution for X(z0) and U(z0) such that Y(z0) is null, then z0 is a zero of the system [22]. Combining the two parts of (4.2) we must satisfy the following requirement:[ziI-A-BuC10][X(z)U(z)]=[00].
Also, the compensator can be described as follows xc(k+1)=Acxc(k)+Bcuc(k),
where uc(k)=C1x(k).
A more general method to introduce r(k) is to add a term Mr(k) to xc(k+1) and also a term Nr(k) to the control equation u(k)=Ccxc(k), as shown in Figure 3. The controller, with these additions, becomes equal to:xc(k+1)=Acxc(k)+Bcuc(k)+Mr(k),u(k)=Ccxc(k)+Nr(k).
Now, considering the controller (4.5), if there exists a transmission zero from r(k) to u(k), then necessarily there exists a transmission zero from r(k) to y(k), unless a pole of the plant cancels the zero. The equation to obtain zi from r(k) to u(k) (we let y(k)=0 because we are considering only the effects of r(k), then uc(k)=0) in (4.5), is the following:[ziI-Ac-MCcN][xcoro]=[00].
Because the coefficient matrix in (4.6) is square, the condition for a nontrivial solution is that the determinant of this matrix must be zero. Thus we havedet[ziI-Ac-MCcN]=0.
Multiplying the second column of the matrix described in (4.7) at right by a nonzero matrix N-1 and then adding to the first column of (4.7) the product of -Cc by the last column, we have:det[ziI-Ac+MN-1Cc-MN-101]=0.
And so, considering zi=z,
det(zI-Ac+MN-1Cc)=0,
where the modified zeros from r(k) to u(k) are the solutions z=zi. It is important to notice that the gain N and the vector M do not only modify the system zeros but also are used to obtain the optimal solution of the tracking problem.
5. Tracking Design
The solution of the tracking problem is based on the design of the matrices of the controller M and N, that minimize the ℋ∞ norm of (Am,Bm,-Cm,1). Weighted frequency is added in the tracking system in order to track signals in a frequency band specified in the design. In tracking design with wheigthed frequency, the goal is to find a global solution that optimize the problem described as follows:min‖Hm(z)V(z)‖∞,
where V(z)=(Av,Bv,Cv,Dv) is a dynamic system designed to specify wheighted frequency in the output. A stable, linear and time invariant system realization Hm=(Am,Bm,-Cm,Dm) is considered as indicated in (2.7). Figure 4 illustrates the structure of inclusion of frequency wheighted in the design of tracking system.
System structure to the tracking with frequency wheighted.
The system (2.4) can be represented by state variables in function of xm(k) and xv(k), as follows:[xm(k+1)xv(k+1)]=[Am0-BvCmAv][xm(k)xv(k)]+[BmBvDm]r(k),yv(k)=[0Cv][xm(k)xv(k)].
In addition, a possible state space realization of Ȟf=Hm(z)V(z) is:[ǍfB̌fČfĎf]=[Am0Bm-BvCmAvBvDm0Cv0].
A metodology for the tracking design problem solution with wheigthed band is proposed in Theorem 5.1 considering the elements of the compensator matrix already fixed.
Theorem 5.1.
Considering the system with output filter (5.3), if there exist a solution for the LMI described in (5.4) and (5.5), then the gain N and vector M that minimize the ℋ∞ norm from r(k) to e(k) can be obtained
[Q11Q12Q13Q12′Q22Q23Q13′Q23′Q33]>0,μ>0.
Proof.
For the tracking design with band weighted, substitute Ǎf, B̌f, Čf, and Ďf in (5.3) by (3.7). It results in the optimization problem described in (5.4) and (5.5). The gain N and vector M are obtained from this process. The matrix Q is partitioned in the form Qij=Qij′, i,j=1,2,3.
The theorem proposes that one can obtain an optimal zero location that guarantee the global optimization of the tracking error ℋ∞ norm. The proposed design is formulated in linear matrix inequalities (LMIs) framework, such that the optimal solution is obtained.
The controller M and N are the optimal solution of (5.4) and (5.5) and minimize the ℋ∞ norm between the reference input signal r(k) and the tracking error signal r(k)-e(k).‖H‖∞2=minμs.t.[Q11Q12Q130Q12′Q22Q230Q13′Q23′Q330000IN′Bu′M′Bv′IQ12BuCc′+Q11A′Q12Ac′+Q11BcC2′Q13Av′-Q11C2′Bv′-Q13Cv′Q22BuCc′+Q12′A′Q22Ac′+Q12′BcC2′Q23Av′-Q12′C2′Bv′-Q23Cv′Q23′BuCc′+Q13′A′Q23′Ac′+Q13′BcC2′Q33Av′-Q13′C2′Bv′-Q33Cv′BuNAQ11+BuCcQ12′AQ12+BuCcQ22AQ13+BuCcQ23MBcC2Q11+AcQ12′BcC2Q12+AcQ22BcC2Q13+AcQ23Bv-BvC2Q11+AvQ13′-BvC2Q12+AvQ23′-BvC2Q13+AvQ33I-CvQ13′-CvQ23′-CvQ33μI0000Q11Q12Q130Q12′Q22Q230Q13′Q23′Q33]>0.
The filter has the goal of adjusting the controllers M and N for the selected frequency band. Then the zero variation described in LMI framework considers the filter dynamics to adjust the operation of the system tracking in this frequency band. However, in the implementation or simulation of the tracking process control, the filter is discarded.
The paper also proposes the statement of weighted disturbance minimization problem that is achieved based on Theorems 3.1 and 5.1 demonstrations and uses the following formulation: ‖H‖∞2=minμs.t.[RI0IS000IBw′Bw′S0RA′+Cj′Bu′+Afr′+R′C1′Bfr′Aj′Cfr′A′+C1′Bfr′A′S+C2′Bj′C2′BwAR+BuCj+Afr+BfrC1RA+BfrC1SBwAjSA+BjC20Cfr0μI000RI0IS]>0,[-Rr-Ir-Ir-SrRA′+Cj′Bu′+Rq+Afr′+RC1′Bfr′Aj′+IqIq+A′+C1′Bfr′A′S+C1′Bj′+SqAR+BuCj+Rq+Afr+BfrC1RA+Iq+BfrC1Aj+IqSA+BjC1+Sq-Rr-Ir-Ir-Sr]<0,[RIIS]>0.
The following examples illustrate the effectiveness of the proposed method.
6. Example 1
Consider a tank temperature control system described in [22], with a zero-order hold. The goal is to design a tracking system for flow control with disturbance attenuation. The model of the delayed system [22] is,
H(s)=e(-λ*s)sa+1.
A sampling period of 0.01 seconds is used in the design. The parameters a=1 and λ=0.005s are adopted in the design.
We found λ=1*0.01-0.5*0.01, and therefore, l=1 and m=0.5.
Using the process to find the 𝒵 transform of a delayed continuous-time function (6.1), the parameters are substituted and we obtainH(z)=(1-e-0.005)z+(e-0.005-e-0.01)/(1-e-0.005)z(z-e-0.01),
or,
H(z)=0.0050z+0.0050z2-0.9900z.
The state-space description of the system is[x1(k+1)x2(k+1)]=[0.99010][x1(k)x2(k)]+[10]u(k)+[10]w(k),y(k)=[0.0050.005][x1(k)x2(k)],
where x(k) is the state vector, u(k) is the control signal and w(k) is a disturbance signal in the system.
The design of the tracking system must include operation for reference signals of low frequencies (smaller than 0.1 rad/s). In such a case, the following filter J(z) was consideredJ(z)=(0.4500z+0.4500)×10-7z2-1.9999z-0.9999.
Using Theorem 3.1 the controller Kc(z) is designed for the system described in (6.4) and showed in (6.6). This controller minimizes the ℋ∞-norm of w(k) to z(k). In this design we obtain a disk of radius r=0.5 and center in q=0.5 is used as a pole placement constraint:Kc(z)=-24.1395z-4.6956z2-0.2981z+0.1392.
The ℋ∞-norm of w(k) to z(k) for the closed-loop system is 0.0342, implying the attenuation of the effect of the disturbance signal in the system. Figure 5 illustrates the frequency response of Z(z)/W(z).
Frequency response of Z(z)/W(z).
To design a tracking system, the proposed zero-variation methodology given by (5.4) and (5.5) were used and the ℋ∞ norm of r(k) to e(k) was minimized considering the signal r(k) with low frequencies (smaller than 0.1 rad/s). The obtained ℋ∞-norm for all frequency spectra was equal to 1.32, while for frequency band specified in the problem, the largest magnitude of frequency response was 3.13 × 10−3). This implies that the tracking system operated adequately in the frequency band specified in the problem.
Figure 6 illustrates the frequency response of E(z)/R(z) and one can verify that the ℋ∞-norm in frequency band follows the characteristcs of a tracking system. The M and N optimum controller values wereM=[0.0088-0.0030],N=31.4617.
In the simulation, an unit step input is considered. The simulation result is illustrated in Figure 7.
Frequency response of E(z)/R(z).
Response of the output z(k) to a unit step input r(k).
In the second simulation we use the following ramp signal: r(kT)=0.7kT with T=0.01s. The simulation results are illustrated in Figure 8.
Output of the system for ramp input.
Finally, an input signal r(kT)=sen(0.1kT) and a disturbance signal w(k) with random amplitudes were simulated. The maximum amplitude of the random signal was equal to 1. The simulation results are illustrated in Figure 9.
Output signal z(k) and input signal r(k) are almost overlapped.
For this example, the zeros of the system where −1; 0.4446+0.8004i and 0.4446-0.8004i. The poles of the system with a feedback controller where 0.5004+0.3694i; 0.5004-0.3694i; 0.1436+0.2002i and 0.1436-0.2002i. It is possible to see that the poles of the system with a feedback controller were allocated according to the circle constraint specification. The pole-zero map is shown in Figure 10.
Pole-zero map of the closed-loop system obtained with Kc(z), M and N.
The example above shows the methodology effectiveness. The disturbance rejection and the minimization of the tracking error for specified frequency band were reached. It was showed that the methodology works properly for ramp, unit step and sinusoidal signals for any frequency in the specified frequency band.
7. Example 2
Consider a discrete form of a continuous model plant that has a zero in the right-half plane. A sampling period of 0.01 seconds is used in design. The state-space description for the system is as follows:
[x1(k+1)x2(k+1)x3(k+1)]=[0.868-0.263-0.4480.0090.999-0.00200.010.899][x1(k)x2(k)x3(k)]+[0.009300]u(k)+[0.009300]w(k),y(k)=[01-90][x1(k)x2(k)x3(k)],
where x(k) is the state vector, u(k) is the control signal and w(k) is a disturbance signal in the system.
The design of the tracking system must include operation for reference signals of low frequencies (down to 5 rad/s). In such a case, the filter F(z) was consideredF(z)=(0.033z2+0.127z+0.031)×10-4z3-2.999z2+2.805z-0.905.
Using Theorem 3.1 the controller Kc(z) is designed for the system described in (7.2) and shown in (7.4). This controller minimizes the ℋ∞-norm of w(k) to z(k). In this design a disk of radius 0.85 and center in −0.1 is used as a pole placement constraintKc(z)=(2.38z2-3.63z+1.44)×104z3+3.75z2+15.77z+9.79.
The ℋ∞-norm of w(k) to z(k) for the closed-loop system was 1.577 ∗ 10−3, implying the attenuation of the effect of the disturbance signal in the system. Figure 11 illustrates the frequency response of Z(z)/W(z).
Frequency response of Z(z)/W(z).
To design a tracking system, the proposed zero-variation methodology (5.4) was used in which the ℋ∞ norm of r(k) to e(k) is minimized considering the signal r(k) with low frequencies (down to 5 rad/s). The obtained ℋ∞-norm for all frequency spectra was equal to 1.8, while for frequency band specified in the problem, the largest magnitude of frequency response was 0.031. This implies that the tracking system operated adequately in the frequency band specified in the problem.
Figure 12 illustrates the frequency response of E(z)/R(z) and one can verify that the ℋ∞-norm in frequency band follows the characteristcs of a tracking system. The M and N optimum parameters values wereM=[8.5*10-71.87194.6],N=2746.
Then, an input signal r(kT)=sen(0,1kT) and a disturbance signal w(k) with random amplitudes were simulated, and it was found that the maximum amplitude of the random signal was equal to 1. The simulation results are illustrated in Figure 13.
Frequency response of E(z)/R(z).
Output signal y(k) and input signal r(k) are almost overlapped.
8. Conclusion
In this manuscript, it is proposed a methodology to solve the tracking and disturbance rejection problem applied to discrete-time systems. Considering Figure 1 the disturbance signal acting in the plant can be attenuated by minimizing the ℋ∞-norm from w(k) to z(k), by using a dynamic feedback compensation. In the tracking process, a zero-variation methodology is used in order to minimize the ℋ∞-norm between the reference signal and tracking error signal, where the tracking error is the diference between the reference signal r(k) and system output signal z(k). In the tracking design with disturbance rejection, the pole placement is used to attenuate the disturbance signal effect, while the zero variation allows the tracking. The zero modification do not interfere in the design of the disturbance rejection. In the tracking process, the frequency band wheighted allows to choose the frequency band on the reference input signal. The tracking method and disturbance rejection are based on LMI framework. Then, when there exists a feasible solution the design can be obtained by convergence polynomial algorithms [23, 25] available in the literature.
Acknowledgments
The authors gratefully acknowledge the partial financial support by FAPESP, CAPES and CNPQ of Brazil.
TeixeiraM. C. M.AssunçãoE.CardimR.da SilvaN. A. P.MachadoE. R. M. D.On complementary root locus of biproper transfer functions20092009147279082562899ZBL1183.93074De la SenM.Pole-placement in discrete systems by using single and multirate sampling1996333572174610.1016/0016-0032(96)00040-31408334ZBL0882.93032SaberiA.SannutiP.StoorvogelA. A.H2 optimal controllers with measurement feedback for discrete-time systems: flexibility in closed-loop pole placement199733328930410.1016/S0005-1098(96)00195-11442549ZBL0874.93040M'SaadM.OrtegaR.LandauI. D.Adaptive controllers for discrete-time systems with arbitrary zeros: an overview198521441342310.1016/0005-1098(85)90077-9798186ZBL0571.93041MessnerW. C.KempfC. J.Zero placement for designing discrete-time repetitive controllers1996445635692-s2.0-003012910410.1016/0967-0661(96)00040-8MendesR. A. T.2007São Paulo, BrazilUnespAndreaC. Q.2002São Paulo, BrazilUNESPAssunçãoE.AndreaC. Q.TeixeiraM. C. M.ℋ2 and ℋ∞-optimal control for the tracking problem with zero variation2007136826882-s2.0-3424833663410.1049/iet-cta:20050506ChenB. M.LinZ.LiuK.Robust and perfect tracking of discrete-time systems20023822932992-s2.0-003646673310.1016/S0005-1098(01)00198-4OliveiraM. C.1999Campinas, BrazilUnicampParkJ. H.On design of dynamic output feedback controller for GCS of large-scale systems with delays in interconnections: LMI optimization approach2005161242343210.1016/j.amc.2003.12.0372112415ZBL1066.93003AssunçãoE.PeresP. L. D.A global optimization approach for the ℋ2-norm model reduction problemProceedings of the 38th IEEE Conference on Decision and Control (CDC '99)December 1999Phoenix, Ariz, USA1857186210.1109/CDC.1999.8309042-s2.0-0033325754AssunçãoE.TeixeiraM. C. M.FariaF. A.da SilvaN. A. P.CardimR.Robust state-derivative feedback LMI-based designs for multivariable linear systems20078081260127010.1080/002071707012838992354341ZBL1133.93022TeixeiraM. C. M.AssunçãoE.AvellarR. G.On relaxed LMI-based designs for fuzzy regulators and fuzzy observers20031156136232-s2.0-1074421978310.1109/TFUZZ.2003.817840TeixeiraM. C. M.AssunçãoE.PalharesR. M.Discussion on: “ℋ∞ output feedback control design for uncertain fuzzy systems with multiple time scales: an LMI approach”20051121671692-s2.0-3364539222610.3166/ejc.11.167-170TeixeiraM. C. M.AssunçãoE.PietrobomH. C.On relaxed LMI-based design for fuzzy regulators and fuzzy observersProceedings of the 6th European Control Conference2001Porto, Portugal120125TeixeiraE.AssunçãoE.AvelarR. G.Design of SPR systems with dynamic compensators and output variable structure control1Proceedings of the International Workshop on Variable Structure Systems2006Alghero, Italy328333FariaF. A.AssunçãoE.TeixeiraM. C. M.CardimR.da SilvaN. A. P.Robust state-derivative pole placement LMI-based designs for linear systems200982111210.1080/002071708019421882523729ZBL1154.93366CardimR.TeixeiraM. C. M.AssunçãoE.CovacicM. R.Variable-structure control design of switched systems with an application to a DC-DC power converter2009569350535132-s2.0-6934908864710.1109/TIE.2009.2026381FariaF. A.AssunçãoE.TeixeiraM. C. M.CardimR.Robust state-derivative feedback LMI-based designs for linear descriptor systems20102010159273622562914ZBL1191.93071GahinetP.NemirovskA.LaubA. J.ChialiM.1995Natick, Mass, USAThe Mathworks Inc.FranklinG. F.PowellJ. D.WorkmanM. L.19902ndNew York, NY, USAAddison WesleyChilaliM.GahinetP.ℋ∞ design with pole placement constraints: an LMI approach199641335836710.1109/9.4866371382985ZBL0857.93048LimaE. L.20004thRio de Janeiro, BrazilIMPAColeção Mateméatica UniversitéariaPalharesR. M.TakahashiR. H. C.PeresP. L. D.ℋ∞ and Ótimo ℋ2 guaranteed costs computation for uncertain linear systems19972821831882-s2.0-0031076687