微分进化算法解决函数优化问题的matlab代码,Matlab微分进化算法及优化函数测试...

微分进化(Difference Evolution,DE)算法是一种优化算法,据称其比GA(遗传算法)等更为优秀。

借鉴网上实现的DE算法,用Matlab实现了对若干函数优化问题的解法,代码如下:

function [] = de_test

clear all;

close all;

clc;

%解:X = [0, 0, ...]

%nVar = 30

%dims: [-30, 30]

function fitness = sphere(vals)

prod = vals .* vals;

fitness = sum(prod, 2);

end

%f(x) = 1/4000 * sum^n_1(x_i)^2 - prod^n_1 * cos(x_i/sqrt(i)) + 1

%f* = 0, x* = [0, 0, ...];

%nVar = 30

%dims: [-600, 600]

function fitness = griewank(vals)

[h w] = size(vals);

p1 = vals.^2;

p1 = 1/4000 .* sum(p1, 2);

t = sqrt([1:w]);

p2 = vals ./ repmat(t, h, 1);

p2 = cos(p2);

p2 = prod(p2, 2);

fitness = p1 - p2 + 1;

end

%解: X = [1, 1, 1, ...]

function fitness = RosenBroek(vals)

[k n] = size(vals);

vals1 = vals(:, 1:n-1);

vals2 = vals(:, 2:n);

tmp = vals2 - vals1 .* vals1;

tmp = 100 * tmp .* tmp;

tmp1 = (vals1 - 1).^2;

fitness = sum(tmp, 2) + sum(tmp1, 2);

end

%许多局部最小值,最优解:X = [0, 0], 固定2个变量

function fitness = Rastrigin (X)

v = X.^2 - 10 .* cos(2 * pi .* X) + 10;

fitness = sum(v, 2);

%fitness = 20 + X(:, 1).^2 + X(:, 2).^2 - 10 .* (cos(2 * pi * X(:, 1)) + cos(2 * pi * X(:, 2)));

end

%参数

popSize = 100; % 群规模

F = 0.5;                    % 突变因子

C = 0.9;                    % 交叉率

iterTimes = 300;% 迭代次数

Run = 4;        % 测试轮次

id = [1:popSize];

%结果

bestFit = 1e5;  %设取最小值

bestGene = [];

for r = 1:Run  %轮次

Func = r;     % 测试适应度函数

switch Func

case 1

Func = @sphere;

Xmin = -100;

Xmax = 100;

nVar = 30;

case 2

Func = @griewank;

Xmin = -600;

Xmax = 600;

nVar = 2;

case 3

nVar = 2;

Func = @RosenBroek;

Xmin = -100;

Xmax = 100;

case 4

nVar = 2;

Func = @Rastrigin

Xmin = -5.12;

Xmax = 5.12;

end;

Func

tic;

%1.初始化种群

X = Xmin + rand(popSize, nVar) * Xmax;

%2.每轮迭代

for i = 1:iterTimes

X0 = X;

%             F = 2 * (1 - (i-1)/iterTimes);

%个体突变得到V

for j = 1:popSize

ids = id;

ids(j) = [];

rids = randperm(popSize - 1, 3);

rids = ids(rids);

V = X(rids(1), :) + F * (X(rids(2), :) - X(rids(3), :));

%对V(j, :)值域检查

ids = find(V < Xmin);

if length(ids) > 0

V(ids) = Xmin;

end;

ids = find(V > Xmax);

if length(ids) > 0

V(ids) = Xmax;

end;

%对每个X和V的配对,进行交叉操作,结果存在U中

jrand = floor(rand() * nVar + 1);  %必交叉项,保证至少一个交叉

for n = 1:nVar

R1 = rand();

if (R1 < C || n == jrand)

U = V;    %保留子代基因

else

U = X(j, :);    %保留父代基因

end

end

%在子代和父代间做选择运算

if Func(U) < Func(X(j, :))

Tr = U;

else

Tr=X(j, :);

end

% 更新种群基因

X(j,:) = Tr;

%计算新的适应度

fitness = Func(X(j,:));

%记录全局最优解

if fitness < bestFit

bestFit = fitness;

bestGene = X(j, :);

end;

end  %结束个体更新

end

bestGene

bestFit

toc

end

end运行结果如下:

Func =

@de_test/sphere

bestGene =

Columns 1 through 15

-7.5570  -11.9099    9.9957  -37.8403  -17.9445   -7.1438  -21.4304  -33.6260  -22.1812  -66.1438    0.9014  -39.6724   -0.3175  -56.4815   26.5422

Columns 16 through 30

6.5446  -31.9653   -9.3640  -37.1629  -23.7325    2.7271   -6.3413  -21.3204  -13.1450   28.7402  -28.8170  -22.6226   10.9031  -16.6128  -14.7637

bestFit =

