《电力系统中的潮流计算程序》由会员分享,可在线阅读,更多相关《电力系统中的潮流计算程序(7页珍藏版)》请在人人文库网上搜索。
1、B2=input(请输入各节点参数形成的矩阵:B2=);%本程序的功能是用牛顿-拉夫逊法进行潮流计算n=input(请输入节点数:n=);nl=input(请输入支路数:nl=);isb=input(请输入平衡母线节电号:isb=);pr=input(请输入误差精度:pr=);B1=input(请输入由支路参数形成的矩阵:B仁);变压器侧为1,否则为0B2=input(请输入各节点参数形成的矩阵:B2=);X=input(请输入由节点号及其对地阻抗形成的矩阵:X=);X=i nput(请输入由节点号及其对地阻抗形成的矩阵:X=);Y=zeros( n);U=zeros(1, n);cta=ze。
2、ros(1, n);V=zeros(1, n);O=zeros(1, n);S1=zeros (n l);for i=1: nif X(i,2)=0;p=X(i,1);Y(p,p)=X(i,2);endendfor i=1: nlif B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5);Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5F2)+B1(i,4)./2;Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,。
3、4)./2;end %求导纳矩阵G=real(Y);B=imag(Y);for i=1: nc