潮流计算的MATLAB源程序.docxVIP

(完整版)潮流计算的MATLAB 源程序 (完整版)潮流计算的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 end Y 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）)＊real(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（i,3））＊（real(Y(i， j)）＊imag(B2(j,3）)+imag（Y(i，j))＊real（B2（j，3））); end end end for i=1:n ％对 PV 节点的处理，注意这时不可再将 h 初始化为 0 if i～=isbB2（i，6)==3 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）)＊real（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(i,3））*(real（Y（i,j)） *imag(B2(j,3)）+imag（Y（i，j))＊real(B2(j，3)))； end end end OrgS ％创建 PVU 用于存储 PV 节点的初始电压PVU=zeros(n-h—1，1）； t=0; for i=1:n if B2（i,6)==3 t=t+1； PVU（t,1)=B2(i,3）; end end PVU %创建 De