群智能优化算法测试函数matlab源代码

global M;

creatematrix(2);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画ackley图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% ackley x from[-5 5]

% x=-5:0.01:5;

% [x,y]=meshgrid(x);

% temp1=x.^2+y.^2;

% temp2=cos(2*pi*x)+cos(2*pi*y);

% z=20+exp(1)-20*exp(-0.2*sqrt(temp1/2))-exp(temp2/2);

% axis([-5,5,-5,5]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的ackley图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%Rotated ackley x from[-5 5

% x=-5:0.01:5;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% temp1=x.^2+y.^2;

% temp2=cos(2*pi*x)+cos(2*pi*y);

% z=20+exp(1)-20*exp(-0.2*sqrt(temp1/2))-exp(temp2/2);

% axis([-5,5,-5,5]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画cigar图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%cigar

x=-5:0.01:5;

[x,y]=meshgrid(x);

z=x.^2+(10^4)*y.^2;

axis([-5,5,-5,5]);

meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的cigar图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%rotated cigar x from [-5 5]

% x=-7:0.01:7;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% z=x.^2+(10^4)*y.^2;

% axis([-7,7,-7,7]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画ellipse图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%ellipse x from [-5 -5]

% x=-5:0.01:5;

% [x,y]=meshgrid(x);

% z=x.^2+10.^4*(y.^2);

% axis([-5,5,-5,5]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画griewangk图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%griewangk

% x=-600:1:600;

% [x,y]=meshgrid(x);

% temp1=(x.^2+y.^2)/4000;

% temp2=cos(x)*cos(y/sqrt(2));

% z=temp1-temp2+1;

% axis([-600,600,-600,600]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的griewangk图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Rotated griewangk

% x=-600:1:600;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% temp1=(x.^2+y.^2)/4000;

% temp2=cos(x)*cos(y/sqrt(2));

% z=temp1-temp2+1;

% axis([-600,600,-600,600]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画nocon_rastrigin图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%nocon_rastrigin

% x=-10:0.05:10;

% [x,y]=meshgrid(x);

% temp1=0;

% temp2=0;

% z=0;

% if abs(x)<1/2

% temp1=x;

% else

% temp1=round(2.*x)/2;

% end

% if abs(y)<1/2

% temp2=y;

% else

% temp2=round(2.*y)/2;

% end

% z=z+temp1.^2-10*cos(2*pi.*temp1)+temp2.^2-10*cos(2*pi.*temp2)+20;

% axis([-10,10,-10,10]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的nocon_rastrigin图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Rotated nocon_rastrigin

% x=-15:0.05:15;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% temp1=0;

% temp2=0;

% z=0;

% if abs(x)<1/2

% temp1=x;

% else

% temp1=round(2.*x)/2;

% end

% if abs(y)<1/2

% temp2=y;

% else

% temp2=round(2.*y)/2;

% end

% z=z+temp1.^2-10*cos(2*pi.*temp1)+temp2.^2-10*cos(2*pi.*temp2)+20;

% axis([-15,15,-15,15]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画rastrigin图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%rastrigin

% x=-10:0.05:10;

% [x,y]=meshgrid(x);

% temp1=x.^2-10*cos(2*pi.*x);

% temp2=y.^2-10*cos(2*pi.*y);

% z=temp1+temp2+20;

% axis([-10,10,-10,10]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的rastrigin图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Rotated rastrigin

% x=-15:0.05:15;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% temp1=x.^2-10*cos(2*pi.*x);

% temp2=y.^2-10*cos(2*pi.*y);

% z=temp1+temp2+20;

% axis([-15,15,-15,15]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画rosenbrock图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%rosenbrock

% x=-10:0.05:10;

% [x,y]=meshgrid(x);

% z=100*(y-x.^2).^2+(x-1).^2;

% axis([-10,10,-10,10]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画schwefel图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%schwefel

% x=-10:0.05:10;

% [x,y]=meshgrid(x);

% z=(x-x.^2).^2+(x-1).^2+(x-y.^2).^2+(y-1).^2;

% axis([-10,10,-10,10]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的schwefel图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Rotated schwefel

% x=-15:0.05:15;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% z=(x-x.^2).^2+(x-1).^2+(x-y.^2).^2+(y-1).^2;

% axis([-15,15,-15,15]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画sphere图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%sphere

% x=-10:0.05:10;

% [x,y]=meshgrid(x);

% z=x.^2+y.^2;

% axis([-10,10,-10,10]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画sumcan图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%sumcan

% x=-0.16:0.0005:0.16;

% [x,y]=meshgrid(x);

% z=1./(10.^(-5)+abs(x)+abs(x+y));

% z=log10(z);

% axis([-0.16,0.16,-0.16,0.16]);

% mesh(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画tablet图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%tablet

% x=-10:0.05:10;

% [x,y]=meshgrid(x);

% z=10.^4*x.^2+y.^2;

% axis([-10,10,-10,10]);

% mesh(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画Weierstrass图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Weierstrass

% x=-1:0.005:1;

% [x,y]=meshgrid(x);

% z=0;

% for k=0:20

% z=z+0.5^k*cos(2*pi*3^k*(x+0.5))+0.5^k*cos(2*pi*3^k*(y+0.5));

% end

% temp=0;

% for k=0:20

% temp=temp+0.5^k*cos(2*pi*3^k*0.5);

% end

% z=z-2*temp;

% axis([-1,1,-1,1]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的Weierstrass图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Rotated Weierstrass

% x=-1.5:0.005:1.5;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% z=0;

% for k=0:20

% z=z+0.5^k*cos(2*pi*3^k*(x+0.5))+0.5^k*cos(2*pi*3^k*(y+0.5));

% end

% temp=0;

% for k=0:20

% temp=temp+0.5^k*cos(2*pi*3^k*0.5);

% end

% z=z-2*temp;

% axis([-1.5,1.5,-1.5,1.5]);

% meshc(x,y,z);

%

xlabel('x');

ylabel('y');

zlabel('f(x,y)');

colormap jet;

% clear;

%saveas(gcf,'temp.jpg')

　　

群智能优化算法-测试函数matlab源码的更多相关文章

  1. GWO（灰狼优化）算法MATLAB源码逐行中文注解（转载）

    以优化SVM算法的参数c和g为例,对GWO算法MATLAB源码进行了逐行中文注解. tic % 计时器 %% 清空环境变量 close all clear clc format compact %% ...

  2. OpenCV学习笔记（27）KAZE 算法原理与源码分析（一）非线性扩散滤波

    http://blog.csdn.net/chenyusiyuan/article/details/8710462 OpenCV学习笔记(27)KAZE 算法原理与源码分析(一)非线性扩散滤波 201 ...

  3. 模拟退火算法SA原理及python、java、php、c++语言代码实现TSP旅行商问题，智能优化算法，随机寻优算法，全局最短路径

    模拟退火算法SA原理及python.java.php.c++语言代码实现TSP旅行商问题,智能优化算法,随机寻优算法,全局最短路径 模拟退火算法(Simulated Annealing,SA)最早的思 ...

  4. iOS下使用SHA1WithRSA算法加签源码

    首先了解一下几个相关概念,以方便后面遇到的问题的解决: RSA算法:1977年由Ron Rivest.Adi Shamirh和LenAdleman发明的,RSA就是取自他们三个人的名字.算法基于一个数 ...

  5. 智能优化算法对TSP问题的求解研究

    要求: TSP 算法(Traveling Salesman Problem)是指给定 n 个城市和各个城市之间的距离,要 求确定一条经过各个城市当且仅当一次的最短路径,它是一种典型的优化组合问题,其最 ...

  6. 百度智能手环方案开源(含源码，原理图，APP，通信协议等)

    分享一个百度智能手环开源项目的设计方案资料. 项目简介 百度云智能手环的开源方案是基于Apache2.0开源协议,开源内容包括硬件设计文档,原理图.ROM.通讯协议在内的全套方案,同时开放APP和云服 ...

  7. java实现插入排序算法 附单元测试源码

    插入排序算法 public class InsertSortTest { /** * @param args */ public static void main(String[] args) { i ...

  8. HashMap与HashTable的哈希算法——JDK1.9源码阅读总结

    下面是HashTable源码中的put方法: 注意上面注释标注的地方: HashTable对于元素在哈希表中的坐标算法是: 将对象自身的哈希值key.hashCode()变为正数:hash & ...

  9. 推荐算法_CIKM-2019-AnalytiCup 冠军源码解读_2

    最近在为机器学习结合推荐算法的优化方法和数据来源想办法.抱着学习的态度继续解读19-AnalytiCup的冠军源码. 第一部分itemcf解读的连接:https://www.cnblogs.com/m ...

随机推荐

  1. php防止sql注入漏洞代码 && 几种常见攻击的正则表达式

    注入漏洞代码和分析                                                                                           ...

  2. Android实现圆形的图片边角

    ImageView没有边角圆化的设置,但是可以通过代码设置实现.一个国外的码农写的方法. 效果: 地址:http://www.curious-creature.org/2012 代码: package ...

  3. 转: 用 Eclipse 平台进行 C/C++ 开发

    http://www.ibm.com/developerworks/cn/linux/opensource/os-ecc/index.html

  4. com.esotericsoftware.kryo.kryoexception java.util.ConcurentModificationException

    近期 有网友看我的"整合Kafka到Spark Streaming--代码演示样例和挑战"文章, 讲 kafka对象 放到 pool 并通过broadcast广播出去: 然后 在开 ...

  5. angular cli中的--flat参数

    添加--flat 参数,表示不创建单独的目录 例如 ng g c news/news-detail 会将文件创建在app/news/news-detail目录下 ng g c news/news-de ...

  6. SQLServer 之 树查询

    一.SqlServer树查询 1.使用公用表表达式(CTE) 很多人可能想要查询整个树形表关联的内容都会通过循环递归来查...事实上在微软在SQL2005或以上版本就能用别的语法进行查询,下面是示例. ...

  7. ORACLE常用监控语句(未完待续)

    --查询日志的切换频率 select  t1.RECID as srecid        ,t2.RECID as erecid        ,t1.FIRST_TIME as stime     ...

  8. php 内置函数JSON处理

    一:生成json数据 <?php $arr = array ('a'=>1,'b'=>2,'c'=>3,'d'=>4,'e'=>5); echo json_enco ...

  9. poj 3345 Bribing FIPA (树形背包dp | 输入坑)

    题目链接:  poj-3345  hdu-2415 题意 有n个国家,你要获取m个国家的支持,获取第i个国家的支持就要给cost[i]的价钱    其中有一些国家是老大和小弟的关系,也就是说,如果你获 ...

  10. PHP高级教程-多维数组

    PHP 多维数组 一个数组中的值可以是另一个数组,另一个数组的值也可以是一个数组.依照这种方式,我们可以创建二维或者三维数组: 实例 <?php // 二维数组: $cars = array ( ...