Regulator rozmyty
Ireneusz Dominik KAP AGH
Rys. 1. Regulator rozmyty
a)
fuzzyfikacja
Rys.2. Rozkład funkcji przynależności dla1-go wejścia SV
ZawWE:
2,5
– 2 =
0,5
SzerPrzedz: 4-2=
2
Przyn_SV[0] =2,5mm in FP niska: 1-(
0,5
/
2
)=3/4
Przyn_SV[1] =2,5mm in FP srednia: 1-3/4=1/4
Rys.2. Rozkład funkcji przynależności dla 2-go wejścia e
b) wyliczanie wartości wyjściowej dla danych: SV =
2,5
[mm], PV= 1,8 [mm], obliczone e =
0,7
[mm],
output MV=?
Fig. 4. Baza reguł
MV = ¾ *1/7*50 + ¾*6/7*60 + ¼*1/7* 80 + ¼*6/7*85 =
3/28*50 + 18/28*60 + 1/28*80 + 6/28*85 = 65 [A]
-0
.5
(A
) D
uz
y_
m
-
- 0
.3
(B
) M
aly
_m
-0
.2÷
0.1
(C
) z
er
o
0.8
(D
) M
aly
_p
(1) 0
!"#$% (2) 2
# &!"% (3) 4
(
10
15
20
25
30
40
50
60
70
75
80
85
+
u (e)
1
- 0.5
Duzy_m
(A)
- 0.3 - 0.2
0.1
0.8
0
0
Maly_m
(B)
zero
(C)
Maly_p
(D)
[mm]
0,7
6/7
1/7
0
2 4
niska (2)
u (SV)
1
zero (1)
[mm]
2,5
srednia (3)
3/4
1/4
error e
se
t
v
a
lu
e
S
V
(* Wprowadzanie parametrow bazy reguł *)
(*
X := 0;*)
Sing[0] := 10;
Sing[1] := 30;
Sing[2] := 70;
(*
X := 1;*)
Sing[3] := 15;
Sing[4] := 40;
Sing[5] := 75;
(*
X := 2;*)
Sing[6] := 20;
Sing[7] := 50;
Sing[8] := 80;
(*
X := 3;*)
Sing[9] := 25;
Sing[10] := 60;
Sing[11] := 85;
(* Położenie granic FP dla wybranych wejść*)
(* Granice FP dla wejścia 1 - dla uchybu*)
We[0] := e_Duzy_m; (* A -0,5*)
We[1] := e_Maly_m; (* B -0,3*)
We[2] := e_Zero_m; (* C -0,2*)
We[3] := e_Zero_p; (* C 0,1*)
We[4] := e_Maly_p; (* D 0,8*)
(* Granice FP dla wejścia 2 - dla wartości zadanej*)
We[5] := INT_TO_REAL(SV_Zero); (* 1 +0*)
We[6] := INT_TO_REAL(SV_Niska); (* 2 +2*)
We[7] := INT_TO_REAL(SV_Srednia); (* 3 +4*)
(* Rozmywanie e *)
if e > We[0] then
if e < We[1] then
ZawWE := e - We[0];
SzerPrzedz := We[1] - We[0];
Przyn_e[1] := ZawWE / SzerPrzedz;
Przyn_e[0] := 1,0 - Przyn_e[1];
X := 0;
elsif e < We[2] then
ZawWE := e - We[1];
SzerPrzedz := We[2] - We[1];
Przyn_e[1] := ZawWE / SzerPrzedz;
Przyn_e[0] := 1,0 - Przyn_e[1];
X := 1;
elsif e <= We[3] then
Przyn_e[0] := 1,0;
Przyn_e[1] := 0,0;
X := 2;
elsif e < We[4] then
ZawWE := e - We[3];
SzerPrzedz := We[4] - We[3];
Przyn_e[1] := ZawWE / SzerPrzedz;
Przyn_e[0] := 1,0 - Przyn_e[1];
X := 2;
else
Przyn_e[0] := 1,0;
Przyn_e[1] := 0,0;
X := 3;
end_if;
else
Przyn_e[0] := 1,0;
Przyn_e[1] := 0,0;
X := 0;
end_if;
(* Rozmywanie SV *)
if SV > We[5] then
if SV < We[6] then
ZawWE := SV - We[5];
SzerPrzedz := We[6] - We[5];
Przyn_SV[1] := ZawWE / SzerPrzedz;
Przyn_SV[0] := 1,0 - Przyn_SV[1];
Y := 0;
elsif SV < We[7] then
ZawWE := SV - We[6];
SzerPrzedz := We[7] - We[6];
Przyn_SV[1] := ZawWE / SzerPrzedz;
Przyn_SV[0] := 1,0 - Przyn_SV[1];
Y := 1;
else
Przyn_SV[1] := 0,0;
Przyn_SV[0] := 1,0;
Y := 2;
end_if;
else
Przyn_SV[0] := 1,0;
Przyn_SV[1] := 0,0;
Y := 0;
end_if;
(* Wnioskowanie i wyostrzanie (wyznaczanie wartości wyjściowej *)
out := Przyn_e[0] * Przyn_SV[0] * INT_TO_REAL(Sing[Y + 3 * X])
+ Przyn_e[0] * Przyn_SV[1] * INT_TO_REAL(Sing[Y + 1 + 3 * X])
+ Przyn_e[1] * Przyn_SV[0] * INT_TO_REAL(Sing[Y + 3* (X + 1)])
+ Przyn_e[1] * Przyn_SV[1] * INT_TO_REAL(Sing[Y + 1 + 3 * (X + 1)]);