33水电站中长期调度优化的粒子群算法MATLAB源码

水电站中长期调度优化的粒子群算法MATLAB源码

水电站中长期优化调度主要是在满足电力系统出力要求以及下游综合用水

要求等的前提下，合理地安排水库各个月份的发电流量，使得调度期内的总发电量最大。即在给定预报入库流量过程线、下游综合用水流量过程线、调度期起始水位和终止水位等条件下，通过总水量的合理运用使调度期内水电站的总发电量最大。

%%参数设置% GreenSim团队一一专业级算法设计&代写程序

% 欢迎访问GreenSim 团队主页宀/greensim

A=8.5;%出力系数，常数

Tt=730*ones(12,1);%第t个时段的小时数

%注意：一年按363天*24小时算，均分为12个月

HtLB=55*ones(12,1);% 第t时段水位约束的下界，单位：米

HtUB=[65;65;65;61;61;61;61;65;65;65;65;65];% 第t 时段水位约束的上界，单位：米

VtLB=zeros(12,1);

VtUB=zeros(12,1);

for i=1:12

VtLB(i)=Ht2Vt(HtLB(i)); VtUB(i)=Ht2Vt(HtUB(i));

end

%注意：蓄水量Vt和水位Ht之间有----- 对应的关系，单位：立方米

NtLB=260000*ones(12,1);% 出力约束的下界，单位：千瓦NtUB=1400000*ones(12,1);% 出力约束的上届，单位：千瓦% 注意：Nt=A*Qt*Ht

QtLB=308*ones(12,1);% 泄流量下界，单位：立方米/秒

QtUB=29200*ones(12,1);% 泄流量上界，单位：立方米/秒

qt=[373;859;1568;2100;3210;5049;1596;1160;925;781;572;1010];% 入库流量，单位：立方米/ 秒

%注意：以上三个量，时间单位相乘时，小时乘以3600转化成秒

%%调用粒子群算法

K=60;

N=80;

w=0.5;

c1=0.3;

c2=0.2; [BESTX,BESTY ,ALLX,ALLY]=PS0(K,N,w,c1,c2,VtLB,VtUB,QtLB,QtUB,NtLB,NtUB,qt,A,Tt)

5

%% X=BESTX{K};

[Vt,Qt,St]=DeCode(X);

disp('最佳蓄水量');

disp(Vt);

disp('最佳平均发电流量');

disp(Qt);

disp('最佳平均弃水流量');

disp(St);

disp('最大总发电量');

disp(-BESTY(K));

fun ctio n

[BESTX,BESTY ,ALLX,ALLY]=PS0(K,N,w,c1,c2,VtLB,VtUB,QtLB,QtUB,NtLB,NtUB,qt,A,Tt) % GreenSim团队一一专业级算法设计&代写程序

% 欢迎访问GreenSim 团队主页宀/greensim

%%输入参数列表

% K 迭代次数

% N 粒子个数

% w 惯性因子

% c1 加速因子，针对历史最优状态

% c2 加速因子，针对全局最优状态

% VtLB 蓄水量约束下界，立方米，N*1 向量

% VtUB 蓄水量约束上界，立方米，N*1 向量、

% QtLB 平均发电流量约束下界，立方米/秒, N*1 向量

% QtUB 平均发电流量约束上界，立方米/秒，N*1 向量

% NtLB 出力约束下界，千瓦，N*1 向量

% NtUB 出力约束上界，千瓦，N*1 向量

% qt 入库流量，立方米/秒，N*1、向量

% A 出力系数

% Tt 一个周期的小时数

%%输出参数列表

% BESTX K X 1细胞结构，每一个元素是M X 1向量，记录每一代的最优个体

% BESTY K X 1矩阵，记录每一代的最优个体的评价函数值

% ALLX K X 1细胞结构，每一个元素是M X N矩阵，记录全部个体

% ALLY K X N矩阵，记录全部个体的评价函数值

%%第一步：粒子状态初始化

M=3*le ngth(VtLB);

farm=zeros(N,M);

for i=N

[Vt,Qt,St]=l nitialize2(VtLB,VtUB,QtLB,QtUB,NtLB,NtUB,qt,A,Tt);

X=E nCode(Vt,Qt,St);

farm(i,:)=X';

end

%粒子历史最优状态初始化

Pfarm=farm;

%粒子群全局最优状态

for i=1:N

X=Pfarm(i,:);

[Vt,Qt,St]=DeCode(X');

SE=0bjFu n(A,Qt,St,Vt,Tt);

Fit ness(i)=SE;

end

MinF it= min( Fit ness);

POS=fi nd(Fit ness==Mi nFit);

Gfarm=Pfarm(POS(1),:);

%输出变量初始化

ALLX=cell(K,1);% 细胞结构，每一个元素是M X N矩阵，记录每一代的个体ALLY=zeros(K,N);%K X N矩阵，记录每一代评价函数值

BESTX=cell(K,1);% 细胞结构，每一个元素是M X 1向量，记录每一代的最优个体

BESTY=zeros(K,1);%K X 1矩阵，记录每一代的最优个体的评价函数值k=1;%迭代计数器初始化

%%第二步：迭代过程

while k<=K

%%粒子状态更新

n ewfarm=farm;

for i=1:N

n ewfarm(i,:)=w*farm(i,:)+c1*ra nd*(Pfarm(i,:)-farm(i,:))+c2*ra nd*(Gfarm-farm(i,:)); X=n

ewfarm(i,:);

[Vt,Qt,St]=DeCode(X');

[Flag,Vt,Qt,St]=Correct(Vt,Qt,St,VtLB,VtUB,QtLB,QtUB,NtLB,NtUB,qt,A,Tt); X=E

nCode(Vt,Qt,St);

n ewfarm(i,:)=X';

end

%%历史最优状态和全局最优状态更新

NEWFIT=zeros(1,N);

for i=1:(N)

aa=n ewfarm(i,:);

[Vt,Qt,St]=DeCode(aa');

FitA=ObjFu n(A,Qt,St,Vt,Tt); bb=Pfarm(i,:);

[Vt,Qt,St]=DeCode(bb');

FitB=ObjFu n(A,Qt,St,Vt,Tt); NEWFIT(i)=FitA;

if FitA

Pfarm(i,:)=A;

end

if FitA

Gfarm=A;

end

end

%%记录