facat
/
threephasese-onlyenquationconstraint
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

准备改成只有一阶导的形式 

Signed-off-by: dmy@lab <dmy@lab.lab>
This commit is contained in:
dmy@lab 2015-05-18 21:50:26 +08:00
parent 82a88a7ec0
commit 224611b051
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
OPF.m
View File

@ -3,10 +3,10 @@ function [JMeasurement,AME_Volt,AME_VAngle,AME_PD,AME_QD,AME_mVolt,AME_mPD,AME_m
tic tic
clc clc
clear clear
lineZ=readLineZ('feeder13\lineParameter.txt'); lineZ=readLineZ('.\..\DistributionNetwork-Power2Current\modified-feeder69\lineParameter.txt');
[ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ... [ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ...
phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY ... phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY ...
cap]=dataRead(lineZ,'feeder13\data1.txt'); cap]=dataRead(lineZ,'.\..\DistributionNetwork-Power2Current\modified-feeder69\data.txt');
% phaseASpotLoadP(phaseASpotLoadP==0)=0.002; % phaseASpotLoadP(phaseASpotLoadP==0)=0.002;
% phaseBSpotLoadP(phaseBSpotLoadP==0)=0.002; % phaseBSpotLoadP(phaseBSpotLoadP==0)=0.002;
% phaseCSpotLoadP(phaseCSpotLoadP==0)=0.002; % phaseCSpotLoadP(phaseCSpotLoadP==0)=0.002;
@ -220,7 +220,7 @@ noPQi3P=zeros(length(noLoadi)*3,1);
noPQi3P(1:3:end)=(noLoadi-1)*3+1; noPQi3P(1:3:end)=(noLoadi-1)*3+1;
noPQi3P(2:3:end)=(noLoadi-1)*3+2; noPQi3P(2:3:end)=(noLoadi-1)*3+2;
noPQi3P(3:3:end)=(noLoadi-1)*3+3; noPQi3P(3:3:end)=(noLoadi-1)*3+3;
noPQi3P=Loadi; % noPQi3P=Loadi;
%量测量 %量测量
sigma=0.03; sigma=0.03;
VoltSigma=(1+normrnd(0,sigma,length(rVoltABCV),1)); VoltSigma=(1+normrnd(0,sigma,length(rVoltABCV),1));
@ -268,8 +268,10 @@ while(maxD>Precision)
%% %%
%% 形成海森阵 %% 形成海森阵
deltdeltF=func_deltdeltF(wVolt,wPD,wQD,ContrlCount,Loadi); deltdeltF=func_deltdeltF(wVolt,wPD,wQD,ContrlCount,Loadi);
% deltdeltF=zeros(size(deltdeltF));
%% 形成ddHy %% 形成ddHy
ddh=func_ddh(Volt,Init_Y,Busnum,Y,UAngel,r,c,Yangle,Loadi,ContrlCount); ddh=func_ddh(Volt,Init_Y,Busnum,Y,UAngel,r,c,Yangle,Loadi,ContrlCount);
% ddh=zeros(size(ddh));
%% 开始构建ddg %% 开始构建ddg
%% 开始构建deltF %% 开始构建deltF
deltF=func_deltF(wVolt,wPD,wQD,mPD3P,PD3P,QD3P,mQD3P,Volt,mVoltABCV,Busnum,Loadi); deltF=func_deltF(wVolt,wPD,wQD,mPD3P,PD3P,QD3P,mQD3P,Volt,mVoltABCV,Busnum,Loadi);