2.m

来自「基于matlab的人工免疫算法」· M 代码 · 共 19 行

c=load('x.txt');
width2=input('请输入窗宽')
step2=input('请输入步长')
m2=length(c)
for i=1:floor((m2-width2)/step2+1)
t(i,1:width2)=c(1,1+(i-1)*step2:width2+(i-1)*step2);
end
b2=max(max(t))
T=t./b2
e=zeros(1,m2)
    for i=1:floor((m2-width2)/step2+1)
        for k=1:nd
        if (norm(T(i,:)-d{k}))^2<Rd{k}
            e(1,i*step2-1)=1;
            break
        end
    end
end
plot(e,'rd')