【正文】 0 0。 0 0。 0 0。0000 0。0000 0]。B=[ 00。 15 0。 0。00 0。000]。Y=G+j*B。X=[]。%X=△X n=m+v+1。%總的節(jié)點(diǎn)數(shù)FX=ones(2*n2,1)。%F(x)矩陣F1=zeros(n1,n1)。%F(x)導(dǎo)數(shù)矩陣 a=0。%記錄迭代次數(shù)EF=zeros(n1,n1)。%最后的節(jié)點(diǎn)電壓矩陣 while max(FX)=10^(5)for i=1:m %PQ節(jié)點(diǎn)FX(i)=P(i)。%△P FX(n+i1)=Q(i)。%△Qfor w=1:n FX(i)= FX(i)E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)F(i)*G(i,w)*F(w)F(i)*B(i,w)*E(w)。%△PFX(n+i1)=FX(n+i1)F(i)*G(i,w)*E(w)+F(i)*B(i,w)*F(w)+E(i)*G(i,w)*F(w)+E(i)*B(i,w)*E(w)。%△Qend end for i=m+1:n1 %PV節(jié)點(diǎn) FX(i)=PP(im)。%△P FX(n+i1)=V(im)^2E(i)^2F(i)^2。%△Qfor w=1:n FX(i)= FX(i)E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)F(i)*G(i,w)*F(w)F(i)*B(i,w)*E(w)。%△Pend endfor i=1:m %PQ節(jié)點(diǎn)for w=1:n1if i~=wF1(i,w)=(G(i,w)*E(i)+B(i,w)*F(i))。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。F1(n+i1,w)=B(i,w)*E(i)G(i,w)*F(i)。F1(n+i1,n+w1)=G(i,w)*E(i)+B(i,w)*F(i)。elseF1(i,w)=G(i,i)*E(i)B(i,i)*F(i)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,w)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,n+w1)=G(i,i)*E(i)+B(i,i)*F(i)。for k=1:nF1(i,w)=F1(i,w)G(i,k)*E(k)+B(i,k)*F(k)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。F1(n+i1,w)=F1(n+i1,w)+G(i,k)*F(k)+B(i,k)*E(k)。F1(n+i1,n+w1)=F1(n+i1,n+w1)G(i,k)*E(k)+B(i,k)*F(k)。endendend end for i=m+1:n1 %PV節(jié)點(diǎn)for w=1:n1if i~=wF1(i,w)=(G(i,w)*E(i)+B(i,w)*F(i))。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。F1(n+i1,w)=0。F1(n+i1,n+w1)=0。elseF1(i,w)=G(i,i)*E(i)B(i,i)*F(i)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,w)=2*E(i)。F1(n+i1,n+w1)=2*F(i)。for k=1:nF1(i,w)=F1(i,w)G(i,k)*E(k)+B(i,k)*F(k)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。endendend end X=inv(F1)*(FX)。for i=1:n1E(i)=E(i)+X(i)。F(i)=F(i)+X(n+i1)。end a=a+1。fprintf(39。第%d次迭代后的節(jié)點(diǎn)電壓分別為:n39。,a)。disp(E+j*F)。fprintf(39。第%d次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:n39。,a)。disp(FX)。enddisp(39。收斂后的節(jié)點(diǎn)電壓用極坐標(biāo)表示為:39。)。EF=E+j*F。for i=1:n1 fprintf(39。%d號節(jié)點(diǎn)電壓的幅值為:39。,i)disp(abs(EF(i)))。fprintf(39。%d號節(jié)點(diǎn)電壓的相角度數(shù)為39。,i)disp(angle(EF(i))*180/pi)。end PPH=0。for i=1:nPPH=PPH+EF(n)*conj(Y(n,i))*conj(EF(i))。end fprintf(39。平衡節(jié)點(diǎn)的功率39。)。disp(PPH)。運(yùn)行結(jié)果：運(yùn)行結(jié)果復(fù)制如下：第1次迭代后的節(jié)點(diǎn)電壓分別為:第1次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:00第2次迭代后的節(jié)點(diǎn)電壓分別為:第2次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:第3次迭代后的節(jié)點(diǎn)電壓分別為:第3次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:第4次迭代后的節(jié)點(diǎn)電壓分別為:第4次迭代后功率偏差△P △Q電壓偏差△V的平方分別為: *收斂后的節(jié)點(diǎn)電壓用極坐標(biāo)表示為: 1號節(jié)點(diǎn)電壓的幅值為:2號節(jié)點(diǎn)電壓的幅值為:3號節(jié)點(diǎn)電壓的幅值為:4號節(jié)點(diǎn)電壓的幅值為:平衡節(jié)點(diǎn)的功率 第三篇：潮流計算電力系統(tǒng)大作業(yè)(C++編寫)程序設(shè)計所涉及二叉樹解釋任意一棵樹或一個森林都能唯一地對應(yīng)一棵二叉樹，由此而編寫本程序。本程序采用類二叉樹為整體結(jié)構(gòu)，二叉樹類下定義節(jié)點(diǎn)類，每一條支路均為樹的一個節(jié)點(diǎn)，支路所有的參數(shù)均作為節(jié)點(diǎn)的屬性，并給節(jié)點(diǎn)加入屬性“支路編號”，并以支路編號為依據(jù)構(gòu)建二叉樹，這就要求提前根據(jù)二叉樹結(jié)構(gòu)給每一個支路編號。支路編號原則：左子樹上所有編號均小于其雙親的編號，右子樹上所有編號均大于其雙親的編號，為了便于查看，本程序在節(jié)點(diǎn)較少時編號從1開始，逐個遞加至支路數(shù)；當(dāng)支路較多時，可不必拘泥于逐個遞加，只要滿足支路編號原則即可。例如習(xí)題34：++104SDST=++j30MVA20+j15MVA程序二叉樹結(jié)構(gòu)示意圖：本二叉樹中節(jié)點(diǎn)1即支路1為題目中節(jié)點(diǎn)1與節(jié)點(diǎn)2之間的部分； 本二叉樹中節(jié)點(diǎn)2即支路2為題目中節(jié)點(diǎn)2與節(jié)點(diǎn)3之間的部分； 對于習(xí)題33：程序二叉樹結(jié)構(gòu)示意圖：本二叉樹中節(jié)點(diǎn)2即支路2為題目中節(jié)點(diǎn)1與節(jié)點(diǎn)2之間的部分； 本二叉樹中節(jié)點(diǎn)1即支路1為題目中節(jié)點(diǎn)2與節(jié)點(diǎn)3之間的部分； ：如下多支路網(wǎng)絡(luò)：對于三節(jié)點(diǎn)網(wǎng)絡(luò)需先進(jìn)行以下網(wǎng)絡(luò)處理轉(zhuǎn)化為標(biāo)準(zhǔn)二叉樹，而后進(jìn)行計算；圖中三角形表明該項阻抗為零，為純導(dǎo)線，并進(jìn)行相應(yīng)參數(shù)補(bǔ)充進(jìn)行計算。程序說明文檔******************************************************************************************************************************************************* 本程序測試使用方法：在E盤根目錄下建立輸入文件：輸入文件名：。將所附算例對應(yīng)輸入文件內(nèi)容復(fù)制粘貼至上述文件中，在VC++。輸出文件在E盤根目錄下輸出文件名為：。******************************************************************************************************************************************************* 程序功能說明：本程序可以計算任意長度線型開始網(wǎng)絡(luò)潮流；支持多電壓等級下的計算；可在除供電節(jié)點(diǎn)外任意節(jié)點(diǎn)引出負(fù)載；*******************************************************************************************************************************************************輸入格式說明：以支路為基本單位，按潮流方向輸入數(shù)據(jù)：以下例示意：//支路個數(shù)1,110,118,0,0,1,0,0//線路等效的支路2,11,110,0,40,30,10,//變壓器等效的支路.//按此方式知道輸入所有的支路.//輸入從上到下的順序?yàn)槌绷髟?//線型開式網(wǎng)絡(luò)中的流動方向/*第一行輸入支路個數(shù)，回車第二行至后輸入各個支路參數(shù)，回車分隔不同支路；各行輸入的支路參數(shù)順序是：支路編號，末端電壓，始端電壓，線路等效電阻，線路等效感抗，線路等效容納，末端輸入有功，末端輸入無功，變比，變壓器有功勵磁損耗，變壓器無功勵磁損耗對于線路等效電路：變壓器有功/無功損耗輸入零，變比輸入1；對于變壓器等效電路：所有的參數(shù)均歸算至高壓側(cè)，Rt,Xt對應(yīng)輸入線路等效電阻/感抗的位置，線路等效容抗為零；對于個節(jié)點(diǎn)的引出負(fù)荷：輸入至以此節(jié)點(diǎn)為末節(jié)點(diǎn)的支路的末端輸入有功/無功部分；*/ ******************************************************************************************************************************************************** 輸出文件格式說明：輸出的內(nèi)容包括（1）支路信息：每個支路元件的始端有功、無功和末端有功、無功；有功損耗無功損耗；電壓損耗；（2）全網(wǎng)信息：全網(wǎng)的總電源有功、總負(fù)荷有功、有功損耗、網(wǎng)損率；（3）迭代信息：每次完整迭代后的所有內(nèi)容；具體在輸出文件中都明確標(biāo)出。********************************************************************************************************************************************************* 程序中變量定義說明：類中定義的變量class line_part{//定義支路類double U[2]。//支路電壓降落：U[0]電壓降落橫分量，U[1]電壓降落縱分量double k。//變壓器變比double val。//支路排序double U_end。//支路末端電壓double U_begin。//支路首段電壓double X[3]。//支路等效阻抗: X[0]電阻，X[1]感抗，X[2]容納double S_end[2]。//支路末端功率：S_end[0]有功，S_end[1]無功double S_begin[2]。//支路首段功率：S_begin[0]有功，S_begin[1]無功double S0[2]