- 1、本文档共19页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
基于matlab的遗传进化算法示例程序.doc
基于matlab的遗传进化算法示例程序
(1)
,
此问题的全局最优值。
一.程序
(1)主函数: main.m
clear all;
clc;
popsize=60; %种群规模
chromlength=34; %二进制编码,编码精度为0.0001,所以串长l为17
pc=0.7; %杂交概率
pm=0.1; %变异概率
t=0; %进化代数初始为0
pop=initpop(popsize,chromlength); %随机产生初始种群
while t500 %迭代次数
t=t+1;
[objvalue]=calobjvalue(pop); %计算目标函数值
fitvalue=calfitvalue(objvalue); %计算群体中每个个体的适应度
[bestindividual,bestfit]=best(pop,fitvalue); %求出群体中适应度最大的个体及其适应度值
x11=decodechrom(bestindividual,1,14); %将二进制数转换为十进制数
x22=decodechrom(bestindividual,15,14);
x1(t)=-5+10*x11/(pow2(14)-1); %将二值域中的数转换为变量域的数
x2(t)=-5+10*x22/(pow2(14)-1);
y(t)=4*x1(t)^2-2.1*x1(t)^4+1/3*x1(t)^6+x1(t)*x2(t)-4*x2(t)^2+4*x2(t)^4; %计算最佳个体的目标函数值
[newpop1]=selection(pop,fitvalue); %选择算子
[newpop2]=crossover(newpop1,pc); %交叉算子
[newpop3]=mutation(newpop2,pm); %变异算子
objvalue1=calobjvalue(newpop3(1,:));
if objvalue1y(t)
newpop3(1,:)=bestindividual; %保留最佳个体
end
pop=newpop3; %产生新种群
end
y; %每代的最佳目标函数值
x1; %每代的最佳目标函数值对应的自变量
x2;
[gy,k]=min(y) %gy为全局最优值,k为最优值对应的进化代数
gx1=x1(k) %全局最优值对应的自变量
gx2=x2(k)
plot(y) %最优值收敛曲线
title(收敛性曲线);
xlabel(进化代数);
ylabel(函数值);
axis([0,500,-1.5,1.5]);
(2)初始种群:initpop.m
function pop=initpop(popsize,chromlength)
pop=round(rand(popsize,chromlength)); %rand随机产生[0,1]区间的一个小数,rand四舍五入取整
end
(3)计算目标函数值::calobjvalue.m
function [objvalue] =calobjvalue( pop )
temp1=decodechrom(pop,1,14);
temp2=decodechrom(pop,15,14);
x1=-5+(10*temp1)/(pow2(14)-1); %将二值域中的数转化为变量域中的数
x2=-5+(10*temp2)/(pow2(14)-1);
objvalue=4*x1.^2-2.1*x1.^4+1/3*x1.^6+x1.*x2-4*x2.^2+4*x2.^4; %计算目标函数
end
a.二进制转换为十进制:decodechrom.m
function temp=decodechrom(pop,spoint,length )
pop1=pop(:,spoint:spoint+length-1); %按变量个数分组转换,spoint为起始点,length为一个变量的长度
temp=decodebinary(pop1);
end
b.求二进制串对应的十进制数:decodebinary.m
function temp =decodebinary( pop)
[px,py]=size(pop); %求pop行数和列数
for i=1:py
pop1(:,i)=2.^(py-i).*pop(:,i);
end
temp=sum(pop1,2); %每一行求和
end
(4)计算个体适应度:calfitvalue.m
function fitvalue
您可能关注的文档
- 厦门大学生科院微生物细胞生物学生物化学与分子生物学考研信息汇总.doc
- 反义疑问句练习题.doc
- 反义疑问句语法解析练习及答案.doc
- 反刍动物酶制剂的应用及体外评定研究.doc
- 反渗透设备预处理的系统说明.doc
- 发面的过程就是培养酵母菌的过程.doc
- 另一种流感.doc
- 只需1分钟,皮肤就变得白嫩嫩,以及晒黑后变白的方法.doc
- 可遗传变异《染色体变异》复习学案.doc
- 可遗传变异复习学案.doc
- 2024年江西省寻乌县九上数学开学复习检测模拟试题【含答案】.doc
- 2024年江西省省宜春市袁州区数学九上开学学业水平测试模拟试题【含答案】.doc
- 《GB/T 44275.2-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第2部分:术语》.pdf
- 中国国家标准 GB/T 44275.2-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第2部分:术语.pdf
- GB/T 44285.1-2024卡及身份识别安全设备 通过移动设备进行身份管理的构件 第1部分:移动电子身份系统的通用系统架构.pdf
- 《GB/T 44285.1-2024卡及身份识别安全设备 通过移动设备进行身份管理的构件 第1部分:移动电子身份系统的通用系统架构》.pdf
- 中国国家标准 GB/T 44285.1-2024卡及身份识别安全设备 通过移动设备进行身份管理的构件 第1部分:移动电子身份系统的通用系统架构.pdf
- GB/T 44275.11-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第11部分:术语制定指南.pdf
- 中国国家标准 GB/T 44275.11-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第11部分:术语制定指南.pdf
- 《GB/T 44275.11-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第11部分:术语制定指南》.pdf
文档评论(0)