找了个潮流计算的程序,按照程序里的例子输入却说串联的矩阵的维度不一致。错误在最后,可以拉到最后看。
能不能帮忙看看是什么地方出错,谢谢啦!
>> %简单潮流计算的牛顿拉夫逊程序,相关的原始数据数据数据输入格式如下:
%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));
end
OrgS=zeros(2*n-2,1);
DetaS=zeros(2*n-2,1); % 将 OrgS、DetaS 初始化
%创建 OrgS,用于存储初始功率参数 h=0;
j=0;
for i=1:n % 对 PQ节点的处理 if i~=isb&B2(i,6)==2
h=h+1; for j=1:n
OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*r
eal(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))* imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));
OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real( B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2