空间金字塔匹配（Spatial Pyramid Matching，SPM）混沌映射系统

5月14号新鲜文章用到的一个混沌映射初始化种群方法

张建华,张斌,赵明伟,等.基于改进型GNDO算法的Stewart平台位姿误差补偿[J/OL].控制工程:1-11[2024-05-14].https://doi.org/10.14107/j.cnki.kzgc.20230674.

在文章中，r，x（t）均为[0,1]之间 的随机数，n为（0,1）的一个随机数，但未给出u的值所以无法还原文章的柱状图

clear;clc
l = 600;%维度
u=0.3;r=rand;n=0.4;
xn =rand(1,l);
for i=1:l
    if 0<=xn(i)&xn(i)<n
        x(i)=mod((xn(i)/n)+u*sin(pi*xn(i))+r,1);
    elseif n<=xn(i)&xn(i)<0.5
        x(i)=mod((xn(i)/(n*(0.5-n)))+u*sin(pi*xn(i))+r,1);
    elseif 0.5<=xn(i)&xn(i)<1-n
        x(i)=mod(((1-xn(i))/(n*(0.5-n)))+u*sin(pi*(1-xn(i)))+r,1);
    else 1-n<=xn(i)&xn(i)<1
        x(i)=mod(((1-xn(i))/n)+u*sin(pi*(1-xn(i)))+r,1);
    end
end
% 绘制散点图
figure
plot(x,'.')
xlabel('维度')
ylabel('混沌值')
% 绘制直方图
figure
histogram(x)
xlabel('混沌值')
ylabel('频数')