Design of decentralized repetitive control of linear MIMO system

August 19, 2017 | Autor: Edi Kurniawan | Categoria: Control Systems Engineering, Control Systems
Share Embed


Descrição do Produto

Design of Decentralized Repetitive Control of Linear MIMO System Edi Kurniawan, Zhenwei Cao and Zhihong Man Faculty of Engineering and Industrial Sciences, Swinburne University of Technology, Australia [email protected], [email protected], and [email protected]

Abstract— This paper proposes decentralized repetitive control (RC) design for multiple inputs multiple outputs (MIMO) system. The design comprises in both continuous and discrete-time domain. A decentralized stabilizing controller is firstly developed in continuous-time to achieve a desired complimentary sensitivity function. Relative Gain Array (RGA) analysis is used to obtain the best pairings of inputs and outputs. Then, robustness analysis is conducted to determine the impact of neglected couplings. Secondly, the RC is designed in discrete-time based on the information of complimentary sensitivity function. A set of low order, stable and causal repetitive control compensators are developed using a new frequency domain method. The proposed RC compensator works for both minimum and non-minimum phase system. Simulation results for linear MIMO model of pick and place robot arm are presented to validate the effectiveness of the design. Keywords— repetitive control; decentralized control; MIMO system

I.

INTRODUCTION

Tracking periodic commands are common tasks found in many control systems such as disk drive, optical disc player, pick-and-place robot. Repetitive Control (RC) gives superior performance for tracking purpose compared to a nonpredictive control scheme [1]. Repetitive controller requires two main designs; an internal model and a RC compensator design. The internal model refers to the internal model principle originated from [2], which states that the reference model needs to be attached to the controller in order to achieve zero tracking error. The RC compensator is part of RC used to stabilize closed-loop system. Some RC designs to track/reject repetitive signal have been proposed in [3-14]. The RC designs in [3-10] aim to track/reject repetitive signal with fixed period, while the designs in [11-14] are intended to handle time-varying repetitive signal. The designs in [3-10] basically formulated a compensator that is used to compensate the dynamic of single input single output (SISO) model, in which the designs end up with either non-causal or unstable compensator. The designs of RC compensator for MIMO system are still few to be found, and some of them were proposed in [15,16]. Jeong and Fabien [15] proposed the design of Phase Cancellation Inverse Matrix (PCI) which operates by compensating the phase lag in the diagonal elements of the plant model. The idea is initiated from Zero Phase Tracking Error Controller (ZPETC) by [17] targeting to cancel the phase response of the plant. Xu [16]

c 978-1-4673-6322-8/13/$31.00 2013 IEEE

proposed an optimization based compensator to mimics the inverse of the MIMO model. The approaches [15,16] are based on the full MIMO design which results of a controller with the same dimension to the plant. This implies that if we have a ݉‫ ݉ݔ‬MIMO system ( ݉ଶ transfer functions), then we need to design ݉ଶ RC compensators. Moreover, the designs also end up with non-causal compensator which needs to be merged with the internal model to make it realizable. This paper proposes the design of RC compensator for MIMO system based on decentralized control, which means that the MIMO system is considered as a set of SISO systems. The design is motivated by the fact that most of MIMO control problems can be treated on decentralized basis [18]. In addition, the decentralized control is easier as it can apply simpler SISO theories [18]. Another advantage of decentralized RC here is only ݉ RC compensators are required. The design of decentralized RC involves both continuous and discrete-time techniques. Firstly, a decentralized stabilizing controller is designed in continuous time to obtain the desired complementary sensitivity function. Before that, Relative Gain Array (RGA) and Robust analysis are performed to determine whether decentralized control is applicable or not. Secondly, a set of RC compensators are designed in discrete time based on the information of the complimentary sensitivity function. The RC design uses an optimization to obtain a low order, stable and causal compensator. This paper is organized as follows. Section 2 presents the design of decentralized stabilizing controller in continuous time, which also covers the RGA and Robustness analysis. Section 3 covers RC compensator design using optimization. Simulation results are given in Section 4. Section 5 draws the conclusion II.

DECENTRALIZED CONTROL

Let a ݉‫ ݉ݔ‬MIMO system is represented by a following transfer function:

‫ܩ‬ሺ‫ݏ‬ሻ ൌ ൥

݃ଵଵ ሺ‫ݏ‬ሻ ‫ڭ‬ ݃௠ଵ ሺ‫ݏ‬ሻ

‫ڮ‬ ‫ڰ‬ ǥ

݃ଵ௠ ሺ‫ݏ‬ሻ ‫ ڭ‬൩ ݃௠௠ ሺ‫ݏ‬ሻ

(1)

The plant model ‫ܩ‬ሺ‫ݏ‬ሻ has ݉ inputs and ݉ outputs, where

427

