From 3b2edd87e0e30d570e985f85be164b7db0ba5792 Mon Sep 17 00:00:00 2001
From: facat <dugg@21cn.com>
Date: Wed, 6 Mar 2013 22:16:56 +0800
Subject: [PATCH] =?UTF-8?q?=E7=94=A8200PD0=E5=92=8C200QD0=E4=BD=9C?=
 =?UTF-8?q?=E4=B8=BA20=E7=82=B9=E5=8F=A0=E5=8A=A0=E4=BA=86=E5=99=AA?=
 =?UTF-8?q?=E9=9F=B3=E7=9A=84=E6=95=B0=E6=8D=AE?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 200PD0.mat                | Bin 0 -> 672 bytes
 200QD0.mat                | Bin 0 -> 670 bytes
 NormalDistributionCurve.m |  14 ++++++++++++++
 OPF.m                     |  12 +++++++-----
 OPF_Init.asv              |   4 ++--
 OPF_Init.m                |   8 ++++----
 6 files changed, 27 insertions(+), 11 deletions(-)
 create mode 100644 200PD0.mat
 create mode 100644 200QD0.mat
 create mode 100644 NormalDistributionCurve.m

diff --git a/200PD0.mat b/200PD0.mat
new file mode 100644
index 0000000000000000000000000000000000000000..abd9e7b5b5adc77140e5bcc1ceb593da6b750fd4
GIT binary patch
literal 672
zcmeZu4DoSvQZUssQ1EpO(M`+DN!3vZ$Vn_o%P-2cQV4Jk_w>_Ia4t$sEJ;mK$j`G<
z2v1E>@J%dIFfdavG`BJ^ure`LFfuSSRv;NLFnap(GcYhnFflMx%rQM5?J=2Ag#Cm6
z^`k{=THak;aAlk6?P<EnssH9iZaQ`$vDkODuF8o$7o0D?*tEsei)pcvMz2V6PTa3;
zJ8vGyo%O6~)6`VuT{bO?+k14TG_j^0_hzq?Iag^h$G$&tMsg?5=7uVvnfs?Cn)!d<
zdDBNItzgEYnww|txm%<t3koO-PV{y0F`98xT1ZLNXj5y-q_r-7+a~I4nyYcbZ&P=J
z#>u=*?IHrAfx^6+Dp5f`B}-afoMcQ=Sh}P~My35kVv>ql+O~hslb+SZn9QBebWQeT
z{A8nBT}20<c>3irthDVl-6d!Lr}vKj$z^sr`ez^Tt@@RCy@YR$X<NauS5HiJ_w^sS
zwlQ_C{O-ivA%=D{o;9xf_i0tuYVKFh_siDo7S|7Q^$l0<7g&ELdgAJ@k3!yV)_Q$i
z;LOeS`R?;~Ocb}<9v1QDz=`KKGBZ^=KKj4UGc(_CV%yQL>l_4%HXfI;y4M|DZyysG
z@u1=1^MLyfbB@{`{m`3ztE}SBri;(S#9r%ce6swhtVH0fzObtI4(Y#zO+VRqeBUj)
z-^owJH*_^~{#^!t(>=eM*du@6%-7!iPHB10r{9&=dd;g}%fB|&^Qc+nvHe2dw}hms
zCHMQeG9y&#-?tx_K22p?)kU$~%G8?|x%a82?YJj;WYVG9t--J5oqj5}#}|EMT=hU^
l_V$Cl0lqWx=VUZ5>C>*Ac_?W`veLqbE^7Z61KsyU00132Cg}hG

literal 0
HcmV?d00001

