21 lines
581 B
Matlab
21 lines
581 B
Matlab
clear
|
|
clc
|
|
data=dlmread('数据.csv');
|
|
%% 数据处理
|
|
%按时段T把每个场景拿出来
|
|
snWind=data(:,1:3:end);
|
|
snSolar=data(:,2:3:end);
|
|
snP=data(:,3:3:end);
|
|
%% 初始化
|
|
% 随机生成一下初始场景序列
|
|
randSNCount=100; %要小于16^288
|
|
selectedIndex=round(1+(16-1)*rand(randSNCount,96));
|
|
selectedWind=zeros(randSNCount,96);
|
|
selectedSolar=zeros(randSNCount,96);
|
|
selectedP=zeros(randSNCount,96);
|
|
for I=1:randSNCount
|
|
onesIndex=sparse(selectedIndex(I,:),1:96,1,16,96);
|
|
selectedWind(I,:)=sum(snWind.*onesIndex);
|
|
selectedSolar(I,:)=sum(snSolar.*onesIndex);
|
|
selectedP(I,:)=sum(snP.*onesIndex);
|
|
end |