前推回代潮流計算程序(共4頁)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上function Flow1(Z,S,Sb,Vb)% Z 各支路阻抗% S 各節(jié)點運算負(fù)荷% Sb 基準(zhǔn)功率% Zb 基準(zhǔn)阻抗% Flow(Z,S,Sb,Vb) 計算網(wǎng)絡(luò)功率及電壓分布% 輸出收斂迭代次數(shù)及功率分布和電壓分布% 數(shù)據(jù)初始化Z(:,1)=Z(:,1)+1;Z1=1,0,1,0;Z=Z1;Z;%使得支路編號與支路末點編號一致m,n=size(Z);Zb=(Vb2)/Sb*1000;Z(:,4)=Z(:,4)/Zb;S(:,2)=S(:,2)/Sb; k=1;V=ones(m,1);dU=ones(m,1);%前推回代n=0;while max(abs(dU

2、)>10-4 Sr=S(:,2); V1=V; for i=m:-1:1 %前推功率 A=(find(Z(:,2)=i);%求以i節(jié)點為始點的各支路 b,c=size(A); if b=0 for j=1:1:b %求以當(dāng)前節(jié)點為始點的各支路功率總和 N=A(j,1); if N=6 dSr=(abs(Sr(N,1)/V(N,1)2)*Z(N,4); %各支路功率損耗 else dSr=0; end Sr(N,2)=Sr(N,1)+dSr; %N支路功率 Sr(i,1)=Sr(i,1)+Sr(N,2); end else Sr(i,2)=Sr(i,1); end end V(1,1)=1

3、; for t=2:m %回代電壓 if t=6 a=Z(t,2); %由于是輻射型網(wǎng)絡(luò)，只有一個始點 dV1=(real(Sr(t,2)*real(Z(t,4)+imag(Sr(t,2)*imag(Z(t,4)/V(a,1); dV2=(real(Sr(t,2)*imag(Z(t,4)-imag(Sr(t,2)*real(Z(t,4)/V(a,1); V(t,1)=(V(a,1)-dV1)2+dV22)0.5 ; else dk=k*(V(5,1)/V1(5,1)-1) k=k-dk if k>=1.1 k=1.1; elseif k<=0.9 k=0.9; else %k=vp

4、a(k,3); k=round(k*100)/100 end k V(6,1)=k*V(5,1); end end dU=V-V1; n=n+1; if n>15 fprintf('迭代15次不收！n'); break; endend%輸出a=1:1:m;%編號fprintf('迭代收斂次數(shù)為：n');nfprintf('網(wǎng)絡(luò)功率分布為：n');Sr(:,2)=Sr(:,2)*Sb;Sr=a',Srfprintf('網(wǎng)絡(luò)電壓分布為：n');%V=V*Vb;V=a',Vclcclear all;Z2= 1, 1

5、, 2, 1.197+j*0.82 ; 2, 2, 3, 1.796+j*1.231; 3, 3, 4, 1.306+j*0.895; 4, 4, 5, 1.851+j*1.268; 5, 5, 6, 1.524+j*1.044; 6, 6, 7, 1.905+j*1.305; 7, 7, 8, 1.197+j*0.82 ; 8, 8, 9, 0.653+j*0.447; 9, 9, 10, 1.143+j*0.783; 10, 4, 11, 2.823+j*1.172; 11, 11, 12, 1.184+j*0.491; 12, 12, 13, 1.002+j*0.416; 13, 13,

6、 14, 0.455+j*0.189; 14, 14, 15, 0.546+j*0.227; 15, 5, 16, 2.55+j*1.058; 16, 6, 17, 1.366+j*0.567; 17, 17, 18, 0.819+j*0.34 ; 18, 18, 19, 1.548+j*0.642; 19, 19, 20, 1.368+j*0.567; 20, 20, 21, 3.552+j*1.474; 21, 7, 22, 1.548+j*0.642; 22, 22, 23, 1.092+j*0.453; 23, 23, 24, 0.91+j*0.378; 24, 24, 25, 0.4

7、55+j*0.189; 25, 25, 26, 0.364+j*0.151; 26, 8, 27, 1.002+j*0.417; 27, 27, 28, 4.403+j*1.215; ;S2= 1, 0; 2, 35.28; 3, 14; 4, 35.28; 5, 14; 6, 35.28; 7, 35.28; 8, 35.28; 9, 14; 10, 14; 11, 56; 12, 35.28; 13, 35.28; 14, 14; 15, 35.28; 16, 35.28; 17, 8.96; 18, 8.96; 19, 35.28; 20, 35.28; 21, 14; 22, 35.28; 23, 8.96; 24, 56 25, 8.96; 26, 35.28; 27, 35.28; 28, 35.28;%28nodes dataJ= 5.441, 2.946 , 2.946 ; 2.946, 8.848 , 5.258 ; 2.946, 5.258 , 8.376 ;c=0.7;Vb=10.5;Sb=100;Z2(:,4