Matlab histogram 画出十二种常见的混沌映射
在群体智能优化算法中,需要随机产生数组,然而往往效果并不好。由于混沌序列具有非线性、遍历性和不可预测性等特点,因此随机数组可用混沌映射来替代。参考 Yu Y , Gao S , Cheng S , et al. CBSO: a memetic brain storm optimization with chaotic local search[J]. Memetic Computing, 2017.相关参数设置可见其论文,不在此赘述。
注:和论文的顺序不一样。
Matlab代码如下
clear;
clc;
%% Logistic
y_1=zeros(1,10^5);
y_1(1)=0.152;
mu=4;
for i = 1 : 10^5-1y_1(i+1)=mu*y_1(i)*(1-y_1(i));
end
%% Tent
y_2=zeros(1,10^5);
y_2(1)= 0.152;
Beta = 0.4;
for i = 1 : 10^5-1 if (y_2(i)<=Beta && y_2(i)>0)y_2(i+1) = y_2(i)/Beta;else y_2(i+1)=(1-y_2(i))/(1-Beta);end
end
%% Cubic
y_3=zeros(1,10^5);
y_3(1)= 0.242;
rho = 2.59;
for i = 1 : 10^5-1 y_3(i+1) = rho*y_3(i)*(1-y_3(i)^2);
end
%% Bernoulli
y_4=zeros(1,10^5);
y_4(1)= 0.152;
lambda = 0.4;
for i = 1 : 10^5-1 if (y_4(i)<=(1-lambda)) && ((y_4(i)>0))y_4(i+1) = y_4(i)/(1-lambda);else y_4(i+1)=(y_4(i)-1+lambda)/lambda;end
end
%% PWLCM
y_5=zeros(1,10^5);
y_5(1)= 0.002;
P = 0.7;
for i = 1 : 10^5-1 if (y_5(i)<P && y_5(i)>0)y_5(i+1) = y_5(i)/P;else y_5(i+1)=(1-y_5(i))/(1-P);end
end
%% Singer
y_6=zeros(1,10^5);
y_6(1)= 0.152;
mu_1 = 1.073;
for i = 1 : 10^5-1 y_6(i+1) = mu_1*(7.86*y_6(i)-23.31*y_6(i)^2+28.75*y_6(i)^3-13.302875*y_6(i)^4);
end
%% Sine
y_7=zeros(1,10^5);
y_7(1)= 0.152;
a = 4;
for i = 1 : 10^5-1 y_7(i+1) = a/4*sin(pi*y_7(i));
end
%% Gaussian
y_8=zeros(1,10^5);
y_8(1)= 0.152;
mu_2 = 1;
for i = 1 : 10^5-1if (y_8(i)==0)y_8(i+1) =0 ;elsey_8(i+1) =(rem(mu_2/y_8(i),1)) ;end
end
%% Chebyshew
y_9=zeros(1,10^5);
y_9(1)= 0.152;
phi = 5;
for i = 1 : 10^5-1y_9(i+1) =cos(phi*acos(y_9(i))) ;
end
%% Circle
y_10=zeros(1,10^5);
y_10(1)= 0.152;
a_1 = 0.5;
b_1 = 2.2;
for i = 1 : 10^5-1y_10(i+1) = y_10(i)+a_1-mod(b_1/(2*pi)*(sin(2*pi*y_10(i))),1);
end
%% Sinusoidal
y_11=zeros(1,10^5);
y_11(1)= 0.74;
a_2 = 2.3;
for i = 1 : 10^5-1y_11(i+1) = a_2*y_11(i)^2*sin(pi*y_11(i));
end
%% ICMIC
y_12=zeros(1,10^5);
y_12(1)= 0.152;
a_3 = 70;
for i = 1 : 10^5-1y_12(i+1) = sin(a_3/y_12(i));
end
%% 画图
subplot(3,4,1)
h1=histogram(y_1,200);
h1.FaceColor=[0 0 1];
xlim([0,1])%设置x轴范围
xlabel('Logistic map')subplot(3,4,2)
h2=histogram(y_2,200);
h2.FaceColor=[0 0 1];
xlim([0,1])%设置x轴范围
xlabel('Tent map')subplot(3,4,3)
h3=histogram(y_3,200);
h3.FaceColor=[0 0 1];
xlim([0,1])%设置x轴范围
xlabel('Cubic map')subplot(3,4,4)
h4=histogram(y_4,200);
h4.FaceColor=[0 0 1];
xlim([0,1])%设置x轴范围
xlabel('Bernoulli map')subplot(3,4,5)
h5=histogram(y_5,200);
h5.FaceColor=[0 0 1];
xlim([0,1])%设置x轴范围
xlabel('PWLCM map')
subplot(3,4,5)
h5=histogram(y_5,200);
h5.FaceColor=[0 0 1];
xlim([0,1])%设置x轴范围
xlabel('PWLCM map')subplot(3,4,6)
h6=histogram(y_6,200);
h6.FaceColor=[0 0 1];
xlim([0,1])%设置x轴范围
xlabel('Singer map')subplot(3,4,7)
h7=histogram(y_7,200);
h7.FaceColor=[0 0 1];
xlim([0,1])%设置x轴范围
xlabel('Sine map')subplot(3,4,8)
h8=histogram(y_8,200);
h8.FaceColor=[0 0 1];
xlim([0,1])%设置x轴范围
xlabel('Gaussian map')subplot(3,4,9)
h9=histogram(y_9,200);
h9.FaceColor=[0 0 1];
xlim([0,1])%设置x轴范围
xlabel('Chebyshew map')subplot(3,4,10)
h10=histogram(y_10,200);
h10.FaceColor=[0 0 1];
xlim([0,1])%设置x轴范围
xlabel('Circle map')subplot(3,4,11)
h11=histogram(y_11,200);
h11.FaceColor=[0 0 1];
xlim([0,1])%设置x轴范围
xlabel('Sinusoidal map')subplot(3,4,12)
h12=histogram(y_12,200);
h12.FaceColor=[0 0 1];
xlim([0,1])%设置x轴范围
xlabel('ICMIC map')
Matlab histogram 画出十二种常见的混沌映射相关推荐
- 十二种常见设计模式代码详解
零:设计模式分类 设计模式有创建型模式.结构型模式与行为型模式 创建型:单例模式.工厂模式(简单工厂,工厂方法,抽象工厂) 结构型:适配器模式.门面模式.装饰器模式.注册树模式.代理模式.管道模式 行 ...
- 薛斯通道的十二种买入法和四种卖出法
薛斯通道包括两组通道指标,分别是长期大通道指标和短期小通道指标.股价实际上是被短期小通道包容着在长期大通道中上下运行.薛斯通道的大.小通道总共有四条轨道线: 内上轨(黄色线或紫色线)代表短 ...
- 关于Authorware的十二种使用技巧
Authorware是美国Macromedia公司(现已被adobe公司收购)开发的一种多媒体制作软件,它是一个图标导向式的多媒体开发工具.今天我们学习一下Authorware的十二种使用技巧,如果你 ...
- gcms基峰有什么用_GC-MS十二种个常见问题和解决办法
GC-MS十二种个常见问题和解决办法 浏览次数:9662发布日期:2017-03-28 1 样品进样一段时间后,突然进样口中压力不上不去了,这是为什么呢? 其实主要是因为我们进样次数过多,导致隔垫密封 ...
- 中国筷子的十二种忌讳(转)
中国人使用筷子用餐是从远古流传下来的,古时又称其为"箸",日常生活当中对筷子 的运用是非常有讲究的.一般我们在使用筷子时,正确的使用方法讲究得是用右手执筷,大拇指和食指捏住筷子的上 ...
- 十二生肖被戏称为“十二缺”,指的是这十二种动物各自都缺少一样东西
十二生肖被戏称为"十二缺",指的是这十二种动物各自都缺少一样东西 那么这些生肖动物都各自缺少什么呢?你知道你的属相缺少的是什么东西吗?下面就一起来看看吧! [鼠无脑] 所谓兽无完兽 ...
- 利用matlab命令画出以下信号的波形,MATLAB实验报告
文档收集于互联网,已重新整理排版.word 版本可编辑,有帮助欢迎下载支持. 1文档来源为:从网络收集整理.word 版本可编辑. 实验一 名称:连续时间信号分析 姓名:王嘉琦 学号:0636 班级: ...
- 安卓判断服务器返回的状态码,关于服务器返回的十四种常见HTTP状态码详解
原标题:关于服务器返回的十四种常见HTTP状态码详解 HTTP状态码 状态码是由3位数字和原因短语组成的(比如最常见的:200 OK),其中第一位数字表示响应类别,响应类别从1到5分为五种 add:其 ...
- 优秀员工的十二种品质 - 马云在上海演讲说的话
马云在上海演讲说的话,值得大家仔细阅读一下. 一: 忠诚--忠心者不被解雇 单位可能开除有能力的员工,但对一个忠心耿耿的人,不会有领导愿意让他走,他会成为单位这个铁打营盘中最长久的战士,而且是最有发展 ...
最新文章
- vault-使用kubernetes作为认证后端
- 什么行业需要java_java开发可以从事什么行业
- 如何设置WebLogic生成的日志文件的权限?
- ubuntu安装samba
- Procrustes Analysis(普氏分析)
- We7 CMS 2.7beta3 版发布!SVN源码同步开放
- 【递归算法01】递归的调用机制
- Facebook 有序队列服务设计原理和高性能浅析
- 常用算法(冒泡、插入、选择、快速)和二叉树详解
- android自定义横向时间轴,Android自定义时间轴的实现过程
- 窗口间的通信(消息互发与数据传递)
- 方舟生存服务器没有响应怎么办,方舟生存进化搜服务器闪退怎么办? Fata error临时解决办法一览...
- 【论文阅读】Hierarchical Multi-modal Contextual Attention Network for Fake News Detection --- 虚假新闻检测,多模态
- 微信鉴权服务器地址,微信开发之微信授权登录
- 少儿编程和机器人编程哪个更好一点
- 自己封装了一个带有比较日期的datepicker组件
- Linux(安装apache,设置端口,基于端口的配置方法)
- 数据结构与算法——算法基础
- 51单片机的PID水温控制器设计
- 黑板、电梯、签到的测试点