3. Đối tượng và phạm vi nghiên cứu của đề tài
3.2.2. Các chương trình con trong chương trình tính phân bố công suất
và điện áp viết theo phương pháp GS
3.2.2.1. Chương trình con tạo ma trận Ybus
Ma trận Ybus được thành lập như sơ đồ thuật toán hình 2.3. Đầu tiên thành lập các tổng dẫn tương hổ Yij giữa các nút (các phần tử nằm ngoài đường chéo chính của Ybus) , sau khi tất cả các tổng dẫn tương hổ được xác định thì các tổng dẫn riêng cũng được xác định Yii = - ∑Yij với j≠i.
Đoạn chương trình con thành lập Ybus như sau: %---phan tu ngoai duong cheo chinh---% Ybus=zeros(n,n);
for e=1:nh %nh: so nhanh
Ybus(linedata(e,2),linedata(e,3))=- 1/(linedata(e,5)+j*linedata(e,6))*Zcb; Ybus(linedata(e,3),linedata(e,2))= Ybus(linedata(e,2),linedata(e,3)); end
%---phan tu tren duong cheo chinh---% for e=1:n %n: so nut
Ybus(e,e)=0;
Đồ Án Chuyên Ngành Ứng dụng Matlab for f=1:n if f~=e Ybus(e,e)=Ybus(e,e)-Ybus(e,f); end end end
3.2.2.2. Chương trình con tính điện áp các nút tải
Công suất P, Q tại các nút tải đã biết nên trong mỗi bước lặp điện áp các nút này được tính toán trước.
for e=2:n
if typebus(e)==3 %xac nhan nut tai% Ut=0; for f=1:n if f~=e Ut=Ut-Ybus(e,f)*U(f); end end Ut=(conj(Snut(e))/conj(U(e))+Ut)/Ybus(e,e); U(e)=Ut; %ket thuc tinh toan cho nut tai thu e end
end
3.2.2.3. Chương trình con tính CSPK và điện áp các nút PV
Sau khi điện áp các nút PQ được hình thành, CSPK các nút P V được xác định và điện áp các nút này cũng được tính toán.
%---tinh CSPK QPV---% for e=2:n
if typebus(e)==2 %xac nhan nut PV% QPV=0; for f=1:n
QPV=QPV+Ybus(e,f)*U(f); %tong Yij*Uj (j gom ca i) end
QPV=-imag(conj(U(e))*QPV); %conj(Ui)*(tong Yij*Uj)
Đồ Án Chuyên Ngành Ứng dụng Matlab
Snut(e)=Pcho(e)+j*QPV; cong suat nut PV thu e end
end
%---tinh dien ap nut PV---% for e=2:n
if typebus(e)==2 %xac nhan nut PV% UPV=0; for f=1:n
if f~=e
UPV=UPV-Ybus(e,f)*U(f); %tong Yij*Uj (j khac i) end end UPV=(conj(Snut(e))/conj(U(e))+UPV)/Ybus(e,e); Ui=imag(UPV);Ur=sqrt(Um(e)^2-Ui^2);U(e)=Ur+j*Ui; end end
3.3. CHƯƠNG TRÌNH CON TÍNH PHÂN BỐ CÔNG SUẤT VÀ ĐIỆN ÁP CHO MẠNG ĐIỆN PHÂN PHỐI VIẾT THEO PHƯƠNG PHÁP ÁP CHO MẠNG ĐIỆN PHÂN PHỐI VIẾT THEO PHƯƠNG PHÁP