generate_choasserials.m

采用logistic构造的混沌模型

function [immigrantdata]=generate_choasserials(immigrant_pop,group_dim,us,ux)
% Logistic choas programme
% 使用平台 - Matlab6.5 / Matlab7.0
% 表达式
% x(n+1) = lambda * x(n) * (1 - x(n)) 
% 利用logistic产生混沌序列

%保存混沌序列系数的数组
choasdata=zeros(group_dim,immigrant_pop);
%保存混沌序列的数组，用于免疫计算。
immigrantdata=zeros(group_dim,immigrant_pop);

%保存设计变量的上下限数组
%us=zeros(group_dim);
%ux=zeros(group_dim);

%产生混沌序列系数因子

a = 0; b = 1;
x = a + (b-a) * rand(1,group_dim);
for i=1:group_dim
    choasdata(i,1)=x(1,i);
end

%choasdata(1,1)=0.234;
choaslanbuda=4.0;
for i=1:group_dim
    for j=2:immigrant_pop
      choasdata(i,j)=choaslanbuda*choasdata(i,j-1)*(1-choasdata(i,j-1));
    end
end
%choasdata(1,:)=choasdata(immigrant_pop,:);
choasdata
figure(1);
subplot(2,1,1); 
plot(choasdata(1,:));
ylabel('\混沌随机数')
xlabel('数据点序号');
subplot(2,1,2);
plot(choasdata(2,:));
ylabel('\混沌随机数')
xlabel('数据点序号');

figure(2);
plot(choasdata(1,:),choasdata(2,:),'ro');
%产生混沌序列，用于免疫计算。
for i=1:group_dim
    for j=1:immigrant_pop
     immigrantdata(i,j)=ux(i)+choasdata(i,j)*(us(i)-ux(i));
    end
end