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));

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