From 77dd6026718f9d9e4e40f8038a5a5dd221bd30d9 Mon Sep 17 00:00:00 2001
From: "dugg@lab-desk" <dugg@lab-desk>
Date: Wed, 5 Nov 2014 12:06:13 +0800
Subject: [PATCH] =?UTF-8?q?=E5=88=A0=E6=8E=89=E4=B8=80=E4=BA=9B=E6=B3=A8?=
 =?UTF-8?q?=E9=87=8A=EF=BC=8C=E5=85=88=E5=90=A7bar=E5=9B=BE=E5=BD=A2?=
 =?UTF-8?q?=E9=9A=90=E6=8E=89=E3=80=82?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
---
 run.m | 34 +++++-----------------------------
 1 file changed, 5 insertions(+), 29 deletions(-)

diff --git a/run.m b/run.m
index d64e157..fc6b79c 100644
--- a/run.m
+++ b/run.m
@@ -78,26 +78,6 @@ while(k<=kmax && maxD> EPS)
     iterPhaseASpotLoadQ=phaseASpotLoadQ+imag(SA);
     iterPhaseBSpotLoadQ=phaseBSpotLoadQ+imag(SB);
     iterPhaseCSpotLoadQ=phaseCSpotLoadQ+imag(SC);
-%    %全部转换为负荷电流
-%    CurpA=conj((iterPhaseASpotLoadP+1j*iterPhaseASpotLoadQ)./VoltpA);
-%    CurpB=conj((iterPhaseBSpotLoadP+1j*iterPhaseBSpotLoadQ)./VoltpB);
-%    CurpC=conj((iterPhaseCSpotLoadP+1j*iterPhaseCSpotLoadQ)./VoltpC);
-%    %转换为序电流
-%    f012=Tp2f*conj([CurpA';CurpB';CurpC']);
-%    %把三序电流分离出来
-%    If0=conj(f012(1,:)');
-%    If1=conj(f012(2,:)');
-%    If2=conj(f012(3,:)');
-%    %试着算一下正序电流
-%     fsY11*V1;
-%    %形成负荷序电流的测量值
-%    mIf0=If0;
-%    mIf1=If1;
-%    mIf1(3)=-mIf1(2);
-%    mIf2=If2;
-%    %计算
-%     fsY11=fsY11+sparse(Balance,Balance,ones(length(Balance),1),busNum,bus
-%     Num);%这里要置0，置1，否则是奇异的
     %%做最小二乘法
     [dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance,busNum, ...
         PQi,PG,QG,QGi,iterPD,iterQD,Vmf1,Vaf1,fsY1amp,fsY1ang,r,c,Vf2,If2,Vf0,If0);%不平衡量
@@ -118,8 +98,6 @@ while(k<=kmax && maxD> EPS)
     If0=conj(f012(1,:)');
     If1=conj(f012(2,:)');
     If2=conj(f012(3,:)');
-    
-    
     If0(Balance)=0;
     If2(Balance)=0;
     %Vf0=fsY0\If0;
@@ -127,8 +105,6 @@ while(k<=kmax && maxD> EPS)
     %Vf2=fsY2\If2;
     Vf2=fsY2Q*(fsY2U\(fsY2L\(fsY2P*(fsY2R\If2))));
     fprintf('迭代时间%f\n',toc);
-    
-    %
 end
 FortiscueToc=toc;
 fprintf('Fortiscue法计算时间 %f\n',FortiscueToc);
@@ -277,7 +253,7 @@ SEIpABC=full(Tf2p*conj([SEIf0';SEIf1(1:3)';SEIf2']));
 %看一下差多少
 VError=(full(abs(VoltpABC))-abs(SEVoltpABC))./abs(VoltpABC)*100;
 VError=reshape(VError,size(VError,1)*size(VError,2),1);
-barPlot( VError, 10,['相对误差%'],['分布密度'],['电压幅值']);
+% barPlot( VError, 10,['相对误差%'],['分布密度'],['电压幅值']);
 % 三相负荷
 rThreeLoad=[ phaseASpotLoadP'+1j*phaseASpotLoadQ';
     phaseBSpotLoadP'+1j*phaseBSpotLoadQ';
@@ -289,7 +265,7 @@ phaseLoadPError=real(rThreeLoad-SEThreeLoad)./real(rThreeLoad)*100;
 phaseLoadQError=imag(rThreeLoad-SEThreeLoad)./imag(rThreeLoad)*100;
 phaseLoadPError=reshape(phaseLoadPError,size(phaseLoadPError,1)*size(phaseLoadPError,2),1);
 phaseLoadQError=reshape(phaseLoadQError,size(phaseLoadQError,1)*size(phaseLoadQError,2),1);
-figure()
-barPlot( phaseLoadPError, 10,['相对误差%'],['分布密度'],['有功负荷误差']);
-figure()
-barPlot( phaseLoadQError, 10,'相对误差%','分布密度','无功负荷误差');
\ No newline at end of file
+% figure()
+% barPlot( phaseLoadPError, 10,['相对误差%'],['分布密度'],['有功负荷误差']);
+% figure()
+% barPlot( phaseLoadQError, 10,'相对误差%','分布密度','无功负荷误差');
\ No newline at end of file