44 lines
1.9 KiB
Mathematica
44 lines
1.9 KiB
Mathematica
|
function OutResult(e,f,T,k,Gap,fx,transI,transJ,filename,nodeNum,PGR,pgNode,QR,pvNode,capI,capK,firstfx,firstT,firstcapK,Tn,Tn1,...
|
|||
|
|
capKn,capKn1,ifconverge,max_hx,G,B,balNode,toc,firstV,capBi)
|
|||
|
|
%<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ܣ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
|||
|
|
%<EFBFBD><EFBFBD>дʱ<EFBFBD>䣺2010<EFBFBD><EFBFBD>10<EFBFBD><EFBFBD>
|
|||
|
|
|
|||
|
|
Y1=G+1j*B;
|
|||
|
|
V1=e+1j*f;
|
|||
|
|
NodeS=zeros(nodeNum,1);
|
|||
|
|
NodeS(pvNode)=V1(pvNode,:).*conj(Y1(pvNode,:)*V1); %<EFBFBD><EFBFBD><EFBFBD><EFBFBD>ƽ<EFBFBD><EFBFBD><EFBFBD>ڵ<EFBFBD><EFBFBD><EFBFBD>PV<EFBFBD>ڵ<EFBFBD><EFBFBD>Ľڵ㸴<EFBFBD><EFBFBD><EFBFBD>ʣ<EFBFBD>Y1(pvNode,:)*V1Ϊ<EFBFBD>ڵ<EFBFBD><EFBFBD><EFBFBD>ע<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
|||
|
|
balanceP=real(NodeS(balNode,1)); %<EFBFBD><EFBFBD>ȡƽ<EFBFBD><EFBFBD><EFBFBD>ڵ<EFBFBD><EFBFBD><EFBFBD><EFBFBD>й<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
|||
|
|
PGR(balNode,1)=balanceP;
|
|||
|
|
balanceQ=imag(NodeS(balNode,1));
|
|||
|
|
QR(balNode,1)=balanceQ;
|
|||
|
|
resultFile=['result_' filename];
|
|||
|
|
fid=fopen(resultFile,'w');
|
|||
|
|
result=[k;Gap;fx;firstfx;toc;max_hx;ifconverge]; %<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>Ϣ
|
|||
|
|
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)']; %<EFBFBD><EFBFBD>ѹ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>Ƚ<EFBFBD><EFBFBD><EFBFBD>
|
|||
|
|
fprintf(fid,'%2.0f %2.0f %5.4f %5.4f %5.4f %5.4f\r\n',trans);
|
|||
|
|
fprintf(fid,'%1.0f\r\n',0);
|
|||
|
|
nn=1:nodeNum;
|
|||
|
|
e=full(e);
|
|||
|
|
f=full(f);
|
|||
|
|
V=sqrt(e.^2+f.^2);
|
|||
|
|
Vresult=[nn;e';f';V';firstV']; %<EFBFBD>ڵ<EFBFBD><EFBFBD><EFBFBD>ѹ<EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
|||
|
|
fprintf(fid,'%3.0f %5.4f %5.4f %5.4f %5.4f\r\n',Vresult);
|
|||
|
|
fprintf(fid,'%1.0f\r\n',0);
|
|||
|
|
PGresult=[pgNode';PGR(pgNode)']; %<EFBFBD>й<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
|||
|
|
PGresult=full(PGresult);
|
|||
|
|
fprintf(fid,'%2.0f %4.4f\r\n',PGresult);
|
|||
|
|
fprintf(fid,'%1.0f\r\n',0);
|
|||
|
|
QRresult=[pvNode';QR(pvNode)']; %<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
|||
|
|
QRresult=full(QRresult);
|
|||
|
|
fprintf(fid,'%2.0f %4.4f\r\n',QRresult);
|
|||
|
|
fprintf(fid,'%1.0f\r\n',0);
|
|||
|
|
capK=capK.*capBi;
|
|||
|
|
capKn=capKn.*capBi;
|
|||
|
|
capKn1=capKn1.*capBi;
|
|||
|
|
Capresult=[capI';capK';firstcapK';capKn';capKn1']; %<EFBFBD><EFBFBD><EFBFBD>ݵ翹<EFBFBD><EFBFBD>Ͷ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
|||
|
|
Capresult=full(Capresult);
|
|||
|
|
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);
|