the relation between inputs and outputs can be formulated as follows: ‫ݕ‬ଵ ሺ‫ݐ‬ሻ ݃ଵଵ ሺ‫ݏ‬ሻ‫ݑ‬ଵ ሺ‫ݐ‬ሻ ൅ ‫ ڮ‬൅ ݃ଵ௠ ሺ‫ݏ‬ሻ‫ݑ‬௠ ሺ‫ݐ‬ሻ ‫ݕ‬ଶ ሺ‫ݐ‬ሻ ݃ଶଵ ሺ‫ݏ‬ሻ‫ݑ‬ଵ ሺ‫ݐ‬ሻ ൅ ݃ଶଶ ሺ‫ݏ‬ሻ‫ݑ‬ଶ ሺ‫ݐ‬ሻǤ Ǥ ൅݃ଶ௠ ሺ‫ݏ‬ሻ‫ݑ‬௠ ሺ‫ݐ‬ሻ ൦ ൪ൌ൦ ൪ (2) ‫ڭ‬ ‫ڭ‬ ‫ݕ‬௠ ሺ‫ݐ‬ሻ ݃௠ଵ ሺ‫ݏ‬ሻ‫ݑ‬ଵ ሺ‫ݐ‬ሻ൅Ǥ Ǥ ൅݃௠௠ ሺ‫ݏ‬ሻ‫ݑ‬௠ ሺ‫ݐ‬ሻ ሺ‫ݐ‬ሻǡ ǥ ǡ ‫ݕ‬௠ ሺ‫ݐ‬ሻ and ‫ݑ‬ଵ ሺ‫ݐ‬ሻǡ ǥ ǡ ‫ݑ‬௠ ሺ‫ݐ‬ሻ are plant outputs where ‫ݕ‬ଵ and inputs respectively. Decentralized control aims to approximate the MIMO system into a set of independent SISO systems. This is different with decoupling control which tries to convert the full MIMO model into a perfect set of independent SISO models. The assumption in the design of decentralized control is ignoring dynamics that result in weak interactions. Each of the system outputs is approximated from the input response that gives dominant contribution. Therefore, the degree of interaction is necessary to be quantified in the decentralized control. RGA introduced by [19] is one of the techniques termed as dominant interaction control method that can be used to determine the best input output pairings for multivariable control [20]. RGA is defined as matrix Ȧ which is formulated as follows: (3) Ȧ ൌ ‫ܩ‬ሺͲሻǤ‫ כ‬ሾ‫ି ܩ‬ଵ ሺͲሻሿ் ିଵ ሺͲሻ

are system dc gain matrix and its where ‫ܩ‬ሺͲሻ and ‫ܩ‬ inverse, notation Ǥ‫ כ‬and ܶ operate as element wise multiplication and transpose of matrix respectively. For simplicity, suppose the RGA suggestsሾ‫ݑ‬௜ ǡ ‫ݕ‬௜ ሿpairings, where ൌ ͳǡ ǥ ǡ ݉ . This means that we only need to consider transfer function ݃௜௜ ሺ‫ݏ‬ሻ in the design. Now, we need to choose complimentary sensitivity function ሺ•ሻ as the desired stabilized closed-loop model of ݃௜௜ ሺ‫ݏ‬ሻǤ (4) ሺ•ሻ ൌ †‹ƒ‰ሾ– ሺ•ሻǡ ǥ ǡ – ሺ•ሻሿ ଵ

(5)

‫ܩ‬ሺ‫ݏ‬ሻ =‫ܩ‬௡ ሺ‫ݏ‬ሻሾ‫ ܫ‬൅ ‫ ୼ܩ‬ሺ‫ݏ‬ሻሿ where Ͳ ‫ڰ‬ Ͳ

Ͳ Ͳ



(6)

݃௠௠ ሺ‫ݏ‬ሻ

‫ ୼ܩ‬ሺ‫ݏ‬ሻ ൌ  ሾ‫ܩ‬ሺ‫ݏ‬ሻ െ ‫ܩ‬௡ ሺ‫ݏ‬ሻሿൣ‫ܩ‬௡ ିଵ ሺ‫ݏ‬ሻ൧

(7)

and ‫ ܫ‬is a ݉‫ ݉ݔ‬identity matrix. Then, robustness check is given by stability condition as follows [18]: (8) ߪതሾ‫ ܩ‬ሺ݆߱ሻሺŒɘሻሿ ൏ ͳ ‫ א ߱׊‬Թ ୼

Eq. (8) states that the maximum singular values of ሾ‫ ୼ܩ‬ሺ݆߱ሻ •ሺŒɘሻሿ has to be less than one.

428





where …୧ ሺ•ሻǡ ݅ ൌ ͳǡ ǥ ǡ ݉ are obtained to satisfy the following equality …୧ ሺ•ሻ݃௜௜ ሺ‫ݏ‬ሻ (10) ൌ – ୧ ሺ•ሻ ͳ ൅ …୧ ሺ•ሻ݃௜௜ ሺ‫ݏ‬ሻ Transfer function …୧ ሺ•ሻ, ݃௜௜ ሺ‫ݏ‬ሻ, and – ୧ ሺ•ሻ can be expressed in the form of numerator and denumerator as follow: ݊ܿ௜ ሺ‫ݏ‬ሻ ݊݃௜௜ ሺ‫ݏ‬ሻ ݊‫ݐ‬௜ ሺ‫ݏ‬ሻ (11) ǡ ‰ ୧୧ ሺ•ሻ ൌ ǡ – ୧ ሺ•ሻ ൌ ݀ܿ௜ ሺ‫ݏ‬ሻ ݀݃௜௜ ሺ‫ݏ‬ሻ ݀‫ݐ‬௜ ሺ‫ݏ‬ሻ The polynomial ݊ܿ௜ ሺ‫ݏ‬ሻ and ݀ܿ௜ ሺ‫ݏ‬ሻ can be obtained by solving the following Diophantine Equation …୧ ሺ•ሻ ൌ

݊ܿ௜ ሺ‫ݏ‬ሻሾ݊݃௜௜ ሺ‫ݏ‬ሻ݊‫ݐ‬௜ ሺ‫ݏ‬ሻ െ ݊݃௜௜ ሺ‫ݏ‬ሻ݀‫ݐ‬௜ ሺ‫ݏ‬ሿ ൅

(12)

