潮流計算的MATLAB源程序

《潮流計算的MATLAB源程序》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。

1、%簡單潮流計算的牛頓拉夫遜程序，相關(guān)的原始數(shù)據(jù)數(shù)據(jù)數(shù)據(jù)輸入格式如下：%B1是支路參數(shù)矩陣，第一列和第二列是節(jié)點(diǎn)編號。節(jié)點(diǎn)編號由小到大編寫%對于含有變壓器的支路，第一列為低壓側(cè)節(jié)點(diǎn)編號，第二列為高壓側(cè)節(jié)點(diǎn)%編號，將變壓器的串聯(lián)阻抗置于低壓側(cè)處理。%第三列為支路的串列阻抗參數(shù)。%第四列為支路的對地導(dǎo)納參數(shù)。%第五列為含變壓器支路的變壓器的變比%第六列為變壓器是否含有變壓器的參數(shù)，其中“1”為含有變壓器，%“0”為不含有變壓器。%B2為節(jié)點(diǎn)參數(shù)矩陣，其中第一列為節(jié)點(diǎn)注入發(fā)電功率參數(shù)；第二列為節(jié)點(diǎn)負(fù)荷功率參數(shù)；第三列為節(jié)點(diǎn)電壓參數(shù)；第六列為節(jié)點(diǎn)類型參數(shù)，其中“1”為平衡節(jié)點(diǎn)，“

2、2”為PQ節(jié)點(diǎn)，“3”為PV節(jié)點(diǎn)參數(shù)。%X為節(jié)點(diǎn)號和對地參數(shù)矩陣。其中第一列為節(jié)點(diǎn)編號，第二列為節(jié)點(diǎn)對地%參數(shù)。n=input('請輸入節(jié)點(diǎn)數(shù):n=');n1=input('請輸入支路數(shù):n1=');isb=input('請輸入平衡節(jié)點(diǎn)號:isb=');pr=input('請輸入誤差精度:pr=');B1=input('請輸入支路參數(shù):B1=');B2=input('請輸入節(jié)點(diǎn)參數(shù):B2=');X=input('節(jié)點(diǎn)號和對地參數(shù):X=');Y=zeros(n);Times=1;???????????????%置迭代次數(shù)為初始值%創(chuàng)建節(jié)點(diǎn)導(dǎo)納矩陣fori=1:n1???i

3、fB1(i,6)==0??????%不含變壓器的支路???????p=B1(i,1);???????q=B1(i,2);???????Y(p,q)=Y(p,q)-1/B1(i,3);???????Y(q,p)=Y(p,q);???????Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);???????Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4);???else???????????????%含有變壓器的支路???????p=B1(i,1);???????q=B1(i,2);???????Y(p,q)=Y(p,q)-1/(B

4、1(i,3)*B1(i,5));???????Y(q,p)=Y(p,q);???????Y(p,p)=Y(p,p)+1/B1(i,3);???????Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3));???endendYOrgS=zeros(2*n-2,1);DetaS=zeros(2*n-2,1);??%將OrgS、DetaS初始化%創(chuàng)建OrgS，用于存儲初始功率參數(shù)h=0;j=0;fori=1:n???????????%對PQ節(jié)點(diǎn)的處理???ifi~=isb&B2(i,6)==2?????????h=h+1;???????forj=1:n???

5、????????OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));???????????OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2

6、(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));???????end???endendfori=1:n??????????%對PV節(jié)點(diǎn)的處理，注意這時不可再將h初始化為0???ifi~=isb&B2(i,6)==3???????h=h+1;???????forj=1:n???????????OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(

7、B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));???????????OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));???????end???endendOrgS%創(chuàng)建PVU用于存