Điều khiển phi tuyến hệ thống phát điện chạy sức gió sử dụng máy phát không đồng bộ nguồn kép trên cơ sở hệ thụ động Euler - Lagrange và Hamilton

Tạp chí Tin học và Điều khiển học, T.28, S.1 (2012), 9–19
ĐIỀU KHIỂN PHI TUYẾN HỆ THỐNG PHÁT ĐIỆN CHẠY SỨC GIÓ
SỬ DỤNG MÁY PHÁT KHÔNG ĐỒNG BỘ NGUỒN KÉP TRÊN CƠ SỞ
HỆ THỤ ĐỘNG EULER - LAGRANGE VÀ HAMILTON
ĐẶNG DANH HOẰNG1, NGUYỄN PHÙNG QUANG2
1Trường Đại học Kỹ thuật Công nghiệp Thái nguyên
2Trường Đại học Bách khoa Hà Nội
Tóm tắt. Việc áp dụng phương pháp điều khiển mới để điều khiển máy phát điện không đồng bộ
nguồn kép (DFIG) trong hệ thống máy phát điện sức gió nhằm đánh giá khả năng ứng dụng vào thực
tiễn có một ý nghĩa hết sức quan trọng. Bài báo trình bày kết quả nghiên cứu vận dụng nguyên lý
tựa theo thụ động (passivity - based) để thiết kế bộ điều khiển bằng phương pháp kết hợp tựa theo
hệ thụ động Euler - Lagrange (EL) và luật Hamiltonian để điều khiển máy phát điện không đồng bộ
nguồn kép. Bộ điều khiển mới khắc phục được sai lệch tĩnh của bộ điều khiển dòng khi chưa kết hợp
với luật Hamiltonian. Ngoài ra với bộ điều khiển này hệ thống làm việc có chất lượng tốt ở các chế
độ trên và dưới tốc độ đồng bộ của máy phát.
Abstract. Applying a novel control method to control Doubly-Fed Induction Generators (DFIG)
in wind power systems for evaluating an applicability plays an important role in practice. In this
paper, the passivity based-principle is applied to design a controller for DFIG, in which the Euler
- Lagrange (EL) based-method and Hamiltonian rules are incorporated. The proposed controller
overcomes static error of the current controller without Hamiltonian rules. The system can also offer
exellent performances above and below synchronous speeds of the DFIG.
Ký hiệu
Ký hiệu Đơn vị Ý nghĩa
R(x) Ma trận suy giảm
J(x) Ma trận liên kết các khối trong cấu trúc hệ thống
G(x) Ma trận biểu diễn quan hệ vào ra
Lr H Điện cảm của rotor
Tr, Ts s Hằng số thời gian rotor và stator
ωr, ω rad/s Vận tốc góc rotor, vận tốc góc cơ học rotor
ψsd, ψsq Wb = Vs Thành phần d và q của từ thông stator
ψrd, ψrq Wb = Vs Thành phần d và q của từ thông rotor
σ Hệ số tản tổng
Lm H Điện cảm hỗ cảm giữa stator và rotor
Chữ viết tắt
10 ĐẶNG DANH HOẰNG, NGUYỄN PHÙNG QUANG
DFIG Máy phát điện không đồng bộ nguồn kép
PĐSG Phát điện sức gió
EL Euler - Lagrange
PBC Điều khiển thụ động
1. MỞ ĐẦU
Ở nước ta, cũng như trên thế giới việc điều khiển máy phát điện không đồng bộ nguồn kép
trong hệ thống máy phát điện sức gió (PĐSG) đang là một vấn đề được quan tâm. Hiện nay
đã có một số tác giả sử dụng các phương pháp điều khiển như tuyến tính hoá chính xác [6],
backstepping [2] và đã có những kết quả nhất định. Bên cạnh đó tác giả cũng đã có một số
công trình nghiên cứu [3, 4]. Bài báo này đưa ra phương pháp điều khiển passivity – based kết
hợp hệ thụ động EL và luật Hamiltonain để điều khiển máy phát điện không đồng bộ nguồn
kép nhằm mục tiêu:
- Đưa ra thuật toán điều khiển mới nhằm khử sai lệch tĩnh ([4] chưa đề cập đến).
- Điều khiển đảm bảo chất lượng làm việc của hệ thống ở tốc độ máy phát trên và dưới đồng
bộ.
Thực tế khi làm việc thông số máy phát có thể bị thay đổi dẫn đến khi làm việc tồn tại
sai lệch tĩnh. Vì vậy việc khắc phục sai lệch tĩnh là cần thiết. Bên cạnh đó trong quá trình
làm việc máy phát điện sức gió phụ thuộc vào tốc độ gió tự nhiên nên thường xuyên bị thay
đổi tốc độ. Do đó việc điều khiển đảm bảo chất lượng làm việc của hệ thống ở các tốc độ máy
phát khác nhau là hết sức quan trọng. Trong nghiên cứu này sẽ tập trung vào việc giải quyết
2 vấn đề trên.
2. NỘI DUNG CHÍNH
2.1. Cơ sở lý luận phương pháp điều khiển tựa theo thụ động
2.1.1. Hệ thụ động EL
Hệ Euler - Lagrange thụ động là hệ mà động học của chúng được mô tả bởi các phương trình
Euler-Lagrange (EL) và bản thân hệ thống không tự sinh ra năng lượng. Như vậy khi nhắc
đến hệ Euler-Lagrange ta hiểu rằng đó là hệ có bản chất thụ động.
Điều khiển tựa theo thụ động (Passivity Based Control - PBC) là thuật toán điều khiển mà
nguyên lý của nó dựa trên đặc điểm thụ động của đối tượng (hệ hở) với mục tiêu làm cho hệ
kín cũng là một hệ thụ động với hàm lưu giữ năng lượng mong muốn.
Xét một hệ động học có n bậc tự do, động học của hệ có thể được mô tả bởi phương trình
Euler - Lagrange có dạng sau [11]:
d
dt
(
∂L
∂x˙
(x,x˙)
)
−
∂L
∂x
(x,x˙) = Q, (2.1)
trong đó:
• x = (x1, x2, .., xn)T và x˙ là vector trạng thái và đạo hàm vector trạng thái của hệ thống,
• L(x, x˙) được gọi là hàm Lagrangian,
THIẾT KẾ BỘ ĐIỀU KHIỂN MÁY PHÁT ĐIỆN KHÔNG ĐỒNG BỘ NGUỒN KÉP 11
•
Q = −
∂F
∂x˙
(x˙) +B.u+Qn (2.2)
là vector lực tác động lên hệ thống.
với: F (x˙)được gọi là hàm tiêu thụ Rayleigh, và thoả mãn:
x˙T
∂F
∂x˙
(x˙) ≥ 0. (2.3)
• Lực tác động điều khiển B.u ∈ Rn với u ∈ Rnlà vector điều khiển và B ∈ Rn×nu là ma
trận hằng,
• Tác động do nhiễu Qn.
Xét một hệ được ký hiệu là Σ có hàm tổng lưu giữ năng lượng H(x, x˙), vector điều khiển đầu
vào, đầu ra u, y và coi như hệ không chịu tác động của nhiễu. Như vậy tốc độ cung cấp năng
lượng cho hệ sẽ là yTu. Hệ trên được gọi là thụ động nếu:
∫ T
0
yTudt︸ ︷︷ ︸
nang luong cap
≥ H(x(T), x˙(T))−H(x(0), x˙(0))︸ ︷︷ ︸
nang luong luu giu
Thật vậy từ (2.1), (2.2), (2.3) sau khi biến đổi ta được:
H [T ]−H [0]︸ ︷︷ ︸
NL luu giu
+
∫ T
0
x˙T
∂F (x˙)
∂x˙
dt︸ ︷︷ ︸
NL tieu hao
=
∫ T
0
x˙TBudt︸ ︷︷ ︸
NL cung cap
(2.4)
Do điều kiện (2.3), nên H [T]−H[0] ≤
∫ T
0 y
TBudt; (y = x˙) suy ra hệ EL là hệ thụ động,
và một tính chất đặc biệt [5, 11] khi phân tích hệ EL thành các hệ EL con cũng như hệ kín
(có bộ điều khiển) đều thoả mãn là thụ động. Đây là một trong những đặc điểm quan trọng
khi thiết kế bộ điều khiển theo phương pháp PBC.
Từ phương trình (2.4) ta có một số nhận xét sau:
• Nếu u = 0 thì năng lượng của hệ không tăng, vì vậy hệ sẽ ổn định tại trạng thái cân
bằng "tầm thường".
• Hệ sẽ vẫn ổn định nếu như đầu ra Bx˙Tbằng không, trong hệ tuyến tính thì hệ thống
được gọi là pha cực tiểu (minimum phase), tức là hệ ổn định Lyapunov.
• Ta thấy rằng tín hiệu suy giảm có thể được phun vào một cách dễ dàng qua các trạng
thái, được tác động trực tiếp bởi tín hiệu điều khiển nếu như các trạng thái đó có thể
đo được.
12 ĐẶNG DANH HOẰNG, NGUYỄN PHÙNG QUANG
2.1.2. Hệ Hamiltonian
Các hệ thống điều khiển theo hệ Hamiltonian viết dưới dạng tường minh [1]{
x˙ = (J(x) − R(x))∂xH(x) + G(x)u
y = GT(x)∂xH(x)
(2.5)
trong đó:
• x - là vecter trạng thái
• u, y - là các biến vào, ra
• H(x) - là hàm Hamiltonian (thường được để biểu diễn hàm năng lượng của hệ thống)
• ∂x - là phép lấy gradient theo x (do người thiết kế)
Mệnh đề.
• 1.1) Hàm H(x) được giới hạn bởi:
H(x) > c khi x* là nhỏ nhất.
• 1.2) Ma trận J(x) = −JT (x).
• 1.3| R(x) = RT (x) ≥ 0.
Chứng minh.
Lấy đạo hàm hàm Hamiltonian ta được:
H˙(x) = (∂H)Tx˙ = (∂H)T(J(x)−R(x))∂H + (∂H)TG(x)u
.
Từ mệnh đề (1.2) ta được :
H˙(x) = −(∂H)TR(x)∂H + yTu
.
Mặt khác từ mệnh đề (1.3) và xét khi u = 0:
H˙(x) ≤ 0
Từ kết quả trên và mệnh đề (1.2), ta thấy hàm Hamiltonian thỏa mãn là hàm Lyapunov. Như
vậy với việc tính hàm Hamiltonian, ta có thể rút ra được các công việc thiết kế bộ điều khiển
cho hệ thống sau này.
2.2. Cấu trúc hệ thống điều khiển
Theo [2, 4, 6], hệ thống gồm 2 phần điều khiển cơ bản như Hình 2.1.
• Điều khiển phía máy phát sử dụng máy điện không đồng bộ nguồn kép (Doubly -fed
induction machines - DFIG).
THIẾT KẾ BỘ ĐIỀU KHIỂN MÁY PHÁT ĐIỆN KHÔNG ĐỒNG BỘ NGUỒN KÉP 13
Hình 2.1. Sơ đồ cấu trúc hệ thống máy phát điện sử dụng DFIG
• Điều khiển phía lưới.
NLPL: Nghịch lưu phía lưới,
NLPMP: Nghịch lưu phía máy phát,
MĐC: Máy đóng cắt,
IE: Thiết bị đo tốc độ bằng khắc vạch xung.
2.3. Áp dụng phương pháp kết hợp hệ thụ động EL và Hamiltonian để thiết kế
bộ điều khiển
2.3.1. Thiết kế bộ điều khiển dòng điện rôto phía máy phát
Thực hiện theo 2 bước:
• Thiết kế bộ điều khiển tựa theo hệ thụ động EL.
• Kết hợp với luật điều khiển theo hàm Hamiltonian để xây dựng bộ điều khiển mới trên
cơ sở của bộ điều khiển EL.
• Bước 1:
Để áp dụng phương pháp ta tách hệ thống phía rotor máy phát điện thành hai phần động
học phần điện (hàm năng lượng He) và động học phần cơ (hàm năng lượng Hm) - Hình 2.2.
Sau đó áp dụng các phương trình động học đó vào phương trình EL, để hệ trở thành thụ động
[11]:
Từ Hình 2.2 ta xây dựng sơ đồ nguyên lý cấu trúc điều khiển theo phương pháp tựa theo
thụ động như Hình 2.3. Ta có thể cụ thể hóa bằng hình 2.4
Áp dụng phương pháp thiết kế là đưa bộ điều khiển vào hệ động học phần điện với tương tác
của hệ động học phần cơ, sao cho hệ kín thoả mãn là thụ động theo phương trình EL, ta được:
uPBCr = u
*
r −D(ω)(ir−i
∗
r ). (2.6)
với
D(ω) =
L2m
4ε
ω2 + d, 0 < ε < Rr, d ≥ 0. (2.7)
14 ĐẶNG DANH HOẰNG, NGUYỄN PHÙNG QUANG














