39 lines
1.9 KiB
Matlab
39 lines
1.9 KiB
Matlab
function [e]=OutResult(G,B,balNode,e,f,T,k,Gap,fx,transI,transJ,filename,nodeNum,...
|
||
PGR,pgNode,QR,pvNode,capI,capK,firstT,firstcapK,Tn,Tn1,capKn,capKn1,firstfx,toc,max_hx,ifconverge,Vfirst,Pgfirst,Qgfirst)
|
||
%程序功能：输出计算结果
|
||
%编写时间：2010年11月
|
||
|
||
mytime=toc;
|
||
Y1=G+1j*B;
|
||
V1=e+1j*f;
|
||
NodeS=zeros(nodeNum,1);
|
||
NodeS(pvNode)=V1(pvNode,:).*conj(Y1(pvNode,:)*V1); %计算平衡节点和PV节点的节点复功率，Y1(pvNode,:)*V1为节点的注入电流
|
||
balanceP=real(NodeS(balNode,1)); %获取平衡节点的有功功率
|
||
PGR(balNode,1)=balanceP;
|
||
balanceQ=imag(NodeS(balNode,1));
|
||
QR(balNode,1)=balanceQ;
|
||
V=sqrt(e.^2+f.^2);
|
||
resultFile=['result_' filename];
|
||
fid=fopen(resultFile,'w');
|
||
result=[k;Gap;fx;firstfx;mytime;max_hx;ifconverge]; %收敛信息
|
||
fprintf(fid,'%2.0f %4.4e %6.6f %6.6f %4.4f %4.4e %2.0f\r\n',result);
|
||
fprintf(fid,'%1.0f\r\n',0);
|
||
trans=[transI';transJ';(1./T)';(1./firstT)';(1./Tn)';(1./Tn1)']; %变压器变比结果
|
||
fprintf(fid,'%2.0f %2.0f %5.4f %5.4f %5.4f %5.4f\r\n',full(trans));
|
||
fprintf(fid,'%1.0f\r\n',0);
|
||
nn=1:nodeNum;
|
||
Vresult=[nn;e';f';V';Vfirst']; %节点电压结果
|
||
fprintf(fid,'%3.0f %5.4f %5.4f %5.4f %5.4f\r\n',full(Vresult));
|
||
fprintf(fid,'%1.0f\r\n',0);
|
||
PGresult=[pgNode';PGR(pgNode)';Pgfirst(pgNode)']; %有功出力结果
|
||
PGresult=full(PGresult);
|
||
fprintf(fid,'%2.0f %4.4f %4.4f\r\n',PGresult);
|
||
fprintf(fid,'%1.0f\r\n',0);
|
||
QRresult=[pvNode';QR(pvNode)';Qgfirst(pvNode)']; %无功出力结果
|
||
QRresult=full(QRresult);
|
||
fprintf(fid,'%2.0f %4.4f %4.4f\r\n',QRresult);
|
||
fprintf(fid,'%1.0f\r\n',0);
|
||
Capresult=[capI';capK';firstcapK';capKn';capKn1']; %电容电抗器投入组数结果
|
||
fprintf(fid,'%2.0f %3.2f %5.4f %3.2f %3.2f\r\n',full(Capresult));
|
||
fprintf(fid,'%1.0f\r\n',0);
|
||
fclose(fid); |