݀ܿ௜ ሺ‫ݏ‬ሻሾ݀݃௜௜ ሺ‫ݏ‬ሻ݊‫ݐ‬௜ ሺ‫ݏ‬ሻሿ ൌ Ͳ The relative degree of complimentary sensitivity function ‫ݐ‬௜ ሺ‫ݏ‬ሻneeds to be chosen carefully so it gives proper stabilizing controller…୧ ሺ•ሻ. III.

DISCRETE RC DESIGN FOR SISO

In this section, we introduce a new RC design using optimization in the frequency domain which is also discussed in [11]. The design is based on the information of complimentary sensitivity function ሺ•ሻ. The general structure of SISO RC system is shown in Fig 1

z −N



In the decentralized control, robustness analysis is needed to check the impact of the neglected dynamics [18]. To perform robustness analysis, the ‫ܩ‬ሺ‫ݏ‬ሻ firstly is represented as the nominal model ‫ܩ‬௡ ሺ‫ݏ‬ሻ with additive uncertainty ‫ ୼ܩ‬ሺ‫ݏ‬ሻ

݃ଵଵ ሺ‫ݏ‬ሻ ‫ܩ‬௡ ሺ‫ݏ‬ሻ ൌ ൥ Ͳ Ͳ

Once robustness check has been performed, then the decentralized stabilizing controller ‫ܥ‬ሺ‫ݏ‬ሻ can be designed based on the chosen complimentary sensitivity function ሺ•ሻ (9) ሺ•ሻ ൌ †‹ƒ‰ሾ… ሺ•ሻǡ ǥ ǡ … ሺ•ሻሿ

Fig. 1. Block diagram of discrete-time RC system

where ‫ܩ‬௥௖ ሺ‫ݖ‬ሻ is the digital repetitive controller, t ሺ‫ݖ‬ሻ is the stabilized plant model, ‫ݎ‬ሺ݇ሻ is the periodic reference signal, ݁ሺ݇ሻ is tracking error, and ‫ݕ‬ሺ݇ሻ is tracking output. The digital repetitive controller has the following transfer function. ܳሺ‫ݖ‬ሻ (13) ‫ܨ‬ሺ‫ݖ‬ሻ ‫ܩ‬௥௖ ሺ‫ݖ‬ሻ ൌ ே ‫ ݖ‬െ ܳሺ‫ݖ‬ሻ where ܰ ൌ

்ೝ ்ೞ

‫ א‬Գ , ܶ௥ is reference signal period, ܶ௦ is the

sampling period, ‫ܨ‬ሺ‫ݖ‬ሻ is the RC compensator, ܳሺ‫ݖ‬ሻ is zero phase low pass filter. The RC system is stable if the following conditions are fulfilled [4,21] : 1) ‫ݐ‬ሺ‫ݖ‬ሻ is stable. గ (14) 2) ห൫ͳ െ ‫ݐ‬ሺ‫ݖ‬ሻ‫ܨ‬ሺ‫ݖ‬ሻ൯ܳሺ‫ݖ‬ሻห ൏ ͳ‫ Ͳ׊‬൏ ߱௜ ൏

2013 IEEE 8th Conference on Industrial Electronics and Applications (ICIEA)

்ೞ

The RC compensator designs mostly aim to mimics the inverse of the SISO plant model, where end up with either non-causal or unstable compensator. Here, we develop a design scheme to obtain a low order, stable and causal compensator. ‫ݍ‬଴ ‫ ݖ‬௠ ൅ ‫ݍ‬ଵ ‫ ݖ‬௠ିଵ ൅ ‫ ڮ‬൅ ‫ݍ‬௠ (15) ǡ ݉ ൐ Ͳ ‫ܨ‬ሺ‫ݖ‬ሻ ൌ ௠ ‫ ݖ‬൅ ‫ݎ‬ଵ ‫ ݖ‬௠ିଵ ൅ ‫ ڮ‬൅ ‫ݎ‬௠ The ‫ݐ‬ሺ‫ݖ‬ሻ is a known stable transfer function, while ܳሺ‫ݖ‬ሻ is a chosen filter to give specific tracking bandwidth. Here, the RC compensator ‫ܨ‬ሺ‫ݖ‬ሻ is simply designed to satisfy stability condition (14) as both ‫ݐ‬ሺ‫ݖ‬ሻ and ܳሺ‫ݖ‬ሻ are already known. Let ‫ܯ‬௧೔ and ߠ௧೔ , ‫ܯ‬ொ೔ andߠொ೔ , ‫ܯ‬ி೔ and ߠி೔ are magnitude and phase response of ‫ݐ‬ሺ‫ݖ‬ሻ , ܳሺ‫ݖ‬ሻ and ‫ܨ‬ሺ‫ݖ‬ሻ at frequency ߱௜ respectively. Denote the magnitude response of ൣ൫ͳ െ ‫ݐ‬ሺ‫ݖ‬ሻ‫ܳݖܨ‬ሺ‫ݖ‬ሻ at frequency ߱݅ as ݄݅ (16) ݄ ൌ หͳ െ ‫ ݁ ܯ‬௝ఏಷ೔ ‫ ݁ ܯ‬௝ఏ೟೔ หห‫ ݁ ܯ‬௝ఏೂ೔ ห ൏ ͳ

The first constraint guarantees that all poles of ‫ܨ‬ሺ‫ݖ‬ሻ are inside the unit circle with a safe minimum distance ߜ to the unit circle. The second constraint guarantees that the closed loop system is stable with a positive margin of ɒ . The optimization problem (20) is a nonlinear optimization to findሺʹ݉ ൅ ͳሻ parameters of ‫ܨ‬ሺ‫ݖ‬ሻ that satisfy both bound and nonlinear constraints.