Hình 2.2. Phân tích DFIG thành động học phần điện và phần cơ
PBC
IR
e
H
m
H








ứ 	

ω
PBC
r
u riri
PBC
I
R : B

 iu chnh dòng 
theo ph

ng pháp PBC
Hình 2.3. Sơ đồ nguyên lý cấu trúc điều khiển theo phương pháp PBC
	




e
H
m
H

 








ω
PBC
r
u
*
r
i ff
fi
fl
ffi

!
"
#
$
fl

%

&
ffi'
(
)*
(+fl
(,
-
fl
.
/
!
ff
fi
fl
ffiffi
0

-
,1

2
(
3
ω


4ω5ψ4


*
r
u
()D ω
r
i
6
7
8
9
:;
<
=
>
< PBC
I
R
r
i
r
i
Hình 2.4. Cấu trúc bộ điều chỉnh véc tơ dòng PBC bao gồm 2 khối chức năng
Như trong [3, 4] hệ phương trình mô tả mô hình dòng rotor của (DFIG) sau khi được tách ra
trên hệ trục toạ độ dq như sau:


dird
dt
= −
1
σ
(
1
Tr
+
1− σ
Ts
)ird + ωrirq +
1− σ
σ
.
1
Ts
ψ
′
sd
−
1−σ
σ
.ωψ
′
sq
+ 1
σLr
urd −
1−σ
σLm
.u
sd
dirq
dt = −
1
σ (
1
Tr
+ 1−σTs )irq − ωrird +
1−σ
σ
1
Ts
ψ
′
sq
+1−σ
σ
ωψ
′
sd
+ 1
σLr
urq −
1−σ
σLm
.usq
(2.8)
Để đặt bài toán điều khiển ta đặt biến ir là biến điều khiển, với giá trị mong muốn là ir∗
được lấy từ bộ điều chỉnh mô men mG và công suất Q.
THIẾT KẾ BỘ ĐIỀU KHIỂN MÁY PHÁT ĐIỆN KHÔNG ĐỒNG BỘ NGUỒN KÉP 15
Bộ điều khiển tựa theo thụ động (PBC) được xây dựng theo nguyên tắc cần phải đưa vào tín
hiệu điều khiển một tín hiệu suy giảm dạng D(ω).ir và D(ω) gọi là hệ số suy giảm, được xác
định theo (2.6).
Tín hiệu điều khiển được xác định:
uPBC
rd
= u∗rd +D(ω).(ird− i
∗
rd),
uPBC
rq
= u∗rq +D(ω).(irq − i
∗
rq),
(2.9)
trong đó:
• uPBC
rd
; uPBC
rq
là điện áp do bộ điều khiển PBC tạo ra (theo d và q).
• urd
∗; urq
∗: điện áp rotor mong muốn của máy phát (theo d và q) được xác định theo
(2.8).
Với phương pháp trên ta được bộ điều khiển dòng điện rotor theo 2 thành phần:
uPBC
rd
= σLr
di∗
rd
dt
+ Lr(
1
Tr
+
1− σ
Ts
)i∗
rd
−σLrωri
∗
rq
− Lr.
1− σ
Ts
(ψ
′
∗
sd − ωTsψ
′
∗
sq)
+Lr
1− σ
Lm
.u∗
sd
−D(ω).(ird − i
∗
rd),
(2.10)
uPBC
rq
= σLr
di∗
rq
dt
+ Lr(
1
Tr
+ 1−σ
Ts
)i∗
rq
+σLrωri
∗
rd
− Lr
1−σ
Ts
(ψ
′
∗
sq
+ ωTsψsd
′
∗)
+Lr
1−σ
Lm
.u∗
sq
−D(ω).(irq − i
∗
rq) .
(2.11)
• Bước 2:
Sau khi có được bộ điều khiển (2.10) và (2.11) tựa theo hệ thụ động EL. Ta tiếp tục thiết kế
bộ điều khiển được tựa theo luật Hamiltonian:
∫ t
0
uT (τ)y(τ)dτ︸ ︷︷ ︸
Nang luong dua tu ben ngoai
= H(x(t))−H(x(0))︸ ︷︷ ︸
Nang luong du tru
+ d(t)︸︷︷︸
Nang luong tieu thu
.
Theo [1], để hệ kín là thụ động, ta chọn luật điều khiển là thêm vào bộ điều chỉnh khâu tích
phân:
uPBC
r(H)
= uPBC
r
−KI
∫
(ir − i
∗
r)dt. (2.12)
với KI > 0.
Thay (2.12) vào (2.10) và (2.11) ta được:
uPBC
rd(H)
= Lr(
1
Tr
+ 1−σ
Ts
)i∗
rd
+ σLr
di∗
rd
dt
−σLrωri
∗
rq
− Lr.
1−σ
Ts
(ψ
′
∗
sd − ωTsψ
′
∗
sq)
+Lr
1−σ
Lm
.u∗
sd
−D(ω)(ird − i
∗
rd
)−KI
∫
(ird − i
∗
rd
)dt,
(2.13)
16 ĐẶNG DANH HOẰNG, NGUYỄN PHÙNG QUANG
uPBC
rq(H)
= Lr(
1
Tr
+ 1−σ
Ts
)i∗
rq
+ σLr
di∗
rq
dt
+σLrωri
∗
rd
− Lr
1−σ
Ts
(ψ
′
∗
sq
+ ωTsψ
′
∗
sd
)
+Lr
1−σ
Lm
.u∗
sq
−D(ω)(irq − i
∗
rq)−KI
∫
(irq − i
∗
rq)dt.
(2.14)
Từ các bộ điều chỉnh dòng dựa trên thụ động ta thấy nó đã đảm bảo tách kênh nhờ bù liên
kết ngang bởi 2 thành phần ωr.ird∗ và ωr.irq∗ và bù được các thông số khác như điện áp lưới,
từ thông stator, tốc độ của rotor và có thành phần tích phân để khử sai lệch tĩnh. Từ (2.13),
(2.14) và n2.4, ta cụ thể hóa bằng cấu trúc điều khiển tổng thể phía máy phát như hình 2.5.
Hình 2.5. Hệ thống điều khiển máy phát (DFIG) trong hệ thống PĐSG sử dụng bộ điều
khiển Passivity-Based
2.3.2. Thiết kế bộ điều khiển phía lưới
Do yêu cầu của phía lưới chỉ là điều khiển ổn định điện áp một chiều UDC cung cấp cho mạch
một chiều trung gian. Vì vậy trong nội dung bài báo này cũng chọn phương pháp thiết kế
tuyến tính đơn giản đó là phương pháp tuyến tính Dead – Beat thông thường [2, 3, 7].
2.3.3. Sơ đồ và kết quả mô phỏng bằng Matlab – Simulink – Plecs
Mô phỏng với máy phát có thông số:
THIẾT KẾ BỘ ĐIỀU KHIỂN MÁY PHÁT ĐIỆN KHÔNG ĐỒNG BỘ NGUỒN KÉP 17
Pđm = 4 KW Uđmr = 366 V Rr = 1,32 Ω
Uđms =230/400(∆/Y) nđm =1450 V/p Lσs = 0,066H
fđm = 50 Hz Rs =1,07 Ω Lσs = 0,0098H
zp = 2 Cosϕđm = 0,85 Lm = 0,1601H
J= 0,032Kgm2

  	




   

 
 

 
 
 
 

 
 
  
 !  
"    
# 
"  
"   
" $
 ! 
!!
" 
%
&&
'(

 

 

!  
!  

  
   
" %

" 
  
  )
  *)

  
 *  !

*   !
   

  !
  *  !
 !* !

!*  !
   !

  !

 
+

$
, 
" 
 


+

$

  
 *  !

*   !
   

  !
" 

- + + $

$% 

 
" 

  !
   

*   !

  
 * !
 !*  !

!*   !
   !

  !

  

   
  )
  *)


