Chng 4 : iu khin m
Hc kì 1 nm hc 2005-2006
Chng 4
IU KHIN M
Khái nim v logic m đc giáo s L.A Zadeh đa ra ln đu tiên nm
1965, ti trng i hc Berkeley, bang California - M. T đó lý thuyt
m đã đc phát trin và ng dng rng rãi.
Nm 1970 ti trng Mary Queen, London – Anh, Ebrahim Mamdani đã
dùng logic m đ điu khin mt máy hi nc mà ông không th điu khin
đc bng k thut c đin. Ti c Hann Zimmermann đã dùng logic m
cho các h ra quyt đnh. Ti Nht logic m đc ng dng vào nhà máy x
lý nc ca Fuji Electronic vào 1983, h thng xe đin ngm ca Hitachi
vào 1987.
Lý thuyt m ra đi M, ng dng đu tiên Anh nhng phát trin mnh
m nht là Nht. Trong lnh vc T đng hoá logic m ngày càng đc
ng dng rng rãi. Nó thc s hu dng vi các đi tng phc tp mà ta
cha bit rõ hàm truyn, logic m có th gii quyt các vn đ mà điu
khin kinh đin không làm đc.
4.1. Khái nim c bn
hiu rõ khái nim “M” là gì ta hãy thc hin phép so sánh sau :
Trong toán hc ph thông ta đã hc khá nhiu v tp hp, ví d nh tp các
s thc R, tp các s nguyên t P={2,3,5, }… Nhng tp hp nh vy đc
gi là tp hp kinh đin hay tp rõ, tính “RÕ” đây đc hiu là vi mt
tp xác đnh S cha n phn t thì ng vi phn t x ta xác đnh đc mt giá
tr y=S(x).
Gi ta xét phát biu thông thng v tc đ mt chic xe môtô : chm,
trung bình, hi nhanh, rt nhanh. Phát biu “CHM” đây không đc ch
rõ là bao nhiêu km/h, nh vy t “CHM” có min giá tr là mt khong
(x) : B
→
[0 1]
trong đó :
μ
F
gi là hàm thuc , B gi là tp nn.
4.1.2. Các thut ng trong logic m
• cao tp m F là giá tr h = Sup
μ
F
(x), trong đó sup
μ
F
(x) ch giá tr nh
nht trong tt c các chn trên ca hàm
μ
F
(x).
• Min xác đnh ca tp m F, ký hiu là S là tp con tho mãn :
S = Supp
μ
F
(x) = { x
∈
B |
μ
F
min tin cy
MX
Chng 4 : iu khin m
Trang 3
4.1.3. Bin ngôn ng
Bin ngôn ng là phn t ch đo trong các h thng dùng logic m. đây
các thành phn ngôn ng ca cùng mt ng cnh đc kt hp li vi nhau.
minh ho v hàm thuc và bin ngôn ng ta xét ví d sau :
Xét tc đ ca mt chic xe môtô, ta có th phát biu xe đang chy:
- Rt chm (VS)
- Chm (S)
- Trung bình (M)
- Nhanh (F)
- Rt nhanh (VF)
Nhng phát biu nh vy gi là bin ngôn ng ca tp m. Gi x là giá tr
ca bin tc đ, ví d x =10km/h, x = 60km/h … Hàm thuc tng ng ca
các bin ngôn ng trên đc ký hiu là :
μ
VS
(x),
μ
S
(x),
μ
M
(x),
μ
(x),
μ
VF
(x) }
Ví d hàm thuc ti giá tr rõ x=65km/h là :
μ
X
(65) = { 0;0;0.75;0.25;0 }
VS S M F VF
0 20 40 60 65 80 100 tc đ
μ 1
0.75
0.25
Hình 4.2:
PGS.TS Nguyn Th Phng Hà
4.1.4. Các phép toán trên tp m
Cho X,Y là hai tp m trên không gian nn B, có các hàm thuc tng ng
là
μ
X
,
μ
Y
Y
(b) =
μ
X
(b) +
μ
Y
(b) -
μ
X
(b).
μ
Y
(b)
- Phép giao hai tp m : X∩Y
+ Theo lut Min
μ
X
∪
Y
(b) = Min{
μ
X
(b) ,
μ
Y
(b) }
+ Theo lut Lukasiewicz
μ
X
4.1.5. Lut hp thành
1. Mnh đ hp thành
Ví d điu khin mc nc trong bn cha, ta quan tâm đn 2 yu t :
+ Mc nc trong bn L = {rt thp, thp, va}
+ Góc m van ng dn G = {đóng, nh, ln}
Ta có th suy din cách thc điu khin nh th này :
Nu mc nc = rt thp Thì góc m van = ln
Nu mc nc = thp Thì góc m van = nh
Nu mc nc = va Thì góc m van = đóng
Trong ví d trên ta thy có cu trúc chung là “Nu A thì B” . Cu trúc này
gi là mnh đ hp thành, A là mnh đ điu kin, C = A
⇒
B là mnh đ kt
lun.
nh lý Mamdani :
“ ph thuc ca kt lun không đc ln hn đ ph thuc điu kin”
Nu h thng có nhiu đu vào và nhiu đu ra thì mnh đ suy din có
dng tng quát nh sau :
If N = n
i
and M = m
i
and … Then R = r
i
and K = k
i
and ….
2. Lut hp thành m
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
),( )1,(
),2( )1,2(
),1( )1,1(
ymxnyxn
ymxyx
ymxyx
RR
RR
RR
μμ
μμ
μμ
=
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
(y) = { l
1
,l
2
,l
3
,…,l
m
} vi l
k
=maxmin{a
i
,r
ik
}.
b. Thut toán xây dng mnh đ hp thành cho h MISO
Lut m cho h MISO có dng :
“If cd
1
= A
1
and cd
2
= A
2
and … Then rs = B”
Các bc xây dng lut hp thành R :
• Ri rc các hàm thuc
μ
i
). T đó suy ra
H = Min{
μ
A1
(c
1
),
μ
A2
(c
2
), …,
μ
An
(c
n
) }
• Lp ma trn R gm các hàm thuc giá tr m đu ra cho tng véct giá tr
m đu vào:
μ
B’
(y) = Min{ H,
μ
B
(y) } hoc
μ
B’
(y) = H.
μ
2
21 yy
+
• Nguyên lý cn trái : chn y’ = y1
• Nguyên lý cn phi : chn y’ = y2
2. Phng pháp trng tâm
im y’ đc xác đnh là hoành đ ca đim trng tâm min đc bao bi
trc hoành và đng
μ
B’
(y).
Công thc xác đnh :
y’ =
∫
∫
S
S
(y)dy
)(
μ
μ
dyyy
trong đó S là min xác đnh ca tp m B’
y1 y2
y
μ H
∑
∫
∑
∫
∑
=
=
=
=
=
=
=
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
=
⎟
⎠
⎞
⎜
⎝
⎛
m
k
k
1
'
1
'
)(
)(
)(
)(
μ
μ
μ
μ
(4.1)
trong đó M
i
=
∫
S
'
)( dyyy
kB
μ
và A
i
=
∫
S
'
)( dyy
kB
y’ =
∑
∑
=
=
m
k
k
m
k
kk
H
Hy
1
1
vi H
k
=
μ
B’k
(y
k
)
ây là công thc gii m theo phng pháp đ cao.
y
m1 m2
a
b
μ
)()( +=
$
. Nu
toàn b các lut ca h thng đc xây dng thì có th mô t toàn b trng
thái ca h trong toàn cc. Trong (4.2) ma trn A(x
k
) và B(x
k
) là nhng ma
trn hng ca h thng trng tâm ca min LX
k
đc xác đnh t các
chng trình nhn dng. T đó rút ra đc :
∑
+= ))()(( uxBxxAwx
kk
k
$
(4.3)
vi w
k
(x)
∈
[0 , 1] là đ tho mãn đã chun hoá ca x* đi vi vùng m LX
k
Lut điu khin tng ng vi (4.2) s là :
R
ck
và đu ra y.
R
1
: If x
1
= BIG and x
2
= MEDIUM Then y
1
= x
1
-3x
2
R
2
: If x
1
= SMALL and x
2
= BIG Then y
2
= 4+2x
1
u vào rõ đo đc là x
1
* = 4 và x
2
* = 60. T hình bên di ta xác đnh
và R
2
là (0.3 ; -176) và (0.35 ; 12). Theo phng
pháp tng trng s trung bình ta có:
77.74
35.03.0
1235.0)176(3.0
−=
+
×
+
−
×
=y
4.2. B điu khin m
4.2.1. Cu trúc mt b điu khin m
Mt b điu khin m gm 3 khâu c bn:
+ Khâu m hoá
+ Thc hin lut hp thành
+ Khâu gii m
Xét b điu khin m MISO sau, vi véct đu vào X =
[
]
T
n
Hình 4.4:
PGS.TS Nguyn Th Phng Hà
4.2.2. Nguyên lý điu khin m ♦ Các bc thit k h thng điu khin m.
+ Giao din đu vào gm các khâu: m hóa và các khâu hiu chnh nh
t l, tích phân, vi phân …
+ Thip b hp thành : s trin khai lut hp thành R
+ Giao din đu ra gm : khâu gii m và các khâu giao din trc tip
vi đi tng.
4.2.3. Thit k b điu khin m
• Các bc thit k:
B1 : nh ngha tt c các bin ngôn ng vào/ra.
B2 : Xác đnh các tp m cho tng bin vào/ra (m hoá).
+ Min giá tr vt lý ca các bin ngôn ng.
+ S lng tp m.
+ Xác đnh hàm thuc.
+ Ri rc hoá tp m.
B3 : Xây dng lut hp thành.
B4 : Chn thit b hp thành.
B5 : Gii m và ti u hoá.
Hình 4.5:
e
μ
B
iv. iu khin Tagaki/Sugeno (TSFC)
4.2.4. Ví d ng dng
Dùng điu khin m đ điu khin h thng bm x nc t đng. H thng
s duy trì đ cao bn nc mt giá tr đt trc nh mô hình bên di.
♦
Mô hình :
Ba b điu khin m (control) s điu khin : bm, van1, van2 sao cho mc
nc 2 bn đt giá tr đt trc (set).
♦
S đ simulink:
PGS.TS Nguyn Th Phng Hà
♦
S đ khi điu khin: Chng 4 : iu khin m
Trang 13
♦
Thit lp h thng điu khin m :
•Xác đnh các ngõ vào/ra :
+ Có 4 ngõ vào gm : sai lch e1, e2; đo hàm sai lch de1, de2
+ Có 3 ngõ ra gm : control1, control2, control3
ZR
OF OS NC CS CF
PM
CS
PB
NC CF CF
DE
Khi controller2
ERROR DB DM ZR IM IB
NB
CF CF NC
NM
CS
ZR
CF CS NC OS OF
PM
OS
PB
NC OF OF
+ Khi “control3”
ây là khi điu tit lu lng cho bn 2, ta đa ra mc u tiên nh sau :
Khi sai lch bn 1 ln thì van2 s điu tit đ sai lch này nh ri mi đn
bn 2.
If error1=NB and de1=DB Then control=CF
If error1=NB and de1=DM Then control=CS
=[0 0]
z (m)
thi gian (s)
PGS.TS Nguyn Th Phng Hà
mc nc đt Z
dat
=[0.5 0.4]
mc nc ban đu Z
init
=[0.8 0] 4.3. Thit k PID m
Có th nói trong lnh vc điu khin, b PID đc xem nh mt gii pháp
đa nng cho các ng dng điu khin Analog cng nh Digital. Vic thit k
b PID kinh đin thng da trên phng pháp Zeigler-Nichols, Offerein,
Reinish … Ngày nay ngi ta thng dùng k thut hiu chnh PID mm
(da trên phm mm), đây chính là c s ca thit k PID m hay PID thích
nghi.
4.3.1. S đ điu khin s dng PID m :
Hình 4.6:
K
D
I
P
++
Các tham s K
P
, K
I
, K
D
đc chnh đnh theo tng b điu khin m riêng
bit da trên sai lch e(t) và đo hàm de(t). Có nhiu phng pháp khác
nhau đ chnh đnh b PID ( xem các phn sau) nh là da trên phim hàm
mc tiêu, chnh đnh trc tip, chnh đnh theo Zhao, Tomizuka và Isaka …
Nguyên tc chung là bt đu vi các tr K
P
, K
I
, K
D
theo Zeigler-Nichols, sau
đó da vào đáp ng và thay đi dn đ tìm ra hng chnh đnh thích hp.
4.3.2. Lut chnh đnh PID:
+ Lân cn a
1
ta cn lut K mnh đ rút ngn thi gian lên, do vy chn: K
P
ln, K
THIT B
CHNH
NH
I T
NG
dt
de
PGS.TS Nguyn Th Phng Hà
+ Lân cn b
1
ta tránh vt l ln nên chn: K
P
nh, K
D
ln, K
I
nh.
+ Lân cn c
1
và d
1
ging nh lân cn a
1
và b
1
h s tích phân
+ K
D
h s vi phân
• S lng bin ngôn ng
ET = {âm nhiu, âm va, âm ít, zero, dng ít, dng va, dng nhiu}
ET = { N3, N2, N1, ZE, P1, P2, P3 }
DET = { âm nhiu, âm va, âm ít, zero, dng ít, dng va, dng nhiu}
DET = { N31, N21, N11, ZE1, P11, P21, P31 }
K
P
/K
D
= { zero, nh, trung bình, ln, rt ln } = {Z, S, M, L,U}
K
I
= {mc 1,mc 2,mc 3,mc 4, mc 5} = {L1,L2,L3,L4,L5}
N3 N2 N1 ZE P1 P2 P3
-12 -8 -4 0 4 8 12
0
C
μ
ET
Chng 4 : iu khin m
Trang 19
DET
K
P
N31 N21 N11 ZE1 P11 P21 P31
N3
U U U U U U U
N2
L L L L L L L
N1
M M M M M M M
ZE
Z Z Z Z Z Z Z
P1
M M M M M M M
P2
L L L L L L L
ET
P3
U U U U U U U
Lut chnh đnh K
D
:
DET
K
N2
L3 L2 L2 L1 L2 L2 L3
N1
L4 L3 L2 L1 L2 L3 L4
ZE
L5 L4 L3 L2 L3 L4 L5
P1
L4 L3 L2 L1 L2 L3 L4
P2
L3 L2 L2 L1 L2 L2 L3
ET
P3
L1 L1 L1 L1 L1 L1 L1
Chng 4 : iu khin m
Trang 21
Biu din lut chnh đnh K
P
trong không gian
3. Chn lut và gii m
+ Chn lut hp thành theo quy tc Max-Min
+ Gii m theo phng pháp trng tâm.
4. Kt qu mô phng
2. H m lai cascade 3. Công tc m
iu khin h thng theo kiu chuyn đi khâu điu khin có tham s đòi
hi thit b điu khin phi cha đng tt c các cu trúc và tham s khác
nhau cho tng trng hp. H thng s t chn khâu điu khin có tham s
phù hp v
i đi tng.
Hình 4.8
B
K
I TNG
B tin
X lý m
Hình 4.9
x
Δ
u
u
+
y
BK M
BK
KINH IN
I T
S dng Simulink kt hp vi toolbox FIS Editor ca Matlab đ mô phng
h thng trên.
áp ng h thng khi không có b m:
Δ
x
+
y
x
DE
E
s
T
K
I
R
1
+
i tng
Δ
u
-
Trong đó x1, x2, x3 ∈ {NB, NS, ZE, PS, PB}
μ
-1 0 1 E
-20 0 20 D
E
NB NS ZE PS PB
-1 0 1
Δ
x
NB NS ZE PS PB
μ
Chng 4 : iu khin m
Trang 25
T hai đ th trên ta thy đc b m đã ci thin rt tt đc tính đng ca
h thng. Th vi nhiu
Δ
u khác nhau ta s thy đáp ng hu nh không
ph thuc vào
Δ
u .
4.5. H m mng nron và ng dng
4.5.1. Mng nron nhân to
Mng nron là s tái to bng k thut nhng chc nng ca h thn kinh
con ngi. Mng nron gm vô s các nron liên kt vi nhau nh hình sau
θ
,
trong đó
θ
là ngng kích hot nron, w
k
là các trng s, f là hàm kích hot.
X
y
B
tng
Hà
m
p
hi tu
y
n
Nhân
Axon
Hình 4.11
Khp ni