ܳሺ‫ݖ‬ሻ is a low pass filter with zero phase for all frequencies. Thus, (16) can be rewritten as: (17) ݄ ൌ หͳ െ ‫ ݁ ܯ‬௝ఏಷ೔ ‫ ݁ ܯ‬௝ఏ೟೔ ห‫ܯ‬

where ݃ଵଵ ሺ‫ݏ‬ሻǡ ݃ଵଶ ሺ‫ݏ‬ሻǡ ݃ଶଵ ሺ‫ݏ‬ሻ , and ݃ଶଶ ሺ‫ݏ‬ሻ are given in the bottom of this page. This is a stable and minimum phase MIMO plant as all zeros and poles of ݀݁‫ݐ‬ሾ‫ܩ‬ሺ‫ݏ‬ሻሿ are on the left half plane. The RGA of ‫ܩ‬ሺ‫ݏ‬ሻ is calculated as follows:



ி೔

௧೔



ி೔

SIMULATION RESULTS

Simulation is now performed to validate the effectiveness of the design. A two-input two -output MIMO model of pick and place robot arm [22] is used in the simulation. ‫ܩ‬ሺ‫ݏ‬ሻ ൌ ൤

݃ଵଵ ሺ‫ݏ‬ሻ ݃ଶଵ ሺ‫ݏ‬ሻ

݃ଵଶ ሺ‫ݏ‬ሻ ൨ ݃ଶଶ ሺ‫ݏ‬ሻ

(21)

ொ೔

௧೔

ொ೔

Eq (17) can be further derived to భ

݄௜ ൌ ൣͳ െ ʹ‫்ܯ‬೔ ܿ‫ݏ݋‬൫ߠ்೔ ൯ ൅ ‫்ܯ‬೔ ଶ ൧మ ‫ܯ‬ொ೔

(18)

where ‫்ܯ‬೔ ൌ ‫ܯ‬ி೔ ‫ܯ‬௧೔ , and ߠ்೔ ൌ ߠி೔ ൅ ߠ௧೔ To be stable, the stability condition (14) states that ݄௜ has to be less than one for all frequencies. Here we only consider harmonic frequencies and define the objective function as follows: ேȀଶ

݄ ்௢௧௔௟ ൌ ෍ ݄௜ ‫߱׊‬௜ ൌ ʹߨ ௜ୀଵ

IV.

݅ ݅ ൌ ͳǡʹǤ Ǥ ǡ ܰȀʹ ܰܶ௦

(19)

Then, the following optimization problem is introduced: ݄ ்௢௧௔௟  ‹

‫ ݐ‬ሺ‫ݏ‬ሻ ܶሺ‫ݏ‬ሻ ൌ ൤ ଵ Ͳ

ሺ௣భǡ ǡǥǡ௣೘ǡ௤బ ǡǥǡ௤೘ ሻ

—„Œ‡…––‘ǣ     ‫݌‬ଵ െͳ ൅ ߜ ͳെߜ Ǥ 1. ൦ Ǥ ൪ ൏ ቎ ቏ ൏ ቎ Ǥ ቏ Ǥ Ǥ Ǥ ‫݌‬௠ െͳ ൅ ߜ ͳെߜ ௜ ݅ ൌ ͳǡʹǤ Ǥ ǡ ܰȀʹ 2. ݄௜ ൏ ͳ െ ߬ǡ ‫߱׊‬௜ ൌ ʹߨ



(22) ͳ Ͳ ͳ Ͳ ͳ Ͳ ିଵ ቃ ቃ Ǥ‫ כ‬൬ቂ ቃ ൰ ൌቂ Ͳ ͳ Ͳ ͳ Ͳ ͳ The RGA value Ȧ shows that the dynamic ݃ଵଵ ሺ‫ݏ‬ሻ and ݃ଶଶ ሺ‫ݏ‬ሻ give dominant interaction to output ‫ݕ‬ଵ ሺ‫ݐ‬ሻ and ‫ݕ‬ଶ ሺ‫ݐ‬ሻ respectively. The RGA suggests that the best pairings are ሾ‫ݑ‬ଵ ǡ ‫ݕ‬ଵ ሿ and ሾ‫ݑ‬ଶ ǡ ‫ݕ‬ଶ ሿ. Robust analysis needs to be performed to check the impact of the neglected dynamics ݃ଵଶ ሺ‫ݏ‬ሻ and ݃ଶଶ ሺ‫ݏ‬ሻ . Let the desired complimentary sensitivity function ܶሺ‫ݏ‬ሻ is Ȧൌቂ

(20)

ே்ೞ

Ͳ ൨ ‫ݐ‬ଶ ሺ‫ݏ‬ሻ

ͶͲͲ ‫ ۍ‬ଶ ሺ‫ ݏ‬൅ ͸‫ ݏ‬൅ ͳ͸ሻሺ‫ ݏ‬൅ ʹͷሻ ൌ‫ێ‬ ‫ێ‬ Ͳ ‫ۏ‬

(23) ‫ې‬ ‫ۑ‬ ͶͲͲ ‫ۑ‬ ሺ‫ ݏ‬ଶ ൅ ͸‫ ݏ‬൅ ͳ͸ሻሺ‫ ݏ‬൅ ʹͷሻ‫ے‬ Ͳ

where ‫݌‬ଵǡ ‫݌‬ଶ ǡ ǥ ǡ ‫݌‬௠ are ݉ real poles of ‫ܨ‬ሺ‫ݖ‬ሻ , ߜ and ߬ are

