initialpop2.m

malab语言实现bspline曲线拟合

function [chrom,v,pbest_chrom,pbest,gbest_chrom,gbest] = initialpop2(NIND,m,p,Q)

midchrom = [ ];
for e = 0:1/(NIND-1):1
    midchrom = [midchrom;datapara(Q,m,e)];
end
chrom = [ ];
v = [];
for i = 1:NIND
    w = midchrom(i,:);
    t = rand(1); 
    ww = floor(t*(p+1)+(1-t)*(m+1)/2);                  % 朝负无穷方向舍入4<=ww<=17,rand:0~1;randn:-1~1
    www = [w,ww];
    chrom = [chrom;www];
    vv(i,:) = rand(1,m-1)*0.2-0.1;                      % 初始化种群微粒速度，值为-0.01-0.01，间隔为0.01                           
    vvv = rand()-0.5;                                  % 初始化ww的速度为-1~1
    v(i,:) = [vv(i,:),vvv];
    pbest_chrom(i,:) = chrom(i,:);
end
[objv,fitnv] = objfun2(m,p,chrom,Q);
pbest = fitnv;
[gbest,index_gbest] = max(fitnv);
gbest_chrom = chrom(index_gbest,:);