gausstest.m

来自「氧乐果控制过程」· M 代码 · 共 53 行

chrom=rand(50,77);
n=7;
m=7;
ns=0;
nmtol=77;
n1=size(chrom)
while (ns<n1)
 x=rand(1,2);   
ns=ns+1;

j=0; 
for se=1:n
    j=j+1;
w11(j)=chrom(ns,se);
end

j=0; 
for se=n+1:n*2
    j=j+1;
b11(j)=chrom(ns,se);
end

j=0; 
for se=n*2+1:n*2+m
    j=j+1;
w12(j)=chrom(ns,se);
end

j=0; 
for se=n*2+m+1:(n+m)*2
    j=j+1;
b12(j)=chrom(ns,se);
end

j=0; 
for se=(n+m)*2+1:nmtol
    j=j+1;
w2(j)=chrom(ns,se);
end

a11=(w11.*x(1)+b11).*(w11.*x(1)+b11);
a12=(w12.*x(2)+b12).*(w12.*x(2)+b12);
out=[];
for se=1:n
 out=[out a11(se)+a12];
end
out=exp(-out);
w2
out
u=out*w2'

end