已阅读5页,还剩5页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
附页:一遗传算法源程序:clc;clear;population;%评价目标函数值for uim=1:popsize vector=population(uim,:); obj(uim)=hanshu(hromlength,vector,phen);end%obj%min(obj)clear uim;objmin=min(obj);for sequ=1:popsize if obj(sequ)=objmin opti=population(sequ,:); endendclear sequ;fmax=22000;%=for gen=1:maxgen%选择操作%将求最小值的函数转化为适应度函数for indivi=1:popsize obj1(indivi)=1/obj(indivi);endclear indivi;%适应度函数累加总合total=0;for indivi=1:popsize total=total+obj1(indivi);endclear indivi;%每条染色体被选中的几率for indivi=1:popsize fitness1(indivi)=obj1(indivi)/total;endclear indivi;%各条染色体被选中的范围for indivi=1:popsize fitness(indivi)=0; for j=1:indivi fitness(indivi)=fitness(indivi)+fitness1(j); endendclear j;fitness;%选择适应度高的个体for ranseti=1:popsize ran=rand; while (ran1|ran0) ran=rand; end ran; if ranfitness(fet-1)&(ran=fitness(fet) newpopulation(ranseti,:)=population(fet,:); end end endendclear ran;newpopulation;%交叉for int=1:2:popsize-1 popmoth=newpopulation(int,:); popfath=newpopulation(int+1,:); popcross(int,:)=popmoth; popcross(int+1,:)=popfath; randnum=rand; if(randnum cpoint1=round(rand*hromlength); cpoint2=round(rand*hromlength); while (cpoint2=cpoint1) cpoint2=round(rand*hromlength); end if cpoint1cpoint2 tem=cpoint1; cpoint1=cpoint2; cpoint2=tem; end cpoint1; cpoint2; for term=cpoint1+1:cpoint2 for ss=1:hromlength if popcross(int,ss)=popfath(term) tem1=popcross(int,ss); popcross(int,ss)=popcross(int,term); popcross(int,term)=tem1; end end clear tem1; end for term=cpoint1+1:cpoint2 for ss=1:hromlength if popcross(int+1,ss)=popmoth(term) tem1=popcross(int+1,ss); popcross(int+1,ss)=popcross(int+1,term); popcross(int+1,term)=tem1; end end clear tem1; end end clear term;endclear randnum;popcross;%变异操作newpop=popcross;for int=1:popsize randnum=rand; if randnum cpoint12=round(rand*hromlength); cpoint22=round(rand*hromlength); if (cpoint12=0) cpoint12=1; end if (cpoint22=0) cpoint22=1; end while (cpoint22=cpoint12) cpoint22=round(rand*hromlength); if cpoint22=0; cpoint22=1; end end temp=newpop(int,cpoint12); newpop(int,cpoint12)=newpop(int,cpoint22); newpop(int,cpoint22)=temp; endendnewpop;clear cpoint12;clear cpoint22;clear randnum;clear int;for ium=1:popsize vector1=newpop(ium,:); obj1(ium)=hanshu(hromlength,vector1,phen);endclear ium;obj1max=max(obj1);for ar=1:popsize if obj1(ar)=obj1max newpop(ar,:)=opti; endend%遗传操作结束二粒子群算法源程序:%-初始格式化-clear all;clc;format long;%-给定初始化条件-c1=1.4962;%学习因子1c2=1.4962;%学习因子2w=0.7298;%惯性权重MaxDT=100;%最大迭代次数D=2;%搜索空间维数(未知数个数)N=40;%初始化群体个体数目eps=10(-6);%设置精度(在已知最小值时候用)%-初始化种群的个体(可以在这里限定位置和速度的范围)-for i=1:N for j=1:D x(i,j)=randn;%随机初始化位置 v(i,j)=randn;%随机初始化速度 endend%-先计算各个粒子的适应度,并初始化Pi和Pg-for i=1:N p(i)=fitness(x(i,:),D); y(i,:)=x(i,:);endpg=x(1,:);%Pg为全局最优for i=2:N if fitness(x(i,:),D)FITNESS(pg,D) pg=x(i,:); endend%-进入主要循环,按照公式依次迭代,直到满足精度要求-for t=1:MaxDT t for i=1:N v(i,:)=w*v(i,:)+c1*rand*(y(i,:)-x(i,:)+c2*rand*(pg-x(i,:); x(i,:)=x(i,:)+v(i,:); if fitness(x(i,:),D)p(i) p(i)=fitness(x(i,:),D); y(i,:)=x(i,:); end if p(i)FITNESS(pg,D) pg=y(i,:); end end Pbest(t)=fitness(pg,D);end%-进入主要循环,按照公式依次迭代,直到满足精度要求-for t=1:MaxDT for i=1:N v(i,:)=w*v(i,:)+c1*rand*(y(i,:)-x(i,:)+c2*rand*(pg-x(i,:); x(i,:)=x(i,:)+v(i,:); if fitness(x(i,:),D)p(i) p(i)=fitness(x(i,:),D); y(i,:)=x(i,:); end if p(i)FITNESS(pg,D) pg=y(i,:); end end Pbest(t)=fitness(pg,D);end%-最后给出计算结果disp(*)disp(函数的
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 山东海事职业学院《食品毒理学》2023-2024学年第一学期期末试卷
- 山东工艺美术学院《电路分析原理》2023-2024学年第一学期期末试卷
- 山东工业职业学院《建筑认识与表达》2023-2024学年第一学期期末试卷
- 文化课程设计总结
- 幼儿服装课程设计
- 安装电气课程设计
- 幼儿芦苇课程设计
- 《挑战性压力源与阻碍性压力源对员工满意度的影响研究》
- 《坎地沙坦酯混合胶束片制备工艺的研究》
- 《我国文化产业税收负担研究》
- 2023年中考英语备考让步状语从句练习题(附答案)
- DL∕T 5028.1-2015 电力工程制图标准 第1部分:一般规则部分
- Unit 6 Section A 课件 人教版2024七年级英语上册
- 江西省物业管理服务收费办法
- 高级临床药学实践概论智慧树知到期末考试答案章节答案2024年沈阳药科大学
- 2024年人教版小学五年级信息技术(下册)期末试卷附答案
- 北京市体育竞赛管理和国际交流中心招聘笔试真题2022
- 社区文化展览活动规划与实施三篇
- JT-T-1134-2017道路客货运运输驾驶员行车操作规范
- 2024年招录考试-法院书记员笔试考试历年典型考题及考点含含答案
- 2024公司与个人汽车租赁合同书
评论
0/150
提交评论