small positive constants. ݃ଵଵ ሺ‫ݏ‬ሻ ൌ ͷ‫ݏ‬ଵଶ

ͷǤ͵݁ െ ݃ଵଶ ሺ‫ݏ‬ሻ ൌ

ͷ‫ݏ‬ଵ଴

ͷǤ͵݁ െ ݃ଶଵ ሺ‫ݏ‬ሻ ൌ

ͲǤͳ͸‫ ݏ‬ଽ ൅ ͳͶǤͷͳ‫ ଼ ݏ‬൅ ͷ͹ͺǤʹ‫ ଻ ݏ‬൅ ͳǤ͵ͻʹ݁Ͷ‫ ଺ ݏ‬൅ ʹǤʹ͸݁ͷ‫ ݏ‬ହ ൅ ʹǤͷͺ݁͸‫ ݏ‬ସ ൅ ʹǤͲͻ݁͹‫ ݏ‬ଷ ൅ ͳǤͳ͹݁ͺ‫ ݏ‬ଶ ൅ ͶǤʹͳ݁ͺ‫ ݏ‬൅ ͹Ǥ͸݁ͺ ൅ ͲǤͲʹ‫ݏ‬ଵଵ ൅ ͲǤͻͳ‫ݏ‬ଵ଴ ൅ ͵ͳǤʹ‫ ݏ‬ଽ ൅ ͹ͳͶǤͳ‫ ଼ ݏ‬൅ ͳǤʹ݁Ͷ‫ ଻ ݏ‬൅ ͳǤͶͷ݁ͷ‫ ଺ ݏ‬൅ ͳǤͶ݁͸‫ ݏ‬ହ ൅ ͳǤͲͳ݁͹‫ ݏ‬ସ ൅ ͷǤ͹݁͹‫ ݏ‬ଷ ൅ ʹǤ͵݁ͺ‫ ݏ‬ଶ ൅ ͷǤͻ݁ͺ‫ ݏ‬൅ ͹Ǥ͸݁ͺ



ͲǤͲͳͶ‫ ݏ‬ଽ

െͲǤͲʹʹ‫ ଻ ݏ‬െ ͵ǤʹͶ‫ ଺ ݏ‬െ ͺͺǤ͵‫ ݏ‬ହ െͳ͵Ͷ͹‫ ݏ‬ସ െ ͳǤͲ͸݁Ͷ‫ ݏ‬ଷ െ ͶǤͷʹ݁Ͷ‫ ݏ‬ଶ ൅ ͲǤ͹ʹ‫ ଼ ݏ‬൅ ʹͲ‫ ଻ ݏ‬൅ ͵͸͵‫ ଺ ݏ‬൅ Ͷ͸Ͷͷ‫ ݏ‬ହ ൅ ͶǤ͵݁Ͷ‫ ݏ‬ସ ൅ ʹǤͻ݁ͷ‫ ݏ‬ଷ ൅ ͳǤͶ݁͸‫ ݏ‬ଶ ൅ ͶǤͳͺ݁͸‫ ݏ‬൅ ͸Ǥ͵ʹ݁͸

െͲǤͳ͸‫ ଻ ݏ‬െ ͺǤ͹‫ ଺ ݏ‬െ ͳͻͶ‫ ݏ‬ହ െʹͶͻͺ‫ ݏ‬ସ െ ͳǤ͹ͺ݁Ͷ‫ ݏ‬ଷ െ ͸Ǥ͸Ͷ݁Ͷ‫ ݏ‬ଶ ͷǤʹ݁ െ ͷ‫ݏ‬ଵ଴ ൅ ͲǤͲͳͶ‫ ݏ‬ଽ ൅ ͲǤ͹ʹ‫ ଼ ݏ‬൅ ʹͲ‫ ଻ ݏ‬൅ ͵͸͵‫ ଺ ݏ‬൅ Ͷ͸Ͷͷ‫ ݏ‬ହ ൅ ͶǤ͵݁Ͷ‫ ݏ‬ସ ൅ ʹǤͻ݁ͷ‫ ݏ‬ଷ ൅ ͳǤͶ݁͸‫ ݏ‬ଶ ൅ ͶǤͳͺ݁͸‫ ݏ‬൅ ͸Ǥ͵ʹ݁͸ ݃ଶଶ ሺ‫ݏ‬ሻ ൌ ͲǤͲʹ͹‫ ݏ‬ଽ ൅ ͶǤͻͷ‫ ଼ ݏ‬൅ ʹ͸Ͷ‫ ଻ ݏ‬൅ ͹͵ͻͶ‫ ଺ ݏ‬൅ ͳǤ͵݁ͷ‫ ݏ‬ହ ൅ ͳǤ͸ͻ݁͸‫ ݏ‬ସ ൅ ͳǤͷ݁͹‫ ݏ‬ଷ ൅ ͻǤͶ݁͹‫ ݏ‬ଶ ൅ ͵Ǥͺ݁ͺ‫ ݏ‬൅ ͹Ǥ͸݁ͺ ͷǤ͵݁ െ ͷ‫ݏ‬ଵଶ ൅ ͲǤͲͳͶ‫ݏ‬ଵଵ ൅ ͲǤͻ‫ ݏ‬ଵ଴ ൅ ͵ͳ‫ ݏ‬ଽ ൅ ͹ͳͶǤͳ‫ ଼ ݏ‬൅ ͳǤͳͻ݁Ͷ‫ ଻ ݏ‬൅ ͳǤͶͺ݁ͷ‫ ଺ ݏ‬൅ ͳǤͶ݁͸‫ ݏ‬ହ ൅ ͳǤͲͶ݁͹‫ ݏ‬ସ ൅ ͷǤ͹݁͹‫ ݏ‬ଷ ൅ ʹǤ͵݁ͺ‫ ݏ‬ଶ ൅ ͷǤͻ݁ͺ‫ ݏ‬൅ ͹Ǥ͸݁ͺ

