基于 MATLAB 的电力系统潮流计算
%简单潮流计算的小程序，相关的原始数据数据数据输入格式如下：
%B1 是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写
%对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点
%编号，将变压器的串联阻抗置于低压侧处理。
%第三列为支路的串列阻抗参数。
%第四列为支路的对地导纳参数。
%第五烈为含变压器支路的变压器的变比
%第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，
%“0”为不含有变压器。
%B2 为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点
%负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中
%“1”为平衡节点，“2”为 PQ 节点，“3”为 PV 节点参数。
%X 为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地
%参数。
n=input(‘请输入节点数:n=‘);
n1=input(‘请输入支路数:n1=‘);
isb=input(‘请输入平衡节点号:isb=‘);
pr=input(‘请输入误差精度:pr=‘);
B1=input(‘请输入支路参数:B1=‘);
B2=input(‘请输入节点参数:B2=‘);
X=input(‘节点号和对地参数:X=‘);
Y=zeros(n);
Times=1; %置迭代次数为初始值
%创建节点导纳矩阵
for i=1:n1
if B1(i,6)==0 %不含变压器的支路
p=B1(i,1);
q=B1(i,2);
Y(p,q)=Y(p,q)-1/B1(i,3);
Y(q,p)=Y(p,q);
Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);
Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4);
else %含有变压器的支路
p=B1(i,1);
q=B1(i,2);
Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5));
Y(q,p)=Y(p,q);
Y(p,p)=Y(p,p)+1/B1(i,3);
Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3));
评论0