diff --git a/200QD0.mat b/200QD0.mat
new file mode 100644
index 0000000000000000000000000000000000000000..e294606a4eef4e798faec8d15a3d50dde2983a18
GIT binary patch
literal 670
zcmeZu4DoSvQZUssQ1EpO(M`+DN!3vZ$Vn_o%P-2cQV4Jk_w>_Ia4t$sEJ;mK$j`G<
z2v1E>@J%dIFfdavG`BJ^urfANFfuSSRv;NLFnap(GcYiSF)=Vy%rQM5>oK`dgzZCp
z<l&raH(DDXE>gKY=j_@3)U97LBu!OvHqSbJ)}*uPl~UnVo=&^8Yiaw|nNHnXtQk>~
z=-C!6sd%;ZU{FM0U~i0L;>1ak8@snWc4x1XIk(edPP%0Bn-fmQSgt%inIxF}I{nw5
znprwCM2~-IJ$3VMt&sbr6IM#9Mw9%Ws5GCw8!*9Vk&4zN*G-;$CsngHHRh<y?wr;j
zG9fp``le&al8!e@M_v5APQ-2Mm+;gyocMQ1&y3F2DU+F(G@ZTrKPKnjaaH4=nh&H(
zHGXFG%9=l&mQXvF-}c|dPgjqbg&e(Z%XR(N##Vm4^LDoyr>D2Q<JJGF^ws%4_s4nN
z4f}lrB|kU>{4F@M{Bc0Zjt3E|Pc`oQ@H^x1sqBs`DkuD0tv5#h+FmYt=fv*SVp5j!
z_qK)G{tLXp5*wfQX>FeMyc%QmuI$=^+tpG*LedK-NZMZG(|4J8fVc5N_w)ECHQ{cy
zF|o;+@7Vo5EPLf2EXU6Dr?qQhCOezCti5^a4whEC<0@yiCu)led;1mqnEF3!*V-Rk
zYR6-lH@^3+5T7P9sl%e@8^5&E)bzkvzrJjCJ<RACcJ%yKzavwAeO$<C?Y8CMDc4(K
zEFn!I$6d75eyZF)zJYuDTlV)`%s+cv_nK$8Qhe_7`u!(`c3k}Z>xtiUd37eUVvXR-
i9RHuLS90zaz8>x}<yIQI$+bNPR~Rt<XD|`kVG95cJ|Y+Z

literal 0
HcmV?d00001

diff --git a/NormalDistributionCurve.m b/NormalDistributionCurve.m
new file mode 100644
index 0000000..1ba077f
--- /dev/null
+++ b/NormalDistributionCurve.m
@@ -0,0 +1,14 @@
+%%
+% 画正态分布曲线
+%%
+
+x=-4:0.1:4;
+y=gaussmf(x,[1 0]);
+plot(x,y)
+hold on
+% line([0 0],[ 0 1])
+% set(gca,'xtick',[],'xticklabel',[])
+% set(gca,'ytick',[],'yticklabel',[])
+axis([-4 4 0 1.2])
+axis normal
+xlabel('gaussmf, P=[2 5]')
\ No newline at end of file
diff --git a/OPF.m b/OPF.m
index f41efec..b59f5ce 100644
--- a/OPF.m
+++ b/OPF.m
@@ -53,10 +53,10 @@ for II=1:1
     %%
     Precision=Precision/1;
     %% 加误差
-    PD0(Loadi)=PD0(Loadi).*(1+normrnd(0,0.5,length(Loadi),1));
-    QD0(Loadi)=QD0(Loadi).*(1+normrnd(0,0.5,length(Loadi),1));
-    load('0PD0.mat');
-    load('0QD0.mat');
+    PD0(Loadi)=PD0(Loadi).*(1+normrnd(0,0.05,length(Loadi),1));
+    QD0(Loadi)=QD0(Loadi).*(1+normrnd(0,0.05,length(Loadi),1));
+    load('200PD0.mat');
+    load('200QD0.mat');
 %     PG0(PGi)=PG0(PGi).*(1+normrnd(0,0.01,length(PGi),1));
 %     QG0(PVi)=QG0(PVi).*(1+normrnd(0,0.01,length(PVi),1));
     %% 读变压器容量
@@ -132,12 +132,14 @@ PD(Loadi)=thesis.MeanPD();
 QD(Loadi)=thesis.MeanQD();
 PG(Balance)=thesis.MeanPG();
 QG(PVi)=thesis.MeanQG();
+thesis.MaxDeviation(RealPG(Balance),RealQG(PVi),RealPD(Loadi),RealQD(Loadi))
+thesis.StatDeviation(RealPG(Balance),RealQG(PVi),RealPD(Loadi),RealQD(Loadi))
 thesis.MaxDeviation(PG0(Balance),QG0(PVi),PD0(Loadi),QD0(Loadi))
 thesis.StatDeviation(PG0(Balance),QG0(PVi),PD0(Loadi),QD0(Loadi))
 % thesis.SqareDeviation(RealPG(Balance(1)),RealQG(PVi(1)),RealPD(Loadi),RealQD(Loadi))
 thesis.SqareDeviation(PG0(Balance(1)),QG0(PVi(1)),PD0(Loadi),QD0(Loadi))
 thesis.PercentOfPass(PG0(Balance(1)),QG0(PVi(1)),PD0(Loadi),QD0(Loadi))
-thesis.MaxBranchDeviation(Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0,Volt0,UAngel0,Volt,UAngel,'E:\算例\柳金Ⅰ926_21671693_2012-09-06\newFIle0.txt',PD0,QD0)
+thesis.MaxBranchDeviation(Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0,Volt0,UAngel0,Volt,UAngel,'E:\算例\柳金Ⅰ926_21671693_2012-09-06\newFIle20.txt',PD0,QD0)
 thesis.StatBranchDeviation(Linei,Linej,Liner,Linex,Lineb,Transfori,Transforj,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0,Volt0,UAngel0,Volt,UAngel)
 toc
 
diff --git a/OPF_Init.asv b/OPF_Init.asv
index c0cfc2d..d4c67dd 100644
--- a/OPF_Init.asv
+++ b/OPF_Init.asv
@@ -20,9 +20,9 @@ wQG=0;
 %randInt=randperm(size(Loadi,1));
 %randPDind=randInt(1:10);
 randPDind=0;
-wPD=1/.5^2*zeros(size(Loadi,1),1);
+wPD=1/.05^2*zeros(size(Loadi,1),1);
 wPD(1:2:end)=0;
-wQD=1/.5^2*ones(size(Loadi,1),1);
+wQD=1/.05^2*ones(size(Loadi,1),1);
 wQD(1:2:end)=0;
 %wD(randPDind)=0;%一些负荷不约束
 %wD(7)=0;
diff --git a/OPF_Init.m b/OPF_Init.m
index 6dfd73d..6ffb2e4 100644
--- a/OPF_Init.m
+++ b/OPF_Init.m
@@ -20,10 +20,10 @@ wQG=0;
 %randInt=randperm(size(Loadi,1));
 %randPDind=randInt(1:10);
 randPDind=0;
-wPD=1/.5^2*ones(size(Loadi,1),1);
-% wPD(1:2:end)=0;
-wQD=1/.5^2*ones(size(Loadi,1),1);
-% wQD(1:2:end)=0;
+wPD=1/.05^2*zeros(size(Loadi,1),1);
+wPD(1:2:end)=0;
+wQD=1/.05^2*zeros(size(Loadi,1),1);
+wQD(1:2:end)=0;
 %wD(randPDind)=0;%一些负荷不约束
 %wD(7)=0;
 % wD(11)=0;