2013 IEEE 8th Conference on Industrial Electronics and Applications (ICIEA)

429

Both ‫ݐ‬ଵ ሺ‫ݏ‬ሻ and ‫ݐ‬ଶ ሺ‫ݏ‬ሻ are chosen to have relative degree 3, and to give tracking bandwidth 3.7 rad/s. The uncertainty model ‫ ୼ܩ‬ሺ‫ݏ‬ሻ is given as follows: Ͳ ‫ ୼ܩ‬ሺ‫ݏ‬ሻ ൌ  ൤ ݃ଶଵ ሺ‫ݏ‬ሻȀ݃ଵଵ ሺ‫ݏ‬ሻ

݃ଵଶ ሺ‫ݏ‬ሻȀ݃ଶଶ ሺ‫ݏ‬ሻ ൨ Ͳ

(24)

Singular values of ‫ ୼ܩ‬ሺ݆߱ሻሺŒɘሻ is shown in Fig.2

ͺǤ͸͵݁ െ Ͷሺ‫ ݖ‬൅ ͵ǤͲͻሻሺ‫ ݖ‬൅ ͲǤʹʹሻ (28) ሺ‫ ݖ‬െ ͲǤͷ͵ሻሺ‫ ݖ‬ଶ െ ͳǤͺͷ‫ ݖ‬൅ ͲǤͺ͸ሻ Both ‫ݐ‬ଵ ሺ‫ݖ‬ሻ and ‫ݐ‬ଶ ሺ‫ݖ‬ሻare non-minimum phase model because they have a zero outside the unit circle. The Q-filterܳሺ‫ݖ‬ሻ, the order of compensator݊, and positive constants ߜ and ߬ are chosen respectively as follow: ܳଵ ൌ ܳଶ ሺ‫ݖ‬ሻ ൌ ͲǤʹͷ‫ ݖ‬൅ ͲǤͷ ൅ ͲǤʹͷ‫ି ݖ‬ଵ , ݊ ൌ ͵, (29) ߜ ൌ ͲǤͲ͹ͷ, and ߬ ൌ ͲǤͲͷ The RC compensator ‫ܨ‬ሺ‫ݖ‬ሻ is obtained by solving optimization problem (20) using an Optimization Toolbox from Matlab. ‫ݐ‬ଵ ሺ‫ݖ‬ሻ ൌ ‫ݐ‬ଶ ሺ‫ݖ‬ሻ ൌ

‫ ܨ‬ሺ‫ݖ‬ሻ ‫ܨ‬ሺ‫ݖ‬ሻ ൌ ൤ ଵ Ͳ

Ͳ ൨ ‫ܨ‬ଶ ሺ‫ݖ‬ሻ

(30)

where ͵ͻͶǤ͹‫ ݖ‬ଷ െ ͻʹ͹Ǥʹ‫ ݖ‬ଶ ൅ ͹ͲͳǤ͸‫ ݖ‬െ ͳ͸͸Ǥ͸ (31) ‫ ݖ‬ଷ ൅ ʹǤͲͺ‫ ݖ‬ଶ ൅ ͳǤʹͺ‫ ݖ‬൅ ͲǤͳͻ In this simulation, two tracking schemes are presented. The first scheme is the first channel ‫ݕ‬ଵ is required to track triangle reference signal, while the second channel ‫ݕ‬ଶ needs to stay idle. The second scheme is both channels are required to track triangle reference signals with different amplitudes. The tracking output and tracking error for the first scheme is shown in Fig 3. ‫ܨ‬ଵ ൌ ‫ܨ‬ଶ ሺ‫ݖ‬ሻ ൌ

Fig. 2. Singular values of ۵ઢ ሺ‫ܒ‬૑ሻ‫܂‬ሺ‫ܒ‬૑ሻ

Fig. 2 shows that maximum singular values of ‫ ୼ܩ‬ሺ݆߱ሻሺŒɘሻ is still less thanͲ݀‫ܤ‬. Thus, the impact of the neglected dynamics is acceptable. The decentralized stabilizing controller ‫ܥ‬ሺ‫ݏ‬ሻ is obtained by solving the following Diophantine Equations ݊ܿଵ ሺ‫ݏ‬ሻሾ݊݃ଵଵ ሺ‫ݏ‬ሻ݊‫ݐ‬ଵ ሺ‫ݏ‬ሻ െ ݊݃ଵଵ ሺ‫ݏ‬ሻ݀‫ݐ‬ଵ ሺ‫ݏ‬ሿ ൅

(25)

݀ܿଵ ሺ‫ݏ‬ሻሾ݀݃ଵଵ ሺ‫ݏ‬ሻ݊‫ݐ‬ଵ ሺ‫ݏ‬ሻሿ ൌ Ͳ ݊ܿଶ ሺ‫ݏ‬ሻሾ݊݃ଶଶ ሺ‫ݏ‬ሻ݊‫ݐ‬ଶ ሺ‫ݏ‬ሻ െ ݊݃ଶଶ ሺ‫ݏ‬ሻ݀‫ݐ‬ଶ ሺ‫ݏ‬ሿ ൅

(a)

(26)

݀ܿଶ ሺ‫ݏ‬ሻሾ݀݃ଶଶ ሺ‫ݏ‬ሻ݊‫ݐ‬ଶ ሺ‫ݏ‬ሻሿ ൌ Ͳ The controller ‫ܥ‬ሺ‫ݏ‬ሻ is ܿଵ ሺ‫ݏ‬ሻ Ͳ (27) ൨ Ͳ ܿଶ ሺ‫ݏ‬ሻ whereܿଵ ሺ‫ݏ‬ሻ and ܿଶ ሺ‫ݏ‬ሻ are shown in the bottom of this page Let the frequency of reference signal and sampling rate be 0.2 Hz and 40 Hz respectively. The discretized model of ‫ݐ‬ଵ ሺ‫ݏ‬ሻ and ‫ݐ‬ଶ ሺ‫ݏ‬ሻ are: ‫ܥ‬ሺ‫ݏ‬ሻ ൌ  ൤

ܿଵ ሺ‫ݏ‬ሻ ൌ

(b)

͵Ǥͻ݁ െ ͺ‫ ݏ‬ଵଶ ൅ ͳǤͳ݁ െ ͷ‫ ݏ‬ଵଵ ൅ ͸Ǥͺ݁ െ Ͷ‫ ݏ‬ଵ଴ ൅ ͲǤͲʹ‫ ݏ‬ଽ ൅ ͲǤͷ͵‫ ଼ ݏ‬൅ ͺǤͻ‫ ଻ ݏ‬൅ ‫ڮ‬ ʹǤͻ݁ െ ͹‫ݏ‬ଵଶ ൅ ͵Ǥ͸݁ െ ͷ‫ ݏ‬ଵଵ ൅ ͳǤͺ݁ െ ͵‫ ݏ‬ଵ଴ ൅ ͲǤͲ͸‫ ݏ‬ଽ ൅ ͳǤͶ‫ ଼ ݏ‬൅ ʹʹǤʹ‫ ଻ ݏ‬൅ ʹͷͺǤ͹‫ ଺ ݏ‬൅ ‫ڮ‬ ͳͲͺǤͶ݁Ͷ‫ ଺ ݏ‬൅ ͳͲͶ͹‫ ݏ‬ହ ൅ ͹ͷͷʹ‫ ݏ‬ସ ൅ ͶǤ͵݁Ͷ‫ ݏ‬ଷ ൅ ͳǤ͹݁ͷ‫ ݏ‬ଶ ൅ ͶǤͶ݁ͷ‫ ݏ‬൅ ͷǤ͹݁ͷ ʹʹ͵Ͳ‫ ݏ‬ହ ൅ ͳǤͶ݁Ͷ‫ ݏ‬ସ ൅ ͸Ǥʹ݁Ͷ‫ ݏ‬ଷ ൅ ͳǤ͹ͷ݁ͷ‫ ݏ‬ଶ ൅ ʹǤ͵͸‫ݏ‬

ܿଶ ሺ‫ݏ‬ሻ ൌ

ʹǤͷ݁ െ ͳʹ‫ݏ‬ଵଶ ൅ ͸Ǥͺ݁ െ ͳͲ‫ݏ‬ଵଵ ൅ ͶǤͶ݁ െ ͺ‫ݏ‬ଵ଴ ൅ ͳǤͷ݁ െ ͸‫ ݏ‬ଽ ൅ ͵Ǥͷ݁ െ ͷ‫ ଼ ݏ‬൅ ͷǤ͹݁ െ Ͷ‫ ଻ ݏ‬൅ ‫ڮ‬ ͵Ǥʹ͹݁ െ ͳʹ‫ݏ‬ଵଶ ൅ ͹݁ െ ͳͲ‫ݏ‬ଵଵ ൅ ͷǤͳ݁ െ ͺ‫ݏ‬ଵ଴ ൅ ͳǤͻͺ݁ െ ͸‫ ݏ‬ଽ ൅ ͶǤͺͺ݁ െ ͷ‫ ଼ ݏ‬൅ ͺǤͶ݁ െ Ͷ‫ ଻ ݏ‬൅ ‫ڮ‬ ͹Ǥʹ݁ െ ͵‫ ଺ ݏ‬൅ ͸Ǥͺ݁ െ ʹ‫ ݏ‬ହ ൅ ͲǤͷ‫ ݏ‬ସ ൅ ʹǤ͹‫ ݏ‬ଷ ൅ ͳͳǤͳ‫ ݏ‬ଶ ൅ ʹͺǤͷ‫ ݏ‬൅ ͵͸Ǥͺ ͳǤͲ͹݁ െ ʹ‫ ଺ ݏ‬൅ ͲǤͳ‫ ݏ‬ହ ൅ ͲǤ͹‫ ݏ‬ସ ൅ ͵ǤͶ‫ ݏ‬ଷ ൅ ͳͲǤͶͺ‫ ݏ‬ଶ ൅ ͳͷǤʹ͸‫ݏ‬

430

2013 IEEE 8th Conference on Industrial Electronics and Applications (ICIEA)

V.

(c) Fig. 3. The first scheme a) Tracking output ࢟૚ ,b) Tracking output ࢟૛ , c).Tracking error ࢋ૚ and ࢋ૛

The first scheme aims to emphasize the effect of neglected dynamics in this decentralized control. On Channel 1, tracking output ‫ݕ‬ଵ is able to follow periodic reference ‫ݎ‬ଵ after 4 repetitions. On channel 2, a small excitation due to coupling effect appears. The tracking output and tracking error for the second scheme is shown in Fig 4.

