📄 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 + -