output.m

程序实现用牛顿－拉夫逊迭代法进行潮流计算。不仅有潮流计算

%*******************************
%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);