jac.m

利用该MATLAB程序可实现对一般电力系统的节点、支路及其损耗的潮流计算。

function J=jac(bus,Y)
%形成雅可比矩阵
global n;
global m;
for i=1:(n-1)  %计算PQ、PV节点的有功
    P1=0;T=0;
    for j=1:n
        T=bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)-bus(j,3))+imag(Y(i,j))*sin(bus(i,3)-bus(j,3)));
        P1=P1+T;
    end
    bus(i,4)=P1;
end
for i=1:n-1  %计算PV、PQ节点的无功
    Q1=0;T=0;
    for j=1:n
        T=bus(i,2)*bus(j,2)*(real(Y(i,j))*sin(bus(i,3)-bus(j,3))-imag(Y(i,j))*cos(bus(i,3)-bus(j,3)));
        Q1=Q1+T;
    end
    bus(i,5)=Q1;
end
for i=1:n-1   %计算H
    for j=1:n-1
        if i~=j
            H(i,j)=-bus(i,2)*bus(j,2)*(real(Y(i,j))*sin(bus(i,3)-bus(j,3))-imag(Y(i,j))*cos(bus(i,3)-bus(j,3)));
            N(i,j)=-bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)-bus(j,3))+imag(Y(i,j))*sin(bus(i,3)-bus(j,3)));
            K(i,j)=bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)-bus(j,3))+imag(Y(i,j))*sin(bus(i,3)-bus(j,3)));
            L(i,j)=H(i,j);
        end
    end
end
for i=1:n-1   %计算H
    for j=1:n-1
        if j==i
           H(i,i)=bus(i,5)+imag(Y(i,i))*bus(i,2)^2;
           N(i,i)=-bus(i,4)-real(Y(i,i))*bus(i,2)^2;
           K(i,i)=N(i,i)+2*real(Y(i,i))*bus(i,2)^2;
           L(i,i)=-H(i,i)+2*imag(Y(i,i))*bus(i,2)^2;
        end
    end
end
N=N(1:n-1,1:m);
K=K(1:m,1:n-1);
L=L(1:m,1:m);
J=[H,N;K,L];