将数学规划原理和常规潮流计算相有机结合_形成的带有最优乘子的牛顿算法_最优潮流_matlab

function [J] = jacobian1(nodeNum,PVI,PQI,aij,Y1,e,f,PQPVI) if isempty(PVI)~=1&&isempty(PQI)~=1 H=zeros(nodeNum-1);N=H; M=zeros(size(PQI,1),nodeNum-1);L=M; R=zeros(size(PVI,1),nodeNum-1);S=R; for i=1:nodeNum-1%对H矩阵进行构成 for j=1:nodeNum-1 if PQPVI(i)==PQPVI(j) H(i,j)=H(i,j)-e(PQPVI(i))*Y1(PQPVI(i),PQPVI(i))*cos(aij(PQPVI(i),PQPVI(i)))-f(PQPVI(i))*Y1(PQPVI(i),PQPVI(i))*sin(aij(PQPVI(i),PQPVI(i))); for hi=1:nodeNum H(i,j)= H(i,j)-(Y1(PQPVI(i),hi)*cos(aij(PQPVI(i),hi))*e(hi)-Y1(PQPVI(i),hi)*sin(aij(PQPVI(i),hi))*f(hi)); end else H(i,j)=-e(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*cos(aij(PQPVI(i),PQPVI(j)))-f(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*sin(aij(PQPVI(i),PQPVI(j))); end end end for i=1:nodeNum-1%对N矩阵进行构成 for j=1:nodeNum-1 if PQPVI(i)==PQPVI(j) N(i,j)=e(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*sin(aij(PQPVI(i),PQPVI(j)))-f(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*cos(aij(PQPVI(i),PQPVI(j))); for ni=1:nodeNum N(i,j)=N(i,j)-(Y1(PQPVI(i),ni)*cos(aij(PQPVI(i),ni))*f(hi)+Y1(PQPVI(i),ni)*sin(aij(PQPVI(i),ni))*e(ni)); end else N(i,j)=e(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*sin(aij(PQPVI(i),PQPVI(j)))-f(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*cos(aij(PQPVI(i),PQPVI(j))); end end end for i=1:size(PQI,1)%对M矩阵 for j=1:nodeNum-1 if PQI(i)==PQPVI(j) M(i,j)=e(PQI(i))*Y1(PQI(i),PQI(i))*sin(aij(PQI(i),PQI(i)))-f(PQI(i))*Y1(PQI(i),PQI(i))*cos(aij(PQI(i),PQI(i))); for mi=1:nodeNum M(i,j)=M(i,j)+(Y1(PQI(i),mi)*cos(aij(PQI(i),mi))*f(mi)-Y1(PQI(i),mi)*sin(aij(PQI(i),mi))*e(mi)); end else M(i,j)=e(PQI(i))*Y1(PQI(i),PQPVI(j))*sin(aij(PQI(i),PQPVI(j)))-f(PQI(i))*Y1(PQI(i),PQPVI(j))*cos(aij(PQI(i),PQPVI(j))); end end end for i=1:size(PQI,1)%对L矩阵进行构成 for j=1:nodeNum-1 if PQI(i)==PQPVI(j) L(i,j)=e(PQI(i))*Y1(PQI(i),PQI(i))*cos(aij(PQI(i),PQI(i)))+f(PQI(i))*Y1(PQI(i),PQI(i))*sin(aij(PQI(i),PQI(i))); for li=1:nodeNum L(i,j)=L(i,j)-(Y1(PQI(i),li)*cos(aij(PQI(i),li))*e(li)-Y1(PQI(i),li)*sin(aij(PQI(i),li))*f(li)); end else L(i,j)=e(PQI(i))*Y1(PQI(i),PQPVI(j))*cos(aij(PQI(i),PQPVI(j)))+f(PQI(i))*Y1(PQI(i),PQPVI(j))*sin(aij(PQI(i),PQPVI(j))); end end end for i=1:size(PVI,1)%对R矩阵 for j=1:nodeNum-1 if PVI(i)==PQPVI(j) R(i,j)=-2*e(PVI(i)); end end end for i=1:size(PVI,1)%对S矩阵 for j=1:nodeNum-1 if PVI(i)==PQPVI(j) S(i,j)=-2*f(PVI(i)); end end end J=[H,N;M,L;R,S];%形成雅克比矩阵J elseif isempty(PVI)==1&&isempty(PQI)~=1%假如没有PV节点 H=zeros(nodeNum-1);N=H; M=zeros(size(PQI,1),nodeNum-1);L=M; for i=1:nodeNum-1%对H矩阵进行构成 for j=1:nodeNum-1 if PQPVI(i)==PQPVI(j) H(i,j)=H(i,j)-e(PQPVI(i))*Y1(PQPVI(i),PQPVI(i))*cos(aij(PQPVI(i),PQPVI(i)))-f(PQPVI(i))*Y1(PQPVI(i),PQPVI(i))*sin(aij(PQPVI(i),PQPVI(i))); for hi=1:nodeNum H(i,j)= H(i,j)-(Y1(PQPVI(i),hi)*cos(aij(PQPVI(i),hi))*e(hi)-Y1(PQPVI(i),hi)*sin(aij(PQPVI(i),hi))*f(hi)); end else H(i,j)=-e(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*cos(aij(PQPVI(i),PQPVI(j)))-f(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*sin(aij(PQPVI(i),PQPVI(j))); end end end for i=1:nodeNum-1%对N矩阵进行构成 for j=1:nodeNum-1 if PQPVI(i)==PQPVI(j) N(i,j)=e(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*sin(aij(PQPVI(i),PQPVI(j)))-f(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*cos(aij(PQPVI(i),PQPVI(j))); for ni=1:nodeNum N(i,j)=N(i,j)-(Y1(PQPVI(i),ni)*cos(aij(PQPVI(i),ni))*f(hi)+Y1(PQPVI(i),ni)*sin(aij(PQPVI(i),ni))*e(ni)); end else N(i,j)=e(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*sin(aij(PQPVI(i),PQPVI(j)))-f(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*cos(aij(PQPVI(i),PQPVI(j))); end end end for i=1:size(PQI,1)%对M矩阵 for j=1:nodeNum-1 if PQI(i)==PQPVI(j) M(i,j)=e(PQI(i))*Y1(PQI(i),PQI(i))*sin(aij(PQI(i),PQI(i)))-f(PQI(i))*Y1(PQI(i),PQI(i))*cos(aij(PQI(i),PQI(i))); for mi=1:nodeNum M(i,j)=M(i,j)+(Y1(PQI(i),mi)*cos(aij(PQI(i),mi))*f(mi)-Y1(PQI(i),mi)*sin(aij(PQI(i),mi))*e(mi)); end else M(i,j)=e(PQI(i))*Y1(PQI(i),PQPVI(j))*sin(aij(PQI(i),PQPVI(j)))-f(PQI(i))*Y1(PQI(i),PQPVI(j))*cos(aij(PQI(i),PQPVI(j))); end end end for i=1:size(PQI,1)%对L矩阵进行构成 for j=1:nodeNum-1 if PQI(i)==PQPVI(j) L(i,j)=e(PQI(i))*Y1(PQI(i),PQI(i))*cos(aij(PQI(i),PQI(i)))+f(PQI(i))*Y1(PQI(i),PQI(i))*sin(aij(PQI(i),PQI(i))); for li=1:nodeNum L(i,j)=L(i,j)-(Y1(PQI(i),li)*cos(aij(PQI(i),li))*e(li)-Y1(PQI(i),li)*sin(aij(PQI(i),li))*f(li)); end else L(i,j)=e(PQI(i))*Y1(PQI(i),PQPVI(j))*cos(aij(PQI(i),PQPVI(j)))+f(PQI(i))*Y1(PQI(i),PQPVI(j))*sin(aij(PQI(i),PQPVI(j))); end end end J=[H,N;M,L];%形成雅克比矩阵J else%假如没有PQ节点 H=zeros(nodeNum-1);N=H; R=zeros(size(PVI,1),nodeNum-1);S=R; for i=1:nodeNum-1%对H矩阵进行构成 for j=1:nodeNum-1 if PQPVI(i)==PQPVI(j) H(i,j)=H(i,j)-e(PQPVI(i))*Y1(PQPVI(i),PQPVI(i))*cos(aij(PQPVI(i),PQPVI(i)))-f(PQPVI(i))*Y1(PQPVI(i),PQPVI(i))*sin(aij(PQPVI(i),PQPVI(i))); for hi=1:nodeNum H(i,j)= H(i,j)-(Y1(PQPVI(i),hi)*cos(aij(PQPVI(i),hi))*e(hi)-Y1(PQPVI(i),hi)*sin(aij(PQPVI(i),hi))*f(hi)); end else H(i,j)=-e(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*cos(aij(PQPVI(i),PQPVI(j)))-f(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*sin(aij(PQPVI(i),PQPVI(j))); end end end for i=1:nodeNum-1%对N矩阵进行构成 for j=1:nodeNum-1 if PQPVI(i)==PQPVI(j) N(i,j)=e(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*sin(aij(PQPVI(i),PQPVI(j)))-f(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*cos(aij(PQPVI(i),PQPVI(j))); for ni=1:nodeNum N(i,j)=N(i,j)-(Y1(PQPVI(i),ni)*cos(aij(PQPVI(i),ni))*f(hi)+Y1(PQPVI(i),ni)*sin(aij(PQPVI(i),ni))*e(ni)); end else N(i,j)=e(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*sin(aij(PQPVI(i),PQPVI(j)))-f(PQPVI(i))*Y1(PQPVI(i),PQPVI(j))*cos(aij(PQPVI(i),PQPVI(j))); end end end for i=1:size(PVI,1)%对R矩阵 for j=1:nodeNum-1 if PVI(i)==PQPVI(j) R(i,j)=-2*e(PVI(i)); end end end for i=1:size(PVI,1)%对S矩阵 for j=1:nodeNum-1 if PVI(i)==PQPVI(j) S(i,j)=-2*f(PVI(i)); end end end J=[H,N;R,S];%形成雅克比矩阵J end end