直流潮流算法程序

x = xlsread（‘data1.xlsx‘‘B6 : D8‘）;
P = xlsread（‘data1.xlsx‘‘C3 : D3‘）;
P = P‘;
% idx= x==0; 
% x（idx） = inf; 
idx = find（isnan（x））; 
x（idx） = inf; 
b = 1 ./ x;
[rc] = size（b）;
B = zeros（r-1）;
for i = 1 : （r-1）
    for j = 1 : （r-1）
        if i ==j
            bi = b（i+1:）;
            B（ij） = sum（bi）;  
        else
            B（ij） = -b（i+1j+1）;
        end
    end
end