【正文】 比矩陣換行k=0。DetaU=inv(Jacbi)*DetaS。for i=1:n %對PQ節(jié)點處理 if B2(i,6)==2 j=j+1。end end for i=1:n %對PV節(jié)點的處理 if B2(i,6)==3 j=j+1。end end B2 %開始循環(huán)********************************************************************** while abs(max(DetaU))pr OrgS=zeros(2*n2,1)。j=0。B2(i,6)==2 h=h+1。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)))。B2(i,6)==3 h=h+1。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)))。for i=1:n if i~=isbamp。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。end end t=0。B2(i,6)==3 h=h+1。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。end end DetaS %創(chuàng)建I i=zeros(n1,1)。for i=1:n if i~=isb h=h+1。end end I %創(chuàng)建Jacbi Jacbi=zeros(2*n2)。k=0。for j=1:n if j~=isb k=k+1。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1))。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)2*imag(I(h,1))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)。end end end end end k=0。for j=1:n if j~=isb k=k+1。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1))。Jacbi(2*h,2*k)=2*real(B2(i,3))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k)=0。end end end end end Jacbi DetaU=zeros(2*n2,1)。DetaU %修正節(jié)點電壓 j=0。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。%迭代次數加1 end Times 一個原始數據的例子 節(jié)點數 5 支路數 5平衡節(jié)點編號 5 精度pr B1(支路參數矩陣)[1 2 + 1 0。2 3 + 1 0。5 3 0 1] B2(節(jié)點參數矩陣)[ 1 0 0 2。 1 0 0 2。0 0 0 1] X(節(jié)點號和對地參數)[1 0。3 0。5 0]電力系統潮流計算——9結點算例PQ法:標號,起始結點,終止結點,支路電阻參數,支路電抗參數,支路對地導納參數 1,2,5, 2,5,9, 3,6,9, 4,3,6, 5,6,8, 6,4,8, 7,5,7, 8,4,7, 9,1,4, : include include define N 9 /*總結點數*/ define M 6 /*PQ結點數*/ define K 9 /*線路數*/ define eps 1e4 void guass(int n,int m,float c[],float b[][N],float x[])/*高斯函數*/ { float a[N][N],y[N]。for(i=0。i){ for(j=i+1。j=tEndnode1。lx=tX。lg=lr/(lr*lr+lx*lx)。g[i][i]+=lg。b[i][i]+=lb+lb1。h[i][j]=h[j][i]=lb1。g[j][i]=lg。b[j][i]=lb。printf(“n=====jie dian dao na ju zhen=====n”)。i。for(i=1。} if(max=eps){ for(i=0。for(i=1。} if(max=eps){ for(i=0。for(i=NM。if(kq==0)val(u,g,b,r,ku,kr,h)。} } } void val(float u[N],float g[N][N],float b[N][N],float r[N],int ku, int kr,float h[N][N]){ float ps=0,pv1=0,pv2=0。float p[N][N]={0}。float s[N][N]。float dq[N][N]={0}。float dSp=0,dSq=0。FILE *fp1。getch()。iprintf(”n=======shu ju bao cun=====n“)。{ fprintf(fp1,”xian lu cao liu:n“)。i第二篇：用matlab電力系統潮流計算題目：潮流計算與matlab教學單位 電氣信息學院姓 名 學 號年 級專 業(yè) 電氣工程及其自動化指導教師職 稱 副教授摘 要電力系統穩(wěn)態(tài)分析包括潮流計算和靜態(tài)安全分析。本位就是運用潮流計算具體分析，并有MATLAB仿真。Flow calculation?？梢哉f,它是電力系統分析中最基本、最重要的計算,是系統安全、經濟分析和實時控制與調度的基礎。潮流計算的結果是電力系統穩(wěn)定計算和故障分析的基礎。同時，為了實時監(jiān)控電力系統的運行狀態(tài)，也需要進行大量而快速的潮流計算。在系統規(guī)劃設計和安排系統的運行方式時，采用離線潮流計算；在電力系統運行狀態(tài)的實時監(jiān)控中，則采用在線潮流計算。它既是對電力系統規(guī)劃設計和運行方式的合理性、可靠性及經濟性進行定量分析的依據，又是電力系統靜態(tài)和暫態(tài)穩(wěn)定計算的基礎。對系統進行潮流分析，本文利用 MATLAB中的SimpowerSystems工具箱設計電力系統，在simulink 環(huán)境下，不僅可以仿真系統的動態(tài)過程，還可以對系統進行穩(wěn)態(tài)潮流分析。在加人測量模塊，并對各元件的參數進行設置后，用measurement和sink中的儀器可以觀察各元件的電壓、電流、功率的大小。63MWVA15MWGGGG120MW10kV236。I%=%=,d1116MWVA4180。xd=====11kW239。I%=%=,d112180。20MWVA20MWGGG4180。xd======89kW237。I%=%=,d112180。Dp0=44kW239。I%=%=,Y,d112180。12MW1180。xd====180。xd========“d在Simulink中按照電力系統原型選擇元件進行建模。在電力系統模型的建立工程中主要涉及到的是：元器件的選擇及其參數的設置；發(fā)電機選型；變壓器選擇；線路的選擇；負荷模型的選擇；母線選擇。運行的具體結果和分析也在詳細設計中詳述。Dp0=44kW239。I%=%===63kW239。Vs(23)%=(13)%= 詳細設計 計算前提首先是發(fā)電機的參數計算，先對5個發(fā)電廠簡化為5臺發(fā)電機來計算。15=60MWQ1=60180。63=252MWQ2=252180。12=36MWQ3=36180。50=50MWQ4=50180。25=50MWQ5=50180。變壓器T1：2Dps180。11023RT1=180。103=(16180。180。10=180。103DS01=Dp0+j變壓器T2：（雙并聯）I0%180。VN189180。180。180。10)21Vs%180。1102=180。10=180。10= 32SN220180。(Dp0+j變壓器T3：（四并聯）I0%180。VN21121180。180。180。10)XT31Vs%180。1102=180。10=180。10= 4SN463180。SN)=(+)MVA100DS03=4180。3523=180。10= 322(10180。[Vs(12)%+Vs(13)%Vs(23)%]=%=180。[Vs(13)%+Vs(23)%Vs(12)%]=%=21Vs1%180。180。VN=180。10=%180。180。(DP06+j變壓器T7：（雙并聯）I0%180。VN150180。180。180。10)21Vs%180。352=180。10=180。10=180。(Dp0+jI0%180。根據教材查得r0=180。l1=180。l1=180。l1=180。40=180。l2=180。l2=180。l2=180。130=180。rl1L30180。180。x0180。180。b40180。180。70=180。l4=180。l4=180。l4=180。60=180。r0180。180。x0180。180。b0180。180。20=180。=(+)MVADSL2=FLRSd=SG2+Sc39。Sd=(+)MVADSL339。=(+)MVADSL4Sf=SG4+Se39。(+)=(+)MVA35Sh=15+DST7+DS07+jDQL5=(+)+=180。(+)=(+)+=180。(+)=()MVA352Sg=SfDST61+SG5DST62DST63Si35=(+)MVA計算每一個FLR的功率分布和電壓分布計算如下： FLR1：DVT2=PR+QX40180。== Vb=115DVT2=+180。16DVL1====VbDVL1=：功率分布：SL2=ZZ+Z*T3*L2*Sd+T3(VbVN)Z+ZL2****180。(+)=()MVAST3=ZZ+Z*L2*L2*Sd+T3(VbVN)Z+ZL2180。(+)=()MVA 電壓分布：Sc1=SL2+DSL2=()+(+)=()180。52==VbDVL2=()==功率分布：FLR3：SL3=ZZ+Z*T4*L3*SeT4(VG3Vd)Z+ZL3****180。(+)=()MVAST4=ZZ+Z*L3*L3*Se+T4(VbVN)Z+ZL3180。(+)+=(+)MVA 電壓分布：Se1=SL3+DSL3=()+(+)=(+)180。14==VdDVL3==功率分布：FLR4：SL4=ZZ+Z*T5*L4*SfT5(VG3Vd)Z+ZL4****180。(+)=(+)MVAST4=ZZ+Z*L4*L4*SfT5(VG3Vd)Z+ZL4180。(+)=(+)MVA 電壓分布：Se1=SL3+DSL3=(+)+(+)=(+)180。24==VdDVL3==FLR5: 這里我們先將f點和發(fā)電機G5當做電源，經過ZT61和ZT63構成兩端供電網絡以g點作為運算負荷進行計算。(+)+()180。(+)()180。+180。+180。+180。節(jié)點編號由小到大編寫。B2為節(jié)點參數矩陣，其中第一列為節(jié)點注入發(fā)電功率參數；第二列為節(jié)點負荷功率參數；第三列為節(jié)點電壓參數；第六列為節(jié)點類型參數，其中“1”為平衡節(jié)點，“2”為PQ節(jié)點，“3”為PV節(jié)點參數。其中第一列為節(jié)點編號，第二列為節(jié)點對地參數。n=11。請輸入節(jié)點數：n=39。nl=11。請輸入支路數：nl=39。isb=1。請輸入平衡母線節(jié)點號：isb=39。pr=。請輸入誤差精度：pr=39。B1=[1+0。+0。+0。+01。+0。+0] B2=[0001。0+002。0+002。0+002。003。00002]。請輸入各節(jié)點參數形成的矩陣： B2=39。Y=zeros(n)。f=zeros(1,n)。sida=zeros(1,n)。% % %for i=1:nl%支路數if B1(i,6)==0%左節(jié)點處于1側p=B1(i,1)。else%左節(jié)點處于K側p=B1(i,2)。endY(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。%非對角元Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。%對角元1側end %求導納矩陣disp(39。)。B=imag(Y)。f(i)=imag(B2(i,3))。%PV節(jié)點電壓給定模值end for i=1:n%給定各節(jié)點注入功率S(i)=B2(i,1)B2(i,2)。%i節(jié)點無功補償量end %===================== P=real(S)。%分解出各節(jié)點注入的有功和無功功率 ICT1=0。N0=2*n。a=0。a=a+1。D(i)=0。%Σ(Gij*ejBij*fj)D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1)。%節(jié)點功率P計算eiΣ(Gij*ejBij*fj)+fiΣ(Gij*fj+Bij*ej)Q1=C(i)*f(i)e(i)*D(i)。,Q39。%電壓模平方%========= 以下針對非PV節(jié)點來求取功率差及Jacobi矩陣元素 =========if B2(i,6)~=3%非PV節(jié)點DP=P(i)P1。%節(jié)點無功功率差%==