例 函数 1012)(iixxf对于适应度函数fitness对其参数w ，1c ，3c 做出不同方式的比较以测试其对函数结果影响。 当22111cc，5.12212cc，2.1w。 (适应函数1012)(iixxf) 程序1 当22111cc，5.12212cc，2.1w。 a)%主函数源程序（main.m） %------基本粒子群算法 （particle swarm optimization） %------名称： 基本粒子群算法 %------初始格式化 clear all; %清除所有变量 clc; %清屏 format long; %将数据显示为长整形科学计数 %------给定初始条条件------------------ N=40; %³初始化群体 个 数 D=10; %初始化群体 维 数 T=100; %初始化群体 最 迭 代 次 数 c11=2; %学习 因 子1 c21=2; %学习 因 子2 c12=1.5; c22=1.5; w=1.2; %惯 性 权 重 eps=10^(-6); %设 置 精 度（在 已 知 最 小 值 的时 候 用 ） %------初始化种 群个 体 （限 定位 置 和 速 度）------------ x=zeros(N,D); %x 是 位 置 ，初始化位 置 空 间 （矩 阵 ） v=zeros(N,D); %v 是 速 度，初始化速 度空 间 （矩 阵 ） for i=1:N for j=1:D x(i,j)=randn; %随 机 初始化位 置 ，randn 返 回 一 个 随 机 变化的符 合 正 态 分 布 的数 v(i,j)=randn; %随 机 初始化速 度 end end %------显示群位置---------------------- figu re(1) for j=1:D if(rem(D,2)>0) su bplot((D+1)/2,2,j) else su bplot(D/2,2,j) end plot(x (:,j),'b*');grid on %’b*’表示颜色是绿的，用*显示在图上 x label('粒子') y label('初始位置') tInfo=strcat('第',char(j+48),'维'); %strcat 使括号里的东西连成字符串 if(j>9) tInfo=strcat('第',char(floor(j/10)+48)); %floor 向负无穷方向取整 char(rem(j,10)+48,'维'); %rem 取余 end title(tInfo) end %------显示种群速度 figu re(2) for j=1:D if(rem(D,2)>0) su bplot((D+1)/2,2,j) else su bplot(D/2,2,j) end plot(v (:,j),'b*');grid on %是不是应该是 v (:,j) x label('粒子') y label('初始速度') tInfo=strcat('第',char(j+48),'维'); if(j>9) tInfo=strcat('第',char(floor(j/10)+48),'维'); char(rem(j,10)+48,'维); end title(tInfo) end figu re(3) %第一个图 su bplot(1,2,1) %-----...
