技术标签: 【啃书】《MATLAB智能算法30个案例分析(第2版)》
%20201229lu注:该代码在matlab2019b成功运行
%% 该函数演示多目标perota优化问题
%清空环境
clc
clear
load data %导入data.mat数据,data.mat中即储存的问题中P,R,C
%% 初始参数
objnum=size(P,1); %类中物品个数
weight=92; %总重量限制
%初始化程序
Dim=5; %粒子维数
xSize=50; %种群个数
MaxIt=200; %迭代次数
c1=0.8; %算法参数
c2=0.8; %算法参数
wmax=1.2; %惯性因子
wmin=0.1; %惯性因子
x=unidrnd(4,xSize,Dim); %粒子初始化
v=zeros(xSize,Dim); %速度初始化
xbest=x; %个体最佳值
gbest=x(1,:); %粒子群最佳位置
% 粒子适应度值
px=zeros(1,xSize); %粒子价值目标
rx=zeros(1,xSize); %粒子体积目标
cx=zeros(1,xSize); %重量约束
% 最优值初始化
pxbest=zeros(1,xSize); %粒子最优价值目标
rxbest=zeros(1,xSize); %粒子最优体积目标
cxbest=zeros(1,xSize); %记录重量,以求约束
% 上一次的值
pxPrior=zeros(1,xSize);%粒子价值目标
rxPrior=zeros(1,xSize);%粒子体积目标
cxPrior=zeros(1,xSize);%记录重量,以求约束
%计算初始目标向量
for i=1:xSize
for j=1:Dim %控制类别
px(i) = px(i)+P(x(i,j),j); %粒子价值
rx(i) = rx(i)+R(x(i,j),j); %粒子体积
cx(i) = cx(i)+C(x(i,j),j); %粒子重量
end
end
% 粒子最优位置
pxbest=px;rxbest=rx;cxbest=cx;
%% 初始筛选非劣解
flj=[];
fljx=[];
fljNum=0;
%两个实数相等精度
tol=1e-7;
for i=1:xSize
flag=0; %支配标志
for j=1:xSize
if j~=i
if ((px(i)<px(j)) && (rx(i)>rx(j))) ||((abs(px(i)-px(j))<tol)...
&& (rx(i)>rx(j)))||((px(i)<px(j)) && (abs(rx(i)-rx(j))<tol)) || (cx(i)>weight)
flag=1;
break;
end
end
end
%判断有无被支配
if flag==0
fljNum=fljNum+1;
% 记录非劣解
flj(fljNum,1)=px(i);flj(fljNum,2)=rx(i);flj(fljNum,3)=cx(i);
% 非劣解位置
fljx(fljNum,:)=x(i,:);
end
end
%% 循环迭代
for iter=1:MaxIt
% 权值更新
w=wmax-(wmax-wmin)*iter/MaxIt;
%从非劣解中选择粒子作为全局最优解
s=size(fljx,1);
index=randi(s,1,1);
gbest=fljx(index,:);
%% 群体更新
for i=1:xSize
%速度更新
v(i,:)=w*v(i,:)+c1*rand(1,1)*(xbest(i,:)-x(i,:))+c2*rand(1,1)*(gbest-x(i,:));
%位置更新
x(i,:)=x(i,:)+v(i,:);
x(i,:) = rem(x(i,:),objnum)/double(objnum);
index1=find(x(i,:)<=0);
if ~isempty(index1)
x(i,index1)=rand(size(index1));
end
x(i,:)=ceil(4*x(i,:));
end
%% 计算个体适应度
pxPrior(:)=0;
rxPrior(:)=0;
cxPrior(:)=0;
for i=1:xSize
for j=1:Dim %控制类别
pxPrior(i) = pxPrior(i)+P(x(i,j),j); %计算粒子i 价值
rxPrior(i) = rxPrior(i)+R(x(i,j),j); %计算粒子i 体积
cxPrior(i) = cxPrior(i)+C(x(i,j),j); %计算粒子i 重量
end
end
%% 更新粒子历史最佳
for i=1:xSize
%现在的支配原有的,替代原有的
if ((px(i)<pxPrior(i)) && (rx(i)>rxPrior(i))) ||((abs(px(i)-pxPrior(i))<tol)...
&& (rx(i)>rxPrior(i)))||((px(i)<pxPrior(i)) && (abs(rx(i)-rxPrior(i))<tol)) || (cx(i)>weight)
xbest(i,:)=x(i,:);%没有记录目标值
pxbest(i)=pxPrior(i);rxbest(i)=rxPrior(i);cxbest(i)=cxPrior(i);
end
%彼此不受支配,随机决定
if ~( ((px(i)<pxPrior(i)) && (rx(i)>rxPrior(i))) ||((abs(px(i)-pxPrior(i))<tol)...
&& (rx(i)>rxPrior(i)))||((px(i)<pxPrior(i)) && (abs(rx(i)-rxPrior(i))<tol)) || (cx(i)>weight) )...
&& ~( ((pxPrior(i)<px(i)) && (rxPrior(i)>rx(i))) ||((abs(pxPrior(i)-px(i))<tol) && (rxPrior(i)>rx(i)))...
||((pxPrior(i)<px(i)) && (abs(rxPrior(i)-rx(i))<tol)) || (cxPrior(i)>weight) )
if rand(1,1)<0.5
xbest(i,:)=x(i,:);
pxbest(i)=pxPrior(i);rxbest(i)=rxPrior(i);cxbest(i)=cxPrior(i);
end
end
end
%% 更新非劣解集合
px=pxPrior;
rx=rxPrior;
cx=cxPrior;
%更新升级非劣解集合
s=size(flj,1);%目前非劣解集合中元素个数
%先将非劣解集合和xbest合并
pppx=zeros(1,s+xSize);
rrrx=zeros(1,s+xSize);
cccx=zeros(1,s+xSize);
pppx(1:xSize)=pxbest;pppx(xSize+1:end)=flj(:,1)';
rrrx(1:xSize)=rxbest;rrrx(xSize+1:end)=flj(:,2)';
cccx(1:xSize)=cxbest;cccx(xSize+1:end)=flj(:,3)';
xxbest=zeros(s+xSize,Dim);
xxbest(1:xSize,:)=xbest;
xxbest(xSize+1:end,:)=fljx;
%筛选非劣解
flj=[];
fljx=[];
k=0;
tol=1e-7;
for i=1:xSize+s
flag=0;%没有被支配
%判断该点是否非劣
for j=1:xSize+s
if j~=i
if ((pppx(i)<pppx(j)) && (rrrx(i)>rrrx(j))) ||((abs(pppx(i)-pppx(j))<tol) ...
&& (rrrx(i)>rrrx(j)))||((pppx(i)<pppx(j)) && (abs(rrrx(i)-rrrx(j))<tol)) ...
|| (cccx(i)>weight) %有一次被支配
flag=1;
break;
end
end
end
%判断有无被支配
if flag==0
k=k+1;
flj(k,1)=pppx(i);flj(k,2)=rrrx(i);flj(k,3)=cccx(i);%记录非劣解
fljx(k,:)=xxbest(i,:);%非劣解位置
end
end
%去掉重复粒子
repflag=0; %重复标志
k=1; %不同非劣解粒子数
flj2=[]; %存储不同非劣解
fljx2=[]; %存储不同非劣解粒子位置
flj2(k,:)=flj(1,:);
fljx2(k,:)=fljx(1,:);
for j=2:size(flj,1)
repflag=0; %重复标志
for i=1:size(flj2,1)
result=(fljx(j,:)==fljx2(i,:));
if length(find(result==1))==Dim
repflag=1;%有重复
end
end
%粒子不同,存储
if repflag==0
k=k+1;
flj2(k,:)=flj(j,:);
fljx2(k,:)=fljx(j,:);
end
end
%非劣解更新
flj=flj2;
fljx=fljx2;
end
%绘制非劣解分布
plot(flj(:,1),flj(:,2),'o')
xlabel('P')
ylabel('R')
title('最终非劣解在目标空间分布')
disp('非劣解flj中三列依次为P,R,C')
中 MATLAB智能算法30个案例分析(第2版)[史峰 - 王辉 - 郁磊][北京航空航天大学出版社][2015-9-1][9787512414112]
文章浏览阅读161次。需要把这个RedisShake安装到k8s集群,然后k8s集群相互检测到redis集群的节点(因为节点直接都是用内Pod内网进行通讯所以必须安装到K8s的Pod内部切记!redis单机之间的迁移比较好解决,但是和redis和k8s-redis集群迁移目前市面的成熟解决方案较少。接下来分享一个较为成熟的解决方案来解决这个问题。文件,写上对应的端口号。和redis版本(很重要)遇到重复的key要求重写或者跳过。主备好pod里面的命令。
文章浏览阅读6.1k次。在国内,我们使用mac电脑,默认都是使用中文的。如果我们想要体验一下其他语言环境下的操作系统,那么如何设置更改Mac的系统语言呢?下面小编就给大家讲解一下苹果电脑系统语言更换的设置方法。1、首先在里面找到。2、点击左下角的号按钮,选择添加新的语言后点击下方的按钮。3、然后系统会提示需要将选择的语言设置为首选语言吗,选择是的话,系统语言就会被更改成我们选择的语言了。4、我们暂时选择不使用,回到设置页,在左边菜单栏中可以看到鼠标已经变成了的样子。5、英文的选项,直接至顶部,那么英文就变成了首选系统语言。6、更改_如何在苹果电脑里设置语言
文章浏览阅读5.6k次,点赞4次,收藏19次。标题 opencv 中常用数据类型的区别以及数据类型的转换8代表8bit 每一个像素点在内存空间占用8bitS–代表—signed int—有符号整形U–代表–unsigned int–无符号整形F–代表–float---------单精度浮点型同时注意 CV_8U 代表的是单通道图像,数据范围是{0-255} CV_32F代表的图像是单通道图像,数据范围是{0.0 - 1.0}imshow在显示图像时均会映射到{0-255}之间 CV_8U显示的就是0-255但要注意,在显示CV_32F类_cv_8u
文章浏览阅读7.7k次,点赞3次,收藏4次。问题1:读取onnx出现报错“error: (-210:Unsupported format or combination of formats) Failed to parse onnx model in function 'ONNXImporter'”出现这个原因网上有挺多种解答,但是首先先检查onnx文件路径是否填写正确。问题2:读取onnx出现报错“error: (-2:Unspecified error) Can't create layer "317" of type "Clip"._error: (-210:unsupported format or combination of formats) failed to parse o
文章浏览阅读1k次,点赞7次,收藏15次。Cortex-M 嵌入式系统 MCU Arm 开发板。_liteos keil
文章浏览阅读1.8k次。对JVM运行参数进行修改是JVM性能调优的重要手段,下面介绍在应用程序开发过程中JVM参数设置的几种方式。方式一java程序运行时指定 -Dproperty=value该参数通常用于设置系统级全局变量值,如配置文件路径,保证该属性在程序中任何地方都可访问。当然,也可以通过在程序中使用System.setProperty进行设置。注意:1、如果-Dproperty=value的value中包含空格,..._idea run property参数
文章浏览阅读486次。因为看到有人对这个概念很模糊,怕自己以后脑子抽了也懵了,所以记录一下。实际上非常简单
文章浏览阅读759次,点赞19次,收藏12次。【考点1】函数的定义这里函数名为func,返回类型为int,并接受两个int类型参数。【考点2】库函数例如,调用C语言的数学库函数需要引入头文件。此代码段会计算4.0的平方根。【考点3】函数的返回值函数可以通过 return 语句返回一个值。这个函数返回两个整数的和。【考点4】函数的声明【考点5】函数的调用此代码调用之前声明的multiply函数,并将返回值存储在result中。【考点6】函数的参数及值传递在C语言中,函数参数默认是值传递。
文章浏览阅读2.4k次,点赞3次,收藏9次。newcommand, renewcommand, newtheorem的使用和总结_latex renewcommand
文章浏览阅读51次。内存泄漏检测工具2007年08月08日1. ccmalloc-Linux和Solaris下对C和C++程序的简单的使用内存泄漏和malloc调试库。 2. Dmalloc-Debug Malloc Library. 3. Electric Fence-Linux分发版中由Bruce Perens编写的malloc()调试库。 4. Leaky..._修改1286400
文章浏览阅读184次。【BILIBILI网课——电子】Verilog硬件描述语言 西安电子科技大学 蔡觉平等主讲硬件安全技术+数字IC串讲——摆渡沧桑程序员宅基地牛客网——verilogHDLBITS——Verilog练习SOC总线的相关知识verilog——FPGA 20个例程专栏(作者: 青青豌豆)字节FPGA实习字节IP实习华为个人中心_fpga学习课程推荐 csdn
文章浏览阅读231次。/client.setCallback(callback)在vscode中只能放在loop函数中而不能放在setup中,而且在arduino中无论是放在哪里速度都差不多。2.stm32代码请到github下载https://github.com/fengyib/Remote-Control-of-Stepper-Motor。首先本项目通过esp8266连接服务器,实现远程通信,在和stm32串口通信实现pca9685的远程控制。1.esp8266代码。//回调函数只能有一个。_服务器下发指令通过sp8266控制stm32