和声搜索算法的matlab源代码使用说明

    本人以前借鉴使用他人的和声搜索算法代码时，感觉在套用到自己的实际问题中，需要修改其中不少地方。而且对他人的代码处理过程，读明白也需要花费不少时间，效率也不高。后来决定将该算法用Matlab语言，采用面向对象的方式重新编写一遍。将大部分不需要修改的地方都封装起来，用户只需要根据自己的实际问题需要，填写基本的几个参数，修改一下问题的适应度函数，就能够快速让和声搜索算法跑起来，然后看到算法运行的效果。这对于那些不熟悉编程的科研工作者和学生，会节省大量的时间，并能帮助他们把精力集中于自己要解决的主要问题上，而不是花费在代码编写上。而如果他们有兴趣，也可以自己进一步修改代码，从而使得代码更符合自己的使用习惯。

  下面会逐步展开介绍使用本人编写的和声搜索算法matlab代码的过程。

一 和声算法的使用

和声搜索算法的原理介绍可以在网上搜索到多篇文章，仔细阅读就可以对算法过程得到清楚的了解，因此这里不在赘述。这里主要说明在使用和声算法时，需要完成的两个关键任务。

第1，把要解决的具体问题转化成一个求取函数极值的问题。

一般是求取函数最大值问题，求取最小值问题是等价的，做个变换就可以成为最大值问题。该函数在和声算法中一般称为适应度函数，有时也叫目标函数。一个函数最大值问题，其数学表达如下：

第2，使用和声算法时，需要预先设定各种必需参数值。最基本的参数有如下5个，还有与适应度函数相关的3个重要参数。

5个基本参数分别为：

合声记忆库的大小HMS，就是和声库中保存的合成向量的数量。
合声记忆库取值的概率HMCR，就是从和声库中取出一个和声的概率。
音调微调的概率PAR，就是对取出的和声进行微调的概率。
音调微调的带宽BW，也就是微调的幅度。
创作的次数Tmax，也就是程序迭代的次数。

3个相关重要参数分别为：
和声向量的维数，即问题的解空间维数。
和声向量在各个维度上的取值下限。
和声向量在各个维度上的取值下限。

这几个参数实际就是确定和声向量的取值范围。

二 和声算法使用举例

下面以使用和声算法搜索一个函数的最大值，来举个例子。该函数的解析式如下：

约束条件为：

该函数的最大值我们已经知道，这里提前说明，其最大值在18.5附近。这样该函数x,y的取值范围在区间[0,10]之间，在三维空间的图形如下：

三 和声算法matlab源代码的修改使用

和声算法matlab源代码的修改使用分为如下几步：

1 下载和声搜索算法代码并解压使用

 该代码已经压缩打包后放在百度网盘上，网址请点击：百度网盘网址，提取码：1234

用户下载后，解压到本地文件夹下，打开matlab软件，进入该文件夹就能使用。代码使用Matlab 2020b版本开发，用户使用的matlab版本最好与这个版本相同或更高些，这样可以保证使用时少出意外问题。

2 先修改适应度函数

   用户根据自己的具体问题，确定了适应度函数以后，只需要在SupportFiles文件夹下的Fitness.m文件中修改fitness函数内容，为自己的适应度函数以后，就完成了使用和声算法需要的第1个修改任务。仍然以上述例子函数为适应度函数，修改后的fitness函数如下：

这里代码中的X(1,1)就是例子函数中的x分量，X(1,2)就是例子函数中的y分量。 

3 参数上的修改

  和声算法的参数修改在主文件mainHS中，在matlab下打开该文件，用户可以按照自己问题的参数需要，修改代码文件中如下图中红框内的和声算法各参数数值即可。红框内是例子函数需要的各参数数值，各参数的功能已经在后边用绿色注释进行了说明，以确保不会填错。

 4 运行

完成以上修改后，点击matlab的运行按钮，就可以正常运行和声搜索算法。在下方的命令窗口内，可以看到迭代信息。如图：

四 和声算法运行结果图示

    算法运行完，会把各轮迭代的最优值画出，最后形成最优适应度变化曲线图。但需要注意的是，因为该算法属于随机搜索，所以不是每次都一定能找到全局最优解，有些时候，可能只能找到次优解。下图是算法运行一次后，得到的最优适应度变化曲线图。

 可以看到，算法在首次发现次优解18.4后，并没有找到全局最优解18.5.。

五 结果的三维图示

   为了方便直观的观察最优解在适应度函数上的位置，这里把最优解用红色星号标注在适应度函数三维图上。可以看到，最优解已经找到了全局最优值所在的主山峰，但可惜，还差一点，没有找到全局最优值。注意，本图只能显示二维函数求取极值的图像，高于二维，就没法画图显示了。

 要做出如上图形，也需要修改一个函数文件，即修改InfoPlot.m文件中的draw3DPoint()函数内容。截图如下：

六 结束语

通过以上操作，研究者只需要根据自己的具体问题，建立数学模型和函数，然后在源代码中简单修改适应度函数和相应的参数，就能立刻运行和声算法程序，跑出结果。而不需要在算法编码和修改上花费太多精力，节省不少时间，我觉得是一件有利于研究的事情。如果用户想加入自己的创意，也可以在原代码基础上进一步修改，提升其功能，那当然更是锦上添花了。