利用python实现PSO算法优化二元函数
作者:爱帕克的喃 发布时间:2022-06-17 20:53:04
标签:python,二元函数,算法
python实现PSO算法优化二元函数,具体代码如下所示:
import numpy as np
import random
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
#----------------------PSO参数设置---------------------------------
class PSO():
def __init__(self,pN,dim,max_iter): #初始化类 设置粒子数量 位置信息维度 最大迭代次数
#self.w = 0.8
self.ws = 0.9
self.we = 0.4
self.c1 = 1.49445
self.c2 = 1.49445
self.r1= 0.6
self.r2= 0.3
self.pN = pN #粒子数量
self.dim = dim #搜索维度
self.max_iter = max_iter #迭代次数
self.X = np.zeros((self.pN,self.dim)) #所有粒子的位置(还要确定取值范围)
self.Xmax = 5
self.Xmin = -5
self.V = np.zeros((self.pN,self.dim)) #所有粒子的速度(还要确定取值范围)
self.Vmax = 1
self.Vmin = -1
self.pbest = np.zeros((self.pN,self.dim)) #个体经历的最佳位置
self.gbest = np.zeros((1,self.dim)) #全局最佳位置
self.p_fit = np.zeros(self.pN) #每个个体的历史最佳适应值
self.fit = 0 #全局最佳适应值
#---------------------目标函数Sphere函数-----------------------------
def function(self,x):
y = np.sin(10*np.pi*x)/x
return y
def Holder_table(self,x,y):
z = -np.abs(np.sin(x) * np.cos(y) * np.exp(np.abs(1 - np.sqrt(x**2 + y**2)/np.pi)))
return z
def fuck(self,x,y):
z = x**2 + y**2 - 10*np.cos(2*np.pi*x) - 10*np.cos(2*np.pi*y) + 20
return z
#---------------------初始化种群----------------------------------
def init_Population(self):
for i in range(self.pN): #遍历所有粒子
for j in range(self.dim): #每一个粒子的纬度
self.X[i][j] = random.uniform(-5,5) #给每一个粒子的位置赋一个初始随机值(在一定范围内)
self.V[i][j] = random.uniform(-1,1) #给每一个粒子的速度给一个初始随机值(在一定范围内)
self.pbest[i] = self.X[i] #把当前粒子位置作为这个粒子的最优位置
tmp = self.fuck(self.X[i][0],self.X[i][1]) #计算这个粒子的适应度值
self.p_fit[i] = tmp #当前粒子的适应度值作为个体最优值
if(tmp > self.fit): #与当前全局最优值做比较并选取更佳的全局最优值
self.fit = tmp
self.gbest = self.X[i]
#---------------------更新粒子位置----------------------------------
def iterator(self):
fitness = []
for t in range(self.max_iter):
w = self.ws - (self.ws - self.we) * (t / self.max_iter)
for i in range(self.pN):
#更新速度
self.V[i] = w*self.V[i] + self.c1*self.r1*(self.pbest[i] - self.X[i]) + self.c2*self.r2*(self.gbest - self.X[i])
if self.V[i][0] > self.Vmax:
self.V[i][0] = self.Vmax
elif self.V[i][0] < self.Vmin:
self.V[i][0] = self.Vmin
if self.V[i][1] > self.Vmax:
self.V[i][1] = self.Vmax
elif self.V[i][1] < self.Vmin:
self.V[i][1] = self.Vmin
#更新位置
self.X[i] = self.X[i] + self.V[i]
if self.X[i][0] > self.Xmax:
self.X[i][0] = self.Xmax
elif self.X[i][0] < self.Xmin:
self.X[i][0] = self.Xmin
if self.X[i][1] > self.Xmax:
self.X[i][1] = self.Xmax
elif self.X[i][1] < self.Xmin:
self.X[i][1] = self.Xmin
for i in range(self.pN): #更新gbest\pbest
temp = self.fuck(self.X[i][0],self.X[i][1])
if(temp > self.p_fit[i]): #更新个体最优
self.pbest[i] = self.X[i]
self.p_fit[i] = temp
if(temp > self.fit): #更新全局最优
self.gbest = self.X[i]
self.fit = temp
fitness.append(self.fit)
print('最优值为:',self.fit)#输出最优值
z1 = self.fit
print('最优位置为:',self.X[i][0],self.X[i][1])
x1 = self.X[i][0]
y1 = self.X[i][1]
return fitness, z1, x1,y1
#----------------------程序执行-----------------------
my_pso = PSO(pN=100,dim=2,max_iter=200)
my_pso.init_Population()
fitness,z1,x1,y1 = my_pso.iterator()
plt.figure(1)
plt.title("Figure1")
plt.xlabel("iterators", size=14)
plt.ylabel("fitness", size=14)
t = np.array([t for t in range(0,200)])
fitness = np.array(fitness)
plt.plot(t,fitness, color='b',linewidth=3)
plt.show()
fig = plt.figure(figsize=(15,10))
ax = Axes3D(fig)
X = np.arange(-5,5,0.1)
Y = np.arange(-5,5,0.1)
X,Y = np.meshgrid(X,Y)
def f(x,y):
return (x**2 + y**2 - 10*np.cos(2*np.pi*x) - 10*np.cos(2*np.pi*y) + 20)
ax.plot_surface(X,Y,f(X,Y),rstride=1,cstride=1,cmap= plt.get_cmap('rainbow'))
ax.scatter(x1, y1, z1,s=400,c='k',marker = '*')
plt.show()
效果图如下
总结
以上所述是小编给大家介绍的利用python实现PSO算法优化二元函数网站的支持!
如果你觉得本文对你有帮助,欢迎转载,烦请注明出处,谢谢!
来源:https://blog.csdn.net/weixin_42553833/article/details/103007620
0
投稿
猜你喜欢
- 官方文档https://developers.weixin.qq.com/miniprogram/dev/framework/open-ab
- rss.asp格式的 下面代码保存为rss.asp 代码如下:<!--#include file="conn.as
- 由于新版站长资讯即将发布,我每天都在想如何防止采集,目前想到的几种办法:1、多做几个列表和内容模板,随机使用,对目前大多数cms来说,这种方
- asp十进制转二进制;二进制转十进制;二进制转十六进制;十六进制转二进制;八进制转二进制'二进制转八进制;八进制转十进制;十六进制转
- 1、使用mysqldump工具将MySql数据库备份mysqldump -u root -p -c --default-character-
- 本文实例讲述了Python实现分割文件及合并文件的方法。分享给大家供大家参考。具体如下:分割文件split.py如下:#!/usr/bin/
- 有时候我们会在页面上显示用户的所在地区,这个原理是:先得到用户的IP,然后去查询将IP转换成一个数值,最后去查这个数值所在的范围,来得到用户
- 格式为: SCRENC [/s] [/f] [/xl] [/l defLanguage ] [/e de
- 在通过拼组sql语句来实现数据插入的应用中,我们很有可能会遇到需要插入大型数据的情况,例如,在oracle中需要插入字节数超过4000的字段
- 一、检测网络信息和系统信息 在Frontpage 2000 的Explorer管理器中选择帮助(Help)|关于Frontpage管理器(A
- 本文实例讲述了Python简单获取网卡名称及其IP地址的方法。分享给大家供大家参考,具体如下:windows上想知道网卡和IP地址,可以使用
- Python-OpenCV环境的配置看上一篇OpenCV环境的配置本篇主要介绍一下OpenCV的基本使用和相关函数的介绍。以下所有操作都基于
- 一、前言在调用函数时,大多数情况下,主调函数和被调函用之间有数据传递关系,这就是有参数的函数形式。函数参数的作用是传递数据给函数使用,函数利
- 从技术上来说,在ASP环境中,读入并管理XML文本的主要方法有三种: 创建MSXML对象,并且将XML文档载入DOM; 使用服务器端嵌入(S
- by yemoo有时在编写网页代码时发现,img底部莫名奇妙多出大约3px的空白,无论怎么调节css都不可以,今天再次遇到此问题,网上看了一
- 第一节:WAP的潜能 这些日子,我们常听到WAP技术,一种手机上网的技术。从技术上讲,移动电话不可能和PC来竞争,移动电话的屏幕只能容下很少
- 内容摘要:现在博客很流行,相信应该上网时间稍微长点的朋友都会在这或者在那的有一个自己的博客。对于一些有一定能力的朋友,可能更喜欢自己去下载一
- 在项目里碰到需要把类似'450000'的数字转换为会计记账所用的格式,'450,000.00',分隔千分位和
- asp无组件上传VBS编写的大家见的多了,这个是纯javascript实现的上传,原来unicode可以解决读取位置的问题,这次真的是纯JS
- A 定义数组有两种方式:DIM和REDIM。DIM定义的是固定个数、数据类型的数组;而REDIM则不同,它可以定义不同类型的数据,也可以定义