和声搜索算法

和声搜索算法，一种智能优化算法，国外学者编写的源码，可以运行，有测试函数和算法的简要介绍。

% Harmony Search Algorithm
%By Sajjad Yazdani
%
% base on:
% [1]:
clc;clear all;close all
%% Problem Prametters
Dim=2; % problem Dimention
Low=[-100 -100]; % Low Boundry of Problem
High=[100 100]; % High Boundry of Problem

Min=1; % Minimaization or maximaiz of Fun? if Min=1 it will be minimaze the function and if Min=0 it will be maximized the function.

%% Harmony Search Parametters

HMS=100;%Harmony Memory Size （Population Number）
bw=0.2;
HMCR=0.95;%[1] Harmony Memory Considering Rate
PAR=0.3;%[1] Pitch Adjustment Rate

MaxItr=10000;% Maximum number of Iteration

%% Initialization
HM=zeros（HMSDim）;
HF=zeros（HMS1）;
for i=1:HMS
    HM（i:）=Low+（High-Low）.*rand（1Dim）;
    HF（i1）=MyFun（HM（i:））;
end

if Min==1
    [WorstFitWorstLoc]=max（HF）;
else
    [WorstFitWorstLoc]=min（HF）;
end


%% Iteration Loop
for Itr=1:MaxItr
    HarmonyIndex=fix（rand（1Dim）*HMS）+1;% Random Selection of Harmony
    Harmony=diag（HM（HarmonyIndex1:Dim））‘;% Extraxt Value of harmony from Memory（Can Be better???）
    CMMask=rand（1Dim）    NHMask=（1-CMMask）;
    PAMask=（rand（1Dim）    CMMask=CMMask.*（1-PAMask）;
    NewHarmony=CMMask.*Harmony+PAMask.*（Harmony+bw*（2*rand（1Dim）-1））+NHMask.*（Low+（High-Low）.*rand（1Dim））;
    OutOfBoundry=（NewHarmony>High）+（NewHarmony    NewHarmony（OutOfBoundry==1）=Harmony（OutOfBoundry==1）;
    NHF=MyFun（NewHarmony）;
    if （NHF        HM（WorstLoc:）=NewHarmony;
        HF（WorstLoc）=NHF;
        [WorstFitWorstLoc]=max（HF）;
    elseif （NHF        HM（WorstLoc:）=NewHarmony;
        HF（WorstLoc）=NHF;
        [WorstFitWorstLoc]=min（HF）;
    end
end

%% Present Best Answer
if Min==1
    [BestFitBestLoc]=min（HF）;
else
    [BestFitBestLoc]=max（HF）;
end
Best=HM（BestLoc:）;

display（Best）
display（BestFit）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    1124324  2009-12-29 02:14  HS\Harmony Search Algorithm .pdf

     文件       1877  2014-04-03 17:30  HS\HS.m

     文件         72  2014-04-03 17:26  HS\MyFun.m

     文件         73  2014-04-03 17:27  HS\Sphere.m

     目录          0  2014-12-05 21:07  HS

----------- ---------  ---------- -----  ----

              1126346                    5