2.0583e+04  (从此结果看,算法效果不佳)

时间已过 2.285354 秒。

Func =

@de_test/griewank

bestGene =

1.0e-08 *

0.1702    0.1424

bestFit =

0

时间已过 3.337347 秒。

Func =

@de_test/RosenBroek

bestGene =

1.0e-08 *

0.1702    0.1424

bestFit =

0

时间已过 1.747765 秒。

Func =

@de_test/Rastrigin

Func =

@de_test/Rastrigin

bestGene =

1.0e-08 *

0.1702    0.1424

bestFit =

0

时间已过 1.403871 秒。

原文:http://blog.csdn.net/miscclp/article/details/38344685

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/529500.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Oracle中获取文件中的数据,操作oracle中的数据文件

收藏于http://dev.yesky.com/296/8090796.shtmlOracle数据库中管理表空间和数据文件数据库的三大组成部分:数据文件&#xff0c;控制文件&#xff0c;Redo日志。表空间分为系统表空间和非系统表空间。SYSTEM表空间是最基本的&#xff0c;必须的&#xff0c;Oracle建议你为用户数…

oracle 存储过程设置回滚点,(转)oracle 存储过程事宜使用断点回滚 -savepoint

学习存储过程中使用断点回滚事务时&#xff0c;发现目前网络上存在一个问题&#xff0c;那就是使用断点回滚后&#xff0c;都忘记了一个很重要的事情&#xff0c;提交事务。虽然使用了断点回滚&#xff0c;但是断点回滚不像rollBack或commit一样结束当前事务&#xff0c;而使用…

php tp5清空数据表并主键,tp5数据库——更新数据

更新数据更新数据表中的数据Db::table(think_user)->where(id, 1)->update([name > thinkphp]);如果数据中包含主键&#xff0c;可以直接使用&#xff1a;Db::table(think_user)->update([name > thinkphp,id>1]);update 方法返回影响数据的条数&#xff0c;…

oracle几个网络,ORACLE网络的几个重点概念

数据库名(DB_NAME)、数据库实例名(INSTANCE_NAME)、操作系统环境变量ORACLE_SID、数据库服务名 (SERVICE_NAME)、数据库域名(DB_DOMAIN)以及全局数据库名(GLOBAL_DB_NAME)是几个使用Oracle数据库容易混 淆的概念。1、数据库名与实例名数据库名(DB_NAME)&#xff1a;是区分数据的…

oracle每季度补丁,Oracle 2020 年第四季度补丁发布

半个月前&#xff0c;也就是 10 月 20 日&#xff0c; Oracle 发布了今年最后一次补丁更新&#xff0c;那么很多人都想打最新的 PSU&#xff0c;理由是有被扫到各种漏洞&#xff0c;有的扫描工具着实太坑&#xff0c;这里就不用说了&#xff0c;前几天看到盖总发布的文章2020年…

php flash上传进度条,PHP_PHP+FLASH实现上传文件进度条相关文件 下载,PHP之所以很难实现上传进度条 - phpStudy...

PHPFLASH实现上传文件进度条相关文件 下载PHP之所以很难实现上传进度条是因为在我们上传文件到服务器的时候&#xff0c;要等到文件全部送到服务器之后&#xff0c;才执行相应的php文件。在这之前&#xff0c;文件数据保存在一个临时文件里面&#xff0c;而php无法获得这个文件…

oracle截取时间函数6,SQL 取日期部分函数

convert(varchar(10),getdate(),120)这个函数就可以得到日期部分了&#xff0c;我觉得还是比较方便的&#xff0c;我觉得SQL在日期方面&#xff0c;没有ORACLE函数坐的好。。Sql Server中的日期与时间函数1. 当前系统日期、时间select getdate()2. dateadd 在向指定日期加上一…

oracle 找不到程序,Oracle 找不到监听程序

