节点电压为极坐标情况下的直流潮流计算

一、 直流潮流的原理

二、例题

三、代码演示
这边仅仅展示核心代码。

% 求得直流潮流中的zita角度
zita = zeros(nb-1,1);
B = zeros(nb-1,nb-1);
for i = 1:nb - 1
    for j = 1:nb - 1
        B(i,j) = initial_B(i,j);
    end
end
disp(B)
B = inv(B)
P = bus(1:nb-1,4);
zita = B * P;
disp(zita)

zita_new = zeros(nb,1);
for i = 1:nb-1
    zita_new(i,1) = zita(i,1);
end
zita_new(nb,1) = 0;

% 记录直流潮流的功率运算
linePlowerFlow = zeros(nl,3);
for i = 1:nl
    I = line(i,1);
    J = line(i,2);
    linePlowerFlow(i,1) = line(i,1);
    linePlowerFlow(i,2) = line(i,2);
    linePlowerFlow(i,3) = ((zita_new(I)-zita_new(J))/line(i,4))* 100;
end