From 27bd5cb77c29b7379c0b91e720aa42ed58d2d4f4 Mon Sep 17 00:00:00 2001
From: facat <dugg@21cn.com>
Date: Mon, 15 Apr 2013 11:14:34 +0800
Subject: [PATCH] =?UTF-8?q?=E4=BF=AE=E6=94=B9=E4=B8=80=E4=BA=9B=E6=B1=82?=
 =?UTF-8?q?=E8=A7=A3=E5=99=A8=E7=9A=84=E5=8F=82=E6=95=B0?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: facat <dugg@21cn.com>
---
 run.m | 19 +++++++++++--------
 1 file changed, 11 insertions(+), 8 deletions(-)

diff --git a/run.m b/run.m
index 7647912..d7cfedb 100644
--- a/run.m
+++ b/run.m
@@ -10,7 +10,7 @@ addpath('.\Powerflow')
 % 电压 相角
 %%
 %% 开始生成量测量
-sigma=0.01;% 标准差
+sigma=0.05;% 标准差
 %% 电压
 %电压幅值
 rVolt=Volt; %幅值
@@ -53,7 +53,7 @@ mPG=rPG.*(normrnd(0,sigma,length(rPG),1)+1);
 mQG=rQG.*(normrnd(0,sigma,length(rQG),1)+1);
 %% 冗余度计算
 stateVarCount=2*length(Volt);
-measurements=length(mVolt)+length(mBranchP)+length(mBranchQ)+length(mPG)+length(mQG);
+measurements=length(mVolt)+length(mBranchI)+length(mBranchP)+length(mBranchQ)+length(mPG)+length(mQG);
 fprintf('冗余度 %f\n',measurements/stateVarCount);
 %% 进入状态估计计算
 % SEVolt=sdpvar(length(Volt),1,'full');
@@ -96,12 +96,15 @@ onlyQG=setdiff(QGi,PDQDi);
 Busnum=length(Volt);
 seOpti=SEOpti();
 seOpti=seOpti.init(Busnum, mVolt,sigma,newwordParameter,zerosInjectionIndex,cmpY,onlyPG,onlyQG,mPG,mQG,Balance,mBranchI,mBranchP,mBranchQ);
-opts = optiset('solver','ipopt');
-opts.maxiter=3500;
-opts.tolrfun=1e-3;
-opts.tolafun=1e-3;
-opts.warnings='warnings';
-opts.display='iter';
+opts = optiset('solver','filtersd');
+opts.maxiter=15500;
+opts.maxtime=3000;
+opts.maxfeval=15000;
+opts.maxnodes=15000;
+opts.tolrfun=1e-5;
+opts.tolafun=1e-5;
+opts.warnings='all';
+opts.display='off';
 x0=[ones(Busnum,1);zeros(Busnum,1)];
 [~,seOpti]=seOpti.equ(x0);
 nlrhs=seOpti.nlrhs();