修复电压幅值计算方差

Signed-off-by: dmy@lab <dmy@lab.lab>
添加评价迭代次数的统计。
2015-06-26 11:20:25 +08:00 · 2015-06-04 11:50:57 +08:00 · 2015-06-03 21:45:35 +08:00 · 2015-06-03 21:31:09 +08:00 · 2015-05-21 14:50:06 +08:00 · 2015-05-20 21:01:52 +08:00
195 changed files with 5080 additions and 296 deletions
--- a/AssignXX.m
+++ b/AssignXX.m
@ -1,5 +1,7 @@
-function [deltZ,deltL,deltW,deltU,deltX,deltY]=AssignXX(XX,ContrlCount,RestraintCount,Busnum)
+function [deltZ,deltL,deltW,deltU,deltX,deltY]=AssignXX(XX,ContrlCount,RestraintCount,Busnum,Balance)
 deltX=XX(1:ContrlCount);
+deltX(Balance)=0;
+deltX(Balance+Busnum)=0;
 k1=ContrlCount+2*Busnum;
 deltY=XX(ContrlCount+1:k1);
 k2=k1;
--- a/FormG.m
+++ b/FormG.m
@ -1,8 +1,14 @@
 function Mat_G=FormG(I1r,I1i)
-% Mat_G=[
-%     I1r;
-%     I1i;
-%     ];
+%t1=PG(PVi);
+%GP=t1;%发电机P
+%GP=[4.5 4.5]';
+%%线路
+%发电机Q
+% t1=Volt'*Volt;
+% t2=real(GB).*sin(AngleIJMat)-imag(GB).*cos(AngleIJMat);
+% t3=t1.*t2;
+% t4=sum(t3,2);%发电机Q
+%GQ=t4;
 Mat_G=[
    I1r.^2+I1i.^2;
    ];
--- a/FormGstate1.m
+++ b/FormGstate1.m
@ -0,0 +1,16 @@
+function Mat_G=FormGstate1(I1r,I1i,Loadi)
+%t1=PG(PVi);
+%GP=t1;%发电机P
+%GP=[4.5 4.5]';
+%%线路
+%发电机Q
+% t1=Volt'*Volt;
+% t2=real(GB).*sin(AngleIJMat)-imag(GB).*cos(AngleIJMat);
+% t3=t1.*t2;
+% t4=sum(t3,2);%发电机Q
+%GQ=t4;
+Mat_G=[
+    I1r(Loadi);
+    I1i(Loadi);
+    ];
+end
--- a/FormH.m
+++ b/FormH.m
@ -1,15 +1,6 @@
 function Mat_H=FormH(fsY11,Loadi,V1r,V1i,I1r,I1i,BalI1r,BalI1i,Balance)
