Bài giảng Giải tích hệ thống điện nâng cao - Chương 2: Ma trận tổng trở - Võ Ngọc Điều

ppt 67 trang ngocly 1770
Bạn đang xem 20 trang mẫu của tài liệu "Bài giảng Giải tích hệ thống điện nâng cao - Chương 2: Ma trận tổng trở - Võ Ngọc Điều", để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên

Tài liệu đính kèm:

  • pptbai_giang_giai_tich_he_thong_dien_nang_cao_chuong_2_ma_tran.ppt

Nội dung text: Bài giảng Giải tích hệ thống điện nâng cao - Chương 2: Ma trận tổng trở - Võ Ngọc Điều

  1. GIẢI TÍCH HỆ THỐNG ĐIỆN NÂNG CAO CHƯƠNG 2: MA TRẬN TỔNG TRỞ Võ Ngọc Điều Bộ Môn Hệ Thống Điện Khoa Điện – Điện tử Trường ĐH Bách Khoa 1
  2. Ma tRận Tổng Trở Nút ❖ Định nghĩa: ❖ Thành lập ma trận trực tiếp - Nghịch đảo ma trận tổng dẫn nút là công việc đòi hỏi nỗ lực n3. - Đối với các mạng có kích cỡ nhỏ hoặc trung bình, việc thành lập ma trận trực tiếp đỡ tốn công sức hơn. - Đối với các mạng lớn, lập trình ma trận thưa với kỹ thuật khử Gauss được ưu tiên lựa chọn. 2
  3. Thành Lập Ma Trận Tổng Trở Nút ❖ Các kỹ thuật lý thuyết graph được dùng để giải thích quá trình xây dựng. 3
  4. Thành Lập Ma Trận Tổng Trở Nút ❖ Xây dựng cơ sở cho mạng và ma trận 4
  5. Thêm Một Đường Dây Mới ❖ Trường hợp có thêm một nút mới, có 2 trường hợp: Thêm nút mới nối nút có sẵn Thêm nút mới nối nút Ref. 5
  6. Thêm Một Đường Dây Mới 1) Thêm một đường dây mới nối với nút có sẵn p hay nối một nút mới q với một nút có sẵn p. p Mạng nguyên thủy I p p Iq q Z pq 0 Reference Dòng điện Iq chảy vào mạng ở nút p sẽ làm tăng điện áp nguyên 0 thủy Vp bởi điện áp: Iq*Zpp 6
  7. Thêm Một Đường Dây Mới * Điện thế tại nút p cũ: Vp = Z p1I1 + Z p2 I2 + + Z pp (I p + Iq ) + + Z pmIm 0 Vp =V p + IqZ pp 0 Vp sẽ lớn hơn giá trị Vp bởi điện áp IqZpp * Điện thế tại nút q mới 0 Vq =Vp + Iq zpq =Vp + IqZ pp + Iq zpq Vq = Z p1I1 + Z p2 I2 + + Z pmIm + Iq (Z pp + z pq ) 0 V p * Giá trị mới bây giờ phải thêm vào Zorig là hàng p, cột p và phần tử trên đường chéo là Zpp + zpq. 7
  8. Thêm Một Đường Dây Mới * Ma trận mới được thành lập 8
  9. Thêm Một Đường Dây Mới * Quy tắc thực hiện: - Giả sử mạng có sẵn có ma trận Znut với kích cỡ mxm, trong đó m là số nút. - Khi có một nút mới q nối vào nút cũ p thì ma trận Znut có kích cỡ là (m+1)x(m+1). - Các phần tử ở vị trí mới như sau: + Cột mới m+1 với các hàng 1 m = với cột p trong ma trận cũ, trong đó p là nút mà nhánh mới nối vào. + Hàng mới m+1 với các cột 1 m = với hàng p trong ma trận cũ, với p là nút mà nhánh mới nối vào. + Phần tử mới trên đường chéo (m+1,m+1) = phần tử thứ p trên đường chéo của ma trận cũ (Zpp) + tổng trở nối từ điểm mới tới điểm p (Zpq), với q là điểm mới. 9
  10. Thêm Một Đường Dây Mới 2) Thêm một đường dây mới nối với nút reference hay nối một nút mới q với nút reference 10
  11. Thêm Một Đường Dây Mới * Quy tắc thực hiện: - Giả sử mạng có sẵn có ma trận Znut với kích cỡ mxm, trong đó m là số nút. - Khi có một nút mới q nối vào nút reference thì ma trận Znut có kích cỡ là (m+1)x(m+1). - Các phần tử ở vị trí mới như sau: + Cột mới thứ m+1 với các hàng 1 m = 0 (= cột của nút reference). + Hàng mới thứ m+1 với các cột 1 m = 0 (= hàng của nút reference). + Phần tử mới trên đường chéo (m+1,m+1) = tổng trở nối từ điểm mới tới điểm reference (Z0q), với q là điểm mới. 11
  12. Thêm Một Đường Dây Mới ❖ Trường hợp không có thêm nút mới, có 2 trường hợp: Nhánh mới nối 2 nút có sẵn Nhánh mới nối nút có sẵn và Ref. 12
  13. Thêm Một Đường Dây Mới 1) Thêm một đường dây mới nối với 2 nút có sẵn p I p − I qp I p Original network with bus p q and the Iqp z I pq reference node q extracted Iq + Iqp q 0 Reference * Trường hợp này dòng điện bơm vào mạng từ nút q là Iq+Iqp thay vì Iq. Tương tụ, dòng điện bơm vào mạng từ nút p là Ip-Iqp thay vì Ip. 13
  14. Thêm Một Đường Dây Mới * Điện áp tại nút thứ i trong hệ thống Vi = Zi1I1 + Zi2I2 + + Zip (I p − Iqp ) + Ziq (Iq + Iqp ) + + ZimIm = Zi1I1 + Zi2I2 + + ZipI p + ZiqIq + + ZimIm + Iqp (Ziq − Zip) * Tương tự, điện áp tại nút thứ p và q trong hệ thống Vp = Z p1I1 + Z p2 I2 + + Z ppI p + Z pqIq + + Z pmIm + Iqp (Z pq − Z pp ) Vq = Zq1I1 + Zq2 I2 + + ZqpI p + ZqqIq + + ZqmIm + Iqp (Zqq − Zqp ) Vp −Vq = z pqIqp == 0 = z pqIqp −Vp +Vq * Khử dòng điện Iqp 0 = (Zq1 − Z p1)I1 + + (Zqp − Z pp )I p + (Zqq − Z pq )Iq + (Z pp + Zqq − 2Z pq + z pq )Iqp 14
  15. Thêm Một Đường Dây Mới Ma trận sau khi thành lập với nhánh mới thêm vào: 15
  16. Thêm Một Đường Dây Mới * Quy tắc thực hiện: - Giả sử mạng có sẵn có ma trận Znut với kích cỡ mxm, trong đó m là số nút. - Khi có một nhánh mới nối giữa 2 nút có sẵn thì ma trận Znut tạm thời có kích cỡ là (m+1)x(m+1). - Các phần tử ở vị trí mới như sau: + Cột mới thứ m+1 với các hàng 1 m = tổng trở cột q (Ziq, i = 1 m) – tổng trở cột p (Zip, i = 1 m). + Hàng mới thứ m+1 với các cột 1 m = tổng trở hàng q (Zqj, j = 1 m) – tổng trở hàng p (Zpj, j = 1 m). + Phần tử mới trên đường chéo (m+1,m+1) = Zll, trong đó Zll = zpq (mới thêm) + Zpp + Zqq - 2Zpq. - Dùng phép khử Kron cho hàng m+1 và cột m+1 16
  17. Thêm Một Đường Dây Mới 2) Thêm một đường dây mới nối giữa 1 nút có sẵn và nút reference (chuẩn) 17
  18. Thêm Một Đường Dây Mới * Quy tắc thực hiện: - Giả sử mạng có sẵn có ma trận Znut với kích cỡ mxm, trong đó m là số nút. - Khi có một nhánh mới nối giữa 1 nút có sẵn p và nút ref. thì ma trận Znut tạm thời có kích cỡ là (m+1)x(m+1). - Các phần tử ở vị trí mới như sau: + Cột mới thứ m+1 với các hàng 1 m = - tổng trở cột p (Zip, i = 1 m). + Hàng mới thứ m+1 với các cột 1 m = - tổng trở hàng p (Zpj, j = 1 m). + Phần tử mới trên đường chéo (m+1,m+1) = Zll, trong đó Zll = zp0 (mới thêm) + Zpp (từ ma trận). - Dùng phép khử Kron cho hàng m+1 và cột m+1 18
  19. Nhắc Lại Phép Khử Kron * Phép khử Kron dùng để khử bỏ một hệ trục (hàng & cột) trong một ma trận trong khi vẫn còn duy trì tác dụng số của hệ trục đó. 19
  20. Các Quy Tắc Xây Dựng Zbus ❖ Quy tắc 1: Thêm một nhánh mới nối giữa một nút mới p và nút reference (chuẩn) - Bắt đầu với ma trận mạng có sẵn (mxm), m là số nút của mạng có sẵn. - Tạo ra một ma trận mới có kích cỡ (m+1)x(m+1) với + Các hàng và cột ngoài đường chéo mới được lắp đầy bằng số 0. + Phần tử trên đường chéo (m+1)x(m+1) được lắp vào tổng trở thành phần zp0. 20
  21. Các Quy Tắc Xây Dựng Zbus ❖ Quy tắc 2: Thêm một nhánh mới nối giữa nút cũ p và nút mới q. - Bắt đầu với ma trận mạng có sẵn (mxm), m là số nút của mạng có sẵn. - Tạo ra một ma trận mới có kích cỡ (m+1)x(m+1) với + Các hàng và cột ngoài đường chéo mới được lắp đầy bằng cách sao chép lại hàng p (với hàng mới) và cột p (với cột mới). + Phần tử trên đường chéo (m+1)x(m+1) được lắp vào tổng trở thành phần zpq (mới thêm) + tổng trở trên đường chéo Zpp. 21
  22. Các Quy Tắc Xây Dựng Zbus ❖ Quy tắc 3: Thêm một nhánh mới nối giữa nút cũ p và nút reference. - Bắt đầu với ma trận mạng có sẵn (mxm), m là số nút của mạng có sẵn. - Tạo ra một ma trận mới có kích cỡ (m+1)x(m+1) với + Các hàng và cột ngoài đường chéo mới được lắp đầy bằng cách sao chép lại trừ hàng p (với hàng mới) và trừ cột p (với cột mới). + Phần tử trên đường chéo (m+1)x(m+1) được lắp vào tổng trở zp0 + Zpp. - Thực hiện phép thử Kron cho hàng và cột thứ m+1 22
  23. Các Quy Tắc Xây Dựng Zbus ❖ Quy tắc 4: Thêm một nhánh mới nối giữa 2 nút cũ p và q. - Bắt đầu với ma trận mạng có sẵn (mxm), m là số nút của mạng có sẵn. - Tạo ra một ma trận mới có kích cỡ (m+1)x(m+1) với + Các hàng và cột ngoài đường chéo mới được lắp đầy bằng cách sao chép lại hàng q trừ hàng p (với hàng mới) và cột q trừ cột p (với cột mới). + Phần tử trên đường chéo (m+1)x(m+1) được lắp vào tổng trở zpq + Zpp + Zqq – 2Zpq. - Thực hiện phép thử Kron cho hàng và cột thứ m+1 23
  24. Các Bước Xây Dựng Zbus 1) Đầu tiên bắt đẩu từ phương trình nối 1 nút mới đầu tiên với nút chuẩn quan tổng trở Za V1= Za  I1 2) Sau đó thêm một nút mới nối với nút chuẩn hoặc nút mới đầu tiên; giả sử nút thứ hai nối với nút chuẩn qua Zb V1 Za 0 I1 = V2 0 Zb I2 3) Thêm những nút khác vào giống như các trường hợp đã nêu cho đến nhánh cuối cùng. 24
  25. Ví dụ ❖ Ví dụ 1: Thực hiện xây dựng Zbus (hay Znut) trực tiếp cho mạng sau: ❖ Thứ tự các đường dây thêm vào: 1-0, 2-0, 1-3, 1-3 và 2- 3 25
  26. Ví dụ 26
  27. Ví dụ 27
  28. Ví dụ ❖ Ví dụ 2: Thực hiện xây dựng Zbus hay Znut trực tiếp cho mạng sau: Z 4 bus j0.2 (5) j0.125 (6) j0.25 (2) j0.4 (3) 1 2 3 j1.25 (1) (4) j1.25 0 28
  29. Ví dụ 1 Trường hợp 1 1 1 V1  =  j1.25 I1  1 j1.25 1 Zbus,1 = j1.25 Thành lập bus 2 với tổng trở của nó j0.25 nối vào bus 1, có Trường hợp 2: 2 1 j0.25 1 j1.25 j1.25 1 2 Zbus,2 = 2 j1.25 j1.5 Số j1.50 ở trên chính là tổng của j1.25 và j0.25 . 29
  30. Ví dụ Kế tiếp, nối nút mới 3 vào nút 2 qua trở kháng j0.4: có Trường hợp 2 2 3 1 j0.25 j0.4 1 j1.25 j1.25 j1.25 1 2 3 Z =2 j1.25 j1.50 j1.50 bus,3 j1.25 3 j1.25 j1.50 j1.90 Số j1.90 ở trên chính là tổng của Z22 (j1.50) của ma trận cũ và Zpq (j0.4) của nhánh mới nối 2, nút 2 và 3. 30
  31. Ví dụ Thêm tổng trở j1.25 nối từ nút 3 vào nút chuẩn, có Trường hợp 3: 1 2 3 p 1 j1.25 j1.25 j1.25 j1.25 2 j1.25 j1.50 j1.50 j1.50 Zbus,4 = 3 j1.25 j1.50 j1.90 j1.90 p j1.25 j1.50 j1.90 j3.15 j0.4 Trong đó j3,15 chính là tổng của j0.25 1 3 2 Z33 + Z pq j1.25 31
  32. Ví dụ Sau khi loại bỏ hàng p cột p: 1 2 3 1 j0.75397 j0.65476 j0.49603 Z = 2 j0.65476 j0.78571 j0.59524 bus,5 3 j0.49603 j0.59524 j0.75397 ( j1.25)( j1.25) Z = j1.25 − = j0.75397 11(new) j3.15 ( j1.50)( j1.50) Z = j1.50 − = j0.78571 22(new) j3.15 ( j1.50)( j1.90) Z = Z = j1.50 − = j0.59524 23(new) 32(new) j3.15 32
  33. Ví dụ Thêm tổng trở Zpq = j0.20 giữa nút cũ 3 và nút mới 4, có Trường hợp 2: 1 2 3 4 1 j0.75397 j0.65476 j0.49603 j0.49603 2 j0.65476 j0.78571 j0.59524 j0.59524 Z = bus,6 j0.49603 j0.59524 j0.75397 j0.75397 3 4 j0.49603 j0.59524 j0.75397 j0.95397 Phần tử nằm trên đường chéo chính 4 j0.2 Z44 là tổng của Z33 và Zpq j0.25 1 3 2 j1.25 33
  34. Ví dụ Cuối cùng thêm Zpq = j0.125 vào giữa 2 nút cũ là 2 và 4: Có Trường hợp 4 Z15 = Z12 − Z14 = j0.65476 − j0.49603 = j0.15873 Z25 = Z22 − Z24 = j0.78571− j0.59524 = j0.19047 Z35 = Z32 − Z34 = j0.59524 − j0.75397 = − j0.15873 Z45 = Z42 − Z44 = j0.59524 − j0.95397 = − j0.35873 Z55 = Z22 + Z44 −2Z24 + Zb = j{(0.78571+ 0.95397) − 2(0.59524)} + j0.125 = j0.67421 4 j0.2 j0.125 1 j0.25 j0.4 3 2 j1.25 34
  35. Ví dụ Ma trận mới: q j0.15873 Zbus,6 j0.19047 − j0.15873 − j0.35873 q j0.15873 j0.19147 − j0.15873 − j0.35873 j0.67421 35
  36. Ví dụ Sau khi loại bỏ hàng q cột q: 1 2 3 4 1 j0.71660 j0.60992 j0.53340 j0.58044 2 j0.60992 j0.73109 j0.64008 j0.69659 Zbus = 3 j0.53340 j0.64008 j0.71660 j0.66951 4 j0.58049 j0.69659 j0.66951 j0.76310 36
  37. Bài Tập Tự Làm ❖ Viết thuật toán xây dựng ma trận Zbus. ❖ Thành lập ma trận Zbus có kể đến phần tử hỗ cảm 37
  38. Tổng Trở Thevenin và Ma Trận Zbus ❖ Để thành lập mối quan hệ giữa các phần tử trong ma trận Zbus và tổng trở tương đương Thevenin, hãy xem xét một ví dụ sau. ❖ Ví dụ: Xem xét một hệ thống có 2 nút. Có thể thấy rằng các điện áp hở mạch của các nút a và b lần lượt là Va và Vb. 38
  39. Tổng Trở Thevenin và Ma Trận Zbus ❖ Ma trận tổng dẫn được thành lập như sau: ❖ Định thức ma trận trên 39
  40. Tổng Trở Thevenin và Ma Trận Zbus ❖ Ma trận Zbus được thành lập dựa trên nghịch đảo ma trận tổng dẫn nút: ❖ Sau khi thực hiện biến đổi 40
  41. Tổng Trở Thevenin và Ma Trận Zbus ❖ Bây giờ xem xét lại hình vẽ, ma trận tổng trở Thevenin nhìn vào mạch tại điểm a là tổ hợp song song của Zaa và Zab + Zbb, tức là: ❖ Tương tự, tổng trở Thevenin có được bằng cách nhìn vào hệ thống tại nút b là tổ hợp song song của Zbb, và Zaa+Zab, tức là: 41
  42. Tổng Trở Thevenin và Ma Trận Zbus ❖ Vì thế, các tổng trở ở điểm quan sát của 2 nút là tổng trở Thevenin của chúng. ❖ Bây giờ xem xét tổng trở Thevenin khi nhìn vao hệ thống giữa 2 nút a và b. Từ hình vẽ, rõ ràng tổng trở Thevenin là tổ hợp song song của Zab và Zaa+Zbb, tức là: 42
  43. Tổng Trở Thevenin và Ma Trận Zbus ❖ Từ giá trị của Zbus, viết lại: ❖ So sánh hai phương trình cuối: 43
  44. Tổng Trở Thevenin và Ma Trận Zbus ❖ Như đã thấy trong ví dụ trong mối quan hệ V = Zbus*I, các điện áp nút Vi, i = 1, n là các điện áp hở mạch. ❖ Giả sử dòng điện bơm vào tại các nút 1, , k-1 và k+1, , n là 0 khi khi xảy ra ngắn mạch ở nút k. Tổng trở Thevenin tại nút k: ❖ Từ các công thức của Zth,a, Zth,b và Zth,ab cho thấy rằng tổng trở điểm quan sát ở mỗi nút chính là tổng trở Thevenin. 44
  45. Tổng Trở Thevenin và Ma Trận Zbus ❖ Bây giờ tìm tổng trở Thevenin giữa 2 nút j và k của một hệ thống. Ký hiệu các điện áp hở mạch là V0 và các dòng điện tương ứng được ký hiệu là I0 sao cho: ❖ Giả sử dòng điện thay đổi một lượng I sao cho các điện áp thay đổi một lượng V. Thì: ❖ So sánh hai phương trình trên, ta có: 45
  46. Tổng Trở Thevenin và Ma Trận Zbus ❖ Bây giờ giả sử rằng các dòng điện thêm vào Ik và Ij được bơm vào lần lượt các nút k và j trong khi dòng điện bơm vào các nút khác vẫn không đổi. ❖ Phương trình của V được viết lại: 46
  47. Tổng Trở Thevenin và Ma Trận Zbus ❖ Do đó, có thể viết được 2 phương trình như sau: 47
  48. Tổng Trở Thevenin và Ma Trận Zbus ❖ Vì Zjk = Zkj, mạng có thể được vẽ lại: 48
  49. Tổng Trở Thevenin và Ma Trận Zbus ❖ Quan sát thấy rằng điện áp hở mạch giữa các nút k và j: ❖ Dòng điện ngắn mạch qua hai nút này: ❖ Trong thời gian ngắn mạch, Vk – Vj = 0. Thay thế các giá trị của Vj và Vk: ❖ Kết quả này giống như đã dẫn ra trong ví dụ. 49
  50. Tính Các Phần Tử Zbus Từ Ybus ❖ Khi toàn bộ ma trận Zbus không cần tính cho tính toán thì các phần tử của Zbus được tính riêng lẻ nếu ma trận Ybus được thừa số hóa. ❖ Xem xét việc trích ra một hàng m từ ma trận Zbus bằng cách nhân với một vector có phần tử hàng m bằng 1 còn các phần tử khác là 0: 50
  51. Tính Các Phần Tử Zbus Từ Ybus ❖ Ma trận trận trích ra thực chất là một vector ❖ Do tích của Ybus và Zbus là ma trận đơn vị nên: 51
  52. Tính Các Phần Tử Zbus Từ Ybus ❖ Nếu ma trận tam giác dưới L và ma trận tam giác trên U của Ybus có sẵn thì: (m) ❖ Từ đây rõ ràng rằng các phần tử của Zbus có thể xác định được bằng cách giải phương trình tuyến tính trên sử dụng các phép khử thuận và thay thế nghịch. 52
  53. Tính Các Phần Tử Zbus Từ Ybus (m) ❖ Giả sử chỉ có vài phần tử trong Zbus cần được tính, giả sử cần tính z33 và z43 trong hệ thống có 4 nút, cách thực hiện như sau. - Sử dụng kết qua phân tích thừa số hóa: 53
  54. Tính Các Phần Tử Zbus Từ Ybus (3) - Thực hiện giải Zbus theo 2 bước: trong đó: 54
  55. Tính Các Phần Tử Zbus Từ Ybus - Thay thế thuận, cho kết quả: - Và thay thế nghịch để có kết quả cột 3 của Zbus: ❖ Nếu tính tất cả các phần tử: 55
  56. Tính Các Phần Tử Zbus Từ Ybus ❖ Ví dụ tự đọc: - Ví dụ 8.6 trang 308 sách của Stevenson. 56
  57. Tính Các Phần Tử Zbus Từ Ybus ❖ Ví dụ tự đọc: - Ví dụ 8.6 trang 308 sách của Stevenson. 57
  58. Tính Các Phần Tử Zbus Từ Ybus ❖ Biến đổi này nhằm bảm toàn công suất, tức là vẫn duy trì được mối quan hệ trong ma trận tổng trở nút. ❖ Công suất phức: ❖ Giả sử cần chuyển đổi dòng điện nút I tới dòng điện mới Inew sử dụng ma trận biến đổi C: 58
  59. Tính Các Phần Tử Zbus Từ Ybus ❖ Việc chuyển đổi này được thực hiện trong trương hợp khi muốn chuyển nút refrence sang một nút khác và yêu cầu tính lại Zbus(new). ❖ Mối quan hệ điện áp và dòng điện trong trường hợp cũ và mới: ❖ Phải tìm điều kiện để thỏa mãn để Vnew và Zbus(new) sao cho công suất vẫn không biến đỗi khi dòng điện thay đổi. ❖ Công suất phức với thay thế V: 59
  60. Tính Các Phần Tử Zbus Từ Ybus ❖ Thay thế I = CInew vào công suất: ❖ Sau khi biến đổi: ❖ Rút ra mối quan hệ giữa Zbus và Zbus(new) 60
  61. Tính Các Phần Tử Zbus Từ Ybus ❖ Thay thế các giá trị mới vào công suất phức: ❖ Từ ban đầu, có: ❖ Mối quan hệ giữa V và Vnew 61
  62. Tính Các Phần Tử Zbus Từ Ybus ❖ Tổ thất công suất được phân ra như sau: ❖ Thực hiện SL + SL*, kết quả: ❖ Khi Zbus đối xứng, có thể viết: ❖ Tổn thất công suất thực 62
  63. Tính Các Phần Tử Zbus Từ Ybus ❖ Xem xét một ví dụ hệ thống có 4 nút, chuyển nút refrence từ nút n sang nút 4. ❖ Quan hệ dòng và áp: ❖ Theo quan hệ Kirchhoff: 63
  64. Tính Các Phần Tử Zbus Từ Ybus ❖ Viết theo quan hệ mới: ❖ Tính Zbus(new) theo ma trận chuyển đổi C: 64
  65. Tính Các Phần Tử Zbus Từ Ybus ❖ Nhân ma trận theo 2 bước: - Bước 1: - Bước 2: 65
  66. Tính Các Phần Tử Zbus Từ Ybus ❖ Xác định được Zbus(new) ❖ Điện áp trong hệ mới: 66
  67. Tính Các Phần Tử Zbus Từ Ybus ❖ Thuật toán: Khi một nút k trong Zbus hiện hữu được chọn làm reference mới, thuậ toán như sau: 67