PSO求解函数最小值的MATLAB例程|MATLAB源代码

本篇文章适合PSO入门 ,进阶的可能会觉得太简单的。

目录

PSO

Particle Swarm Optimization,粒子群优化算法,通过模拟鸟群或鱼群的行为来寻找最优解。在计算时通过对一群粒子的位置和速度 进行迭代,是全局搜索算法。

相应的还有蚁群、遗传算法、模拟退火等智能算法。后面更新

例程

作用

这是一段PSO的例程。只有一个m文件,在里面填写好需要求解的函数后,就能得到PSO求解的函数最小值和对应的自变量取值。

在现在的函数示例中,设置的待求函数为F = x^2 + x - 6,从详情页能看到,求出来是当x=-0.25时,函数得到最小值-6.25

运行结果

设置如下的目标函数:
f ( x ) = x 2 + x − 6 f(x) = x^2 + x - 6 f(x)=x2+x−6

在代码编辑窗口里面如下:

代码

matlab 复制代码
% PSO求解函数最小值
% 2024-7-29/Ver1
clear;clc;close all;
rng(0);
[xm, fv] = PSO(@fitness, 1, 1.5, 1.5, 0.5, 100000, 1);
fprintf('满足目标函数取最小值时的自变量为:%f\n',xm);
fprintf('目标函数的最小值为:%f\n',fv);


function F = fitness(x) %设置目标函数
F = x^2 + x - 6;
end

function [xm, fv] = PSO(fitness, pop_size, c1, c2, w, epochs, chromosome_size) %PSO计算函数
% 完整代码下载链接:https://gf.bilibili.com/item/detail/1105924012

end

函数解释

PSO这个函数形式如下:
xm, fv = PSO(fitness, pop_size, c1, c2, w, epochs, chromosome_size)

其中,fitness是目标函数,所以在使用的时候这里填写目标函数"@fitness",pop_size和chromosome_size是每次计算的未知数个数(群体鬼墨),c1和c2是加速常数,w是惯性权重()使微粒保持运动的惯性,使其有扩展搜索空间的趋势,有能力探索新的区域)。epochs是迭代步数。

例程修改tips

  • 加大PSO函数的"epochs"这一个量(也就是倒数第二个输入量),可以加大迭代次数,计算时间增加、精度提升
  • 修改pop_size和chromosome_size,可以更改每次由PSO计算的量的个数
相关推荐
心祭丫丫2 小时前
【每天一个 Python 小案例17】Python 实现通讯录
开发语言·人工智能·笔记·python·学习
露天赏雪2 小时前
第3篇:Sentinel熔断降级实战:大模型接口抖动时的服务兜底方案
java·开发语言·人工智能·spring boot·spring·servlet·sentinel
yeccheng4 小时前
MySQL 索引超详解:从底层 B+ 树到面试高频考点
java·开发语言·jvm·mysql·面试
联犀4 小时前
纯 Go 的 CAD 开源库 go-cad:直接解析 DWG/DXF,现在也能导出 PDF
开发语言·golang·pdf
zjxtxdy5 小时前
C++面向对象程序开发-类与对象
开发语言·c++
Evand J7 小时前
【PMSM例程】(14)永磁同步电机无迹卡尔曼滤波(UKF)定子电阻在线辨识与补偿。电阻变化跟踪及预测误差对比仿真。附完整代码的下载链接
matlab·电机·滤波·卡尔曼滤波·永磁同步电机·无迹卡尔曼滤波·pmsm
qq_406981037 小时前
c++拷贝构造函数
开发语言·c++
爱喝雪碧的可乐8 小时前
MATLAB
开发语言·算法·matlab
H.莓飛9 小时前
【C++】【STL】vector_OJ题
linux·服务器·开发语言·数据结构·c++·visual studio
严肃小白兔10 小时前
MATLAB 优秀开源项目与工具箱推荐:涵盖通信、图像处理、脑机接口等九大领域
图像处理·matlab·信息与通信·开源项目