From e7e390c20aa4404216198a9e6710e395d55bc38b Mon Sep 17 00:00:00 2001 From: facat Date: Sat, 17 Aug 2013 17:50:36 +0800 Subject: [PATCH] =?UTF-8?q?=E5=AE=8C=E6=88=90=E4=BA=86=E5=90=AB=E6=9C=89?= =?UTF-8?q?=E7=94=B5=E5=8E=8B=E5=B9=85=E5=80=BC=EF=BC=8C=E7=BA=BF=E8=B7=AF?= =?UTF-8?q?=E5=8A=9F=E7=8E=87=EF=BC=8C=E5=8F=98=E5=8E=8B=E5=99=A8=E5=8A=9F?= =?UTF-8?q?=E7=8E=87=E7=9A=84=E7=8A=B6=E6=80=81=E4=BC=B0=E8=AE=A1=E7=A8=8B?= =?UTF-8?q?=E5=BA=8F=E3=80=82?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: facat --- run.m | 15 ++++++++------- 1 file changed, 8 insertions(+), 7 deletions(-) diff --git a/run.m b/run.m index 2040b8f..ab9b9fd 100644 --- a/run.m +++ b/run.m @@ -208,17 +208,18 @@ while maxD>1e-3; SEBranchI=BranchI( SEVolt.*exp(1j*SEVAngel),lineI,lineJ,lineR,lineX );%复数支路电流 SEBranchP=BranchP( SEVolt.*exp(1j*SEVAngel),SEBranchI,lineI,lineB2 ); SEBranchQ=BranchQ( SEVolt.*exp(1j*SEVAngel),SEBranchI,lineI,lineB2 ); - SETransP=TransPower( newwordParameter,cmpY,rVolt,rVAngel ); - h=[SEVolt;SEBranchP;SEBranchQ]; - z=[mVolt;mBranchP;mBranchP]; + SETransP=TransPower( newwordParameter,SEVolt,SEVAngel ); + SETransQ=TransReactivePower( newwordParameter,SEVolt,SEVAngel ); + h=[SEVolt;SEBranchP;SEBranchQ;SETransP;SETransQ]; + z=[mVolt;mBranchP;mBranchQ;mTransP;mTransQ]; W=sparse(1:length(h),1:length(h),1/sigma.^2,length(h),length(h)); G=H'*W*H; g=-H'*W*(z-h); % 平衡节点相角恒定; - % G(length(mVolt)+Balance,:)=0; - % G(:,length(mVolt)+Balance)=0; - % G=G+sparse(length(mVolt)+Balance,length(mVolt)+Balance,1,length(mVolt)*2,length(mVolt)*2); - % g(length(mVolt)+Balance)=0; + G(length(mVolt)+Balance,:)=0; + G(:,length(mVolt)+Balance)=0; + G=G+sparse(length(mVolt)+Balance,length(mVolt)+Balance,1,length(mVolt)*2,length(mVolt)*2); + g(length(mVolt)+Balance)=0; dX=G\-g; maxD=max(abs(dX)) % 更新变量