潮流计算牛拉法

潮流计算程序 %5.解修正方程 deltaVS=Jaco\-deltaPQ; %上半部是δ，上半部是△V/V deltaVS=deltaVS.*[ones(BusNum,1);abs(V)]; %求△V。(abs求绝对值，正值) %修正V和S的部分 angV=angle(V)+deltaVS(1:BusNum); %修正角度，1到n是δ modV=abs(V)+deltaVS(BusNum+1:BusNum*2); %修正幅值，n到n+1是V V=modV.*cos(angV)+i*modV.*sin(angV); %修正后的节点电