神剑山庄资源网 Design By www.hcban.com
测试函数主要是用来评估优化算法特性的,这里我用python3绘制了部分测试函数的图像。具体的测试函数可以结合维基百科来了解。想要显示某个测试函数的图片把代码结尾对应的注释去掉即可,具体代码如下:
import numpy as np import matplotlib.pyplot as plt from mpl_toolkits.mplot3d import Axes3D def draw_pic(X, Y, Z, z_max, title, z_min=0): fig = plt.figure() ax = Axes3D(fig) ax.plot_surface(X, Y, Z, rstride=1, cstride=1, cmap=plt.cm.hot) # ax.contourf(X, Y, Z, zdir='z', offset=-2, cmap=plt.cm.hot) ax.set_zlim(z_min, z_max) ax.set_title(title) # plt.savefig("./myProject/Algorithm/pic/%s.png" % title) # 保存图片 plt.show() def get_X_AND_Y(X_min, X_max, Y_min, Y_max): X = np.arange(X_min, X_max, 0.1) Y = np.arange(Y_min, Y_max, 0.1) X, Y = np.meshgrid(X, Y) return X, Y # rastrigin测试函数 def Rastrigin(X_min = -5.52, X_max = 5.12, Y_min = -5.12, Y_max = 5.12): A = 10 X, Y = get_X_AND_Y(X_min, X_max, Y_min, Y_max) Z = 2 * A + X ** 2 - A * np.cos(2 * np.pi * X) + Y ** 2 - A * np.cos(2 * np.pi * Y) return X, Y, Z, 100, "Rastrigin function" # Ackley测试函数 def Ackley(X_min = -5, X_max = 5, Y_min = -5, Y_max = 5): X, Y = get_X_AND_Y(X_min, X_max, Y_min, Y_max) Z = -20 * np.exp(-0.2 * np.sqrt(0.5 * (X**2 + Y**2))) - np.exp(0.5 * (np.cos(2 * np.pi * X) + np.cos(2 * np.pi * Y))) + np.e + 20 return X, Y, Z, 15, "Ackley function" # Sphere测试函数 def Sphere(X_min = -3, X_max = 3, Y_min = -3, Y_max = 3): X, Y = get_X_AND_Y(X_min, X_max, Y_min, Y_max) Z = X**2 + Y**2 return X, Y, Z, 20, "Sphere function" # beale测试函数 def Beale(X_min = -4.5, X_max = 4.5, Y_min = -4.5, Y_max = 4.5): X, Y = get_X_AND_Y(X_min, X_max, Y_min, Y_max) Z = np.power(1.5 - X + X * Y, 2) + np.power(2.25 - X + X * (Y ** 2), 2) + np.power(2.625 - X + X * (Y ** 3), 2) return X, Y, Z, 150000, "Beale function" # Booth测试函数 def Booth(X_min = -10, X_max = 10, Y_min = -10, Y_max = 10): X, Y = get_X_AND_Y(X_min, X_max, Y_min, Y_max) Z = np.power(X + 2*Y - 7, 2) + np.power(2 * X + Y - 5, 2) return X, Y, Z, 2500, "Booth function" # Bukin测试函数 def Bukin(X_min = -15, X_max = -5, Y_min = -3, Y_max = 3): X, Y = get_X_AND_Y(X_min, X_max, Y_min, Y_max) Z = 100 * np.sqrt(np.abs(Y - 0.01 * X**2)) + 0.01 * np.abs(X + 10) return X, Y, Z, 200, "Bukin function" # Three-hump camel测试函数 def three_humpCamel(X_min = -5, X_max = 5, Y_min = -5, Y_max = 5): X, Y = get_X_AND_Y(X_min, X_max, Y_min, Y_max) Z = 2 * X**2 - 1.05 * X**4 + (1/6) * X**6 + X*Y + Y*2 return X, Y, Z, 2000, "three-hump camel function" # Hölder table测试函数 def Holder_table(X_min = -10, X_max = 10, Y_min = -10, Y_max = 10): X, Y = get_X_AND_Y(X_min, X_max, Y_min, Y_max) Z = -np.abs(np.sin(X) * np.cos(Y) * np.exp(np.abs(1 - np.sqrt(X**2 + Y**2)/np.pi))) return X, Y, Z, 0, "Hölder table function", -20 z_min = None # X, Y, Z, z_max, title = Rastrigin() # X, Y, Z, z_max, title = Ackley() # X, Y, Z, z_max, title = Sphere() # X, Y, Z, z_max, title = Beale() X, Y, Z, z_max, title = Booth() # X, Y, Z, z_max, title = Bukin() # X, Y, Z, z_max, title = three_humpCamel() # X, Y, Z, z_max, title, z_min = Holder_table() draw_pic(X, Y, Z, z_max, title, z_min)
以下是上述代码绘制的测试函数的图像:
感觉图像的颜色还不是很好看,等之后优化了来改
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持。
标签:
python,测试函数
神剑山庄资源网 Design By www.hcban.com
神剑山庄资源网
免责声明:本站文章均来自网站采集或用户投稿,网站不提供任何软件下载或自行开发的软件!
如有用户或公司发现本站内容信息存在侵权行为,请邮件告知! 858582#qq.com
神剑山庄资源网 Design By www.hcban.com
暂无python绘制评估优化算法性能的测试函数的评论...
稳了!魔兽国服回归的3条重磅消息!官宣时间再确认!
昨天有一位朋友在大神群里分享,自己亚服账号被封号之后居然弹出了国服的封号信息对话框。
这里面让他访问的是一个国服的战网网址,com.cn和后面的zh都非常明白地表明这就是国服战网。
而他在复制这个网址并且进行登录之后,确实是网易的网址,也就是我们熟悉的停服之后国服发布的暴雪游戏产品运营到期开放退款的说明。这是一件比较奇怪的事情,因为以前都没有出现这样的情况,现在突然提示跳转到国服战网的网址,是不是说明了简体中文客户端已经开始进行更新了呢?
更新日志
2024年10月01日
2024年10月01日
- 孙悦.1996-伙伴【正大国际】【WAV+CUE】
- 纪钧瀚《钢琴阅读时光 雨中书店聆听轻音乐》[FLAC/分轨][399.62MB]
- 证声音乐图书馆《走向自然 疗心爵士乐》[320K/MP3][87.4MB]
- 证声音乐图书馆《走向自然 疗心爵士乐》[FLAC/分轨][184.94MB]
- 陈慧娴.2018-Priscilla-Ism演唱会3CD(2024环球红馆40复刻系列)【环球】【WAV+CUE】
- 郑秀文.1999-我应该得到(国)【华纳】【WAV+CUE】
- 陈家慧.2011-钢琴酒吧2CD【龙吟唱片】【WAV+CUE】
- 证声音乐图书馆《雨季 蓝调吉他 Rainy Blues》[320K/MP3][45.01MB]
- 证声音乐图书馆《雨季 蓝调吉他 Rainy Blues》[FLAC/分轨][109.13MB]
- 赞多《序章》[320K/MP3][45.54MB]
- 许巍.2004-每一刻都是崭新的【步升大风】【WAV+CUE】
- 群星.2024-四方馆影视原声带【韶愔音乐】【FLAC分轨】
- 陈雷.1997-安锁咧【金圆唱片】【WAV+CUE】
- 关淑怡.2013-MY.FAVORITE.SK.3CD【环球】【WAV+CUE】
- Sweety.2006-花言乔语【丰华】【WAV+CUE】