数据挖掘上机报告 - 图文(5)

function r = entropy(z) s=sum(z); r=0;

for ri=1:length(z)

r=r-z(ri)/s.*log2(z(ri)/s); end;

验证：

第四次上机训练报告

代码：

function [z]=my_fun(x) z=x*sin(10*pi*x)+2; z=-z; end

>> options = gaoptimset('Generations',200,'PopulationSize',20,'PlotFcns',@gaplotbestf,'StallGenLimit', 200); >> options.Generations = 2000; %最大迭代数设为2000

>> [X,FVAL,EXITFLAG,OUTPUT] =ga(@my_fun, 1 ,[], [],[],[],[],[2,-1],[],options); Optimization terminated: average change in the fitness value less than options.TolFun. 运行结果如图