MATLAB实现电力系统七节点潮流计算代码解析

%% 利用牛顿-拉夫逊法进行潮流计算
%% author：张永辉
%% 程序说明：为了编程方便，将平衡节点1改为节点7，其他节点号依次减1，即原节点2为现节点1.
% B1矩阵 ：支路参数矩阵 每一列分别代表：起始节点  终止节点  支路阻抗 支路对地导纳  
% B2矩阵 ：节点参数矩阵 每一列分别代表：节点注入功率  负荷消耗功率  节点初始电压 节点初始角度 节点类型  
% T1矩阵：发电机注入功率                     T2矩阵：负荷消耗的功率
% Node_U_It矩阵：最终得到的电压幅值矩阵      angle_It矩阵： 最终得到的电压相角矩阵
% Node_S_It矩阵：最终得到的节点注入功率矩阵  Branch_It矩阵：最终得到的支路功率损耗矩阵
disp('【 节点数:】');
%[n1]=xlsread('C:\Users\lenovo\Desktop\input.xlsx','A3:A3')%节点数
[n1]=xlsread('input.xlsx','A3:A3')%节点数
disp('【 支路数:】');
[n]=xlsread('input.xlsx','B3:B3')%支路数
disp('【 精度:】');
Accuracy=xlsread('input.xlsx','B4:B4')%精度
[branch]=xlsread('input.xlsx','E4:K10');
[node]=xlsread('input.xlsx','M4:S10');
Data_B1=branch;%支路参数
Data_B2=node;%节点参数
T1=zeros(n,2);
T2=zeros(n1,3);
i=sqrt(-1);
format short
for j=1:n
    T1(j,1)=Data_B1(j,3)+Data_B1(j,4)*1i;
    T1(j,2)=Data_B1(j,5)*1i;
end
for j=1:n1
    T2(j,1)=Data_B2(j,1)+Data_B2(j,2)*1i;
    T2(j,2)=Data_B2(j,3)+Data_B2(j,4)*1i;

end
B1=zeros(n,6);
B2=zeros(n1,5);
for j=1:n
    B1(j,1)=Data_B1(j,1);
    B1(j,2)=Data_B1(j,2);
    B1(j,3)=T1(j,1);
    B1(j,4)=T1(j,2);
    B1(j,5)=Data_B1(j,6);
    B1(j,6)=Data_B1(j,7);
end
for j=1:n1
    B2(j,1)=T2(j,1);
    B2(j,2)=T2(j,2); 
    B2(j,3)=Data_B2(j,5);
    B2(j,4)=Data_B2(j,6);
    B2(j,5)=Data_B2(j,7);
end
disp('【 支路参数矩阵:】');
B1  %显示支路参数矩阵
disp('【 节点参数矩阵:】');
B2  %显示节点参数矩阵
% 以上为从excel中导入初值的程序
Y=zeros(n1);
for i=1:n

        p=B1(i,1);
        q=B1(i,2);
        Y(p,q)=Y(p,q)-1/B1(i,3);
        Y(q,p)=Y(p,q);