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

水电站中长期调度优化的粒子群算法 MATLAB源码 水电站中长期优化调度主要是在满足电力系统出力要求以及下游综合用水 要求等的前提下，合理地安排水库各个月份的发电流量，使得调度期内的总发电 量最大。即在给定预报入库 流量过程线、下游 综合用水流量过程线、 调度期起始 水位和终止水位等条件下，通 过总水量的合理运用使调度期内水电站的总发电 量最大。 % 参数设置 % GreenSim 团队专业级算法设计%出力系数，常数 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=PSO(K,N,w,c1,c2,VtLB,VtUB,QtLB,QtUB,NtLB,NtUB,qt,A,Tt); % X=BESTXK; Vt,Qt,St=DeCode(X); disp(最佳蓄水量); disp(Vt); disp(最佳平均发电流量); disp(Qt); disp(最佳平均弃水流量); disp(St); disp(最大总发电量); disp(-BESTY(K); function BESTX,BESTY,ALLX,ALLY=PSO(K,N,w,c1,c2,VtLB,VtUB,QtLB,QtUB,NtLB,NtUB,qt,A,Tt) % GreenSim 团队专业级算法设计 farm=zeros(N,M); for i=N Vt,Qt,St=Initialize2(VtLB,VtUB,QtLB,QtUB,NtLB,NtUB,qt,A,Tt); X=EnCode(Vt,Qt,St); farm(i,:)=X; end %粒子历史最优状态初始化 欢迎访问 GreenSim 团队主页/greensim 邮箱： 第 3 页 Pfarm=farm; %粒子群全局最优状态 for i=1:N X=Pfarm(i,:); Vt,Qt,St=DeCode(X); SE=ObjFun(A,Qt,St,Vt,Tt); Fitness(i)=SE; end MinFit=min(Fitness); POS=find(Fitness=MinFit); Gfarm=Pfarm(POS(1),:); %输出变量初始化 ALLX=cell(K,1);%细胞结构，每一个元素是 MN 矩阵，记录每一代的个体 ALLY=zeros(K,N);%KN 矩阵，记录每一代评价函数值 BESTX=cell(K,1);%细胞结构，每一个元素是 M1 向量，记录每一代的最优个体 BESTY=zeros(K,1);%K1 矩阵，记录每一代的最优个体的评价函数值 k=1;%迭代计数器初始化 % 第二步：迭代过程 while k=K % 粒子状态更新 newfarm=farm; for i=1:N newfarm(i,:)=w*farm(i,:)+c1*rand*(Pfarm(i,:)-farm(i,:)+c2*rand*(Gfarm-farm(i,:); X=newfarm(i,:); Vt,Qt,St=DeCode(X); Flag,Vt,Qt,St=Correct(Vt,Qt,St,VtLB,VtUB,QtLB,QtUB,NtLB,NtUB,qt,A,Tt); X=EnCode(Vt,Qt,St); newfarm(i,:)=X; end % 历史最优状态和全局最优状态更新 NEWFIT=zeros(1,N); for i=1:(N) aa=newfarm(i,:); Vt,Qt,St=DeCode(aa); FitA=ObjFun(A,Qt,St,Vt,Tt); bb=Pfarm(i,:); Vt,Qt,St=DeCode(bb); FitB=ObjFun(A,Qt,St,Vt,Tt); NEWFIT(i)=FitA; if FitAFitB Pfarm(i,:)=A; end if FitAGfarm Gfarm=A; end end % 记录 ALLXk=newfar