output_iterate.m

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

%*******************************
%Filename:output_iterate.m
%Author:Hweel_Zheng(郑奕辉)     
%First created:2008.08.23      
%Last mended:2008.08.25
%******************************
%output_iterate
%输出迭代过程中的相关数据
fid=fopen(resultfile,'a');
fprintf(fid,'----------------------第%d次迭代结果------------------\r\n',...
    iloop);
fprintf(fid,'--------------------------------------第%d次迭代的',iloop);
fprintf(fid,'雅可比矩阵Jac-----------------------------------\r\n');
for i=1:iP+iQ
    fprintf(fid,'第%d行:',i);
    for jl=1:iP+iQ
        fprintf(fid,'%10.6f',Jac(i,jl));
    end
    fprintf(fid,'\r\n');
end
fprintf(fid,'\r\n-------------------------------------第%d次迭代的',iloop);
fprintf(fid,'功率偏差△P和△Q-----------------------------------\r\n');
for i=1:iP
    fprintf(fid,'△P(%d):',i);
    fprintf(fid,'%10.6f \r\n',deltaPQ(i));
end
for i=1:iQ
    fprintf(fid,'△Q(%d):',i);
    fprintf(fid,'%10.6f \r\n',deltaPQ(i+iP));
end
fprintf(fid,'\r\n-------------------------------------第%d次迭代的',iloop);
fprintf(fid,'节点相角和电压偏差△δ和△U------------------------------\r\n');
for i=1:iP
    fprintf(fid,'△δ(%d):',i);
    fprintf(fid,'%10.6f \r\n',DctaU(i));
end
for i=1:iQ
    fprintf(fid,'△U(%d):',i);
    fprintf(fid,'%10.6f \r\n',DctaU(i+iP));
end
fprintf(fid,'\r\n-------------------------------------第%d次迭代的',iloop);
fprintf(fid,'节点相角和电压δ和U------------------------------\r\n');
for i=1:iP
    fprintf(fid,'δ(%d):',i);
    fprintf(fid,'%10.6f \r\n',ctaU(i));
end
for i=1:iQ
    fprintf(fid,'U(%d):',i);
    fprintf(fid,'%10.6f \r\n',ctaU(i+iP));
end
fclose(fid);