$&
Hình 2.6. Sơ đồ mô phỏng hệ thống máy phát điện sức gió sử dụng DFIG
(a) (b)
Hình 2.7. Đáp ứng dòng điện ird và irq theo giá trị đặt với bộ điều khiển tựa theo hệ thụ
động EL (a) và với bộ điều khiển tựa theo hệ thụ động EL và Hamiltonian (b).
Kết quả mô phỏng làm việc của hệ thống ở tốc độ định mức (dưới đồng bộ n = 1450v/ph)
như Hình 2.8. Kết quả mô phỏng ở tốc độ dưới đồng bộ (n = 1350 v/ph); công suất phát (mô
men có bước nhảy từ M = -5Nm lên M= -10Nm, công suất vô công Q = 500Var) như Hình
2.9. Kết quả mô phỏng ở tốc độ trên đồng bộ (n = 1650 v/ph); công suất phát (mô men có
bước nhảy từ M = -5Nm lên M = -10Nm, công suất vô công Q = 500Var) như Hình 2.10.
18 ĐẶNG DANH HOẰNG, NGUYỄN PHÙNG QUANG
(a) (b)
Hình 2.8. Đáp ứng dòng ird và irq theo i∗rd và i
∗
rd khi công suất đặt thay đổi.
(a) (b)
Hình 2.9. Kết quả mô phỏng dòng ird và irq theo giá trị đặt (a) và kết quả mô phỏng mômen
theo giá trị đặt (b).
Hình 2.10. Kết quả mô phỏng dòng ird và irq theo giá trị đặt (a) và kết quả mô phỏng
mômen theo giá trị đặt (b).
3. KẾT LUẬN
- Bài báo đã đưa ra được thuật toán điều khiển mới bổ sung cho [4].
- Kết quả mô phỏng cho thấy tính đúng đắn của thuật toán như: Bộ điều khiển có khả năng
THIẾT KẾ BỘ ĐIỀU KHIỂN MÁY PHÁT ĐIỆN KHÔNG ĐỒNG BỘ NGUỒN KÉP 19
khử được sai lệch tĩnh, điều chỉnh được các dòng điện thành phần ird, irq bám theo giá trị
đặt i∗rd, i
∗
rq (Hình 2.7) và đảm bảo chất lượng ở các chế độ làm việc trên và dưới tốc độ đồng
bộ (Hình 2.8, 2.9, 2.10). Như vậy với kết quả mô phỏng cho thấy hệ thống điều khiển đã giải
quyết được các vấn đề nghiên cứu đặt ra.
- Nghiên cứu này mở ra một phương pháp thiết kế phi tuyến mới để điều khiển quá trình làm
việc của máy phát điện không đồng bộ nguồn kép trong hệ thống máy phát điện sức gió, từ
đó là cơ sở để ứng dụng vào thực tiễn.
TÀI LIỆU THAM KHẢO
[1] Arnau D‘oria-Cerezo, “Modeling, simulation and control of a doubly-fed induction machine con-
trolled by a back-to-back converter”, PhD Thesis (2006).
[2] Cao Xuân Tuyển, “Tổng hợp các thuật toán phi tuyến trên cơ sở phương pháp backstepping để
điều khiển máy điện dị bộ nguồn kép trong hệ thống máy phát điện sức gió”, Luận án tiến sĩ kỹ
thuật, Đại học Bách khoa Hà nội (2008).
[3] Đặng Danh Hoằng, Hoà đồng bộ máy phát điện lên lưới bằng phương pháp điều khiển passivity-
based, Tạp chí Khoa học Công nghệ, Đại học Thái nguyên, (2010).
[4] Đặng Danh Hoằng, Nguyễn Phùng Quang, Thiết kế bộ điều khiển dựa trên thụ động “Passivity
- based” để điều khiển máy phát điện không đồng bộ nguồn kép, Tạp chí Khoa học Công nghệ
các trường đại học kỹ thuật (76) (2010).
[5] Levent U. godere, Marwan A. Simaan, Charles W. Brice, Passivity – based control of saturated
Induction Motors, IEEE (1997).
[6] Lan, Ph. Ng, “Linear and nonlinear control approach of doubly - fed induction generator in wind
power generation”, P.h.D thesis, TU-Dresden (2006).
[7] Ng. Ph. Quang, Andreas Dittrich, Truyền động điện thông minh, NXB Khoa học và Kĩ thuật,
Hà Nội, 2002.
[8] Ng. Ph. Quang, Matlab & Simulink dành cho kỹ sư điều khiển tự động, NXB Khoa học và
Kĩ thuật, Hà Nội, 2004.
[9] Ng. Ph. Quang, Điều khiển tự động truyền động điện xoay chiều ba pha, NXB Giáo dục,
Hà Nội, 1996.
[10] N. D. Phước, P. X. Minh, H. T. Trung, Lý thuyết điều khiển phi tuyến, NXB Khoa học và Kĩ
thuật, Hà Nội, 2003.
[11] R. Ortega, A. Loria, P. J. Nicklasson, H. Sira-Ramírez, Passivity-based Control of Euler
Lagrange Systems: Mechanical, Electrical and Electromechanical Applications, Springer-
Verlay, London-Berlin-Heidelberg (1998).
Nhận bài ngày 20 - 12 - 2011