From 2dd7ffb1837162caaa534f9a28a20257343bc00c Mon Sep 17 00:00:00 2001 From: "dmy@lab" Date: Thu, 19 Mar 2015 13:40:41 +0800 Subject: [PATCH] =?UTF-8?q?=E9=80=9A=E8=BF=87=E6=BD=AE=E6=B5=81=E6=96=B9?= =?UTF-8?q?=E7=A8=8B=E6=8A=8APD,QD=E5=8F=8D=E6=B1=82=E5=87=BA=E6=9D=A5?= =?UTF-8?q?=E3=80=82?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: dmy@lab --- Run_YALMIP2.m | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) diff --git a/Run_YALMIP2.m b/Run_YALMIP2.m index 7486365..d5b4310 100644 --- a/Run_YALMIP2.m +++ b/Run_YALMIP2.m @@ -112,8 +112,11 @@ for I=1:1 SEVolt=x(1:length(Volt)); SEVAngel=x(length(Volt)+1:end); fprintf('×î´óÆ«²î\n') - PD=x(1:length(PDi)); - QD=x(length(PDi)+1:length(PDi)+length(QDi)); + AngleIJ=sparse(r,c,SEVAngel(r)-SEVAngel(c)-Angle,Busnum,Busnum); + zeroP=diag(SEVolt)*Y.*cos(AngleIJ)*SEVolt; + zeroQ=diag(SEVolt)*Y.*sin(AngleIJ)*SEVolt; + PD=-zeroP(PDi); + QD=-zeroQ(QDi); maxDev1=MaxDeviation(sigma,mVolt,SEVolt,rVAngel,SEVAngel,rPD(PDi),rQD(QDi),PD,QD) fprintf('ͳ¼ÆÆ«²î\n') statDev1=StatDeviation(sigma,rVolt,SEVolt,rVAngel,SEVAngel,rPD(PDi),rQD(QDi),PD,QD)