run_nr.m

来自「这个程序是采用模块化的方式对潮流进行求解」· M 代码 · 共 22 行

function run_NR()
global G B V delta P Q  n m nl myf bus line l J 
N=10000; 
eps1=1.0e-10;
eps2=1.0e-10;
for l=1:N
    [delP,delQ] = dPQ(Y,bus);
     J=form_jac(bus,Y);       % 求取jacci矩阵 (雅克比矩阵) 
     b = [delP;delQ];
    dx=J\b;           %矩阵左除，即b除以A
    ddelta= dx(1:n-1,1);
    dV    = dx(n:n+m-1,1).*V(1:m,1);
    
    delta(1:n-1,1) = delta(1:n-1,1) - ddelta(1:n-1,1);
    V(1:m,1) = V(1:m,1) - dV(1:m,1);
    
    printData(l,J,delP,delQ,ddelta,dV);
    
    if(norm(ddelta)<eps1 & norm(dV)<eps2)
        break;
    end
end