資源描述：
《潮流計算地MATLAB源程序》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、實用標準文案%簡單潮流計算的牛頓拉夫遜程序，相關(guān)的原始數(shù)據(jù)數(shù)據(jù)數(shù)據(jù)輸入格式如下：%B1是支路參數(shù)矩陣，第一列和第二列是節(jié)點編號。節(jié)點編號由小到大編寫%對于含有變壓器的支路，第一列為低壓側(cè)節(jié)點編號，第二列為高壓側(cè)節(jié)點%編號，將變壓器的串聯(lián)阻抗置于低壓側(cè)處理。%第三列為支路的串列阻抗參數(shù)。%第四列為支路的對地導(dǎo)納參數(shù)。%第五列為含變壓器支路的變壓器的變比%第六列為變壓器是否含有變壓器的參數(shù)，其中“1”為含有變壓器，%“0”為不含有變壓器。%B2為節(jié)點參數(shù)矩陣，其中第一列為節(jié)點注入發(fā)電功率參數(shù)；第二列為節(jié)點負荷功率
2、參數(shù)；第三列為節(jié)點電壓參數(shù)；第六列為節(jié)點類型參數(shù)，其中“1”為平衡節(jié)點，“2”為PQ節(jié)點，“3”為PV節(jié)點參數(shù)。%X為節(jié)點號和對地參數(shù)矩陣。其中第一列為節(jié)點編號，第二列為節(jié)點對地%參數(shù)。n=input('請輸入節(jié)點數(shù):n=');n1=input('請輸入支路數(shù):n1=');isb=input('請輸入平衡節(jié)點號:isb=');pr=input('請輸入誤差精度:pr=');B1=input('請輸入支路參數(shù):B1=');B2=input('請輸入節(jié)點參數(shù):B2=');X=input('節(jié)點號和對地參數(shù):X='
3、);Y=zeros(n);Times=1;???????????????%置迭代次數(shù)為初始值%創(chuàng)建節(jié)點導(dǎo)納矩陣fori=1:n1???ifB1(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
4、*B1(i,4);???else???????????????%含有變壓器的支路???????p=B1(i,1);???????q=B1(i,2);???????Y(p,q)=Y(p,q)-1/(B1(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));???end文檔大全實用標準文案endYOrgS=zeros(2*n-2,1);DetaS=z
5、eros(2*n-2,1);??%將OrgS、DetaS初始化%創(chuàng)建OrgS，用于存儲初始功率參數(shù)h=0;j=0;fori=1:n???????????%對PQ節(jié)點的處理???ifi~=isb&B2(i,6)==2?????????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(B2(i,3
6、))*(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???endendfori
7、=1:n??????????%對PV節(jié)點的處理，注意這時不可再將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(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(
8、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???end
