95号资源-源程序：论文可在知网下载《计及需求响应和电能交互的多主体综合能源系统 主从博弈优化调度策略 》本人博客有解读

clc clear; Max_Dt=20;%最大迭代次数300 D=100;%搜索空间维数（未知数个数） （1-24为售电价格，25-48为售热价格，49-96电需求响应价格，96-144为热需求响应价格） N=10;%粒子个数600 w_max=0.9; w_min=0.4; v_max=0.2; s=1; %%%%%粒子初始化%%%%% x=zeros(N,D); for i=1:N for j=1:9 v(i,j)=0.0; x(i,j)=0.35+rand()*0.03; end for j=10:12 v(i,j)=0.0; x(i,j)=0.35+rand()*0.45; end for j=13:16 v(i,j)=0.0; x(i,j)=0.35+rand()*0.85; end for j=17:20 v(i,j)=0.0; x(i,j)=0.35+rand()*0.45; end for j=21:22 v(i,j)=0.0; x(i,j)=0.35+rand()*0.85; end for j=23:24 v(i,j)=0.0; x(i,j)=0.35+rand()*0.45; end for j=25:48 v(i,j)=0.0; x(i,j)=0.2+rand()*0.3; end for j=1+48:9+48 v(i,j)=0.0; x(i,j)=0.1+rand()*0.1; end for j=10+48:12+48 v(i,j)=0.0; x(i,j)=0.1+rand()*0.3; end for j=13+48:16+48 v(i,j)=0.0; x(i,j)=0.1+rand()*0.55; end for j=17+48:20+48 v(i,j)=0.0; x(i,j)=0.1+rand()*0.3; end for j=21+48:22+48 v(i,j)=0.0; x(i,j)=0.1+rand()*0.55; end for j=23+48:24+48 v(i,j)=0.0; x(i,j)=0.1+rand()*0.3; end for j=1+72:9+72 v(i,j)=0.0; x(i,j)=0.05+rand()*0.05; end for j=10+72:12+72 v(i,j)=0.0; x(i,j)=0.05+rand()*0.25; end for j=13+72:16+72 v(i,j)=0.0; x(i,j)=0.05+rand()*0.45; end for j=17+72:20+72 v(i,j)=0.0; x(i,j)=0.05+rand()*0.25; end for j=21+72:22+72 v(i,j)=0.0; x(i,j)=0.05+rand()*0.45; end for j=23+72:24+72 v(i,j)=0.0; x(i,j)=0.05+rand()*0.25; end for j=97:144 v(i,j)=0.0; x(i,j)=0.05+rand()*0.15; end end for i=1:N p(i)=fitness(x(i,:),s); y(i,:)=x(i,:);%每个粒子的个体寻优值 end Pbest=fitness(x(1,:),s); pg=x(1,:);%Pg为全局最优 for i=2:N if fitness(x(i,:),s)<fitness(pg,s) Pbest=fitness(x(i,:),s); pg=x(i,:);%全局最优更新 end end %进入主循环***************************************** for t=1:Max_Dt t for i=1:N w=w_max-(w_max-w_min)*t/Max_Dt;%惯性权重更新 c1=(0.5-2.5)*t/Max_Dt+2.5; %认知 c2=(2.5-0.5)*t/Max_Dt+0.5; %社会认识 w=0.7; v(i,:)=w*v(i,:)+c1*rand()*(y(i,:)-x(i,:))+c2*rand()*(pg-x(i,:)); for m=1:D if(v(i,m)>v_max) v(i,m)=v_max; elseif(v(i,m)<-v_max) v(i,m)=-v_max; end end x(i,:)=x(i,:)+v(i,:); %对粒子边界处理***************************** for n=1:9 if x(i,n)>0.38 x(i,n)=0.38; v(i,n)=-v(i,n); elseif x(i,n)<0.35 x(i,n)=0.35; v(i,n)=-v(i,n); end end for n=10:12 if x(i,n)>0.8 x(i,n)=0.8; v(i,n)=-v(i,n); elseif x(i,n)<0.35 x(i,n)=0.35; v(i,n)=-v(i,n); end end for n=13:16 if x(i,n)>1.2 x(i,n)=1.2; v(i,n)=-v(i,n); elseif x(i,n)<0.35 x(i,n)=0.35; v(i,n)=-v(i,n); end end for n=17:20 if x(i,n)>0.8 x(i,n)=0.8; v(i,n)=-v(i,n); elseif x(i,n)<0.35 x(i,n)=0.35; v(i,n)=-v(i,n); end end for n=21:22 if x(i,n)>1.2 x(i,n)=1.2; v(i,n)=-v(i,n); elseif x(i,n)<0.35 x(i,n)=0.35; v(i,n)=-v(i,n); end end for n=23:24 if x(i,n)>0.8 x(i,n)=0.8; v(i,n)=-v(i,n); elseif x(i,n)<0.35 x(i,n)=0.35; v(i,n)=-v(i,n); end end for n=25:48 if x(i,n)>0.5 x(i,n)=0.5; v(i,n)=-v(i,n); elseif x(i,n)<0.2 x(i,n)=0.2; v(i,n)=-v(i,n); end end for n=1+48:9+48 if x(i,n)>0.2 x(i,n)=0.2; v(i,n)=-v(i,n); elseif x(i,n)<0.1 x(i,n)=0.1; v(i,n)=-v(i,n); end end for n=10+48:12+48 if x(i,n)>0.4 x(i,n)=0.4; v(i,n)=-v(i,n); elseif x(i,n)<0.1 x(i,n)=0.1; v(i,n)=-v(i,n); end end for n=13+48:16+48 if x(i,n)>0.65 x(i,n)=0.65; v(i,n)=-v(i,n); elseif x(i,n)<0.1 x(i,n)=0.1; v(i,n)=-v(i,n); end end for n=17+48:20+48 if x(i,n)>0.4 x(i,n)=0.4; v(i,n)=-v(i,n); elseif x(i,n)<0.1 x(i,n)=0.1; v(i,n)=-v(i,n); end end for n=21+48:22+48 if x(i,n)>0.65 x(i,n)=0.65; v(i,n)=-v(i,n); elseif x(i,n)<0.1 x(i,n)=0.1; v(i,n)=-v(i,n); end end for n=23+48:24+48 if x(i,n)>0.4 x(i,n)=0.4; v(i,n)=-v(i,n); elseif x(i,n)<0.1 x(i,n)=0.1; v(i,n)=-v(i,n); end end for n=1+72:9+72 if x(i,n)>0.1 x(i,n)=0.1; v(i,n)=-v(i,n); elseif x(i,n)<0.05 x(i,n)=0.05; v(i,n)=-v(i,n); end end for n=10+72:12+72 if x(i,n)>0.3 x(i,n)=0.3; v(i,n)=-v(i,n); elseif x(i,n)<0.05 x(i,n)=0.05; v(i,n)=-v(i,n); end end for n=13+72:16+72 if x(i,n)>0.5 x(i,n)=0.5; v(i,n)=-v(i,n); elseif x(i,n)<0.05 x(i,n)=0.05; v(i,n)=-v(i,n); end end for n=17+72:20+72 if x(i,n)>0.3 x(i,n)=0.3; v(i,n)=-v(i,n); elseif x(i,n)<0.05 x(i,n)=0.