资源简介
用粒子群算法来求解函数的最优值,标准粒子群算法。
代码片段和文件信息
function [pso F] = pso_2D()
% FUNCTION PSO --------USE Particle Swarm Optimization Algorithm
%global present;
% close all;
pop_size = 10; % pop_size 种群大小
part_size =2; % part_size 粒子大小 ** =n-D
gbest = zeros(1part_size+1); % gbest 当前搜索到的最小的值
max_gen = 80; % max_gen 最大迭代次数
region=zeros(part_size1); % 设定搜索空间范围
region=[-33;-33]; % **每一维设定不同范围
rand(‘state‘sum(100*clock)); % 重置随机数发生器状态
arr_present = ini_pos(pop_sizepart_size); % present 当前位置随机初始化rand()的范围为0~1
v=ini_v(pop_sizepart_size); % 初始化当前速度
pbest = zeros(pop_sizepart_size+1); % pbest 粒子以前搜索到的最优值,最后一列包括这些值的适应度
w_max = 0.9; % w_max 权系数最大值
w_min = 0.4;
v_max = 0.2; % **最大速度为粒子的范围宽度
c1 = 2; % 学习因子
c2 = 2; % 学习因子
best_record = zeros(1max_gen); % best_record记录最好的粒子的适应度。
% ————————————————————————
% 计算原始种群的适应度及初始化
% ————————————————————————
arr_present(:end)=ini_fit(arr_presentpop_sizepart_size);
% for k=1:pop_size
% present(kend) = fitness(present(k1:part_size)); %计算原始种群的适应度
% end
pbest = arr_present; %初始化各个粒子最优值
[best_value best_index] = min(arr_present(:end)); %初始化全局最优,即适应度为全局最小的值,根据需要也可以选取为最大值
gbest = arr_present(best_index:);
%v = zeros(pop_size1); % v 速度
% ————————————————————————
% 迭代
% ————————————————————————
% global m;
% m = moviein(1000); %生成帧矩阵
x=[-3:0.01:3];
y=[-3:0.01:3];
z=@(xy) 3*(1-x).^2.*exp(-(x.^2) - (y+1).^2) ...
- 10*(x/5 - x.^3 - y.^5).*exp(-x.^2-y.^2) ...
- 1/3*exp(-(x+1).^2 - y.^2);
for i=1:max_gen
grid on;
% plot3(xyz);
% subplot(121)ezmesh(z)hold ongrid onplot3(arr_present(:1)arr_present(:2)arr_present(:3)‘*‘)hold off;
% subplot(122)ezmesh(z)view([14590])hold ongrid onplot3(arr_present(:1)arr_present(:2)arr_present(:3)‘*‘)hold off;
ezmesh(z) hold ongrid onplot3(arr_present(:1)arr_present(:2)arr_present(:3)‘*‘)hold off;
drawnow
F(i)=getframe;
% ezmesh(z)
% % view([-3790])
% hold on;
% grid on;
% % plot(-0.08980.7126‘ro‘);
% plot3(arr_present(:1)arr_present(:2)arr_present(:3)‘*‘); %改为三维
% axis([-2*pi2*pi-pipi-5010]);
% hold off;
pause(0.01);
% m(:i) = getframe; %添加图形
w = w_max-(w_max-w_min)*i/max_gen;
% fprintf(‘# %i 代开始!\n‘i);
% 确定是否对打散已经收敛的粒子群——————————————————————————————
reset = 0; % reset = 1时设置为粒子群过分收敛时将其打散,如果=1则不打散
if
相关资源
- 计算可见光通信OOK调制方式的BER mat
- matlab编的一个用于对统计数据进行威
- 最优化 一维不精确搜索wolfe matlab
- Matlab 广义高斯分布建模
- 遗传算法 车辆调度Matlab仿真
- matlab编译的bm3d去噪代码
- 利用Gabor滤波器实现的MATLAB掌纹识别系
- SPS并联机器人可达工作空间绘制的M
- puma机器人运动学正解MATLAB算法
- 追赶法求解三对角方程的MATLAB程序代
- 自适应遗传算法matlab代码
- matlab 在Hsv空间里画出颜色直方图程序
- 人脸识别算法matlab
- 层次聚类matlab程序
- rgb2lab程序matlab实现RGB到LAB转换
- 二维图像卷积matlab程序
- BP神经网络时间序列预测matlab代码ma
- 语音信号处理matlab GUI
- 关于四元数法捷联惯导解算算法的m
- matlab均衡器67280
- matlab快速傅里叶变换进行数字图像处
- 电池储能仿真
- matlab实现水平集方法
- matlab 圆形识别并标注
- 车牌字符分割垂直投影切分matlab程序
- LFM及匹配滤波的MATLAB程序
- 多光谱图像融合 matlab.zip
- matlab中仿真自适应信号处理LMSNewton算
- matlab中仿真自适应信号处理LMS算法
- 四步相移法程序
评论
共有 条评论