%*******************************
%Filename:output.m
%Author:Hweel_Zheng(鄭奕輝)     
%First created:2008.08.23      
%Last mended:2008.08.25
%******************************
fid=fopen(resultfile,'a');
fprintf(fid,'\r\n');
fprintf('---------------牛頓－拉夫遜潮流計算結果如下---------------\n');
fprintf(fid,'---------------牛頓－拉夫遜潮流計算結果如下---------------\r\n');
fprintf('共進行了%d次迭代\n',iloop);
fprintf('節點計算結果數據如下\n'); %tempbus
fprintf(fid,'節點計算結果數據如下\r\n');
fprintf('(bus#) ( volt ) ( ang )    (  p )      (  q ) (bus type)\n');
fprintf(fid,'(bus#) ( volt ) ( ang )    (  p )     (  q ) (bus type)\r\n');
for i=1:nb
    fprintf('%3d %10.6f %10.6f %10.6f %10.6fj %3d\n',tempbus(i,1),...
        tempbus(i,2),tempbus(i,3),tempbus(i,4),tempbus(i,5),...
        tempbus(i,6));
    fprintf(fid,'%3d %10.6f %10.6f %10.6f %10.6fj %3d\r\n',tempbus(i,1),...
        tempbus(i,2),tempbus(i,3),tempbus(i,4),tempbus(i,5),...
        tempbus(i,6));
end
fprintf('線路計算結果數據如下\n'); %flow
fprintf(fid,'線路計算結果數據如下\r\n');
fprintf('節點I 節點J    線路功率S(I,J)         線路功率S(J,I)'); 
fprintf(fid,'節點I 節點J    線路功率S(I,J)         線路功率S(J,I)');
fprintf('            線路損耗dS\n');
fprintf(fid,'            線路損耗dS\r\n');
for i=1:nl
    fprintf('%3d %4d %10.6f %10.6fj %10.6f %10.6fj '...
        ,flow(i,1),flow(i,2),real(flow(i,3)),imag(flow(i,3)),...
        real(flow(i,4)),imag(flow(i,4)));
    fprintf(fid,'%3d %4d %10.6f %10.6fj %10.6f %10.6fj '...
        ,flow(i,1),flow(i,2),real(flow(i,3)),imag(flow(i,3)),...
        real(flow(i,4)),imag(flow(i,4)));
    fprintf('%10.6f %10.6fj\n',real(flow(i,5)),imag(flow(i,5)));
    fprintf(fid,'%10.6f %10.6fj\r\n',real(flow(i,5)),imag(flow(i,5)));
end
fclose(fid);