Trang 33 - Số liệu ban đầu
basemva = 10; accuracy = 0.00001; maxiter = 500;
% No Loai U goc Phu tai May phat Qbu Udm % P(MW) Q(MVAr) P(MW) Q(MVAr) Qmin Qmax kV busdata=[1 1 1.0 0.0 0.0 0.0 0.0 0.0 0 0 0 22 2 0 1.0 0.0 2.0 1.5 0.0 0.0 0 0 0 22 3 0 1.0 0.0 2.0 1.5 0.0 0.0 0 0 0 22 4 0 1.0 0.0 2.0 1.5 0.0 0.0 0 0 0 22 5 0 1.0 0.0 2.0 1.5 0.0 0.0 0 0 0 22 6 0 1.0 0.0 2.0 1.5 0.0 0.0 0 0 0 22 7 0 1.0 0.0 2.0 1.5 0.0 0.0 0 0 0 22 8 0 1.0 0.0 2.0 1.5 0.0 0.0 0 0 0 22 9 0 1.0 0.0 2.0 1.5 0.0 0.0 0 0 0 22 10 0 1.0 0.0 2.0 1.5 0.0 0.0 0 0 0 22 11 0 1.0 0.0 1.0 0.75 0.0 0.0 0 0 0 22 12 0 1.0 0.0 0.5 0.375 0.0 0.0 0 0 0 22 13 0 1.0 0.0 0.5 0.375 0.0 0.0 0 0 0 22 14 0 1.0 0.0 0.5 0.375 0.0 0.0 0 0 0 22 15 0 1.0 0.0 0.5 0.375 0.0 0.0 0 0 0 22 16 0 1.0 0.0 0.5 0.375 0.0 0.0 0 0 0 22 ];
% dau cuoi R(dvtd) X(dvtd) Y/2(dvtd) tap linedata=[1 2 0.005372 0.01446 0.0 1 2 3 0.005372 0.01446 0.0 1 2 4 0.005372 0.01446 0.0 1 2 5 0.005372 0.01446 0.0 1 4 6 0.005372 0.01446 0.0 1 6 7 0.005372 0.01446 0.0 1 6 8 0.005372 0.01446 0.0 1 6 9 0.005372 0.01446 0.0 1 8 10 0.005372 0.01446 0.0 1 10 11 0.005372 0.01446 0.0 1 11 12 0.005372 0.01446 0.0 1 12 13 0.005372 0.01446 0.0 1 13 14 0.005372 0.01446 0.0 1 14 15 0.005372 0.01446 0.0 1 15 16 0.005372 0.01446 0.0 1 ];
Trang 34
- Khai báo số tcsc : so_tcsc, số nút được điều chỉnh : sonut_dc
- Thành lập ma trận cột vpv2 kích thước sonutpvx1 chứa bình phương điện áp nút vv2 của các nút PV , (bustype =1)
-Thành lập ma trận cột v kích thước (sonut-sonutpv)x1 chứa bình phương điện áp nút vv2 của các nút PQ, (bustype=0 )
- Thành lập ma trận cột ql kích thước (sonut-sonutpv)x1 chúa qload của các nút PQ (bustype=0)
- Các ma trận Ac , A1 , A1cong , A1tru , A1phay vẫn không thay đổi -Tính các ma trận :
lamda_Act_vpv2 = lamda*Ac'*vpv2
lamdaAc_At = -(lamda*A1cong'+A1tru'); - Khai báo các ma trận zeros
- Thành lập ma trận hệ thống : sys_mat = [ A zero1 zero2 zero4 A1 zero3
2*R 2*X lamdaAc_At ];
- Thành lập ma trân Achuyen có kích thước so_tcscxsonhanh suy ra từ Ma trận A ứng với các nút điều chỉnh .
- Thành lập ma trận hang stt_pq chứa số thứ tự nút PQ có bustype = 0, số nút PQ là sonut_pq .
- Thành lập ma trân hàng stt_nutdc kích thước 1xso_tcsc chứa chỉ số trong mảng stt_pq ứng vơi phần tử trong mảng stt_pq có nội dung bằng với phần tử trong mảng nut_dieuchinh .
- Thành lập các ma trân p2q2, l, m , k - Thành lập các ma trận y1, y2, yy như sau: y1 = [ pl ql lamda_Act_vpv2 ]; y2 = [ -Aphay*l -A1phay*m -k ]; yy = y1+y2+y_chuyen;
Trang 35 - Gỉai ẩn x từ Phương trình : x = inv(sys_mat_svc)*yy; Tổng số ẩn là N = 2*sonhanh+sonut-sonutpv_svc trong đó : p = x(1:sonhanh,1); q = x(sonhanh+1:2*sonhanh,1); v = x(2*sonhanh+1:N,1);
Tiến hành phép lặp lần lượt qua các bước sau :
- Tính ma trận cột vpv2 kích thước sonutpvx1 chứa bình phương điện áp nút vv2 của các nút PV , (bustype =1)
- Tính ma trận cột v kích thước (sonut-sonutpv)x1 chứa bình phương điện áp nút vv2 của các nút PQ, (bustype=0 )
- Ma trận cột ql kích thước (sonut-sonutpv)x1 chúa qload của các nút PQ (bustype=0)
- Tính các ma trận :
lamda_Act_vpv2 = lamda*Ac'*vpv2
lamdaAc_At = -(lamda*A1cong'+A1tru'); - Khai báo các ma trận zeros
- Thành lập ma trận hệ thống: sys_mat = [ A zero1 zero2 zero4 A1 zero3
2*R 2*X lamdaAc_At ];
- Thành lập ma trân Achuyen có kích thước so_tcscxsonhanh suy ra từ Ma trận A ứng với các nút điều chỉnh .
- Thành lập ma trận hang stt_pq chứa số thứ tự nút PQ có bustype = 0, số nút PQ là sonut_pq.
- Thành lập ma trân hàng stt_nutdc kích thước 1xso_tcsc chứa chỉ số trong mảng stt_pq ứng vơi phần tử trong mảng stt_pq có nội dung bằng với phần tử trong mảng nut_dieuchinh.
- Tính các ma trân p2q2, l như đã biết
Tính ma trận cột mottrenV2 kích thước (sonut-1)x1 mottrenV2(ii-1,1)=1/vv(ii)^2
với ii chay từ 2 đến sonut Tính các ma trân l như đã biết
l = p2q2*R*tru_Atru'*mottrenV2;
- Tính ma trận cột mottrenV2 kích thước (sonut-1)x1 mottrenV2(ii-1,1)=1/vv(ii)^2
Trang 36 Với ii chay từ 2 đến sonut
- Thành lập ma trận X_tcsc suy ra từ X vối các phần tử nhánh có đặt tcsc thì điều chỉnh lại như sau :
- Thành lập ma trận R2X2_tcsc suy ra từ R2X2 với các phần tử nhánh có đặt tcsc thì điều chỉnh lại như sau :
- Tính ma trận m và k như sau :
m = p2q2*X_tcsc*tru_Atru'*mottrenV2;
k = p2q2 * R2X2_tcsc * tru_Atru' * mottrenV2; - Thành lập các ma trận y1, y2, yy như sau :
y1 = [ pl ql lamda_Act_vpv2 ]; y2 = [ -Aphay*l -A1phay*m -k ]; yy = y1+y2+y_chuyen; - Gỉai ẩn x từ Phuong trình : x = inv(sys_mat_tcscvc)*yy; Tổng số ẩn là N = 2*sonhanh+sonut-sonutpv_tcsc trong đó : p = x(1:sonhanh,1); q = x(sonhanh+1:2*sonhanh,1); v = x(2*sonhanh+1:N,1);
- Qúa trình lặp tiếp tục cho đến khi sai số điện áp nhỏ hơn một độ chính xác yêu cầu - Khai báo mảng Xcbu :
for ii=1:so_tcsc
Xcbu(nhanh_tcsc(ii))=Xc(ii,1); end
- Sau khi hội tụ thực hiện các việc sau : + In dòng công suất tác dụng trên nhánh p. + In dòng công suất phản kháng trên nhánh q . + In tổn thất công suất tác dụng l và tổng tổn thất. + In tổn thất công suất tác dụng m và tổng tổn thất.
Trang 37 - In trị số dung kháng tụ bù dọc Xc
- Tính và in công suất nguồn Pnguon = p(1,1)+l(1,1); Qnguon = q(1,1)+m(1,1); Hay tính kiểm tra công suất nguồn :
Pnguon_kt = 0; Qnguon_kt =0; for ii =2 : sonut Pnguon_kt=Pnguon_kt-pload(ii); Qnguon_kt=Qnguon_kt-qload(ii); end Pnguon_kt = Pnguon_kt+deltaP; Qnguon_kt = Qnguon_kt+deltaQ;