Oracle服务器文件 app\Administrator\product\11.2.0\dbhome_1\NETWORK\ADMIN\listener.oraSID_LIST_LISTENER (SID_LIST (SID_DESC (SID_NAME CLRExtProc)(ORACLE_HOME D:\app\Administrator\product\11.2.0\dbhome_1)(PROGRAM extproc)(ENVS "EXTPROC_DLLSONLY:D:\…

oracle中的存储过程教程,oracle 存储过程

-- 查看系统定时任务SELECT * FROM DBA_JOBS-- 新建定时任务DECLAREjobno NUMBER;BEGINSYS.DBMS_JOB.SUBMIT(JOB > jobno, /*自动生成JOB_ID*/WHAT > qmcb_ls_data;, /*需要执行的存储过程名称或SQL语句*/NEXT_DATE > SYSDATE, /*初次执行时间-立即执行*/INTERVAL >…

Oracle segment啥意思,关于oracle数据库段segment的小结

段(segment)是一种在数据库中消耗物理存储空间的任何实体(一个段可能存在于多个数据文件中&#xff0c;因为物理的数据文件是组成逻辑表空间的基本物理存储单位)今天碰到一个高水位问题&#xff1a;一个分区表&#xff0c;删除某个分区中的一些数据后&#xff0c;表空间并没有回…

Linux在文件第一行添加字符,使用shell命令给文件中每一行的前面、后面、指定列添加字符...

shell给一个文件中的每一行开头插入字符的方法&#xff1a;awk {print "X"$0} fileNameshell给一个文件中的每一行结尾插入字符的方法&#xff1a;awk {print $0"X"} fileNameshell给一个文件中的每一行的指定列插入字符的方法&#xff1a;awk $O$O" …

路由器 刷 linux系统版本,在Linux下用tftp刷写路由器固件

(以Buffalo WHR-G300N V2路由器为例)以发行版Ubuntu为例(如果你在Windows下&#xff0c;可下载Ubuntu的ISO文件&#xff0c;再用wubi安装Ubuntu&#xff0c;可免去重新分区的麻烦)&#xff0c;下载Buffalo WHR-G300N V2路由器的FTP版固件文件&#xff0c;放到Ubuntu的/home目录…

linux没有etho网卡,关于怎么解决CENTOS7没有ETH0网卡这个问题

CentOS7系统安装完毕之后&#xff0c;输入ifconfig命令发现没有eth0,不符合我们的习惯。而且也无法远程ssh连接。1.进入目录/etc/sysconfig/network-scripts/2.将文件ifcfg-ens33重命名为ifcfg-eth0;(注意&#xff1a;修改需要切换至root用户)3.编辑ifcfg-eth0使用命令&#xf…

linux添加物理卷编辑文件夹,Red hat Linux下的逻辑卷管理器LVM-上

【IT168 专稿】Red hat 下的LVM 上 LVM是Logical Volume Manager(逻辑卷管理器)的简写&#xff0c;它为主机提供了更高层次的磁盘存储管理能力。LVM可以帮助系统管理员为应用与用户方便地分配存储空间。在LVM管理下的逻辑卷可以按需改变大小或添加移除。另外&#xff0c;LVM可以…

linux的nfs端口号,#Linux NFS服务 固定端口及防火墙配置#

#Linux NFS服务 固定端口及防火墙配置#1.在Linux上正常安装NFS服务2.修改/etc/service,添加以下内容(端口号必须在1024以下&#xff0c;且未被占用)# Local servicesmountd 1011/tcp #rpc.mountdmountd 1011/udp #rpc.mountdrquotad 1012/tcp #rpc.rquotadrquotad 1012/udp #rp…

嵌入式linux添加环境变量,Linux环境变量

前言&#xff1a;1. 环境变量是操作系统环境设置的变量&#xff0c;适用于整个系统的用户进程&#xff1b;2. 环境变量可以在命令中设置&#xff0c;但是用户注销的时候将会丢失这些设置值&#xff1b;3. 若要重复适用&#xff0c;则最好在.profile中定义&#xff1b;环境变量的…

eval函数linux,Python中的eval()、exec()及其相关函数

刚好前些天有人提到eval()与exec()这两个函数&#xff0c;所以就翻了下Python的文档。这里就来简单说一下这两个函数以及与它们相关的几个函数&#xff0c;如globals()、locals()和compile()&#xff1a;1. eval函数函数的作用&#xff1a;计算指定表达式的值。也就是说它要执行…

linux更改语言脚本,Linux shell脚本入门——shell语言脚本【CentOS】

认识脚本是使用一种特定的描述性语言&#xff0c;依据一定的格式编写的可执行文件。脚本语言又被称为扩建的语言, 或者动态语言, 是一种编程语言, 用来控制软件应用程序, 脚本通常是以文本 (ASCⅡ) 保存, 只是在被调用时进行解释或者编译。学习shell脚本的用途对于一个合格的系…

linux可用的ftp,linuxunix下有很多可用的ftp服务器

一、ftp服务器简介linux/unix下有很多可用的ftp服务器&#xff0c;根据这些ftp服务器的可配置性大概可以分为3类&#xff1a;弱、中等、高。功能比较简单的ftp服务器软件有“ftpd”和“oftpd”。前者与ftp客户端工具“ftp”类似&#xff0c;只实现了标准的ftp功能&#xff0c;支…

linux视图版怎么输入命令,分享在Linux命令下操作MySQL视图实例代码

视图VIEW命令简介&#xff1a;VIEW视图就是存储下来的SELECT语句数据1。创建视图命令格式&#xff1a;创建[或更换][ALGORITHM {UNDEFINED | MERGE | 不是Temptable}][DEFINER {user | 当前用户 }][SQL SECURITY {DEFINER | INVOKER}]VIEW view_name [(column_list)]AS selec…