【优化求解】基于matlab改进的粒子群算法求解混合储能系统容量优化问题【含Matlab源码 1823期】.zip

%% 清空环境 clc clear %% 参数初始化 %粒子群算法中的两个参数 c1 = 1.49445; c2 = 1.49445; maxgen=400; %进化次数 sizepop=100; %种群规模 Vmax=10; %速度最大值 Vmin=-10; %速度最小值 popmax=[49060 9847457]; %粒子最大值 popmin=[1400 295198]; %粒子最小值 %% 产生初始粒子和速度 for i=1:sizepop %随机产生一个种群 pop(i,1)=(49060-1400)*rand+1400; %初始种群 pop(i,2)=(9847457-295198)*rand+295198; V(i,:)=0.5*rands(1,2); %初始化速度 %计算适应度 fitness(i)=funm(pop(i,:)); %染色体的适应度 %调用funm end %% 个体极值和群体极值 [bestfitness bestindex]=min(fitness); %找最小适应度 zbest=pop(bestindex,:); %全局最佳 gbest=pop; %个体最佳 fitnessgbest=fitness; %个体最佳适应度值 fitnesszbest=bestfitness; %全局最佳适应度值 %% 迭代寻优 for i=1:maxgen for j=1:sizepop %计算权重 weight=1; weight=0.9+(0.9-0.4)*exp(-20*(i^6)/(maxgen^6)); %计算学习因子 %速度更新 V(j,:) = weight*V(j,:) + c1*rand*(gbest(j,:) - pop(j,:)) + c2*rand*(zbest - pop(j,:)); V(j,find(V(j,:)>Vmax))=Vmax; V(j,find(V(j,:)<Vmin))=Vmin; %种群更新 pop(j,:)=pop(j,:)+V(j,:); for k=1:2 if pop(j,k)>popmax(k) pop(j,k)=popmax(k); end end for k=1:2 if pop(j,k)<popmin(k) pop(j,k)=popmin(k); end end %适应度值 fitness(j)=funm(pop(j,:)); end for j=1:sizepop %个体最优更新 if fitness(j) < fitnessgbest(j) gbest(j,:) = pop(j,:); fitnessgbest(j) = fitness(j); end %群体最优更新 if fitness(j) < fitnesszbest zbest = pop(j,:); fitnesszbest = fitness(j); end end yy(i)=fitnesszbest; %保存下每一代的最优值 end %% 结果分析 plot(yy,'k') title('最优个体适应度','fontsize',12); xlabel('进化代数','fontsize',12);ylabel('适应度','fontsize',12); zbest fitnesszbest