function [xv,fv] = myGA(fitness, a, b, NP, NG, Pc, Pm, eps)
% 用遗传算法求解一维无约束优化问题
%
% 待优化的目标函数 fitness
% 自变量下界 a
% 自变量上界 b
% 种群个体数 NP
% 最大进化代数 NG
% 杂交概率 Pc
% 变异概率 Pm
% 自变量离散精度 eps
% 目标变量取最大值时自变量的值: xm
% 目标函数的最大值 fv
%
% Example:
% function F = fitness(x)
% F = x^3-60*x^2+900*x+100;
% -------------------------------
% [xv,fv] = myGA(@fitness,0, 30, 50, 100, 0.9, 0.04, 0.01);
% --------------------------------------------------
% xv = 10
% fv = 4100
%
% 本程序在《精通MATLAB最优化计算》页315程序的基础上修改 L = ceil(log2((b-a) / eps + 1)); %编码长度
x = zeros(NP, L); %种群
nx = zeros(size(x)); %滚动数组
fx = zeros(NP, 1); %适应度
for i = 1:NP
x(i,:) = Initial(L);
end fv = -inf; for k = 1 : NG
for i = 1 : NP
fx(i) = fitness(Dec(a, b, x(i, :), L));
if (fx(i) > fv)
xv = Dec(a, b, x(i, :), L);
fv = fx(i);
end
end sumfx = sum(fx);
Px = fx / sumfx; PPx = zeros(NP, 1);
PPx(1) = Px(1); %概率叠加
for i = 2 : NP
PPx(i) = PPx(i - 1) + Px(i);
end selFather = 0;
for i = 1 : NP
sita = rand();
for j = 1 : NP
if (sita <= PPx(j))
selFather = j; %使用轮盘赌法进行选择父亲
break;
end
end selMother = floor(rand() * NP) + 1; %母亲随机选择
posCut = floor(rand() * (L - 1)) + 1; %交叉点 r1 = rand();
if (r1 <= Pc)
nx(i, 1 : posCut) = x(selFather, 1:posCut);
nx(i, (posCut + 1) : L) = x(selMother, (posCut + 1) : L);
r2 = rand();
if (r2 <= Pm)
posMut = floor(rand() * L) + 1;
nx(i, posMut) = ~nx(i, posMut);
end
else
nx(i, :) = x(selFather, :);
end
end x = nx;
end %--------------------------------------------------------
% 初始化种群
function result = Initial(length)
result = zeros(size(length()));
for i = 1 : length
r = rand();
result(i) = round(r);
end %----------------------------------------------------------
% 编码转换
function y = Dec(a, b, x, L)
base = 2 .^ ((L - 1) : -1: 0);
y = dot(base, x);
y = a + y * (b - a) / (2 ^ L - 1);

matlab遗传算法的更多相关文章

  1. matlab遗传算法工具箱

    转自http://blog.sina.com.cn/s/blog_5ebcc0240101pnrj.html matlab遗传算法工具箱函数及实例讲解 (2014-01-10 13:03:57)   ...

  2. Matlab遗传算法优化问题求解的演示样例代码

    代码例如以下: function m_main() clear clc Max_gen = 100;% 执行代数 pop_size = 100;%种群大小 chromsome = 10;%染色体的长度 ...

  3. 遗传算法MATLAB实现(2):一元函数优化举例

    遗传算法提供了一种求解非线性.多模型.多目标等复杂系统优化问题的通用框架. 先从例子开始,慢慢再总结理论... [例]利用遗传算法计算函数f(x)=x*cos(5*pi*x)+3.5在区间[-1,2. ...

  4. matlab中文论坛视频谷普教程MATLAB压缩包介绍

    matlab中文论坛视频谷普教程MATLAB压缩包介绍 我也正在学习这个软件 ,看到这个教程就在这里分享了,希望大家喜欢!Matlab 初学者视频教学1. Matlab视频:Matlab中文论坛为新手 ...

  5. 史上最全的Matlab资源电子书教程和视频下载合集【超级推荐】

    收藏吧,网上搜集的,费了老大劲了,推荐给有需要的人,^_^.   MATLAB课件2007北京交通大学.zip 4.87 MB   A Guide to MATLAB for Beginners an ...

  6. 一些matlab教程资源收藏,使用matlab编程的人还是挺多的

    Matlab教程专题资源免费下载整理合集收藏 <MATLAB从入门到精通>高清文字版[PDF] 103.9MB 简体中文 <矩阵实验室>(Mathworks.Matlab.R2 ...

  7. geatpy - 遗传和进化算法相关算子的库函数(python)

    Geatpy The Genetic and Evolutionary Algorithm Toolbox for Python Introduction Website (including doc ...

  8. Python遗传和进化算法框架(一)Geatpy快速入门

    https://blog.csdn.net/qq_33353186/article/details/82014986 Geatpy是一个高性能的Python遗传算法库以及开放式进化算法框架,由华南理工 ...

  9. MATLAB中如何使用遗传算法

    matlab有遗传算法工具箱. 核心函数:   (1)function [pop]=initializega(num,bounds,eevalFN,eevalOps,options)--初始种群的生成 ...

随机推荐

  1. resin config 中文(resin.xml)

    <!-- - Resin 3.1 配置文件. --> <resin xmlns="http://caucho.com/ns/resin" xmlns:resin= ...

  2. Jdbc 事务

    package com.j1; import java.sql.Connection; import java.sql.SQLException; import com.mysql.jdbc.Prep ...

  3. Beacon浅析

    作者:hongbosun 一.Beacon简介 Beacon是基于BLE技术实现的物理设备.BLE(全称Bluetooth Low Energy)是蓝牙4.0技术规范的一部分.它起源于Nokia的Wi ...

  4. centos静默式安装Oracle11g

    1.   Centos及Oracle版本 Centos:CentOS release 6.4 (Final) Oracle:linux.x64_Oracle_11gR2_database 2.   硬 ...

  5. 爆炸!iOS资源大礼包(持续更新...)

    今天为大家整理了一些关于iOS学习的干货,献给正在奋斗的你们,首先声明一下,在整理的过程中参考了大量的博客和文章,知识的分享终究会增值,在此表示感谢,希望这篇文章给大家带来帮助. 基础部分: C语言教 ...

  6. canvas ---1

    Canvas1 (关键词:canvas)   canvas :就是html5中提供的一个标签,只是用来展示绘图的内容 canvas  标签的默认宽高:300*150 如果给canvas来设置高度和宽度 ...

  7. easyui之datagrid(不定时补充)

    1,datagrid之formatter formatter格式化函数有三个参数: value:字段值(一般为后台传递给前台的值): row:当前行数据: index:当前行索引. return值是显 ...

  8. MySQL 使用索引扫描来做排序

    MySQL有两种方式可以生成有序的结果:通过排序操作:或者按照索引顺序扫描:如果EXPLAIN 出来的结果的type列的值为“index”,则说明MySQL使用了索引扫描来做排序(不要和Extra列的 ...

  9. 3月23日html(五) 格式与布局练习:360浏览器布局

    <meta http-equiv="Content-Type" content="text/html; charset=utf-8" /> < ...

  10. PHP学习系列(1)——字符串处理函数(1)

    从09年开始使用JAVA到现在差不多5年多了,然后由于即将要去的公司(研究僧终于要毕业了!)是使用PHP,Python作为后台开发语言的,所以要开始一段双P的学习旅程.用过PHP的都说这是很简单的一门 ...