【正文】 )。請輸入節(jié)點數(shù):n=39。通過對電力系統(tǒng)潮流分布的分析和計算，可進一步對系統(tǒng)運行的安全性，經(jīng)濟性進行分析、評估，提出改進措施。ylabel(39。subplot(2,2,2)。ylabel(39。)。)。Sjy(i)=Sj(q,p)。,39。39。disp(V)。沒有達到精度要求的個數(shù)：39。%消去運算end%用當前行K2列元素減去當前行k列元素乘以第k行K2列元素J(k3,k)=0。J(p,N)=DV。J(p,N)=DV。J(m,q)=X1。% X3=dQ/de J(p,N)=DQ節(jié)點無功功率差J(m,q)=X1。%節(jié)點功率Q計算fiΣ(Gij*ejBij*fj)eiΣ(Gij*fj+Bij*ej)%求i節(jié)點有功和無功功率P39。%i節(jié)點注入功率SGSLB(i,i)=B(i,i)+B2(i,5)。q=B1(i,2)。0+002。+01。%input(39。DS(i)=Si(p,q)+Sj(q,p)。ZF=[39。,num2str(SSi(p,q))]。各條支路的首端功率Si為(順序同您輸入B1時一致)：39。)。for k=1:nV(k)=sqrt(e(k)^2+f(k)^2)。%當前行第3列元素已消為0endendend %====上面是用線性變換方式高斯消去法將Jacobi矩陣化成單位矩陣=====for k=3:2:N01L=(k+1)./2。J(m,N)=DP。J(m,N)=DP。J(m,N)=DP。% X1=dP/de J(m,N)=DP節(jié)點有功功率差J(p,q)=X4。%電壓模平方%========= 以下針對非PV節(jié)點來求取功率差及Jacobi矩陣元素 =========if B2(i,6)~=3%非PV節(jié)點DP=P(i)P1。%分解出各節(jié)點注入的有功和無功功率ICT1=0。endY(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。00002。B1=[1+0。n=10。q=B1(i,2)。SSj(q,p)=Sj(q,p)。)=39。disp(39。39。disp(ICT2)。%消去運算（當前行3列元素消為0）end%用當前行K2列元素減去當前行3列元素乘以第三行K2列元素J(k3,k)=0。J(m,q)=X1。J(m,q)=X1。J(p,q)=X4。q=q+1。的計算值V2=e(i)^2+f(i)^2。Q=imag(S)。q=B1(i,1)。00002。+01。%input(39。35=(+)()180。180。Se=SG3+Sd3025DST4DS04jQL3DSL3=(+)MVAFLR4： DST5P2+Q2502+312=(RT5+jXT5)=(+)=(+)MVA22VN110P2++=(RL4+jXL4)=(+j24)=(+)MVA22VN110Se39。l5=180。70=14W BL3=2180。104S DQ1L1=2BL1V2N==r0180。180。180。180。180。1102XT1=180。tan()=45MVar發(fā)電機G2：P2=4180。50MW“236。63MW“236。常規(guī)潮流計算的任務是根據(jù)給定的運行條件和網(wǎng)路結構確定整個系統(tǒng)的運行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡中的功率分布以及功率損耗等。潮流計算的結果是電力系統(tǒng)穩(wěn)定計算和故障分析的基礎。xd=====11kW239。xd====180。63=252MWQ2=252180。10=180。10= 32SN220180。10= 4SN463180。10=%180。10=180。l2=180。b40180。180。=(+)MVADSL4Sf=SG4+Se39。52==VbDVL2=()==功率分布：FLR3：SL3=ZZ+Z*T4*L3*SeT4(VG3Vd)Z+ZL3****180。(+)()180。請輸入節(jié)點數(shù)：n=39。+0。003。endY(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。%分解出各節(jié)點注入的有功和無功功率 ICT1=0。%電壓模平方%========= 以下針對非PV節(jié)點來求取功率差及Jacobi矩陣元素 =========if B2(i,6)~=3%非PV節(jié)點DP=P(i)P1。% X1=dP/de J(m,N)=DP節(jié)點有功功率差J(p,q)=X4。J(m,N)=DP。J(m,N)=DP。J(m,N)=DP。%當前行第3列元素已消為0endendend %====上面是用線性變換方式高斯消去法將Jacobi矩陣化成單位矩陣=====for k=3:2:N01L=(k+1)./2。for k=1:nV(k)=sqrt(e(k)^2+f(k)^2)。)。各條支路的首端功率Si為(順序同您輸入B1時一致)：39。,num2str(SSi(p,q))]。ZF=[39。DS(i)=Si(p,q)+Sj(q,p)。%input(39。+01。003。%非對角元Y(q,p)=Y(p,q)。IT2=1。%節(jié)點有功功率差DQ=Q(i)Q1。J(m,q)=X2。%擴展列△PJ(m,q)=X2。q=q+1。q=q+1。e(L)=e(L)J(k,N)。%計算各節(jié)點電壓的模值sida(k)=atan(f(k)./e(k))*180./pi。disp(39。)。disp(ZF)。S(39。disp(DS(i))。請輸入節(jié)點數(shù)：n=39。+0。00002。else%左節(jié)點處于K側p=B1(i,2)。%i節(jié)點無功補償量end %===================== P=real(S)。,Q39。J(m,N)=DP。q=q+1。% PV節(jié)點電壓誤差m=p+1。% PV節(jié)點電壓誤差m=p+1。%當前行第k列元素已消為0endelse%是第三行k=3%====================== 第三行k=3的前代運算 ========================for k3=k1:N0%從第四行到2n行（最后一行）for k2=k1:N1%從第四列到2n+1列（即擴展列）J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。)。%顯示各節(jié)點的電壓大小V的模值 disp(39。)。,num2str(q),39。enddisp(Sj(q,p))。for i=1:nlp=B1(i,1)。end如果源程序的運行結果需要作圖可用下面的程序figure(1)。電壓角度39。bar(Q)。支路首端注入有功39。同時潮流分布也是電力系統(tǒng)規(guī)劃設計的一項基礎工作。)。X=input(39。請輸入節(jié)點參數(shù):B2=39。n=input(39。電力系統(tǒng)分析的潮流計算是電力系統(tǒng)分析的一個重要的部分。)。grid on。)。39。各條支路的功率損耗DS為(順序同您輸入B1時一致)：39。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))))。,num2str(p),39。%顯示各節(jié)點的注入功率disp(39。)。disp(39。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)。J(p,q)=X5。J(p,q)=X5。%擴展列△Qm=p+1。m=p+1。%節(jié)點功率P計算eiΣ(Gij*ejBij*fj)+fiΣ(Gij*fj+Bij*ej)Q1=C(i)*f(i)e(i)*D(i)。%PV節(jié)點電壓給定模值end for i=1:n%給定各節(jié)點注入功率S(i)=B2(i,1)B2(i,2)。% % %for i=1:nl%支路數(shù)if B1(i,6)==0%左節(jié)點處于1側p=B1(i,1)。00002。B1=[1+0。n=12。q=B1(i,2)。SSj(q,p)=Sj(q,p)。)=39。disp(39。39。disp(ICT2)。%消去運算（當前行3列元素消為0）end%用當前行K2列元素減去當前行3列元素乘以第三行K2列元素J(k3,k)=0。J(m,q)=X1。J(m,q)=X1。J(p,q)=X4。q=q+1。的計算值V2=e(i)^2+f(i)^2。Q=imag(S)。q=B1(i,1)。0+002。)。end%本程序的功能是用牛頓——拉夫遜法進行10節(jié)點潮流計算 %本程序的功能是用牛頓——拉夫遜法進行潮流計算 clear。for i=1:nlp=B1(i,1)。enddisp(Sj(q,p))。,num2str(q),39。)。%顯示各節(jié)點的電壓大小V的模值 disp(39。)。%當前行第k列元素已消為0endelse%是第三行k=3%====================== 第三行k=3的前代運算 ========================for k3=k1:N0%從第四行到2n行（最后一行）for k2=k1:N1%從第四列到2n+1列（即擴展列）J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。% PV節(jié)點電壓誤差m=p+1。% PV節(jié)點電壓誤差m=p+1。q=q+1。J(m,N)=DP。,Q39。%i節(jié)點無功補償量end %===================== P=real(S)。else%左節(jié)點處于K側p=B1(i,2)。0+002。B1=[1+0。n=11。(+)+()180。(+)=()MVA 電壓分布：Sc1=SL2+DSL2=()+(+)=()180。Sd=(+)MVADSL339。r0180。180。40=180。352=180。VN=180。10=180。10=180。180。15=60MWQ1=60180。12MW1180。I%=%=,d1116MWVA4180?？梢哉f,它是電力系統(tǒng)分析中最基本、最重要的計算,是系統(tǒng)安全、經(jīng)濟分析和實時控制與調度的基礎。在電力系統(tǒng)運行方式和規(guī)劃方案的研究中，都需要進行潮流計算以比較運行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟性。Dps=50kW237。25MW”236。tan()=156MVar發(fā)電機G3：P3=3180。10=180。10DS02=2180。103I0%180。VN=180。180。130==x0180。l3=2180。20==180。80DST5DS05jQL4DSL4=(+)MVAFLR5: 152DST7=2180。VN=()180。35==(+)：ST631=ST63+DST63=(+)180。)。+01。0+002。%非對角元Y(q,p)=Y(p,q)。IT2=1。%節(jié)點有功功率差DQ=Q(i)Q1。J(m,q)=X2。%擴展列△PJ(m,q)=X2。q=q+1。q=q+1。e(L)=e(L)J(k,N)。%計算各節(jié)點電壓的模值sida(k)=atan(f(k)./e(k))*180./pi。disp(39。)。disp(ZF)。S(39。disp(DS(i))。請輸入節(jié)點數(shù)：n=39。+0。0+002]。%非對角元Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。N0=2*n。%節(jié)點無功功率差%=============== 以上為除平衡節(jié)點外其它節(jié)點的功率計算 ================= %================= 求取Jacobi矩陣 ===================for j1=1:nif j1~=isbamp。% X4=dQ/df X2=dp/dfelseif j1==iamp。endendelse%=============== 下面是針對PV節(jié)點來求取Jacobi矩陣的元素===========DP=P(i)P1。J(p,q)=X6。J(p,q)=X6。%修改節(jié)點電壓實部k1=k+1。%計算各節(jié)點電壓的角度E(k)=e(k)+f(k)*j。各節(jié)點的電壓相角sida為(節(jié)點號從小到大排列)：39。for i=1:nlp=B1(i,1)。disp(39。,num2str(q),39。DDS(i)=DS(i)。)。+01。003。q=B1(i,1)。Q=imag(S)。的計算值V2=e(i)^2+f(i)^2。q=q+1。J(p,q)=X4。J(m,q)=X1。J(m,q)=X1。%消去運算（當前行3列元素消為0）end%用當前行K2列元素減去當前行3列元素乘以第三行K2列元素J(k3,k)=0。disp(ICT2)。39。disp(39。)=39。SSj(q,p)=Sj(q,p)。q=B1(i,2)。subplot(1,2,1)。)。xlabel(39。)。整個計算過程的模型建立并不是十分復雜，但計算過程十分繁瑣、計算量相當?shù)拇螅矣捎谥?jié)太多很容易算錯。n1=input(39。節(jié)點號和對地參數(shù):X=39。B2=input(39。其中第一列為節(jié)點編號，第二列為節(jié)點對地 %參數(shù)。電力系統(tǒng)的穩(wěn)定運行狀況即是正常運行狀況