-%% 
-%QDcos=textread('D:\Project\×îÐ¡»¯³±Á÷\×îÐ¡³±Á÷ËãÀý\ÏÉº£919PDQDglys.txt');
-%QD(QD~=0)=PD(QD~=0)./tan(QDcos);
-%QD(QD_NON_ZERO_IND)=QD_NON_ZERO;
-%% 
-%PD(Loadi)=QD(Loadi)./tan(acos(0.98));
-% AngleIJ=sparse(r,c,UAngel(r)-UAngel(c)-Angle,Busnum,Busnum);
-% dP=PG-PD-diag(Volt)*Y.*cos(AngleIJ)*Volt;
-% dQ=QG-QD-diag(Volt)*Y.*sin(AngleIJ)*Volt;
-% 
-% Mat_H=[dP;dQ;];
+
+%%%%一下是学姐给的公式
 H=[
    real(fsY11),-imag(fsY11);
    imag(fsY11),real(fsY11);
@ -21,4 +12,5 @@ Mat_H=H*[V1r;V1i];
 Mat_H=Mat_H+sparse(Loadi,1,-I1r,busNum*2,1)+sparse(Loadi+busNum,1,-I1i,busNum*2,1);
 Mat_H(Balance)=Mat_H(Balance)-BalI1r;
 Mat_H(Balance+busNum)=Mat_H(Balance+busNum)-BalI1i;
+
 end
--- a/FormLw.m
+++ b/FormLw.m
@ -1,12 +1,5 @@
-function Lw=FormLw(Loadi,Mat_G,Init_U)
-% VoltU=sparse(1.2*ones(Busnum,1));
-% PDU=sparse(2*ones(length(Loadi),1));
-% QDU=sparse(2*ones(length(Loadi),1));
-% t1=[PDU;
-%     QDU;
-%     VoltU];
-% t2=Mat_G+Init_U-t1;
-% Lw=t2;
-upper=2*sparse(ones(length(Loadi)*1,1));
+function Lw=FormLw(Loadi,Mat_G,Init_U,I1measurement)
+
+upper=abs(I1measurement).^2.*1.1.*sparse(ones(length(Loadi)*1,1));
 Lw=Mat_G+Init_U-upper;
 end
--- a/FormLwstate1.m
+++ b/FormLwstate1.m
@ -0,0 +1,46 @@
+function Lw=FormLwstate1(Loadi,Mat_G,Init_U,I1measurement,dI_F,flag,guessIf1,rIf1,upper)
+
+
+
+% Loadi=[1,4,5,6,8,10];
+
+% lower=-0.2*sparse(ones(length(Loadi)*2,1));
+% upper=ones(length(Loadi)*2,1);
+% upper(pIr)=1.4*Ir(pIr);
+% upper(nIr)=0.6*Ir(nIr);
+% upper(pIi+length(Ir))=1.4*Ii(pIi);
+% upper(nIi+length(Ir))=0.6*Ii(nIi);
+
+% upper(pIr)=abs(real(dI_F(pIr,flag)))+Ir(pIr);
+% upper(nIr)=abs(real(dI_F(nIr,flag)))+Ir(nIr);
+% upper(pIi+length(Ir))=abs(imag(dI_F(pIi,flag)))+Ii(pIi);
+% upper(nIi+length(Ir))=abs(imag(dI_F(nIi,flag)))+Ii(nIi);
+
+
+
+
+
+% if flag~=2
+% upper(pIr)=0.8;
+% upper(nIr)=0.8;
+% upper(pIi+length(Ir))=0.8;
+% upper(nIi+length(Ir))=0.8;
+% end
+
+
+% if real(rIf1(Loadi))>upper( union( pIr,nIr) )
+%     warning('not bounded');
+% end
+% 
+% if imag(rIf1(Loadi))>upper(union( pIr,nIr)+length(Ir))
+%     warning('not bounded');
+% end
+
+%Ì«Ð¡µÄÊý¶¼Òª·Å¿íÒ»Ð©
+% tooSmall=find(abs(Ir)<0.05);
+% upper(tooSmall)=0.8;
+% tooSmall=find(abs(Ii)<0.05);
+% upper(tooSmall+length(Ir))=0.8;
+
+Lw=Mat_G+Init_U-upper;
+end
--- a/FormLz.m
+++ b/FormLz.m
@ -1,14 +1,5 @@
-function Lz=FormLz(Loadi,Mat_G,Init_L)
-% VoltL=sparse(0.8*ones(Busnum,1));
-% PDL=sparse(-2*ones(length(Loadi),1));
-% QDL=sparse(-2*ones(length(Loadi),1));
-% t1=[PDL;
-%     QDL;
-%     VoltL
-%     ];
-% t2=Mat_G-Init_L-t1;
-% Lz=t2;
-lower=0*sparse(ones(length(Loadi)*1,1));
+function Lz=FormLz(Loadi,Mat_G,Init_L,I1measurement)
+
+lower=abs(I1measurement).^2.*0.9.*sparse(ones(length(Loadi)*1,1));
 Lz=Mat_G-Init_L-lower;
-% Lz=-[]
 end
--- a/FormLzstate1.m
+++ b/FormLzstate1.m
@ -0,0 +1,51 @@
+function Lz=FormLzstate1(Loadi,Mat_G,Init_L,I1measurement,dI_F,flag,guessIf1,rIf1,lower)
+
+
+% Loadi=[1,4,5,6,8,10];
+
+
+% lower(pIr)=0.6*Ir(pIr);
+% lower(nIr)=1.4*Ir(nIr);
+% lower(pIi+length(Ir))=0.6*Ii(pIi);
+% lower(nIi+length(Ir))=1.4*Ii(nIi);
+
+% lower(pIr)=-abs(real(dI_F(pIr,flag)))+Ir(pIr);
+% lower(nIr)=-abs(real(dI_F(nIr,flag)))+Ir(nIr);
+% lower(pIi+length(Ir))=-abs(imag(dI_F(pIi,flag)))+Ii(pIi);
+% lower(nIi+length(Ir))=-abs(imag(dI_F(nIi,flag)))+Ii(nIi);
+
+% lower(pIr)=-abs(real(dI_F(pIr,flag)))+Ir(pIr);
+% lower(nIr)=-abs(real(dI_F(nIr,flag)))+Ir(nIr);
+% lower(pIi+length(Ir))=-abs(imag(dI_F(pIi,flag)))+Ii(pIi);
+% lower(nIi+length(Ir))=-abs(imag(dI_F(nIi,flag)))+Ii(nIi);
+
+
+% lower(nIr(1))=-18;
+% lower(nIr(2))=-18;
+% lower(nIr(3))=-18;
+% lower(nIr(4))=-18;
+
+
+% if flag~=2
+% lower(pIr)=-18;
+% lower(nIr)=-18;
+% lower(pIi+length(Ir))=-18;
+% lower(nIi+length(Ir))=-18;
+% end
+
+% if real(rIf1(Loadi))<lower( union( pIr,nIr) )
+%     warning('not bounded');
+% end
+% 
+% if imag(rIf1(Loadi))<lower(union( pIr,nIr)+length(Ir))
+%     warning('not bounded');
+% end
+
+%Ì«Ğ¡µÄÊı¶¼Òª·Å¿íÒ»Ğ©
+% tooSmall=find(abs(Ir)<0.05);
+% lower(tooSmall)=-0.8;
+% tooSmall=find(abs(Ii)<0.05);
+% lower(tooSmall+length(Ir))=-0.8;
+
+Lz=Mat_G-Init_L-lower;
+end
--- a/IPMLoop.m
+++ b/IPMLoop.m
@ -1,61 +1,134 @@
-function [ V1r,V1i,I1r,I1i ] = IPMLoop(measurement,BalI1r,BalI1i,busNum,Loadi,fsY1,Balance,Vref )
+function [ V1r,V1i,I1r,I1i,isConverged,calcuTime,KK ] = IPMLoop(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,BalI1r,BalI1i,busNum,Loadi,fsY1,Balance,Vref,dI_F,flag,guessIf1,rIf1,noLoadi )
+
 %把每个序的循环写在这个函数中。其实也就是内点法循环。
-V1r=1*ones(busNum,1);
+V1r=Vref*ones(busNum,1);
 V1i=0*ones(busNum,1);
-I1r=0.01*ones(length(Loadi),1);%注入电流，相当于放电机电流，与负荷电流负荷相反
-I1i=0.01*ones(length(Loadi),1);%注入电流，相当于放电机电流，与负荷电流负荷相反
+% I1r=0.1*ones(length(Loadi),1);%注入电流，相当于放电机电流，与负荷电流负荷相反
+% I1i=0.1*ones(length(Loadi),1);%注入电流，相当于放电机电流，与负荷电流负荷相反
+I1r=real(guessIf1);
+I1i=imag(guessIf1);
 KK=0;
 plotGap=zeros(1,60);
+
+% ContrlCount=busNum*2+length(Loadi)*2;
+noBoundedLoadi=setdiff(Loadi,noLoadi);
+% noBoundedLoadi=[];
+boundedLoadi=setdiff(Loadi,noBoundedLoadi);
 %初始化
 %状态量为 SEPD SEQD SEVmf1 SEVaf1
-RestraintCount=length(Loadi)*1;
+% RestraintCount=length(Loadi)*2;
+RestraintCount=length(boundedLoadi)*2;
 ContrlCount=busNum*2+length(Loadi)*2;
 CenterA=0.1;
-Init_Z=sparse(ones(RestraintCount,1));
+Init_Z=1*sparse(ones(RestraintCount,1));
 Init_W=sparse(-1*ones(RestraintCount,1));
 Init_L=1*sparse(ones(RestraintCount,1));
 Init_U=1*sparse(ones(RestraintCount,1));
 Init_Y=sparse(2*busNum,1);%等式约束乘子
 Gap=(Init_L'*Init_Z-Init_U'*Init_W);
-% PG=sparse(busNum,1);
-% PG(Balance)=0.1105;
-% QG=sparse(busNum,1);
-% QG(Balance)=0.0984;
-% SEPD(2)=0.1105;
-% SEQD(2)=0.0984;
-while Gap>1e-5 && KK<20
-    KK=KK+1;
+kmax=15;
+isSetBound=0;
+largerBound=0;
+realBound=1;
+
+
+Ir=real(guessIf1(ismember(Loadi,boundedLoadi)) );
+pIr=find(Ir>0);
+nIr=find(Ir<0);
+Ii=imag(guessIf1(ismember(Loadi,boundedLoadi)));
+pIi=find(Ii>0);
+nIi=find(Ii<0);
+lower=ones(length(boundedLoadi)*2,1);
+if realBound==1
+    isSetBound=1;
+    lower(pIr)=(0.7)*Ir(pIr);
+    lower(nIr)=(1.3)*Ir(nIr);
+    lower(pIi+length(Ir))=(0.7)*Ii(pIi);
+    lower(nIi+length(Ir))=(1.3)*Ii(nIi);
+
+% lower(pIr)=-0.1;
+%     lower(nIr)=-0.1;
+%     lower(pIi+length(Ir))=-0.1;
+%     lower(nIi+length(Ir))=-0.1;
+
+end
+
+Ir=real(guessIf1(ismember(Loadi,boundedLoadi)));
+pIr=find(Ir>0);
+nIr=find(Ir<0);
+Ii=imag(guessIf1(ismember(Loadi,boundedLoadi)));
+pIi=find(Ii>0);
+nIi=find(Ii<0);
+upper=ones(length(boundedLoadi)*2,1);
+%     if abs(Gap)<1 || isSetBound==1
+if realBound==1
+    isSetBound=1;
+    upper(pIr)=(1.3)*Ir(pIr);
+    upper(nIr)=(0.7)*Ir(nIr);
+    upper(pIi+length(Ir))=(1.3)*Ii(pIi);
+    upper(nIi+length(Ir))=(0.7)*Ii(nIi);
+% upper(pIr)=0.1;
+%     upper(nIr)=0.1;
+%     upper(pIi+length(Ir))=0.1;
+%     upper(nIi+length(Ir))=0.1;
+end
+tic
+while(abs(Gap)>1e-5)
+    if KK>=kmax
+        break;
+    end
+    plotGap(KK+1)=Gap;
    Init_u=Gap/2/RestraintCount*CenterA;
+    %% 开始计算OPF
+    %% 形成等式约束的雅克比
    deltH=func_deltH(busNum,fsY1,Loadi,Balance);
-    deltG=func_deltG(busNum,Loadi,I1r,I1i);
-%     deltG=0;
+    %% 形成不等式约束的雅克比
+%      deltG=func_deltGstate1(busNum,Loadi,I1r,I1i);
+deltG=func_deltGstate1(busNum,Loadi,boundedLoadi,I1r,I1i);
+    %%
 %     L_1Z=diag(Init_Z./Init_L);
 %     U_1W=diag(Init_W./Init_U);
-    deltdeltF=func_deltdeltF(busNum,fsY1,Loadi);
+    %% 形成海森阵
+    deltdeltF=func_deltdeltF(busNum,fsY1,Loadi,wV1r,wV1i,wI1r,wI1i,V1measurement,V1r,V1i);
 %     deltdeltF=0;
-%     ddh=func_ddh(SEVmf1,Init_Y,busNum,fsY1amp,SEVaf1,r,c,fsY1ang,Loadi,ContrlCount);
+    %% 形成ddHy
+%     ddh=func_ddh(busNum,Loadi,Init_Z,Init_W);
    ddh=0;
-    ddg=func_ddg(busNum,Loadi,Init_Z,Init_W);
-    deltF=func_deltF(measurement,busNum,fsY1,Loadi,V1r,V1i,I1r,I1i);
+    %%  开始构建ddg
+    ddg=0;
+    %%  开始构建deltF
+    deltF=func_deltF(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,busNum,fsY1,Loadi,V1r,V1i,I1r,I1i);
 %     deltF=0;
+    %%
    Luu=Init_U.*Init_W+Init_u*ones(RestraintCount,1);
    Lul=Init_L.*Init_Z-Init_u*ones(RestraintCount,1);
-    Mat_G=FormG(I1r,I1i);
+%     Mat_G=FormGstate1(I1r,I1i);
+    Mat_G=FormGstate1(I1r,I1i,ismember(Loadi,boundedLoadi) );
    Mat_H=FormH(fsY1,Loadi,V1r,V1i,I1r,I1i,BalI1r,BalI1i,Balance);
    Ly=Mat_H;
-    Lz=FormLz(Loadi,Mat_G,Init_L);
-    Lw=FormLw(Loadi,Mat_G,Init_U);
+
+    Lz=FormLzstate1(Loadi,Mat_G,Init_L,I1measurement,dI_F,flag,guessIf1,rIf1,lower);
+    Lw=FormLwstate1(Loadi,Mat_G,Init_U,I1measurement,dI_F,flag,guessIf1,rIf1,upper);
    Lx=FormLx(deltF,deltH,Init_Y,deltG,Init_Z,Init_W);
-%     YY=FormYY(Lul,Lz,Ly,Luu,Lw,Lx);
+    YY=FormYY(Lul,Lz,Ly,Luu,Lw,Lx);
    %% 开始解方程
-    plotGap(KK)=Gap;
-    fprintf('迭代次数 %d Gap %f\n',KK,plotGap(KK));
-    XX=SolveIt(deltF,deltG,Init_L,Init_Z,Init_U,Init_W,deltdeltF,ddh,ddg,deltH,Init_Y,Ly,Lz,ContrlCount,Lw,Lul,Luu,Lx,Balance,busNum,Loadi);
-    [deltZ,deltL,deltW,deltU,deltX,deltY]=AssignXX(XX,ContrlCount,RestraintCount,busNum);
+    XX=SolveIt(deltF,deltG,Init_L,Init_Z,Init_U,Init_W,deltdeltF,ddh,ddg,deltH,Init_Y,Ly,Lz,ContrlCount,Lw,Lul,Luu,RestraintCount,Lx,Balance,busNum);
+    %%取各分量
+    [deltZ,deltL,deltW,deltU,deltX,deltY]=AssignXX(XX,ContrlCount,RestraintCount,busNum,Balance);
    [Init_Z,Init_L,Init_W,Init_U,Init_Y,V1r,V1i,I1r,I1i]=Modification(Init_Z,Init_L,Init_W,Init_U,Init_Y,deltZ,deltL,deltW,deltU,deltX,deltY,V1r,V1i,I1r,I1i,ContrlCount,Balance,busNum,Loadi,Vref);
-    Gap=(Init_L'*Init_Z-Init_U'*Init_W);
+    Gap=(Init_L'*Init_Z-Init_U'*Init_W);+max(abs(deltX));
+%     Gap=max([max(abs(Mat_H)),max(Lz),max(Lw) ]);
+    fprintf('Gap %f :%d\n',full(Gap),KK+1);
+    KK=KK+1;
+end
+calcuTime=0;
+% toc
+calcuTime=toc;
+% f=sum(([real(I1measurement);imag(I1measurement)]-[I1r;I1i]).^2)+sum((real(V1measurement)-V1r(Loadi)).^2)+sum((imag(V1measurement)-V1i(Loadi)).^2);
+if abs(Gap)<1e-5 && KK<kmax
+    isConverged=1;
+else
+    isConverged=0;
 end
-f=sum(([real(measurement);imag(measurement)]-[-I1r;-I1i]).^2);
-
 end

--- a/Modification.m
+++ b/Modification.m
@ -1,22 +1,27 @@
-function [Init_Z,Init_L,Init_W,Init_U,Init_Y,V1r,V1i,I1r,I1i]=Modification(Init_Z,Init_L,Init_W,Init_U,Init_Y,deltZ,deltL,deltW,deltU,deltX,deltY, ... 
-    V1r,V1i,I1r,I1i,ContrlCount,Balance,busNum,Loadi,Vref)
+function [Init_Z,Init_L,Init_W,Init_U,Init_Y,V1r,V1i,I1r,I1i]=Modification(Init_Z,Init_L,Init_W,Init_U,Init_Y,deltZ,deltL,deltW,deltU,deltX,deltY,V1r,V1i,I1r,I1i,ContrlCount,Balance,Busnum,Loadi,Vref)
 AlphaP=FormAlphaP(Init_L,deltL,Init_U,deltU);
+% fprintf('AlphaP %f\n',full(AlphaP));
 AlphaD=FormAlphaD(Init_Z,deltZ,Init_W,deltW);
+% fprintf('AlphaD %f\n',full(AlphaD));
+
 Init_Z=Init_Z+AlphaD*deltZ;
 Init_L=Init_L+AlphaP*deltL;
 Init_W=Init_W+AlphaD*deltW;
 Init_U=Init_U+AlphaP*deltU;
 Init_Y=Init_Y+AlphaD*deltY;
-t=deltX(1:busNum*2);
-V1r=V1r+AlphaP*t(1:busNum);
+%PG(PGi)=PG(PGi)+deltX(1:size(PGi,1));
+% PG(PGi)=PG(PGi)+AlphaP*deltX(1:size(PGi,1));
+V1r=V1r+AlphaP*deltX(1:Busnum);
 V1r(Balance)=Vref;
-V1i=V1i+AlphaP*t(busNum+1:end);
+V1i=V1i+AlphaP*deltX(Busnum+1:2*Busnum);
 V1i(Balance)=0;
-t=deltX(busNum*2+1:ContrlCount);
+%QG(PVi)=QG(PVi)+deltX(size(PGi,1)+1:size(PVi,1)+size(PGi,1) );
+% QG(PVi)=QG(PVi)+AlphaP*deltX(size(PGi,1)+1:size(PVi,1)+size(PGi,1) );
+t=deltX(Busnum*2+1:ContrlCount);
 I1r=I1r+AlphaP*t(1:length(Loadi));
 I1i=I1i+AlphaP*t(length(Loadi)+1:end);
-%balVolt=1;%Volt(Balance);
-%Volt=Volt+AlphaP*t(1:Busnum);
-%Volt(Balance)=balVolt;
-%UAngel=UAngel+AlphaP*t(Busnum+1:2*Busnum);
+
+%Volt=Volt+AlphaP*t(2:2:2*Busnum);ÔÝÊ±¸ÄÒ»ÏÂ20111227
+%UAngel=UAngel+AlphaP*t(1:2:2*Busnum);ÔÝÊ±¸ÄÒ»ÏÂ20111227
+
 end
--- a/SolveIt.m
+++ b/SolveIt.m
@ -1,31 +1,26 @@
-function XX=SolveIt(deltF,deltG,Init_L,Init_Z,Init_U,Init_W,deltdeltF,ddh,ddg,deltH,Init_Y,Ly,Lz,ContrlCount,Lw,Lul,Luu,Lx,Balance,Busnum,Loadi)
+function XX=SolveIt(deltF,deltG,Init_L,Init_Z,Init_U,Init_W,deltdeltF,ddh,ddg,deltH,Init_Y,Ly,Lz,ContrlCount,Lw,Lul,Luu,ConstraintCount,Lx,Balance,Busnum)
 LxComa=FormLxComa(deltF,deltG,deltH,Init_L,Luu,Lul,Init_Z,Init_Y,Lz,Init_U,Init_W,Lw,Lx);
-H=-deltdeltF+ddh+ddg;
-t1=diag(Init_Z./Init_L-Init_W./Init_U);
+H=-deltdeltF+ddh+ddg;%+ddg*(Init_Z'+Init_W');
+t1=diag(Init_L.\Init_Z-Init_U.\Init_W);
 t2=-deltG*( t1 )*deltG';
 aa=[
    (H+t2),deltH;
-    deltH',zeros(length(Init_Y));
+    deltH',sparse(length(Init_Y),length(Init_Y));
    ];
 yy=[LxComa;-Ly];
-%% 平衡节点电压实部

+% t=size(PVi,1)+size(PGi,1);%电压不变
+% yy(t+Balance)=0;
 aa(Balance,:)=0;
 aa(:,Balance)=0;
 aa=aa+sparse(Balance,Balance,ones(length(Balance),1),ContrlCount+2*Busnum,ContrlCount+2*Busnum);
-deltG(Balance,:)=0;
-%% 平衡节点电压虚部
-t=Busnum*1;
+
+%暂时改一下
+t=Busnum;%相角不变
 aa(t+Balance,:)=0;
 aa(:,t+Balance)=0;
-aa=aa+sparse(t+Balance,t+Balance,ones(length(Balance),1),ContrlCount+2*Busnum,ContrlCount+2*Busnum);
-deltG(t+Balance,:)=0;
-%%
+aa(t+Balance,t+Balance)=1;
 dxdy=aa\yy;
-%% KLU
-%spy(aa)
-%dxdy = klu(aa,'\',full(yy));
-%%
 dX=dxdy(1:ContrlCount);
 dY=dxdy(ContrlCount+1:ContrlCount+2*Busnum);
 dL=Lz+deltG'*dX;
--- a/barPlot.m
+++ b/barPlot.m
@ -0,0 +1,13 @@
+function barPlot( data, split_number,xWord,yWord,figTitle)
+ymin=min(data);
+ymax=max(data);
+% split_number=10;
+x=linspace(ymin,ymax,split_number); %将最大最小区间分成split_number个等分点(19等分),然后分别计算各个区间的个数
+yy=hist(data,x); %计算各个区间的个数
+% yy=yy/(sum(yy)*(ymax-ymin)/split_number); %计算各个区间的个数,除以总面积，总面积计算的方式为：所有小分割的面积的和即:
+bar(x,yy) %画出概率密度分布图
+xlabel(xWord);
+ylabel(yWord);
+title(figTitle);
+end
+
--- a/calIWi.m
+++ b/calIWi.m
@ -0,0 +1,12 @@
+function [ output_args ] = calIWi( Vr,Vi,sigmaP,sigmaQ ,P,Q)
+dIrdP=(Vi./(Vr.^2+Vi.^2)).^2.*(sigmaP.^2);
+dIrdQ=(Vr./(Vr.^2+Vi.^2)).^2.*(sigmaQ.^2);
+dIrdVr=-(Q.*(Vr.^2+Vi.^2)-2*(Q.*Vr-P.*Vi).*Vr)./(Vr.^2+Vi.^2);
+dIrdVi=-(P.*(Vr.^2+Vi.^2)-2*(Q.*Vr-P.*Vi).*Vi)./(Vr.^2+Vi.^2);
+% dIrdVr=(dIrdVr).^2*0.1;
+% dIrdVi=(dIrdVi).^2*0.1;
+dIrdVr=zeros(size(dIrdVr));
+dIrdVi=zeros(size(dIrdVi));
+output_args=max(abs([dIrdP,dIrdQ,dIrdVr,dIrdVi]),[],2);
+end
+
--- a/calIWr.m
+++ b/calIWr.m
@ -0,0 +1,13 @@
+function [ output_args ] = calIWr(Vr,Vi,sigmaP,sigmaQ,P,Q )
+%计算电流实部的方差
+dIrdP=(Vr./(Vr.^2+Vi.^2)).^2.*(sigmaP.^2);
+dIrdQ=(Vi./(Vr.^2+Vi.^2)).^2.*(sigmaQ.^2);
+dIrdVr=(P.*(Vr.^2+Vi.^2)-2*(P.*Vr+Q.*Vi).*Vr)./(Vr.^2+Vi.^2);
+% dIrdVr=(dIrdVr).^2.*0.1;
+dIrdVi=(Q.*(Vr.^2+Vi.^2)-2*(P.*Vr+Q.*Vi).*Vi)./(Vr.^2+Vi.^2);
+% dIrdVi=(dIrdVi).^2.*0.1;
+dIrdVr=zeros(size(dIrdVr));
+dIrdVi=zeros(size(dIrdVi));
+output_args=max(abs([dIrdP,dIrdQ,dIrdVr,dIrdVi]),[],2);
+end
+
--- a/feeder104/Cap
+++ b/feeder104/Cap
--- a/feeder104/Distributed
+++ b/feeder104/Distributed
--- a/feeder104/IEEE
+++ b/feeder104/IEEE
--- a/feeder104/Line
+++ b/feeder104/Line
--- a/feeder104/Regulator
+++ b/feeder104/Regulator
--- a/feeder104/Spot
+++ b/feeder104/Spot
--- a/feeder104/Transformer
+++ b/feeder104/Transformer
--- a/feeder104/UG
+++ b/feeder104/UG
--- a/feeder104/config.xls
+++ b/feeder104/config.xls
--- a/feeder104/data
+++ b/feeder104/data
@ -0,0 +1,13 @@
+650
+0
+602	632	633	500
+601	650	632	2000
+601	632	671	2000
+601	671	680	1000
+0
+0
+633 160 120 120 110 90 90
+671 385 385 385 220 220 220
+0
+
+
--- a/feeder104/data
+++ b/feeder104/data
@ -0,0 +1,16 @@
+650
+0
+600	632	645	500
+602	632	633	500
+601	633	634	0
+600	645	646	300
+601	650	632	2000
+602	684	652	800
+601	632	671	2000
+600	671	684	300
+602	671	680	1000
+600	671	692	0
+602	684	611	300
+601	692	675	500
+0
+
--- a/feeder104/data.xls
+++ b/feeder104/data.xls
--- a/feeder104/data1.txt
+++ b/feeder104/data1.txt
@ -0,0 +1,190 @@
+650 4.16
+0
+600	632	645	50
+602	632	633	50
+601	633	634	60
+600	645	646	30
+601	650	632	200
+602	684	652	80
+601	632	671	200
+600	671	684	30
+600	671	692	30
+602	684	611	30
+601	692	675	50
+600	2632	2645	50
+602	2632	2633	50
+601	2633	2634	60
+600	2645	2646	30
+601	650	2632	200
+602	2684	2652	80
+601	2632	2671	200
+600	2671	2684	30
+600	2671	2692	30
+602	2684	2611	30
+601	2692	2675	50
+600	3632	3645	50
+602	3632	3633	50
+601	3633	3634	60
+600	3645	3646	30
+601	650	3632	200
+602	3684	3652	80
+601	3632	3671	200
+600	3671	3684	30
+600	3671	3692	30
+602	3684	3611	30
+601	3692	3675	50
+600	4632	4645	50
+602	4632	4633	50
+601	4633	4634	60
+600	4645	4646	30
+601	650	4632	200
+602	4684	4652	80
+601	4632	4671	200
+600	4671	4684	30
+600	4671	4692	30
+602	4684	4611	30
+601	4692	4675	50
+600	5632	5645	50
+602	5632	5633	50
+601	5633	5634	60
+600	5645	5646	30
+601	650	5632	200
+602	5684	5652	80
+601	5632	5671	200
+600	5671	5684	30
+600	5671	5692	30
+602	5684	5611	30
+601	5692	5675	50
+600	6632	6645	50
+602	6632	6633	50
+601	6633	6634	60
+600	6645	6646	30
+601	650	6632	200
+602	6684	6652	80
+601	6632	6671	200
+600	6671	6684	30
+600	6671	6692	30
+602	6684	6611	30
+601	6692	6675	50
+600	7632	7645	50
+602	7632	7633	50
+601	7633	7634	60
+600	7645	7646	30
+601	650	7632	200
+602	7684	7652	80
+601	7632	7671	200
+600	7671	7684	30
+600	7671	7692	30
+602	7684	7611	30
+601	7692	7675	50
+600	8632	8645	50
+602	8632	8633	50
+601	8633	8634	60
+600	8645	8646	30
+601	650	8632	200
+602	8684	8652	80
+601	8632	8671	200
+600	8671	8684	30
+600	8671	8692	30
+602	8684	8611	30
+601	8692	8675	50
+0
+692 0 0 0
+0
+634	120	90	120	90	160	110
+645	40	30	50	40	60	50
+646	100	50	70	40	60	40
+652	20	10	60	30	40	40
+671	385	220	385	220	385	220
+675	68	60	290	212	285	190
+692	260	200	40	40	70	40
+611	80	20	60	40	30	20
+2634	120	90	120	90	160	110
+2645	40	30	50	40	60	50
+2646	100	50	70	40	60	40
+2652	20	10	60	30	40	40
+2671	385	220	385	220	385	220
+2675	68	60	290	212	285	190
+2692	260	200	40	40	70	40
+2611	80	20	60	40	30	20
+3634	120	90	120	90	160	110
+3645	40	30	50	40	60	50
+3646	100	50	70	40	60	40
+3652	20	10	60	30	40	40
+3671	385	220	385	220	385	220
+3675	68	60	290	212	285	190
+3692	260	200	40	40	70	40
+3611	80	20	60	40	30	20
+4634	120	90	120	90	160	110
+4645	40	30	50	40	60	50
+4646	100	50	70	40	60	40
+4652	20	10	60	30	40	40
+4671	385	220	385	220	385	220
+4675	68	60	290	212	285	190
+4692	260	200	40	40	70	40
+4611	80	20	60	40	30	20
+5634	120	90	120	90	160	110
+5645	40	30	50	40	60	50
+5646	100	50	70	40	60	40
+5652	20	10	60	30	40	40
+5671	385	220	385	220	385	220
+5675	68	60	290	212	285	190
+5692	260	200	40	40	70	40
+5611	80	20	60	40	30	20
+6634	120	90	120	90	160	110
+6645	40	30	50	40	60	50
+6646	100	50	70	40	60	40
+6652	20	10	60	30	40	40
+6671	385	220	385	220	385	220
+6675	68	60	290	212	285	190
+6692	260	200	40	40	70	40
+6611	80	20	60	40	30	20
+7634	120	90	120	90	160	110
+7645	40	30	50	40	60	50
+7646	100	50	70	40	60	40
+7652	20	10	60	30	40	40
+7671	385	220	385	220	385	220
+7675	68	60	290	212	285	190
+7692	260	200	40	40	70	40
+7611	80	20	60	40	30	20
+8634	120	90	120	90	160	110
+8645	40	30	50	40	60	50
+8646	100	50	70	40	60	40
+8652	20	10	60	30	40	40
+8671	385	220	385	220	385	220
+8675	68	60	290	212	285	190
+8692	260	200	40	40	70	40
+8611	80	20	60	40	30	20
+0
+600	632	645	500
+602	632	633	500
+601	633	634	600
+600	645	646	300
+601	650	632	2000
+602	684	652	800
+601	632	671	2000
+600	671	684	300
+602	671	680	1000
+600	671	692	300
+602	684	611	300
+601	692	675	500
+0
+632 0 0 0
+0
+634 160 110 120 90 120 90
+645 123 102 170 125 133 110
+646 172 150 230 132 0 211 160
+652 128 86 30 10 80 70
+671 385 220 385 220 385 220
+675 485 190 68 60 290 212
+692 142 140 153 121 170 151
+611 156 142 185 60 170 80
+0
+634 160 110 120 90 120 90
+645 123 102 170 125 133 110
+646 17 15 23 13 21 16
+652 128 86 30 10 80 70
+675 85 10 68 60 20 12
+611 52 42 85 60 70 40
+680 12 4 10 6 7 4
+
--- a/feeder104/data1三相不平衡.txt
+++ b/feeder104/data1三相不平衡.txt
@ -0,0 +1,14 @@
+650 4.16
+0
+601	650	632	2000
+602	632	631	1800
+602	631	630	1700
+0
+632 0 0 0
+0
+630 90 70 80 40 30 10
+631 40 10 60 30 100 80
+632 12 12 120 120 11 11
+0
+
+
--- a/feeder104/data2
+++ b/feeder104/data2
@ -0,0 +1,13 @@
+650
+0
+602	632	633	500
+601	650	632	2000
+601	632	671	2000
+601	671	680	1000
+0
+0
+633 160 120 120 110 90 90
+671 385 385 385 220 220 220
+0
+
+
--- a/feeder104/data2.txt
+++ b/feeder104/data2.txt
@ -0,0 +1,12 @@
+650 4.16
+0
+602	632	633	500
+601	650	632	2000
+601	632	671	2000
+601	671	680	1000
+0
+0
+671 385 385 385 220 220 220
+0
+
+
--- a/feeder104/feeder13.zip
+++ b/feeder104/feeder13.zip
--- a/feeder104/ie13.TCW
+++ b/feeder104/ie13.TCW
--- a/feeder104/lineParameter
+++ b/feeder104/lineParameter
@ -0,0 +1,27 @@
+601
+0.3465  0.1560  0.1580
+0.1560   0.3375   0.1535
+0.1580   0.1535   0.3414
+0
+1.0179   0.5017   0.4236
+0.5017   1.0478   0.3849
+0.4236   0.3849    1.0348
+0
+6.2998   -1.9958   -1.2595
+-1.9958   5.9597   -0.7417
+-1.2595   -0.7417   5.6386
+0
+602
+0.7526  0.1580   0.1560
+0.1580  0.7475   0.1535
+0.1560  0.1535   0.7436
+0
+1.1814   0.4236   0.5017
+0.4236   1.1983   0.3849
+0.5017   0.3849   1.2112
+0
+5.6990   -1.0817   -1.6905
+-1.0817   5.1795   -0.6588
+-1.6905  -0.6588   5.4246
+0
+
--- a/feeder104/lineParameter.txt
+++ b/feeder104/lineParameter.txt
@ -0,0 +1,40 @@
+600
+0.3465  0.1560  0.1580
+0.1560   0.3375   0.1535
+0.1580   0.1535   0.3414
+0
+1.0179   0.5017   0.4236
+0.5017   1.0478   0.3849
+0.4236   0.3849    1.0348
+0
+6.2998   -1.9958   -1.2595
+-1.9958   5.9597   -0.7417
+-1.2595   -0.7417   5.6386
+0
+601
+0.3465  0.1560  0.1580
+0.1560   0.3375   0.1535
+0.1580   0.1535   0.3414
+0
+1.0179   0.5017   0.4236
+0.5017   1.0478   0.3849
+0.4236   0.3849    1.0348
+0
+6.2998   -1.9958   -1.2595
+-1.9958   5.9597   -0.7417
+-1.2595   -0.7417   5.6386
+0
+602
+0.7526  0.1580   0.1560
+0.1580  0.7475   0.1535
+0.1560  0.1535   0.7436
+0
+1.1814   0.4236   0.5017
+0.4236   1.1983   0.3849
+0.5017   0.3849   1.2112
+0
+5.6990   -1.0817   -1.6905
+-1.0817   5.1795   -0.6588
+-1.6905  -0.6588   5.4246
+0
+
--- a/feeder104/note.txt
+++ b/feeder104/note.txt
@ -0,0 +1,3 @@
+line data
+XFM-1 -1
+Switch -2
--- a/feeder13/data
+++ b/feeder13/data
@ -1,16 +1,16 @@
 650
 0
-603	632	645	500
+600	632	645	500
 602	632	633	500
-1	633	634	0
-603	645	646	300
+601	633	634	0
+600	645	646	300
 601	650	632	2000
-607	684	652	800
+602	684	652	800
 601	632	671	2000
-604	671	684	300
-601	671	680	1000
-2	671	692	0
-605	684	611	300
-606	692	675	500
+600	671	684	300
+602	671	680	1000
+600	671	692	0
+602	684	611	300
+601	692	675	500
 0

--- a/feeder13/data1.txt
+++ b/feeder13/data1.txt
@ -1,12 +1,57 @@
 650 4.16
 0
+600	632	645	50
+602	632	633	50
+601	633	634	60
+600	645	646	30
+601	650	632	200
+602	684	652	80
+601	632	671	200
+600	671	684	30
+600	671	692	30
+602	684	611	30
+601	692	675	50
+0
+692 0 0 0
+0
+634	120	90	120	90	160	110
+645	40	30	50	40	60	50
+646	100	50	70	40	60	40
+652	20	10	60	30	40	40
+671	385	220	385	220	385	220
+675	68	60	290	212	285	190
+692	260	200	40	40	70	40
+611	80	20	60	40	30	20
+0
+600	632	645	500
+602	632	633	500
+601	633	634	600
+600	645	646	300
 601	650	632	2000
-602	632	631	2000
+602	684	652	800
+601	632	671	2000
+600	671	684	300
+602	671	680	1000
+600	671	692	300
+602	684	611	300
+601	692	675	500
 0
-632 30 0 0
+632 0 0 0
 0
-632 120 120 120 110 90 90
-631 90 100 110 80 90 90
+634 160 110 120 90 120 90
+645 123 102 170 125 133 110
+646 172 150 230 132 0 211 160
+652 128 86 30 10 80 70
+671 385 220 385 220 385 220
+675 485 190 68 60 290 212
+692 142 140 153 121 170 151
+611 156 142 185 60 170 80
 0
-
+634 160 110 120 90 120 90
+645 123 102 170 125 133 110
+646 17 15 23 13 21 16
+652 128 86 30 10 80 70
+675 85 10 68 60 20 12
+611 52 42 85 60 70 40
+680 12 4 10 6 7 4

--- a/feeder13/data1三相不平衡.txt
+++ b/feeder13/data1三相不平衡.txt
@ -0,0 +1,14 @@
+650 4.16
+0
+601	650	632	2000
+602	632	631	1800
+602	631	630	1700
+0
+632 0 0 0
+0
+630 90 70 80 40 30 10
+631 40 10 60 30 100 80
+632 12 12 120 120 11 11
+0
+
+
--- a/feeder26/Cap
+++ b/feeder26/Cap
--- a/feeder26/Distributed
+++ b/feeder26/Distributed
--- a/feeder26/IEEE
+++ b/feeder26/IEEE
--- a/feeder26/Line
+++ b/feeder26/Line
--- a/feeder26/Regulator
+++ b/feeder26/Regulator
--- a/feeder26/Spot
+++ b/feeder26/Spot
--- a/feeder26/Transformer
+++ b/feeder26/Transformer
--- a/feeder26/UG
+++ b/feeder26/UG
--- a/feeder26/config.xls
+++ b/feeder26/config.xls
--- a/feeder26/data
+++ b/feeder26/data
@ -0,0 +1,13 @@
+650
+0
+602	632	633	500
+601	650	632	2000
+601	632	671	2000
+601	671	680	1000
+0
+0
+633 160 120 120 110 90 90
+671 385 385 385 220 220 220
+0
+
+
--- a/feeder26/data
+++ b/feeder26/data
@ -0,0 +1,16 @@
+650
+0
+600	632	645	500
+602	632	633	500
+601	633	634	0
+600	645	646	300
+601	650	632	2000
+602	684	652	800
+601	632	671	2000
+600	671	684	300
+602	671	680	1000
+600	671	692	0
+602	684	611	300
+601	692	675	500
+0
+
--- a/feeder26/data.xls
+++ b/feeder26/data.xls
--- a/feeder26/data1.txt
+++ b/feeder26/data1.txt
@ -0,0 +1,76 @@
+650 4.16
+0
+600	632	645	50
+602	632	633	50
+601	633	634	60
+600	645	646	30
+601	650	632	200
+602	684	652	80
+601	632	671	200
+600	671	684	30
+600	671	692	30
+602	684	611	30
+601	692	675	50
+600	2632	2645	50
+602	2632	2633	50
+601	2633	2634	60
+600	2645	2646	30
+601	650	2632	200
+602	2684	2652	80
+601	2632	2671	200
+600	2671	2684	30
+600	2671	2692	30
+602	2684	2611	30
+601	2692	2675	50
+0
+692 0 0 0
+0
+634	120	90	120	90	160	110
+645	40	30	50	40	60	50
+646	100	50	70	40	60	40
+652	20	10	60	30	40	40
+671	385	220	385	220	385	220
+675	68	60	290	212	285	190
+692	260	200	40	40	70	40
+611	80	20	60	40	30	20
+2634	120	90	120	90	160	110
+2645	40	30	50	40	60	50
+2646	100	50	70	40	60	40
+2652	20	10	60	30	40	40
+2671	385	220	385	220	385	220
+2675	68	60	290	212	285	190
+2692	260	200	40	40	70	40
+2611	80	20	60	40	30	20
+0
+600	632	645	500
+602	632	633	500
+601	633	634	600
+600	645	646	300
+601	650	632	2000
+602	684	652	800
+601	632	671	2000
+600	671	684	300
+602	671	680	1000
+600	671	692	300
+602	684	611	300
+601	692	675	500
+0
+632 0 0 0
+0
+634 160 110 120 90 120 90
+645 123 102 170 125 133 110
+646 172 150 230 132 0 211 160
+652 128 86 30 10 80 70
+671 385 220 385 220 385 220
+675 485 190 68 60 290 212
+692 142 140 153 121 170 151
+611 156 142 185 60 170 80
+0
+634 160 110 120 90 120 90
+645 123 102 170 125 133 110
+646 17 15 23 13 21 16
+652 128 86 30 10 80 70
+675 85 10 68 60 20 12
+611 52 42 85 60 70 40
+680 12 4 10 6 7 4
+
--- a/feeder26/data1三相不平衡.txt
+++ b/feeder26/data1三相不平衡.txt
@ -0,0 +1,14 @@
+650 4.16
+0
+601	650	632	2000
+602	632	631	1800
+602	631	630	1700
+0
+632 0 0 0
+0
+630 90 70 80 40 30 10
+631 40 10 60 30 100 80
+632 12 12 120 120 11 11
+0
+
+
--- a/feeder26/data2
+++ b/feeder26/data2
@ -0,0 +1,13 @@
+650
+0
+602	632	633	500
+601	650	632	2000
+601	632	671	2000
+601	671	680	1000
+0
+0
+633 160 120 120 110 90 90
+671 385 385 385 220 220 220
+0
+
+
--- a/feeder26/data2.txt
+++ b/feeder26/data2.txt
@ -0,0 +1,12 @@
+650 4.16
+0
+602	632	633	500
+601	650	632	2000
+601	632	671	2000
+601	671	680	1000
+0
+0
+671 385 385 385 220 220 220
+0
+
+
--- a/feeder26/feeder13.zip
+++ b/feeder26/feeder13.zip
--- a/feeder26/ie13.TCW
+++ b/feeder26/ie13.TCW
--- a/feeder26/lineParameter
+++ b/feeder26/lineParameter
@ -0,0 +1,27 @@
+601
+0.3465  0.1560  0.1580
+0.1560   0.3375   0.1535
+0.1580   0.1535   0.3414
+0
+1.0179   0.5017   0.4236
+0.5017   1.0478   0.3849
+0.4236   0.3849    1.0348
+0
+6.2998   -1.9958   -1.2595
+-1.9958   5.9597   -0.7417
+-1.2595   -0.7417   5.6386
+0
+602
+0.7526  0.1580   0.1560
+0.1580  0.7475   0.1535
+0.1560  0.1535   0.7436
+0
+1.1814   0.4236   0.5017
+0.4236   1.1983   0.3849
+0.5017   0.3849   1.2112
+0
+5.6990   -1.0817   -1.6905
+-1.0817   5.1795   -0.6588
+-1.6905  -0.6588   5.4246
+0
+
--- a/feeder26/lineParameter.txt
+++ b/feeder26/lineParameter.txt
@ -0,0 +1,40 @@
+600
+0.3465  0.1560  0.1580
+0.1560   0.3375   0.1535
+0.1580   0.1535   0.3414
+0
+1.0179   0.5017   0.4236
+0.5017   1.0478   0.3849
+0.4236   0.3849    1.0348
+0
+6.2998   -1.9958   -1.2595
+-1.9958   5.9597   -0.7417
+-1.2595   -0.7417   5.6386
+0
+601
+0.3465  0.1560  0.1580
+0.1560   0.3375   0.1535
+0.1580   0.1535   0.3414
+0
+1.0179   0.5017   0.4236
+0.5017   1.0478   0.3849
+0.4236   0.3849    1.0348
+0
+6.2998   -1.9958   -1.2595
+-1.9958   5.9597   -0.7417
+-1.2595   -0.7417   5.6386
+0
+602
+0.7526  0.1580   0.1560
+0.1580  0.7475   0.1535
+0.1560  0.1535   0.7436
+0
+1.1814   0.4236   0.5017
+0.4236   1.1983   0.3849
+0.5017   0.3849   1.2112
+0
+5.6990   -1.0817   -1.6905
+-1.0817   5.1795   -0.6588
+-1.6905  -0.6588   5.4246
+0
+
--- a/feeder26/note.txt
+++ b/feeder26/note.txt
@ -0,0 +1,3 @@
+line data
+XFM-1 -1
+Switch -2
--- a/feeder39/Cap
+++ b/feeder39/Cap
--- a/feeder39/Distributed
+++ b/feeder39/Distributed
--- a/feeder39/IEEE
+++ b/feeder39/IEEE
--- a/feeder39/Line
+++ b/feeder39/Line
--- a/feeder39/Regulator
+++ b/feeder39/Regulator
--- a/feeder39/Spot
+++ b/feeder39/Spot
--- a/feeder39/Transformer
+++ b/feeder39/Transformer
--- a/feeder39/UG
+++ b/feeder39/UG
--- a/feeder39/config.xls
+++ b/feeder39/config.xls
--- a/feeder39/data
+++ b/feeder39/data
@ -0,0 +1,13 @@
+650
+0
+602	632	633	500
+601	650	632	2000
+601	632	671	2000
+601	671	680	1000
+0
+0
+633 160 120 120 110 90 90
+671 385 385 385 220 220 220
+0
+
+
--- a/feeder39/data
+++ b/feeder39/data
@ -0,0 +1,16 @@
+650
+0
+600	632	645	500
+602	632	633	500
+601	633	634	0
+600	645	646	300
+601	650	632	2000
+602	684	652	800
+601	632	671	2000
+600	671	684	300
+602	671	680	1000
+600	671	692	0
+602	684	611	300
+601	692	675	500
+0
+
--- a/feeder39/data.xls
+++ b/feeder39/data.xls
--- a/feeder39/data1.txt
+++ b/feeder39/data1.txt
@ -0,0 +1,95 @@
+650 4.16
+0
+600	632	645	50
+602	632	633	50
+601	633	634	60
+600	645	646	30
+601	650	632	200
+602	684	652	80
+601	632	671	200
+600	671	684	30
+600	671	692	30
+602	684	611	30
+601	692	675	50
+600	2632	2645	50
+602	2632	2633	50
+601	2633	2634	60
+600	2645	2646	30
+601	650	2632	200
+602	2684	2652	80
+601	2632	2671	200
+600	2671	2684	30
+600	2671	2692	30
+602	2684	2611	30
+601	2692	2675	50
+600	3632	3645	50
+602	3632	3633	50
+601	3633	3634	60
+600	3645	3646	30
+601	650	3632	200
+602	3684	3652	80
+601	3632	3671	200
+600	3671	3684	30
+600	3671	3692	30
+602	3684	3611	30
+601	3692	3675	50
+0
+692 0 0 0
+0
+634	120	90	120	90	160	110
+645	40	30	50	40	60	50
+646	100	50	70	40	60	40
+652	20	10	60	30	40	40
+671	385	220	385	220	385	220
+675	68	60	290	212	285	190
+692	260	200	40	40	70	40
+611	80	20	60	40	30	20
+2634	120	90	120	90	160	110
+2645	40	30	50	40	60	50
+2646	100	50	70	40	60	40
+2652	20	10	60	30	40	40
+2671	385	220	385	220	385	220
+2675	68	60	290	212	285	190
+2692	260	200	40	40	70	40
+2611	80	20	60	40	30	20
+3634	120	90	120	90	160	110
+3645	40	30	50	40	60	50
+3646	100	50	70	40	60	40
+3652	20	10	60	30	40	40
+3671	385	220	385	220	385	220
+3675	68	60	290	212	285	190
+3692	260	200	40	40	70	40
+3611	80	20	60	40	30	20
+0
+600	632	645	500
+602	632	633	500
+601	633	634	600
+600	645	646	300
+601	650	632	2000
+602	684	652	800
+601	632	671	2000
+600	671	684	300
+602	671	680	1000
+600	671	692	300
+602	684	611	300
+601	692	675	500
+0
+632 0 0 0
+0
+634 160 110 120 90 120 90
+645 123 102 170 125 133 110
+646 172 150 230 132 0 211 160
+652 128 86 30 10 80 70
+671 385 220 385 220 385 220
+675 485 190 68 60 290 212
+692 142 140 153 121 170 151
+611 156 142 185 60 170 80
+0
+634 160 110 120 90 120 90
+645 123 102 170 125 133 110
+646 17 15 23 13 21 16
+652 128 86 30 10 80 70
+675 85 10 68 60 20 12
+611 52 42 85 60 70 40
+680 12 4 10 6 7 4
+
--- a/feeder39/data1三相不平衡.txt
+++ b/feeder39/data1三相不平衡.txt
@ -0,0 +1,14 @@
+650 4.16
+0
+601	650	632	2000
+602	632	631	1800
+602	631	630	1700
+0
+632 0 0 0
+0
+630 90 70 80 40 30 10
+631 40 10 60 30 100 80
+632 12 12 120 120 11 11
+0
+
+
--- a/feeder39/data2
+++ b/feeder39/data2
@ -0,0 +1,13 @@
+650
+0
+602	632	633	500
+601	650	632	2000
+601	632	671	2000
+601	671	680	1000
+0
+0
+633 160 120 120 110 90 90
+671 385 385 385 220 220 220
+0
+
+
--- a/feeder39/data2.txt
+++ b/feeder39/data2.txt
@ -0,0 +1,12 @@
+650 4.16
+0
+602	632	633	500
+601	650	632	2000
+601	632	671	2000
+601	671	680	1000
+0
+0
+671 385 385 385 220 220 220
+0
+
+
--- a/feeder39/feeder13.zip
+++ b/feeder39/feeder13.zip
--- a/feeder39/ie13.TCW
+++ b/feeder39/ie13.TCW
--- a/feeder39/lineParameter
+++ b/feeder39/lineParameter
@ -0,0 +1,27 @@
+601
+0.3465  0.1560  0.1580
+0.1560   0.3375   0.1535
+0.1580   0.1535   0.3414
+0
+1.0179   0.5017   0.4236
+0.5017   1.0478   0.3849
+0.4236   0.3849    1.0348
+0
+6.2998   -1.9958   -1.2595
+-1.9958   5.9597   -0.7417
+-1.2595   -0.7417   5.6386
+0
+602
+0.7526  0.1580   0.1560
+0.1580  0.7475   0.1535
+0.1560  0.1535   0.7436
+0
+1.1814   0.4236   0.5017
+0.4236   1.1983   0.3849
+0.5017   0.3849   1.2112
+0
+5.6990   -1.0817   -1.6905
+-1.0817   5.1795   -0.6588
+-1.6905  -0.6588   5.4246
+0
+
--- a/feeder39/lineParameter.txt
+++ b/feeder39/lineParameter.txt
@ -0,0 +1,40 @@
+600
+0.3465  0.1560  0.1580
+0.1560   0.3375   0.1535
+0.1580   0.1535   0.3414
+0
+1.0179   0.5017   0.4236
+0.5017   1.0478   0.3849
+0.4236   0.3849    1.0348
+0
+6.2998   -1.9958   -1.2595
+-1.9958   5.9597   -0.7417
+-1.2595   -0.7417   5.6386
+0
+601
+0.3465  0.1560  0.1580
+0.1560   0.3375   0.1535
+0.1580   0.1535   0.3414
+0
+1.0179   0.5017   0.4236
+0.5017   1.0478   0.3849
+0.4236   0.3849    1.0348
+0
+6.2998   -1.9958   -1.2595
+-1.9958   5.9597   -0.7417
+-1.2595   -0.7417   5.6386
+0
+602
+0.7526  0.1580   0.1560
+0.1580  0.7475   0.1535
+0.1560  0.1535   0.7436
+0
+1.1814   0.4236   0.5017
+0.4236   1.1983   0.3849
+0.5017   0.3849   1.2112
+0
+5.6990   -1.0817   -1.6905
+-1.0817   5.1795   -0.6588
+-1.6905  -0.6588   5.4246
+0
+
--- a/feeder39/note.txt
+++ b/feeder39/note.txt
@ -0,0 +1,3 @@
+line data
+XFM-1 -1
+Switch -2
--- a/feeder52/Cap
+++ b/feeder52/Cap
--- a/feeder52/Distributed
+++ b/feeder52/Distributed
--- a/feeder52/IEEE
+++ b/feeder52/IEEE
--- a/feeder52/Line
+++ b/feeder52/Line
--- a/feeder52/Regulator
+++ b/feeder52/Regulator
--- a/feeder52/Spot
+++ b/feeder52/Spot
--- a/feeder52/Transformer
+++ b/feeder52/Transformer
--- a/feeder52/UG
+++ b/feeder52/UG
--- a/feeder52/config.xls
+++ b/feeder52/config.xls
--- a/feeder52/data
+++ b/feeder52/data
@ -0,0 +1,13 @@
+650
+0
+602	632	633	500
+601	650	632	2000
+601	632	671	2000
+601	671	680	1000
+0
+0
+633 160 120 120 110 90 90
+671 385 385 385 220 220 220
+0
+
+
--- a/feeder52/data
+++ b/feeder52/data
@ -0,0 +1,16 @@
+650
+0
+600	632	645	500
+602	632	633	500
+601	633	634	0
+600	645	646	300
+601	650	632	2000
+602	684	652	800
+601	632	671	2000
+600	671	684	300
+602	671	680	1000
+600	671	692	0
+602	684	611	300
+601	692	675	500
+0
+
--- a/feeder52/data.xls
+++ b/feeder52/data.xls
--- a/feeder52/data1.txt
+++ b/feeder52/data1.txt
@ -0,0 +1,114 @@
+650 4.16
+0
+600	632	645	50
+602	632	633	50
+601	633	634	60
+600	645	646	30
+601	650	632	200
+602	684	652	80
+601	632	671	200
+600	671	684	30
+600	671	692	30
+602	684	611	30
+601	692	675	50
+600	2632	2645	50
+602	2632	2633	50
+601	2633	2634	60
+600	2645	2646	30
+601	650	2632	200
+602	2684	2652	80
+601	2632	2671	200
+600	2671	2684	30
+600	2671	2692	30
+602	2684	2611	30
+601	2692	2675	50
+600	3632	3645	50
+602	3632	3633	50
+601	3633	3634	60
+600	3645	3646	30
+601	650	3632	200
+602	3684	3652	80
+601	3632	3671	200
+600	3671	3684	30
+600	3671	3692	30
+602	3684	3611	30
+601	3692	3675	50
+600	4632	4645	50
+602	4632	4633	50
+601	4633	4634	60
+600	4645	4646	30
+601	650	4632	200
+602	4684	4652	80
+601	4632	4671	200
+600	4671	4684	30
+600	4671	4692	30
+602	4684	4611	30
+601	4692	4675	50
+0
+692 0 0 0
+0
+634	120	90	120	90	160	110
+645	40	30	50	40	60	50
+646	100	50	70	40	60	40
+652	20	10	60	30	40	40
+671	385	220	385	220	385	220
+675	68	60	290	212	285	190
+692	260	200	40	40	70	40
+611	80	20	60	40	30	20
+2634	120	90	120	90	160	110
+2645	40	30	50	40	60	50
+2646	100	50	70	40	60	40
+2652	20	10	60	30	40	40
+2671	385	220	385	220	385	220
+2675	68	60	290	212	285	190
+2692	260	200	40	40	70	40
+2611	80	20	60	40	30	20
+3634	120	90	120	90	160	110
+3645	40	30	50	40	60	50
+3646	100	50	70	40	60	40
+3652	20	10	60	30	40	40
+3671	385	220	385	220	385	220
+3675	68	60	290	212	285	190
+3692	260	200	40	40	70	40
+3611	80	20	60	40	30	20
+4634	120	90	120	90	160	110
+4645	40	30	50	40	60	50
+4646	100	50	70	40	60	40
+4652	20	10	60	30	40	40
+4671	385	220	385	220	385	220
+4675	68	60	290	212	285	190
+4692	260	200	40	40	70	40
+4611	80	20	60	40	30	20
+0
+600	632	645	500
+602	632	633	500
+601	633	634	600
+600	645	646	300
+601	650	632	2000
+602	684	652	800
+601	632	671	2000
+600	671	684	300
+602	671	680	1000
+600	671	692	300
+602	684	611	300
+601	692	675	500
+0
+632 0 0 0
+0
+634 160 110 120 90 120 90
+645 123 102 170 125 133 110
+646 172 150 230 132 0 211 160
+652 128 86 30 10 80 70
+671 385 220 385 220 385 220
+675 485 190 68 60 290 212
+692 142 140 153 121 170 151
+611 156 142 185 60 170 80
+0
+634 160 110 120 90 120 90
+645 123 102 170 125 133 110
+646 17 15 23 13 21 16
+652 128 86 30 10 80 70
+675 85 10 68 60 20 12
+611 52 42 85 60 70 40
+680 12 4 10 6 7 4
+
--- a/feeder52/data1三相不平衡.txt
+++ b/feeder52/data1三相不平衡.txt
@ -0,0 +1,14 @@
+650 4.16
+0
+601	650	632	2000
+602	632	631	1800
+602	631	630	1700
+0
+632 0 0 0
+0
+630 90 70 80 40 30 10
+631 40 10 60 30 100 80
+632 12 12 120 120 11 11
+0
+
+
--- a/feeder52/data2
+++ b/feeder52/data2
@ -0,0 +1,13 @@
+650
+0
+602	632	633	500
+601	650	632	2000
+601	632	671	2000
+601	671	680	1000
+0
+0
+633 160 120 120 110 90 90
+671 385 385 385 220 220 220
+0
+
+
--- a/feeder52/data2.txt
+++ b/feeder52/data2.txt
@ -0,0 +1,12 @@
+650 4.16
+0
+602	632	633	500
+601	650	632	2000
+601	632	671	2000
+601	671	680	1000
+0
+0
+671 385 385 385 220 220 220
+0
+
+
--- a/feeder52/feeder13.zip
+++ b/feeder52/feeder13.zip
--- a/feeder52/ie13.TCW
+++ b/feeder52/ie13.TCW
--- a/feeder52/lineParameter
+++ b/feeder52/lineParameter
@ -0,0 +1,27 @@
+601
+0.3465  0.1560  0.1580
+0.1560   0.3375   0.1535
+0.1580   0.1535   0.3414
+0
+1.0179   0.5017   0.4236
+0.5017   1.0478   0.3849
+0.4236   0.3849    1.0348
+0
+6.2998   -1.9958   -1.2595
+-1.9958   5.9597   -0.7417
+-1.2595   -0.7417   5.6386
+0
+602
+0.7526  0.1580   0.1560
+0.1580  0.7475   0.1535
+0.1560  0.1535   0.7436
+0
+1.1814   0.4236   0.5017
+0.4236   1.1983   0.3849
+0.5017   0.3849   1.2112
+0
+5.6990   -1.0817   -1.6905
+-1.0817   5.1795   -0.6588
+-1.6905  -0.6588   5.4246
+0
+
--- a/feeder52/lineParameter.txt
+++ b/feeder52/lineParameter.txt
@ -0,0 +1,40 @@
+600
+0.3465  0.1560  0.1580
+0.1560   0.3375   0.1535
+0.1580   0.1535   0.3414
+0
+1.0179   0.5017   0.4236
+0.5017   1.0478   0.3849
+0.4236   0.3849    1.0348
+0
+6.2998   -1.9958   -1.2595
+-1.9958   5.9597   -0.7417
+-1.2595   -0.7417   5.6386
+0
+601
+0.3465  0.1560  0.1580
+0.1560   0.3375   0.1535
+0.1580   0.1535   0.3414
+0
+1.0179   0.5017   0.4236
+0.5017   1.0478   0.3849
+0.4236   0.3849    1.0348
+0
+6.2998   -1.9958   -1.2595
+-1.9958   5.9597   -0.7417
+-1.2595   -0.7417   5.6386
+0
+602
+0.7526  0.1580   0.1560
+0.1580  0.7475   0.1535
+0.1560  0.1535   0.7436
+0
+1.1814   0.4236   0.5017
+0.4236   1.1983   0.3849
+0.5017   0.3849   1.2112
+0
+5.6990   -1.0817   -1.6905
+-1.0817   5.1795   -0.6588
+-1.6905  -0.6588   5.4246
+0
+
--- a/Show More
+++ b/Show More
Author	SHA1	Message	Date
dmy@lab	116d31e395	修复电压幅值计算方差 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-06-26 11:20:25 +08:00
dmy@lab	a740b4e4a6	添加评价迭代次数的统计。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-06-04 11:50:57 +08:00
dmy@lab	5f1a375fa7	修复一个小bug Signed-off-by: dmy@lab <dmy@lab.lab>	2015-06-03 21:45:35 +08:00
dmy@lab	010dc10f6c	用了Popoviciu's inequality估计电压幅值方程，结果和之前的没什么变化。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-06-03 21:31:09 +08:00
dmy@lab	fddb3316bd	加三相不平衡度计算 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-05-21 14:50:06 +08:00
dmy@lab	6b78384380	加入138算例，加入一些并行代码。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-05-20 21:01:52 +08:00
dmy@lab	f7d3e5747c	R/X=8 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-05-16 22:31:18 +08:00
dmy@lab	9385ef2269	增加修改的69节点算例 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-05-16 15:49:16 +08:00
dmy@lab	fa1a29c86b	改用误差传递的方差确定方式 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-05-08 17:36:59 +08:00
dmy@lab	3dd10280c7	修改为只只有不存在实时数据的时候才加到不等式中，收敛速度快很多。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-05-05 22:35:55 +08:00
dmy@lab	eb7782ead1	修复一个收敛判据的bug，也不知道是不是bug，哈哈 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-05-04 21:08:14 +08:00
dmy@lab	932716ceb1	加到104节点，删掉一些没用的注释 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-05-04 20:55:55 +08:00
dmy@lab	b3892896ec	这个是100%最好的参数 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-30 14:27:39 +08:00
dmy@lab	b57e14a499	1.加大系统算例 2.这个是0%最好的参数。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-30 14:24:41 +08:00
dmy@lab	04e95a448e	加计算时间的统计 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-26 14:39:53 +08:00
dmy@lab	286caf83a4	加一个小的理论验证代码片段 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-07 10:39:04 +08:00
dmy@lab	e7f32f1be9	考虑量测配置 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-06 16:03:00 +08:00
dmy@lab	562b209ec8	1.改用新算例。 2.考虑负荷曲线不准确情况。 3.修改统计中的大循环，修复统计错误的问题。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-06 11:51:26 +08:00
dmy@lab	a7141734c0	修复一个大错字的千年错误！！！ Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-02 17:38:43 +08:00
dmy@lab	9b8769071b	AME里面用的应该是真实值，搞错了。已修复。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-02 17:25:48 +08:00
dmy@lab	fa03721b18	之前搞错了，其实可以用量测量做负荷曲线。之前是把三序顺序搞错了。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-02 17:15:39 +08:00
dmy@lab	b6a1f8058c	不能用测量值当负荷曲线，不然上下界不准。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-02 16:42:51 +08:00
dmy@lab	0ad8514b6a	编制了新的节点文件。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-02 15:55:54 +08:00
dmy@lab	4543616833	换一个更简单的计算上下界的方法。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-01 22:08:07 +08:00
dmy@lab	92c16fac04	修复了计算AME的一个错误。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-01 16:58:55 +08:00
dmy@lab	0e591ca2f9	采用新的电压幅值和相角Absolute Mean Error Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-01 16:56:45 +08:00
dmy@lab	28c5270e87	加入了确定不等式上下界的代码。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-01 16:12:37 +08:00
dmy@lab	69ba28f2a9	按我硕士论文中的推导，用了新的电压幅值权重计算。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-01 10:35:11 +08:00
dmy@lab	361909bcde	如果用量测值确定不等式上下界容易出现不收敛的问题。 Signed-off-by: dmy@lab <dmy@lab.lab>	2015-04-01 10:06:11 +08:00
dugg@lab-desk	37bbf6fc86	三相功率太平衡了，不收敛。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-12-10 15:23:20 +08:00
dugg@lab-desk	7d05bcee3a	改了迭代的初值 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-12-10 15:14:11 +08:00
dugg@lab-desk	0b8015e9d0	修复了目标函数二次导数中一个系数的bug。不过这个系数的取值还需要考虑一下。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-12-04 16:33:34 +08:00
dugg@lab-desk	1c2b53d2b5	修复了目标函数求导的错误。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-12-04 11:19:08 +08:00
dugg@lab-desk	5357f5edcf	小修改了一下三序测试 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-12-03 22:05:02 +08:00
dugg@lab-desk	3f0363fadd	1.加一个文件做不平衡度的测试 2.果然，如果三相比较平衡的话，零序和负序就会是一个比较小的数，上下界取得不好就不容易收敛。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-12-02 10:41:34 +08:00
dugg@lab-desk	cf4c3d14c5	还在调试 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-11-17 22:22:39 +08:00
dugg@lab-desk	fb4b9cd9a0	小修改 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-11-11 13:56:12 +08:00
dugg@lab-desk	8210c96601	先把目标函数注释掉了。奇怪，没有用真实平衡节点电流居然这么好收敛？！！ Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-11-11 11:51:07 +08:00
dugg@lab-desk	00a4cfad44	有点问题 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-11-11 11:37:18 +08:00
dugg@lab-desk	3db4877bb7	改了一下目标函数，又不收敛了。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-11-10 22:33:39 +08:00
dugg@lab-desk	d411972751	如果虚拟电压方差太小，就会不容易收敛。要限制方差太小。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-11-10 10:21:23 +08:00
facat@ipso	fb1bde0e22	用虚拟电压做量测量居然不收敛了。 Signed-off-by: facat@ipso <facat@ipso>	2014-11-10 00:19:51 +08:00
facat@ipso	0f4d416ca0	把测试文件的代码改好了。 Signed-off-by: facat@ipso <facat@ipso>	2014-11-09 23:28:36 +08:00
facat@ipso	3a1655112a	加入一个test.m的文件做测试用。 Signed-off-by: facat@ipso <facat@ipso>	2014-11-09 23:04:41 +08:00
dugg@lab-desk	25b73ebad2	加了一些量测评价的代码。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-11-09 22:25:47 +08:00
dugg@lab-desk	92469a3004	换了一种计算序量测量的方式。用最大误差作为量测误差。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-11-09 21:52:19 +08:00
dugg@lab-desk	77dd602671	删掉一些注释，先吧bar图形隐掉。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-11-05 12:06:13 +08:00
dugg@lab-desk	3875907756	修改了数据文件，不要电容了。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-11-04 22:14:13 +08:00
dugg@lab-desk	0c9c2e2eb7	1.添加了画柱状图的函数 2.加了结果展示代码 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-11-03 11:07:25 +08:00
dugg@lab-desk	971d2e3455	修复了最后算原始电压时的bug，例如SEIf0=I0r+1jI0r;写成了SEIf0=I0r+1jI0i; Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-30 11:45:43 +08:00
dugg@lab-desk	ebc7882194	删掉了注释掉的牛顿法求解方法 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-29 21:58:02 +08:00
dugg@lab-desk	4c2f16c8c7	先不考虑电容 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-29 21:56:38 +08:00
dugg@lab-desk	7409723aae	添加把结果转换会三相值 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-29 21:26:39 +08:00
dugg@lab-desk	098fe0e7c1	l<=I<=u 的约束方式调通了。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-29 21:09:38 +08:00
dugg@lab-desk	c1efa7e6a8	修复了目标函数一次导数的一个bug。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-29 17:45:58 +08:00
dugg@lab-desk	d81d1fd26e	改一种电流约束的形式，还没改完，正在进行中……	2014-10-28 22:43:23 +08:00
dugg@lab-desk	420ac03a92	改一种电流约束的形式，还没改完，正在进行中…… Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-19 20:52:21 +08:00
dugg@lab-desk	401b698c37	把方程加到目标函数中。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-19 17:02:55 +08:00
dugg@lab-desk	a3b1c5a30e	如果不置0置1的话，潮流方程又不平衡了，还是该回来。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-19 15:58:55 +08:00
dugg@lab-desk	4544a6d009	由于是直角坐标下写的，平衡节点不要置0置1，否则有问题。现在好了，只把电压放目标函数中，电压误差就为0了。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-19 15:43:08 +08:00
dugg@lab-desk	3c173fde9e	准备把不等式先去掉。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-19 15:15:03 +08:00
dugg@lab-desk	058671a01b	把电压放目标函数中，电流移除目标函数。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-19 14:50:04 +08:00
dugg@lab-desk	ec67f3a775	加了3序。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-19 12:08:01 +08:00
dugg@lab-desk	607d118eca	验证了一下，好像是对的。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-19 11:52:15 +08:00
dugg@lab-desk	17412d8ec5	把原来的OPF重新加入这里面，收敛了，但总觉得有地方不太对。 Signed-off-by: dugg@lab-desk <dugg@lab-desk>	2014-10-19 11:42:44 +08:00