powerresult.m

来自「这个程序是采用模块化的方式对潮流进行求解」· M 代码 · 共 27 行

function [P,Q]=PowerResult()
global P Q n m G B V delta myf nodenum ;


for i = m+1:n
    for j =1:n
        Q(i) = Q(i) + V(i)*V(j)* ( G(i,j)* sin(delta(i) - delta(j)) -  B(i,j)* cos(delta(i) - delta(j)));
    end
end

for j =1:n
    P(n) = P(n) + V(i)*V(j)* (  G(i,j)* cos(delta(n) - delta(j)) + B(i,j)* sin(delta(n) - delta(j)));
end


fprintf(myf,'---------------牛顿－拉夫逊法潮流计算结果----------\n');
fprintf(myf,'节点计算结果：\n');
fprintf(myf,'节点          节点电压           节点相角（角度）            节点注入功率\n');
for i = 1:n
    for k=1:n
    if (nodenum(k,2)==i)
    t= nodenum(k,1);
    fprintf(myf,'%4d   %12.6f      %12.6f        %12.6f+j%12.6f\n',i,V(t),delta(t)*360/2/pi,P(t),Q(t));
    end
end

end