群智能优化算法测试函数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. TreeListControl 不同类别的行使用不同的数据模板

  2. isPostback 的原理及作用（很easy）

    1.IsPostBack用来推断表单是否是回发. (不是第一次请求),是点击表单的提交button回发过来的.是否是回发与get请求还是Post请求无关.可是普通情况下回发都是Post请求. 一般Ge ...

  3. OAuth 2.0 Authorization Framework RFC

    Internet Engineering Task Force (IETF) D. Hardt, Ed.Request for Comments: 6749 MicrosoftObsoletes: 5 ...

  4. [Node.js]31. Level 7: Redis coming for Node.js, Simple Redis Commands

    Let's start practicing using the redis key-value store from our node application. First require the  ...

  5. 8. 利用反射机制， ListArray，intent来实现多Activity的切换

    package com.example.thenewboston; import android.app.ListActivity; import android.content.Intent; im ...

  6. js的正则匹配 和 blur

    <script type="text/javascript" src="http://code.jquery.com/jquery-1.10.2.min.js&qu ...

  7. 能说明你的Javascript技术很烂的五个原因

    Javascript在互联网上名声很臭,但你又很难再找到一个像它这样如此动态.如此被广泛使用.如此根植于我们的生活中的另外一种语言.它的低学习门槛让很多人都称它为学前脚本语言,它另外一个让人嘲笑的东西 ...

  8. nl命令（转）

    原文:http://www.cnblogs.com/peida/archive/2012/11/01/2749048.html nl命令在linux系统中用来计算文件中行号.nl 可以将输出的文件内容 ...

  9. 一个简单的knockout.js 和easyui的绑定

    <!DOCTYPE html><html><head><meta charset="UTF-8"><title>Basi ...

  10. 29、java中阻塞队列

    阻塞队列与普通队列的区别在于,当队列是空的时,从队列中获取元素的操作将会被阻塞,或者当队列是满时,往队列里添加元素的操作会被阻塞.试图从空的阻塞队列中获取元素的线程将会被阻塞,直到其他的线程往空的队列 ...