From 0f95f9c04b514dd972738de7bae887d8ab19a2cb Mon Sep 17 00:00:00 2001 From: facat Date: Wed, 27 Mar 2013 15:55:34 +0800 Subject: [PATCH] =?UTF-8?q?=E5=8A=A0=E5=85=A5=E4=BA=86=E7=94=BB=E7=94=B5?= =?UTF-8?q?=E5=8E=8B=E5=88=86=E5=B8=83=E7=9A=84=E4=BB=A3=E7=A0=81?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: facat --- OPF.m | 3 ++- OPF_Init.asv | 4 ++-- OPF_Init.m | 8 ++++---- case1V.mat | Bin 0 -> 1274 bytes case2V.mat | Bin 0 -> 1276 bytes case3V.mat | Bin 0 -> 1250 bytes caseM.mat | Bin 0 -> 1053 bytes caseR.mat | Bin 0 -> 1016 bytes plotVolt.asv | 26 ++++++++++++++++++++++++++ plotVolt.m | 45 +++++++++++++++++++++++++++++++++++++++++++++ 10 files changed, 79 insertions(+), 7 deletions(-) create mode 100644 case1V.mat create mode 100644 case2V.mat create mode 100644 case3V.mat create mode 100644 caseM.mat create mode 100644 caseR.mat create mode 100644 plotVolt.asv create mode 100644 plotVolt.m diff --git a/OPF.m b/OPF.m index df3f0d6..2daf42d 100644 --- a/OPF.m +++ b/OPF.m @@ -7,7 +7,8 @@ clear thesis=ForThesis(1,45); for II=1:1 [kmax,Precision,UAngel,Volt,Busnum,PVi,PVu,Balance,Y,Angle,P0,Q0,r,c,GB,Linei,Linej,Transfori,Transforj,GenU,GenL,GenC,PG,QG,PD,QD,CenterA,PGi,PVQU,PVQL,Liner,Linex,Lineb,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0]= ... - pf('E:\算例\柳金Ⅰ926_21671693_2012-09-06\newFIle4 - 计算用.txt'); + pf('E:/算例/柳金Ⅰ926_21671693_2012-09-06/newFIle20.txt'); +% pf('E:\算例\柳金Ⅰ926_21671693_2012-09-06\newFIle16.txt'); %pf('D:\Project\青秀降损项目\最小化潮流\最小潮流算例\原始\津头站津视922（3-1）_0.5_120%.txt'); %pf('D:\Project\最小化潮流\最小潮流算例\仙海919.txt'); %pf('c:/file31.txt'); diff --git a/OPF_Init.asv b/OPF_Init.asv index d4c67dd..453a5f4 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/.05^2*zeros(size(Loadi,1),1); +wPD=1/.05^2*ones(size(Loadi,1),1); wPD(1:2:end)=0; -wQD=1/.05^2*ones(size(Loadi,1),1); +wQD=1/.05^2*zeros(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 6ffb2e4..4ccea57 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/.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; +wPD=1/.05^2*ones(size(Loadi,1),1); +% wPD(1:2:end)=0; +wQD=1/.05^2*ones(size(Loadi,1),1); +% wQD(1:2:end)=0; %wD(randPDind)=0;%一些负荷不约束 %wD(7)=0; % wD(11)=0; diff --git a/case1V.mat b/case1V.mat new file mode 100644 index 0000000000000000000000000000000000000000..886f9f76ed65fb034856a4aca1c068205c2949fa GIT binary patch literal 1274 zcmeZu4DoSvQZUssQ1EpO(M`+DN!3vZ$Vn_o%P-2cQV4Jk_w>_Ia4t$sEJ;mK$j`G< z@XgOt@J%dIFfvszu(UF=v@$eMFfuSSRv;NLFnap(GcYg|u`n=H%z1mcJNQbd2>XZm z4M{pH7Idgc@TV?)x|QqnY;JDub*E>SmR7F&z2=6@Pfond=5oJx?Nj6<+=9zN-EI|9_Y5uUgs%&flz+(V&*G@_5L1*6YvS zTYutu!N6VO$ZvMQ?LcYLmxByuhs6w9)+VsrJ|Gd-`Zj@S*#_n~W@7_(?hP$-H*l|M z;N5*tDvq7+MuYB+Mq3L8@oxv>3Ru)1G#c+{(B9E#UD2Rk!N{2Zfzj?qlRU=(d8Pz; zjzs$o1%W072`2@SLkcoV4m|7*OzaN@JQRf%DM%$Lh&^&(Yk$ab#FT>o>ayxnH z^AU}86DDc;-~JP(sd{VS!o?3KPVqV6Za;l_z^4rlGdE~2I;N@g(r=lH=f%T?RxX{n zTbITQy!zU#{I|?+Y+d2B~}HfAE)(NAn_A;;(=Gp^yXNx&b2q0WUB99`^*5}qUG!c? z+i?C-$(~b3s(y8wPruP&-ZG7Mx8VFv?Vss__C48B{|u*0;#-&Ou-;zvuHW0|$`_XZ z`PQiYzEh(zK)-jB%bS(!TA40!Mf+Hr+ZjUS-OO3dHB8l zmo9hp#FoF}`F7P{+tW#UpC(5o9W4=3{hELBR^78;CHWp%t>fiww_U&g`W;vl`6_R2 zhuxW3E?3rlKU68jzpGs_-Lhr#`Q*Cv|5I-BK5pR$PdU@e88)X~ zBlD}?dHsJqlZx+cKKdukYtt&XuYOkofrkbv0}3au*}7Nyl~9{4Uhzj<&lr z>+s*73x7|1x3m7P^{%;RXS&zl+Hq%t(wWWLQ~xdN{rBbmsf{ZS?DhW9d47NU$D^%0 xnZ}PiT#CauAFc7e-^zJL+q?G7yQZ2&AIc(UWyRfn@^r(J`g@_Ia4t$sEJ;mK$j`G< z@XgOt@J%dIFfvszu(UF=urfARFfuSSRv;NLFnap(GcYieurM%G%z1mcJNQbd2>S>5 ze$|e$gEdYHE}KF$I8?ZGDjH+wb=y4I*|A&WboTUiZf}hqh3F=YE=>`s7XnO)eH}qU zZbqq(nAMNCw^%f+_wT#E>$E4%@Lq4w!?;Ls zgO6m|YO{0CtA4qJZ7*r46~0FzFlWkjRt9!|3lVrv+pA2iAC z;eCH#cLKk5qNqiK@)|bv0ygmq{%Z}|9}e)|ImrCu0cYF>F^-C6{vS;YEf3iL2{1GY zupCrq5{O`5D8Tkmfbq4W&>{t?Bn7cY4m=+>FbX9odpK})Zs0U?7CNCI8KfY3NkR6> zf(D)wEjlOKbQiS9oM07wa`5+r^Lzf}8-FgpvDSRPfpAOo!IOsmi=NzRjrccb&LrFR z)vuMBG(T!suiT}%tmSEv{S4XujB5!IQ`|h;_iQ;SdhycA{eMHWRNpm8#f8YIn*H`= zRh;pDA*Z~?fgKZn9^fvtQk|~qR>!@uuKJK)v;9PE>8hq{ni@OBu3gUP`BgoOed4?m zlWUJfE!)_CdhOqUVRolWi{;|U%Z*L#H;6h=#TZKPI|9` zWpv%X=xmU(X}zJmkwf!Hwd~yW{qM8pZE`cwyZ+*|!1j_WdcLl2q%Q0?Z2c~@*~|WB z=dbR*)HgbB()Y%$*y6nBQ*U~=YEnteqww8=&GS>f-sN5Ha{1-CT%H9fme+r*H0qcn zt7jf6?J|X@eBRGn?19H0?pzV;aB{6epP!qEYGv^0Ye!syKONYZx$Atp(C+go^)r41 zE7|DI%)dxs*_qP}&S&SJ)qaVdt$S|iubcj*64yi3UT^&osHFXQ)uq1r3wnQln0%6d z{9eDztuF2HyF+67b5?s#nkss&_D0m)51A7CuJ_)2-LmJ_&9i@uUbW7=E+X@$xG<;Q zsa!>M?zZyesU7`MS98BwYd&r55s`jvh27sr3Qo0xCyakAYwc~{xd7t3$uPQK*gwWC@2==F=!?{w~L ze*8lIy|D3N(f_Zk179$HWI4BdTgzu>#aAZBf0=D8unYUP=a_Ia4t$sEJ;mK$j`G< z@XgOt@J%dIFfvszu(UEVw=%I%FfuSSRv;NLFnap(GcYhju`n=H%z1mcJ2)g%hW*27 z=|n~41F@oJ6C&a?^(#+1a>U$wbLbPB;Do13y={&-xQ89BusCofMcCtJgNuNv7Y5c|9yd{26%Ml;t6t!q!7vmSpo-*^+> z76<8w1vMZs5&1$oJ+TqfDZ;MWgYK2JIb<))fc%KCm#wZ4lmbko_KG{Ddf#}f z*PeNTyI#J1@|)~A0&hOVhMw|WuO^>uo1;2I{I>nlH( zm)CW#@m{SLL!2^WtYGb(cA?A6s-Kgcyxv6rsaoqgCFJemJNA1tGy5-J_x-Y1#rR&8s`{i! zOD?SP5_?u&y68ZKMD4qspZ`7+oH9Y>^diqY>o?we^ysC4*Ey^2Ut=Bl*wuf9FL-~0 zwRC%s&Gk2b?+Ng|QtG?iywJp}YdUY#^o{nf1AczZ{Ns0Xy78m@8qJlQ@_D~i+a}39 zT;61+?^SVV1(#&ja^G7!?|%BXKcy!7OL(r_i;`~(T;w&+8_Q%72qKR~S#7bomM2Uzy@b^DbNraNfW28vERn zFT_jFRkSRY%Xy)7Gs*8{C!5}PyKnaM=iT1e@qDw8XHQE+IJZG}oWtFMh3hQ%OJDzC zF8{~pb!~5b))da@(-$uZ`11ul2*}?g{aHk-$iJuJnc}nXFDDl=CR(VzuFm%QGI5R3 z>B*Ub*(cx4oN?m#pM7x!k5v61-ue6N_n}IG^N&NrDz-=$f0$+c>aNepf2*pk-q(7Z ze7Jn_sXfY5&iKg0t@Szd|G=7Er9#gyndf|bXmPGT>vnqeCwV0+mme<{=WA8>L<&t?rfup$F^NMb$Ui}>TOkBvsXL^ XOV%;m{Al6%q$wxl&hz}8bsT8`gQ80x literal 0 HcmV?d00001 diff --git a/caseM.mat b/caseM.mat new file mode 100644 index 0000000000000000000000000000000000000000..5d96c580bbafbcf767712c7f805a558c29a4d23a GIT binary patch literal 1053 zcmeZu4DoSvQZUssQ1EpO(M`+DN!3vZ$Vn_o%P-2cQV4Jk_w>_Ia4t$sEJ;mK$j`G< z@XgOt@J%dIFfvszu(UEUv@$hQFfuSSRv;NLFnap(GcYhrWoBTgnBzN{Kj3nz#Buk< zf~rhXPBH4%mwf)cdu0_6!>_P#g+g-P!A(-E(hpqzEh{|a6);_P3X|2L4F`BFZN#(E zd+YrUUFK2e)_l>&u==uI#SBsFPaF1$m4DrRyY|~(+w-68&i}5xe?O_w!9ZDr!{vc; zQ~l2el@*S=A8vmxzl(dy?h@5`Kk5T_Oi-Djcygh^`WNOi+^WxJGnE6x2y%fgzzUjENj(VFP3(Bc4mkwU|c3tq|giR5@PQNYmaeb*lT-v-;XM*;aRk_O2tD zwQRrqT+$oazimwWxim`gm$j6}nV#Azxst0{!VkPJ&tCrRp^fjgrehPt7lcc%{dpnp z{_)oGBgUp>Rr4epzxiozy1M@8Py6NmKF)g={Cs-Ye16KLm2y?brLDqEBwp@rT=_F- zXLw&*X6cmcJJ;y?NH0Dj`p-gr;~_cM8Ns1pmH(NpN!#x8`zHBq^(3~bRSWLrx=3u- zn(^}TpT(z=JJh$|)yc1$`0^wF{~g*+F^_5wt>^!o*Jb$C)#3Sz{jM=t{3d@ z$i>{A@b|&qL_3M^x{3uR$6DSWw7k3TiP<0PMRNS_nWK+tPnvaL^?nbG# b|2yz6{aIZ>-uLoUtF!g5|L{&Ue`*N;k=6h? literal 0 HcmV?d00001 diff --git a/caseR.mat b/caseR.mat new file mode 100644 index 0000000000000000000000000000000000000000..d20c4f25216b3e492dbd2020d2c1c7240edfcc46 GIT binary patch literal 1016 zcmeZu4DoSvQZUssQ1EpO(M`+DN!3vZ$Vn_o%P-2cQV4Jk_w>_Ia4t$sEJ;mK$j`G< z@XgOt@J%dIFfvszu(UEZurfAOFfuSSRv;NLFnap(GcYg|Ff%Yz%yFH>8PJm`aopZZ ztLvbI=8+PgqX{ypUF_N6(S5FByIXX3n{Bc>78&WCdtBwIhBlAvY(e4FEtbo8iu~Tq ziaERN$qMu9n^LFQ&$^laZ1(-{^Pm5>JOA14{^K)+G0S<6-)#JLr~Y5lDK=;{`;!H6B@xQf-O8GSikVT(_`B! zbS&-I+VqmpX!m;$T8@0lP33YcKCgY5PdEKZ)U$N)rq6#4r-#lKlo#5o?7b?*#OPk? zd*>(NXS}k^o!`$b@V2ODUfd}4JgEA|<-XKkEROm-cNa@(r=}Fg$1Fcwx%|k7?O8v9 z-#t&g6Vv`{f7dOa7mw$doivvG=x_Mbv*^XL=Ue4ZAH2F<`1@o>_sFHV4ax3A|)iY@(FFq~PW@3D{_S4=7g|w8F87s>V{^@(>@4&Hlw-HRU9Yh=%R$hQCI=AXM& zCtOUpzdlg%LR9U}Y4%-#o@GrM5>>4c7tO`L%NHxEB?`T|8-0}X(R0IDD^B*+%A7nA zd+DD~yQL^wKZk*0z0VuTtqgVHe^2hPwAg=l`7P7)lLIr05)ULF>3ux!vw0#SVe8L-?y#*PDj*Op`%l0xAmR{a29J%29M%M1d2T~`k j%3yvLQL{GMiYMxd?6Zv=Z#MNB{5-w3@*ji2_YWrky_MiI literal 0 HcmV?d00001 diff --git a/plotVolt.asv b/plotVolt.asv new file mode 100644 index 0000000..84e68f8 --- /dev/null +++ b/plotVolt.asv @@ -0,0 +1,26 @@ +%% 画电压，为了写论文用 +load('case1V.mat'); +CV1=Volt; +load('case2V.mat'); +CV2=Volt; +load('case3V.mat'); +CV3=Volt; +load('caseR.mat'); +CVR=Volt; +load('caseM.mat'); +CVM=Volt; +subplot(3,2,1) +hist(CVR) +title('真实值') + + +subplot(2,2,3) +hist(CV1) +title('Case 1') +subplot(2,2,4) +hist(CV2) +title('Case 2') +subplot(2,2,5) +hist(CV3) +title('Case 3') + diff --git a/plotVolt.m b/plotVolt.m new file mode 100644 index 0000000..d8503c3 --- /dev/null +++ b/plotVolt.m @@ -0,0 +1,45 @@ +%% 画电压，为了写论文用 +load('case1V.mat'); +CV1=Volt; +load('case2V.mat'); +CV2=Volt; +load('case3V.mat'); +CV3=Volt; +load('caseR.mat'); +CVR=Volt; +load('caseM.mat'); +CVM=Volt; +subplot(3,2,1) +hist(CVR) +title('真实值') +subplot(3,2,2) +hist(CVM) +title('测量值') +subplot(3,2,3) +hist(CV1) +title('Case 1') +subplot(3,2,4) +hist(CV2) +title('Case 2') +subplot(3,2,5) +hist(CV3) +title('Case 3') + + +% figure() +% hold on +% +% plot(CVR) +% title('真实值') +% +% plot(CVM) +% title('测量值') +% +% plot(CV1) +% title('Case 1') +% +% plot(CV2) +% title('Case 2') +% +% plot(CV3) +% title('Case 3')