This paper proposes decentralized RC design for MIMO system. The decentralized stabilizing controller is firstly developed in continuous-time to achieve the desired complimentary sensitivity function. RGA analysis is used to obtain the best pairing of inputs and outputs. Then, robustness analysis is conducted to determine the impact of neglected couplings. Secondly, the RC is designed in discrete-time based on the information of complimentary sensitivity. The RC design uses an optimization to obtain a low order, stable and causal compensator. Simulation results have been presented to validate the effectiveness of the design. REFERENCES [1]

[2]

[3]

[4]

[5] (a) [6]

[7]

[8]

(b)

[9]

[10]

[11]

[12] (c) Fig. 4. The second scheme a) Tracking output࢟૚ ,b) Tracking output ࢟૛ , c).Tracking error ࢋ૚ and ࢋ૛ [13]

Fig. 4 shows good tracking performance for both channel 1 and 2. The simulation results verify the effectiveness of decentralized RC.

CONCLUSION

[14]

[15]

G. Hillerström and J. Sternby, "Repetitive Control Theory and Applications - A Survey," Proceedings of the 13th IFAC World Congress, vol. D, pp. 1-6, 1996. B. A. Francis and W. M. Wonham, "The internal model principle for linear multivariable regulators," Applied Mathematics & Optimization, vol. 2, pp.170-194, 1975. M. Tomizuka, T.-C.Tsao, and K.-K. Chew, "Analysis and Synthesis of Discrete-Time Repetitive Controllers," Journal of Dynamic Systems, Measurement, and Control, vol. 111, pp. 353-358, 1989. K. K. Chew and M. Tomizuka, "Digital control of repetitive errors in disk drive systems," IEEE Control Systems Magazine, vol. 10, pp. 1620, 1990. C. Cosner, G. Anwar, and M. Tomizuka, "Plug in repetitive control for industrial robotic manipulators,". Proceedings of IEEE International Conference on Robotics and Automation, vol. 3, pp. 1970-1975, 1990. G. F. Ledwich and A. Bolton, "Repetitive and periodic controller design," IEE Proceedings D on Control Theory and Applications, vol. 140, pp. 19-24, 1993. G. Hillerstrom and J. Sternby, "Application of Repetitive Control to a Peristaltic Pump," Proceedings of American Control Conference, pp. 136-141, 1993. B. Panomruttanarug and R. W. Longman, "Repetitive controller design using optimization in the frequency domain," Proceedings of AIAA/AAS Astrodynamics Specialist Conference, pp. 1215-1236, 2004. B. Zhang, D. Wang, K. Zhou, and Y. Wang, "Linear phase lead compensation repetitive control of a CVCF PWM inverter," IEEE Transactions on Industrial Electronics, vol. 55, pp. 1595-1602, 2008. X. H. Wu, S. K. Panda, and J. X. Xu, "Design of a Plug-In Repetitive Control Scheme for Eliminating Supply-Side Current Harmonics of Three-Phase PWM Boost Rectifiers Under Generalized Supply Voltage Conditions," IEEE Transactions on Power Electronics, vol. 25, pp. 1800-1810, 2010. E. Kurniawan, Z. Cao, and Z. Man, "Robust design of repetitive control system," Proceedings of 37th Annual Conference on IEEE Industrial Electronics Society, 2011, pp. 722-727. E. Kurniawan, Z. Cao, and Z. Man, "Adaptive Repetitive Control of System Subject to Periodic Disturbance with Time-Varying Frequency," Proceedings of First International Conference on Informatics and Computational Intelligence , 2011, pp. 185-190. Z. Cao and S. S. Narasimhulu, "Digital PLL-based adaptive repetitive control," Proceedings of 1st International Symposium on Systems and Control in Aerospace and Astronautics, Vols 1and 2, pp. 1468-1471, 2006. Z. Cao and G.F. Ledwich, "Tracking variable periodic signals with fixed sampling rate," Proceedings of the 40th IEEE Conference on Decision and Control , 2001, vol.5, pp. 4885-4890. D. Jeong and B. C. Fabien, "A discrete time repetitive control system for MIMO plants," Proceedings of the American Control Conference, 1999, vol.6, pp. 4295-4299.

2013 IEEE 8th Conference on Industrial Electronics and Applications (ICIEA)

431

[16] K. Xu, "Multi-Input Multi-Output Repetitive Control Theory And Taylor Series Based Repetitive Control Design," Phd thesis, Graduate School of Arts and Sciences –Columbia University, 2012. [17] M. Tomizuka, "Zero phase error tracking algorithm for digital control," Journal of Dynamic Systems, Measurement and Control, Transactions of the ASME, vol. 190, pp. 65-68, 1987. [18] G. C. Goodwin, S. F. Graebe, and M. E. Salgado, Control System Design , Prentice Hall, 2000. [19] E. Bristol, "On a new measure of interaction for multivariable process control," IEEE Transactions on Automatic Control, vol. 11, pp. 133134, 1966. [20] M. T. Tham, "Multivariable Control: An Introduction to Decoupling Control " Dept.Of Chemical and Process Engineering, University of Newcastle upon Tyne, 1999. [21] K. K. Chew and M. Tomizuka, "Steady-state and stochastic performance of a modified discrete-time prototype repetitive controller” Journal of Dynamic Systems, Measurement and Control, Transactions of the ASME, vol. 112, pp. 35-41, 1990. [22] L. Wang, S. Chai, E. Rogers, and C. T. Freeman, "Multivariable Repetitive-Predictive Controllers Using Frequency Decomposition," IEEE Transactions on Control Systems Technology, vol. 20, pp. 15971604, 2012.

432

2013 IEEE 8th Conference on Industrial Electronics and Applications (ICIEA)

Lihat lebih banyak...

Comentários

Copyright © 2017 DADOSPDF Inc.