【正文】 rend of the network and the network of the ponents of the power loss, and getting electric standard is to use the power flow calculation and analysis, the specific have MATLAB words: Power system。常規(guī)潮流計算的任務是根據(jù)給定的運行條件和網(wǎng)路結(jié)構(gòu)確定整個系統(tǒng)的運行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡中的功率分布以及功率損耗等。因此，潮流計算是電力系統(tǒng)中應用最廣泛、最基本和最重要的一種電氣運算。潮流計算經(jīng)歷了一個由手工到應用數(shù)字電子計算機的發(fā)展過程，現(xiàn)在的潮流算法都以計算機的應用為前提用計算機進行潮流計算主要步驟在于編制計算機程序，這是一項非常復雜的工作。10MWYN,d114180。63MW“236。10MWVA35kV32km25MW110kV80km25MW110kV70km110kVYN,d112180。I%=%===73kW237。Dps=121kW10kV237。50MW“236。所建立的模型和建立的方法在詳細設計中詳述。30km35kV0km31236。I0%=%=(12)80MW239。tan()=45MVar發(fā)電機G2：P2=4180。tan()=27MVar發(fā)電機G4：P4=1180。tan()=其次是變電站的參數(shù)計算，我們還是對7個變電站簡化為7臺變壓器來計算。10=180。1102XT1=180。SN=(+)MVA 100RT2XT221Dps180。10=180。180。180。SN)=(+)MVA100 1Dps180。10=180。180。180。(Dp0+j變壓器T4：（雙并聯(lián)）1RT1= XT4=XT1==2DS01=(+)MVART4=變壓器T5：RT5=4RT3==4XT3==1DS03=(+)MVA4163180。10)變壓器T6：（兩個三繞組變壓器并聯(lián)）RT61=RT62=RT631180。VNXT61=180。180。10= 2SNDS06=2180。3523=180。103=(10180。180。103DS07=2180。10S/km 6線路L1：線路L2：線路L3：（雙回路）線路L4：線路L5：（雙回路）RL1=r0180。40=16WB6L1=b0180。104S DQ1L1=2BL1V2N==r0180。130=52WB6L2=b0180。104S DQ1L2=2BL2V2N==12180。70==2180。70=14W BL3=2180。106180。60==x0180。106180。l5=180。l5=180。l5=2180。104S1DQL5=BL3VN2=FLR1：P2102DST1=2(RT1+jXT1)=(+)=(+)MVA2VN110Sa=10MW+DST1+DS01+jDQL1=(+)MVAP2++=(RL1+jXL1)=(+j16)=(+)MVAVN21102DST2P2+Q2402+452=(RT2+jXT2)=(+)=(+)MVAVN21102FLR2Sb=SG120DST2=60+=(+)MVASc=SbSa25jQL1DSL1=(+)MVA：DST3P2+Q22522+1562=(RT3+jXT3)=(+)=(+)MVA22VN110P++(R+jX)=(+j52)=(+)MVAL2L222VN1102222Sc39。Se=SG3+Sd3025DST4DS04jQL3DSL3=(+)MVAFLR4： DST5P2+Q2502+312=(RT5+jXT5)=(+)=(+)MVA22VN110P2++=(RL4+jXL4)=(+j24)=(+)MVA22VN110Se39。(+j4)=(+)MVA352Si=Sh+DSL5+DS06+jDQL5+5=(+)MVA 152+=180。+45180。VN=()180。180。VN=()180。VN=()180。+180。35=(+)()180。==37DVT63==180。4==ViDVL5==4設計圖及源程序相關的原始數(shù)據(jù)輸入格式如下：B1是支路參數(shù)矩陣，第一列和第二列是節(jié)點編號。X為節(jié)點號和對地參數(shù)矩陣。%input(39。%input(39。%input(39。%input(39。+01。+0。+01。00002。00002。% B1矩陣：支路首端號；末端號；支路阻抗；支路對地電納 %支路的變比；支路首端處于K側(cè)為1，1側(cè)為0 % B2矩陣：該節(jié)點發(fā)電機功率；該節(jié)點負荷功率；節(jié)點電壓初始值 %PV節(jié)點電壓V的給定值；節(jié)點所接的無功補償設備的容量 %節(jié)點分類標號：1為平衡節(jié)點（應為1號節(jié)點）；2為PQ節(jié)點； %3為PV節(jié)點；%input(39。e=zeros(1,n)。S1=zeros(nl)。q=B1(i,1)。%對角元K側(cè)Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。disp(Y)%G=real(Y)。V(i)=B2(i,4)。Q=imag(S)。N=N0+1。for i=1:nif i~=isb%非平衡節(jié)點C(i)=0。%Σ(Gij*fj+Bij*ej)endP1=C(i)*e(i)+f(i)*D(i)。的計算值V2=e(i)^2+f(i)^2。j1~=i%非平衡節(jié)點amp。% X2=dp/df X3=dQ/deX4=X1。J(p,N)=DQ。q=q+1。j1~=isb %非平衡節(jié)點amp。% dQ/deX4=C(i)+G(i,i)*e(i)+B(i,i)*f(i)。J(p,N)=DQ。J(p,q)=X4。% PV節(jié)點有功誤差DV=V(i)^2V2。% dP/deX2=B(i,j1)*e(i)G(i,j1)*f(i)。q=2*j11。J(m,q)=X1。% PV節(jié)點有功誤差J(m,q)=X2。% dP/deX2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。q=2*j11。J(m,q)=X1。% PV節(jié)點有功誤差J(m,q)=X2。% 用K行K列對角元素去除K行K列后的非對角元素進行規(guī)格化endJ(k,k)=1。%當前行第k列元素已消為0endif k==N0%若已到最后一行break。%消去運算（當前行3列元素消為0）end%用當前行K2列元素減去當前行3列元素乘以第三行K2列元素J(k3,k)=0。f(L)=f(L)J(k1,N)。%不滿足要求的節(jié)點數(shù)ICT1=ICT1+1。disp(ICT1)。disp(ICT2)。%將各節(jié)點電壓用復數(shù)表示 end %=============== 計算各輸出量 =========================== disp(39。%顯示各節(jié)點的實際電壓標幺值E用復數(shù)表示 disp(39。各節(jié)點的電壓大小V為(節(jié)點號從小到大排列)：39。39。)。%計算各節(jié)點的注入電流的共軛值endS(p)=E(p)*C(p)。disp(S)。disp(39。q=B1(i,2)。Siz(i)=Si(p,q)。S(39。)=39。39。)。Sjy(i)=Sj(q,p)。SSj(q,p)=Sj(q,p)。,39。disp(ZF)。end disp(39。q=B1(i,2)。ZF=[39。,num2str(q),39。disp(39。n=10。nl=10。isb=1。pr=。B1=[1+0。+0。+01。00002。00002。請輸入各節(jié)點參數(shù)形成的矩陣： B2=39。f=zeros(1,n)。% % %for i=1:nl%支路數(shù)if B1(i,6)==0%左節(jié)點處于1側(cè)p=B1(i,1)。endY(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。%對角元1側(cè)end %求導納矩陣disp(39。B=imag(Y)。%PV節(jié)點電壓給定模值end for i=1:n%給定各節(jié)點注入功率S(i)=B2(i,1)B2(i,2)。%分解出各節(jié)點注入的有功和無功功率ICT1=0。a=0。D(i)=0。%節(jié)點功率P計算eiΣ(Gij*ejBij*fj)+fiΣ(Gij*fj+Bij*ej)Q1=C(i)*f(i)e(i)*D(i)。%電壓模平方%========= 以下針對非PV節(jié)點來求取功率差及Jacobi矩陣元素 =========if B2(i,6)~=3%非PV節(jié)點DP=P(i)P1。非對角元X1=G(i,j1)*e(i)B(i,j1)*f(i)。% X1=dP/de X4=dQ/dfp=2*i1。m=p+1。% X1=dP/de J(m,N)=DP節(jié)點有功功率差J(p,q)=X4。對角元X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。% dQ/dfp=2*i1。%擴展列△Qm=p+1。J(m,N)=DP。% PV節(jié)點電壓誤差for j1=1:nif j1~=isbamp。% dP/dfX5=0。J(p,q)=X5。J(m,N)=DP。elseif j1==iamp。% dP/dfX5=2*e(i)。J(p,q)=X5。J(m,N)=DP。endendendendend %========= 以上為求雅可比矩陣的各個元素及擴展列的功率差或電壓差 =====================for k=3:N0% N0=2*n（從第三行開始，第一、二行是平衡節(jié)點）k1=k+1。% 對角元規(guī)格化K行K列對角元素賦1%==================== 回代運算=======================================if k~=3% 不是第三行k 3k4=k1。end%================== 前代運算==================================for k3=k1:N0% 從k+1行到2*n最后一行for k2=k1:N1% 從k+1列到擴展列消去k+1行后各行下三角元素J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。%當前行第3列元素已消為0endendend %====上面是用線性變換方式高斯消去法將Jacobi矩陣化成單位矩陣=====for k=3:2:N01L=(k+1)./2。%修改節(jié)點電壓虛部end%修改節(jié)點電壓for k=3:N0DET=abs(J(k,N))。%迭代次數(shù) end %用高斯消去法解”w=J*V“ disp(39。disp(39。for k=1:nV(k)=sqrt(e(k)^2+f(k)^2)。各節(jié)點的實際電壓標幺值E為(節(jié)點號從小到大排列)：39。39。)。)。disp(sida)。%計算各節(jié)點的功率 S = 電壓 X 注入電流的共軛值 end disp(39。%顯示各節(jié)點的注入功率disp(39。各條支路的首端功率Si為(順序同您輸入B1時一致)：39。if B1(i,6)==0Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))...conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。enddisp(Si(p,q))。,num2str(p),39。,num2str(SSi(p,q))]。)。for i=1:nlp=B1(i,1)。elseSj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))...conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。ZF=[39。,num2str(p),39。disp(39。各條支路的功率損耗DS為(順序同您輸入B1時一致)：39。DS(i)=Si(p,q)+Sj(q,p)。DS(39。)=39。39。%input(39。%input(39。%input(39。%input(39。+01。+0。+01。0+002。0+002。00002。)。V=zeros(1,n)。q=B1(i,2)。%非對角元Y(q,p)=Y(p,q)。導納矩陣 Y=39。%分解出導納陣的實部和虛部for i=1:n%給定各節(jié)點初始電壓的實部和虛部e(i)=real(B2(i,3))。%i節(jié)點注入功率SGSLB(i,i)=B(i,i)+B2(i,5)。IT2=1。%迭代次數(shù)ICTa；不滿足收斂要求的節(jié)點數(shù)IT2 while IT2~=0% N0=2*n 雅可比矩陣的階數(shù)；N=N0+1擴展列IT2=0。for j1=1:nC(i)=C(i)+G(i,j1)*e(j1)B(i,j1)*f(j1)。%節(jié)點功率Q計算fiΣ(Gij*ejBij*fj)eiΣ(Gij*fj+Bij*ej)%求i節(jié)點有功和無功功率P39。%節(jié)點有功功率差DQ=Q(i)Q1。% dP/de=dQ/dfX2=B(i,j1)*e(i)G(i,j1)*f(i)。q=2*j11。% X3=dQ/de J(p,N)=DQ節(jié)點無功功率差J(m,q)=X1。J(m,q)=X2。% dP/deX2=D(i)+B(i,i)*e(i)G(i,i)*f(i