Đề tài: Hệ mờ-nơron nhận dạng và điều khiển điều tốc tuốc bin thủy lực

Hội nghị toàn quốc về Điều khiển và Tự động hoá - VCCA-2011 Hệ mờ-nơron nhận dạng và điều khiển điều tốc tuốc bin thuỷ lực Fuzzy-neural System For Identification And Control The Hydro Turbin Governors Nguyễn Tăng Cường, Nguyễn Chỉ Sáng, Phạm Văn Đa e-Mail: dapv.ima@gmail.com Tóm tắt Bài báo này trình bày việc áp dụng hệ mờ - nơron để nhận dạng và điều khiển hệ chuyển động phi tuyến cũng như cho hệ điều tốc tuốc bin thuỷ lực. Kết quả mô phỏng cho khả năng ứng dụng tốt, ổn định và thời gian quá độ nhỏ. Abstract: This paper describes the fuzzy-neural system applying to identification and controls the non-linear motion system as well as the hydro turbine governors. Simulation results show the feasibility of using, stability and the improved of overshoot time. Ký hiệu Ký hiệu Đơn vị Ý nghĩa f hàm phi tuyến μAj (xi ) hàm liên thuộc x Véc tơ trạng thái Chữ viết tắt SISO Single Input – Single Output MISO Multi Input – Single Output MIMO Multi Input – Multi Output NFCs Neuro Fuzzy Controllers T-S Takagi-Sugeno 1. Phần mở đầu Kết hợp mạng nơron vào bộ điều khiển mờ, tạo thành một bộ điều khiển mới gọi là bộ điều khiển mờ - nơron (NFCs) hay còn gọi là hệ mờ-nơron. Sự kết hợp này chính là sự tích hợp chặt chẽ những kiến thức chuyên gia trong hệ mờ và lặp lại nó trong suốt quá trình học của mạng nơron. Hệ mờ-nơron sẽ sử dụng quá trình học của mạng để chỉnh lại các tham số hàm liên thuộc của tập mờ (như tâm, độ rộng) với mong muốn đầu ra của hệ thỏa mãn yêu cầu cho trước [2]. Hệ chuyển động phi tuyến còn có thể được điều khiển bằng hệ mờ thích nghi [1], để đạt được sai lệch tiệm cậm về 0 cho một lớp đối tượng. Kết quả này cũng áp dụng hiệu quả cho nhận dạng và điều khiển hệ điều tốc tuốc bin thuỷ lực, chỉ ra khả năng ứng dụng của phương pháp, sự ổn định làm việc cũng như thời gian đáp ứng tốt. 2. Hệ NFCs với mô hình Takagi-Sugeno (T-S) áp dụng cho nhận dạng hệ chuyển động. Việc áp dụng mô hình mờ T-S đã được nghiên cứu trong nhiều công trình khác nhau, thường sử dụng 1 mạng nơron 5 lớp [2],[4]. Xét một cấu trúc hệ NFCs gồm 5 lớp, 2 lớp ngoài vào và ra thực hiện nhiệm vụ nhận tín hiệu vào và tính giá trị ra của hệ. 03 lớp trong thực hiện chức năng mờ hoá, thực hiện các luật mờ và giải mờ. Để nhận dạng hệ thống, ta dùng hệ mờ-nơron với luật mờ T-S [2],[4],[9]. Điểm khác nhau cơ bản giữa hai mô hình mờ T-S và mô hình mờ Mamdani là mệnh đề kết luận của các luật hợp thành. Mệnh đề kết luận trong mô hình mờ T-S luôn là những hàm giá trị thực, nó thay thế cho các tập mờ. Giả sử luật mờ thứ j trong mô hình mờ T-S có dạng: Nếu 1 là 1 và x2 là 2 và … và xn là Aj Thì y = fj = a0 +aj x +a2 x2 +...+an xn (1) Với: xi là biến đầu vào ; y là biến đầu ra ; Aj là biến ngôn ngữ của mệnh đề điều kiện với hàm liên thuộc μ j (x ) ; aj là hệ số ; j=1…m ; i=1…n. i Xét một hệ chuyển động SISO có phương trình động lực học phi tuyến tổng quát được mô tả bởi hệ phương trình trạng thái (2) như sau: ⎧ dt = F(x)+G(x)u (2) ⎪y = x Trong đó: x = (x ,x2,...,xn ) là một vector các biến trạng thái của hệ; F(x) và G(x) là hai hàm phi tuyến phụ thuộc vào các biến trạng thái x của hệ; u là tín hiệu điều khiển tác động vào hệ; y là tín hiệu ra của hệ. Nếu hệ đã nêu là hệ thống ổn định, tức khi nó bắt đầu ở một vị trí nào đó, nó sẽ tiếp tục làm việc ở lân cận vị trí này, giá trị đầu ra không vượt quá giới hạn cho phép trong suốt thời gian sau đó và quỹ đạo chuyển động trơn, liên tục, thì chúng ta có thể dùng các công cụ toán học để đưa hệ trên về dạng chuẩn (3) . ⎧&1 = x2 x = x ⎨........... (3) xn = F(x)+G(x)u ⎩y = x Cụ thể hoá, ta xét hệ có 04 biến trạng thái và được mô tả bởi phương trình (4) [3],[4],[9], trong đó Hội nghị toàn quốc về Điều khiển và Tự động hoá - VCCA-2011 x = (x ,x2,x3,x4)là các biến trạng thái của hệ; u là đầu vào và y là đầu ra của hệ. F(x) = F(x ,x2,x3,x4 ),G(x) =G(x ,x2,x3,x4)là hai hàm phi tuyến phụ thuộc vào các biến trạng thái của hệ. Hệ mờ-nơron dùng để nhận dạng hai hàm F(x) và G(x) có cấu trúc giống nhau, các trọng số trong mạng được cập nhật là khác nhau và là hệ MISO có 4 đầu vào và một đầu ra. Chọn mỗi đầu vào có hai hàm liên thuộc, thì sẽ có 2n luật hợp thành. Với n =4 là số biến ⎧x1 = x2 ⎪x2 = x3 ⎨x3 = x4 x4 = F(x) + G(x)u ⎪y = x1 đầu vào ta được 16 luật hợp thành từ R1 đến R16 theo luật mờ T-S (1) như các phương trình từ (5) đến (20) (4) sau: R : Nếu x là A và x2 là A1 và x3 là A và x4 là A Thì : y = f1 = a0 +a1x +a2x2 +a3x3 +a4x4 (5) 1 1 1 2 2 1 1 2 2 3 3 4 4 Thì : y = f2 = a0 +a2x +a2 x2 +a3 x3 + a4 x4 (6) R3 : Nếu x là A1 và x2 là A và x3 là A2 và x4 là A Thì : y = f3 = a0 +a3x +a2x2 +a3 x3 + a4 x4 (7) 1 1 2 2 4 1 1 2 2 3 3 4 4 Thì : y = f4 = a0 +a4x +a2 x2 +a3 x3 + a4 x4 (8) R5 : Nếu x là A1 và x2 là A2 và x3 là A và x4 là A Thì : y = f5 = a0 +a5x +a2 x2 +a3 x3 +a4 x4 (9) R6 : Nếu 1 là A1 và x2 là A2 và x3 là 1 và x4 là A2 Thì : y = f6 = a0 +a6x +a2 x2 +a3 x3 +a4 x4 (10) R7 : Nếu x là A và x2 là A2 và x3 là A2 và x4 là A Thì : y = f7 = a0 +a7x +a2 x2 +a3 x3 +a4 x4 (11) 1 2 2 2 8 1 1 2 2 3 3 4 4 Thì : y = f8 = a0 +a8x +a2 x2 +a3 x3 +a4 x4 (12) 2 1 1 1 9 1 1 2 2 3 3 4 4 Thì : y = f9 = a0 +a9x +a2 x2 +a3 x3 +a4 x4 (13) 2 1 1 2 10 1 1 2 2 3 3 4 4 Thì : y = f10 = a00 +a10x +a20x2 +a30x3 + a40x4 (14) 2 1 2 1 11 1 1 2 2 3 3 4 4 Thì : y = f11 = a01 +a11x +a21x2 +a31x3 +a41x4 (15) 12 : Nếu 1 là A2 và x2 là 1 và x3 là A2 và x4 là A2 Thì : y = f12 = a02 +a12x +a22x2 +a32x3 + a42x4 (16) 2 2 1 1 13 1 1 2 2 3 3 4 4 Thì : y = f13 = a03 +a13x +a23x2 +a33x3 +a43x4 (17) 2 2 1 2 14 1 1 2 2 3 3 4 4 2 Thì : y = f14 = a04 +a14x +a24x2 +a34x3 + a44x4 (18) 2 2 2 1 15 1 1 2 2 3 3 4 4 Thì : y = f15 = a05 +a15x +a25x2 +a35x3 +a45x4 (19) R 6 : Nếu x là A2 và x2 là A2 và x3 là A2 và x4 là A2 Thì : y = f16 = a06 +a16x +a26x2 +a36x3 + a46x4 (20) A1 x1 A2 A1 x2 A2 y* x3 A3 A3 A4 x4 A2 H.1 Cấu trúc hệ NFCs theo luật mờ T-S, với n=4, luật hợp thành 2 Trong đó x ; x2;x3;x4 là các biến đầu vào; y là biến đầu ra ; Aij là biến ngôn ngữ của mệnh đề điều kiện với hàm liên thuộc μ j (x ); aj là hệ số. Sơ đồ i cấu trúc mạng NFCs thực hiện các luật mờ từ (5) đến (20) có cấu trúc như hình H.1. + Lớp 1: làm nhiệm vụ nhận tín hiệu vào. + Lớp 2: tất cả các nút trong lớp này thực hiện chức năng mờ hoá, tạo hàm liên thuộc. Đầu ra của chúng chỉ rõ độ thoả mãn của biến đầu vào xi với mỗi biến Ai . Hàm liên thuộc chọn là hàm Gaus có dạng: −⎜ xi −mij ⎟2 j (x ) = e ⎝ i ⎠ (21) i Với mj; σj lần lượt là tâm và độ rộng của hàm liên thuộc. + Lớp 3: thực hiện phép toán trong mệnh đề điều kiện, theo luật tích đại số, để tìm các giá trị hàm liên thuộc tích tương ứng với các trường hợp, từ 1 đến n. Trường hợp cụ thể, n= 4 ta có các phương trình từ 1 đến 16 như sau. μ1 = μA1 (x1 ).μA1 (x2 ).μA1 (x3 ).μA1 (x4 ) μ2 = μA1 (x1 ).μA1 (x2 ).μA1 (x3 ).μA2 (x4 ) μ3 = μA1 (x1 ).μA1 (x2 ).μA2 (x3 ).μA1 (x4 ) μ4 = μA1 (x1 ).μA1 (x2 ).μA2 (x3 ).μA2 (x4 ) μ5 = μA1 (x1 ).μA2 (x2 ).μA3 (x3 ).μA4 (x4 ) μ6 = μA1 (x1).μA2 (x2 ).μA3 (x3 ).μA4 (x4 ) μ7 = μA1 (x1 ).μA2 (x2 ).μA2 (x3 ).μA1 (x4 ) μ8 = μ 1 (x1).μA2 (x2 ).μA3 (x3 ).μA4 (x4 ) μ9 = μ 1 ( 1).μ 2 (x2 ).μ 3 (x3).μ 4 (x4 ) μ10 = μ 1 ( 1).μ 2 (x2 ).μ 3 (x3).μ 4 (x4 ) μ11 = μ 1 ( 1).μ 2 (x2 ).μ 3 (x3).μ 4 (x4 ) μ12 = μ 1 ( 1).μ 2 (x2 ).μ 3 (x3).μ 4 (x4 ) 3 Hội nghị toàn quốc về Điều khiển và Tự động hoá - VCCA-2011 μ13 = μ 1 ( 1).μ 2 (x2).μ 3 (x3).μ 4 (x4) μ14 = μA2 ( 1).μA2 (x2).μ 1 (x3).μA2 (x4) μ15 = μ 1 ( 1).μ 2 (x2).μ 3 (x3).μ 4 (x4) μ16 = μA2 ( 1).μA2 (x2).μA2 (x3).μA2 (x4) + Lớp 4: thực hiện phép toán theo mệnh đề kết luận y* = μ1 f1 +μ2 f2 +...+μn fn (23) 3. Hệ NFCs với phương pháp tuyến tính hoá chính xác cho điều khiển hệ chuyển động phi tuyến Phương pháp tuyến tính hoá chính xác được áp dụng để điều khiển hệ phi tuyến đạt được độ chính xác cao trong việc điều khiển đầu ra bám theo quỹ đạo đặt mong muốn. Trở lại hệ chuyển động được mô tả toán dạng (3), hai hàm F(x);G(x) được nhận dạng dựa trên hệ mờ-nơron sử dụng mô hình mờ T-S như phần 2 (hình 1). Điều kiện để hệ thống điều khiển được là tất cả các biến trạng thái trong hệ bị chặn và G(x) ≠ 0 [2]. Biến đổi (3) ta được: y(n) = F(x)+G(x)u (24) u = G(x)(y(n) −F(x)) (25) Gọi y là ước lượng đầu ra, ta giả thiết y = y và y(n) = y(n) , với sai số được xác định bởi mô hình mờ T-S, bị chặn và tiệm cận về 0 [3]. Khi đó: μ = μ ;μ2 = μ2 ;...;μn = μn ; Sum Sum Sum Với :μSum = ∑μi (22) i=1 + Lớp 5: tính giá trị đầu ra của hệ theo (23). Như vậy quá trình cập nhật trọng số của mạng, đặc biệt là lớp 2, thì các hàm liên thuộc của tập mờ luôn được điều chỉnh để đảm bảo các luật hợp thành được thực hiện tại lớp 3 thoả mãn mô hình mờ T-S và sai lệch của hệ tiệm cận về 0. u = G(x)(y(n) − F(x)) (26) Gọi lượng ra mong muốn là ym. Khi đó sai lệch em = ym − y . Để em→0 sau một khoảng thời gian hữu hạn thì em là nghiệm của phương trình vi phân tuyến tính [2]: (n) (n−1) m 1 m n−1 m n m Với các hệ số α1,α2,...,αn được chọn sao cho tất cả các nghiệm của phương trình đặc trưng sn +α1sn−1 +...+αn−1s+αn = 0có phần thực nằm bên trái trục ảo. Như vậy việc chọn các hệ số αi,i =1→ n có thể thực hiện theo phương pháp gán điểm cực của Luenberger. Thay biểu thức sai lệch vào (27) ta được (28) và (29). Thay (29) vào (26) ta tìm được luật điều khiển (30). Trường hợp xét n=4, ta được cấu trúc nhận dạng và điều khiển hệ theo phương pháp tuyến tính hoá chính xác thích nghi như hình H.2. Hệ chuyển động mô hình ym Bộ điều khiển với luật (30) )(3) && ) u 1 1 1 1 ) s s s s Hệ mờ-nơron F(x) G(x) Hệ mờ-nơron (-) x4 Hệ chuyển động x3 thực tế x2 y x1 H.2 Mô hình cấu trúc hệ nhận dạng và điều khiển với n=4 4 (y(n) − y(n) )+α1(y(n−1) − y(n−1) )+...+αn−1(ym − y)+αn (ym − y) = 0 (28) y(n) = y(n) +α1(y(n−1) − y(n−1) )+...+αn−1(ym − y)+αn (ym − y) (29) u=G(x)(y(n) + 1(y(n−1) −y(n−1))+...+αn−1(ym −y)+αn(ym −y)−F(x)) (30) 4. Mô phỏng hệ NFCs cho nhận dạng và điều khiển điều tốc tuốc bin thuỷ lực Hệ thống điều khiển điều tốc (governor) trong nhà máy thuỷ điện có chức năng rất quan trọng đối với sự hoạt động ổn định về chất lượng và hiệu suất của toàn hệ thống. Trong tài liệu [5], hệ governor được mô tả như tại hình H.3, trong đó hàm truyền của các hệ thống chuyển động thực tế đã được nghiên cứu và đưa ra và được áp dụng trong nhiều công trình nghiên cứu như [2], [7], [9], [11], ta có thể cụ thể hoá các khối chức năng trong hình H.3 như hình H.4 để thuận lợi cho việc mô phỏng áp dụng hệ NFCs cho nhận dạng và điều khiển điều tốc tuốc bin thuỷ lực. Tần số đặt Governor Hệ thống chuyển động thực tế (plant) Van Thuỷ lực, Tuốc bin Servomotor Cơ khí Máy phát Tải Tần số Lưới Lưới H.3 Khối chức năng hệ thống nhà máy thuỷ điện Tần số Đặt ∑ G(p) - Governor 1 1+T p Servomotor Cánh hướng 1−Tw.p 1+0.5Tw.p Tuốc bin 1 1+Tm.p Máy phát ∑ -ΔP 1 Δw Kd+T .p Lưới H.4 Sơ đồ khối hệ thống điều khiển tổ máy thuỷ điện độc lập H.5 Sai lệch tần số của hệ khi tải ngẫu nhiên 95%, 80%, 120% Pm H.6 Sai lệch tần số của hệ khi tải ngẫu nhiên ổn định tại 20%, 95% Pm 5 ... - tailieumienphi.vn