printdata.m

这个程序是采用模块化的方式对潮流进行求解

 function printData(l,J,delP,delQ,ddelta,dV)
 global l J V delta n m myf  
 
fprintf(myf,'--------------第%d次迭代的结果----------\n',l);
%雅可比矩阵J
fprintf(myf,'-------第%d次迭代的雅可比矩阵J-------\n',l);
for i=1:n+m-1
    for j=1:n+m-1
        fprintf(myf, '       %12.6f', J(i,j));
    end
    fprintf(myf, '\n');
end
fprintf(myf, '\n');

%功率偏差dP和dQ
fprintf(myf,'-------------第%d次迭代的功率偏差dP和dQ----------\n',l);
for i = 1:n-1
    fprintf(myf,'dP%d   %12.6e\n',i,delP(i));
end
for i = 1:m
    fprintf(myf,'dQ%d   %12.6e\n',i,delQ(i));
end
fprintf(myf, '\n');

%节点相角和电压的偏差dx
fprintf(myf,'-------------第%d次迭代的节点相角和电压的偏差dx----------\n',l);
for i = 1:n-1
    fprintf(myf,'dang%d   %12.6e\n',i,ddelta(i));
end
for i = 1:m
    fprintf(myf,'dU%d     %12.6e\n',i,dV(i));
end
fprintf(myf, '\n');

%节点相角和电压的偏差dx
fprintf(myf,'-------------第%d次迭代的节点相角delta(弧度为单位）和电压U----------\n',l);
for i = 1:n-1
    fprintf(myf,'ang%d   %12.6f\n',i,delta(i));
end
for i = 1:m
    fprintf(myf,'U%d     %12.6f\n',i,V(i));
end
fprintf(myf, '\n');