update_position.m

malab语言实现bspline曲线拟合

function chrom = update_position(chrom,v,NIND,m)

% 更新粒子的位置

for i = 1 : NIND
    chrom(i,2:m) = chrom(i,2:m) + v(i,1:m-1);
    chrom(i,m+2) = round(chrom(i,m+2) + v(i,m));
    
    % ===============保持chrom元素的规范============
    if chrom(i,2) < 0
        chrom(i,2) = chrom(i,3)/2;
    end
    for j = 3:m
        if (chrom(i,j) > 1) | (chrom(i,j) < 0)
            chrom(i,j) = (chrom(i,j-1) + chrom(i,j+1))/2;
        end
    end
    % =====================================
    
    for j = 1:m
        if chrom(i,j) > chrom(i,j)
           break;
        end
    end
    
    if chrom(i,m+2) > 17
        chrom(i,m+2) = 17;
    elseif chrom(i,m+2) < 4
        chrom(i,m+2) = 4;
    end
    
end