facat
/
ipmthreephasese
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-04-05 20:22:38 +08:00
parent e26a2f9ef0
commit 98502be311
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
OPF.m
View File

@ -225,7 +225,8 @@ wQD=1./(abs(mQD3P*sigma).^2);
%% 没有量测量 %% 没有量测量
% noLoadi=[1,5,6,10,11]; % noLoadi=[1,5,6,10,11];
noLoadi=Loadi; noLoadi=[1,8,9,10,12];
% noLoadi=Loadi;
% noLoadi=[1,11]; % noLoadi=[1,11];
% noLoadi=[]; % noLoadi=[];
noPQi3P=zeros(length(noLoadi)*3,1); noPQi3P=zeros(length(noLoadi)*3,1);
@ -234,9 +235,9 @@ noPQi3P(2:3:end)=(noLoadi-1)*3+2;
noPQi3P(3:3:end)=(noLoadi-1)*3+3; noPQi3P(3:3:end)=(noLoadi-1)*3+3;
% noLoadi=noPQi3P; % noLoadi=noPQi3P;
wVolt(setdiff(1:length(wVolt),Loadi))=0;%只有负荷处才有电压量测。 wVolt(setdiff(1:length(wVolt),Loadi))=0;%只有负荷处才有电压量测。
wVolt( noLoadi)=0; wVolt( noPQi3P)=0;
wPD(ismember( Loadi,noLoadi))=0; wPD(ismember( Loadi,noPQi3P))=0;
wQD(ismember( Loadi,noLoadi))=0; wQD(ismember( Loadi,noPQi3P))=0;
%% %%
% RestraintCount=size(Loadi,1)*2+length(rVoltABCV); %约束条件数 % RestraintCount=size(Loadi,1)*2+length(rVoltABCV); %约束条件数
RestraintCount=size(Loadi,1)*2; %约束条件数 RestraintCount=size(Loadi,1)*2; %约束条件数