matlab潮流计算（2020年十二月整理）.docVIP

. PAGE PAGE 1 . 附录1 使用牛顿拉夫逊法进行潮流计算的Matlab程序代码 %牛拉法计算潮流程序 % %B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 %5、支路的变比；6、支路首端处于K侧为1，1侧为0 %B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 %4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 %6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；3为PV节点； % clearall; formatlong; n=input(请输入节点数:nodes=); nl=input(请输入支路数:lines=); isb=input(请输入平衡母线节点号:balance=); pr=input(请输入误差精度:precision=); B1=input(请输入由各支路参数形成的矩阵:B1=); B2=input(请输入各节点参数形成的矩阵:B2=); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); % fori=1:nl %支路数 ifB1(i,6)==0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else%左节点处于K侧 p=B1(i