?? compensation2.m
字號:
0 0 0 0 0 0 0 0 0 0 E(7) F(7) -E(6) -F(6) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 -F(7) E(7) F(6) -E(6) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 E(8) F(8) -E(7) -F(7) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 -F(8) E(8) F(7) -E(7) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 E(9) F(9) -E(8) -F(8) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 -F(9) E(9) F(8) -E(8) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 E(10) F(10) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -E(2) -F(2) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 -F(10) E(10) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 F(2) -E(2) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 E(11) F(11) 0 0 0 0 0 0 0 0 0 0 0 0 -E(4) -F(4) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 -F(11) E(11) 0 0 0 0 0 0 0 0 0 0 0 0 F(4) -E(4) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 E(12) F(12) -E(11) -F(11) 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -F(12) E(12) F(11) -E(11) 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 E(13) F(13) -E(12) -F(12) 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -F(13) E(13) F(12) -E(12) 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 E(14) F(14) -E(13) -F(13) 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -F(14) E(14) F(13) -E(13) 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 E(15) F(15) -E(14) -F(14) 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -F(15) E(15) F(14) -E(14) 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 E(16) F(16) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -E(6) -F(6) 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 -F(16) E(16) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 F(6) -E(6) 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 E(17) F(17) -E(16) -F(16) 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -F(17) E(17) F(16) -E(16) 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 E(18) F(18) -E(17) -F(17) 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -F(18) E(18) F(17) -E(17) 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 E(19) F(19) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -E(9) -F(9)
0 0 0 0 0 0 0 0 0 0 0 0 0 0 -F(19) E(19) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 F(9) -E(9)];
Jacobi=B*Y*C ; %雅各比矩陣
epQPcorrection=epQP; %Q、P的修正量
Jacobizz=inv(Jacobi);
epEFcorrection=Jacobizz*epQPcorrection; %E、F的修正量
epQP=epQPcorrection;
epEF=epEFcorrection;
k=k+1;
for j=2:numPQ+1
E(j-1)=E(j-1)+epEF(2*(j-1)-1)*E(j-1);
F(j-1)=F(j-1)+epEF(2*(j-1));
if Rmin==0&&j==20
EEcopy(flag)=E(j-1);
end
if Rmin==1&&j==2
EEEcopy(flag)=E(j-1);
end
end
end
clear epEFcorrection;
clear epQPcorrection;
clear Jacobizz;
end %對應次外層while循環
PQ=PQcopy;
if flag==1
linepower=zeros(numPQ+1); %線路功率分布
for j=1:numPQ+1
for m=1:numPQ+1
if j~=m
if j==1
linepower(j,m)=(E1+i*F1)*((E1-i*F1)-(E(m-1)-i*F(m-1)))*(YBre(1,m)-i*YBim(1,m))*0.1;
elseif m==1
linepower(j,m)=(E(j-1)+i*F(j-1))*((E(j-1)-i*F(j-1))-(E1-i*F1))*(YBre(j,1)-i*YBim(j,1))*0.1;
else
linepower(j,m)=(E(j-1)+i*F(j-1))*((E(j-1)-i*F(j-1))-(E(m-1)-i*F(m-1)))*(YBre(j,m)-i*YBim(j,m))*0.1;
end
end
end
end
linepower=linepower;
lineloss=0;
for j=1:numPQ+1
for m=1:numPQ+1
lineloss=lineloss+linepower(j,m);
end
end
lineloss=lineloss %總功率損耗
end
end %對應最外層for循環
fprintf(flowoutput,'\n');
if Rmin==0
kk=0; %存放補償后對應E(19)的最大值
Kmax=0; %存放補償后對應E(19)最大時的補償節點號
for j=1:4
if kk<EEcopy(j)
kk=EEcopy(j);
Kmax=j;
end
end
fprintf(flowoutput,'當選擇兩個補償點時由電壓恢復情況知,補償最優組合為組合%g\n',Kmax);
if kk>=0.9
fprintf(flowoutput,'并且由結果知:最大負荷情況下選擇該補償組合時電壓能滿足供電要求(電壓標幺值達0.9)\n');
Rmin=2;
else
fprintf(flowoutput,'并且由結果知:當選擇該組合進行無功補償時網絡中節點電壓仍不能滿足供電要求(電壓標幺值達0.9)\n');
break;
end
else
if Kmax==1&&EEEcopy(1)<1.12
fprintf(flowoutput,'由結果知:最小負荷情況下選擇上面補償組合時不會出現過補償現象,故本次畢業設計考慮選擇該組合中的兩節點作為所給配網饋線上的無功集中補償點\n');
else
fprintf(flowoutput,'由結果知:最小負荷情況下選擇上面補償組合時會出現過補償現象,故上面補償組合是不可取的');
end
end
fprintf(flowoutput,'\n');
end %對應最外層while循環
fprintf(flowoutput,'\n');
fclose(flowoutput);
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -