From d758ad50050c393ac471d5a9194795eda8c02bc2 Mon Sep 17 00:00:00 2001 From: "facat@lab.com" Date: Mon, 24 Nov 2014 11:50:06 +0800 Subject: [PATCH] =?UTF-8?q?=E5=86=99=E5=AE=8C=E4=BA=86?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: facat@lab.com --- main.m | 14 ++++++++++++-- 1 file changed, 12 insertions(+), 2 deletions(-) diff --git a/main.m b/main.m index 0e8063a..763cb63 100644 --- a/main.m +++ b/main.m @@ -7,8 +7,8 @@ snWind=data(:,1:3:end); snSolar=data(:,2:3:end); snP=data(:,3:3:end); snP(snP==-100)=1; -randSNCount=100; %要小于16^288 。解的场景数量 -gaCount=100;%粒子数 +randSNCount=500; %要小于16^288 。解的场景数量 +gaCount=20;%粒子数 %% 初始化 ga=repmat(struct(),gaCount,1); for I=1:gaCount @@ -42,5 +42,15 @@ for LOOP=1:100 ga(I).fitness=fitness; end end +%利用优化后的结果给出新的场景 +newSnWind=zeros(randSNCount,96); +newSnSolar=zeros(randSNCount,96); +newSnP=zeros(randSNCount,96); +for I=1:randSNCount + onesIndex=sparse(bestIndex(I,:),1:96,1,16,96); + newSnWind(I,:)=sum(onesIndex.*snWind); + newSnSolar(I,:)=sum(onesIndex.*snSolar); + newSnP(I,:)=sum(onesIndex.*snP); +end