matlab一维变量优化问题.txt

MATLAB绘制出权值wv和阀值bv确定的误差曲面

一维变量优化问题2007/06/02 23:59计算下列函数的最大值
f(x)=x+10*sin(5x)+7*cos(4x)    x[0,9]
%编写目标函数,并存为 gaDemo1Eeval.m
%参数说明
%eval：个体的适应度； 
%sol：当前个体，n+1个元素的行向量。
function [sol, eval] =gaDemo1Eeval(sol,options)
x=sol(1);
eval = x + 10*sin(5*x)+7*cos(4*x);
%另存m文件
%遗传算法求最大值
clc
fplot('x + 10*sin(5*x)+7*cos(4*x)',[0 9])
% 生成初始种群，大小为10 
initPop=initializega(10,[0 9],'gademo1eeval1');
plot (initPop(:,1),initPop(:,2),'b*')
% 调用遗传函数
% 1次遗传迭代
[x endPop] = ga([0 9],'gademo1eeval1',[],initPop,[1e-5 ...
1 1],'maxGenTerm',1,'normGeomSelect',[0.08],['arithXover'], ... 
[20], 'nonUnifMutation',[2 1 3]);
plot (endPop(:,1),endPop(:,2),'bo')
% 25次遗传迭代
[x endPop bpop trace] = ga([0 9],'gademo1eeval1',[],initPop, [1e-6 1 1],'maxGenTerm',25, 'normGeomSelect',[0.08], ['arithXover'],[2],'nonUnifMutation',[2 25 3]);
plot (endPop(:,1),endPop(:,2),'y*')
figure(2)
plot(trace(:,1),trace(:,3),'y-')
hold on
plot(trace(:,1),trace(:,2),'r-')
xlabel('Generation'); ylabel('Fittness');
legend('解的变化','种群平均值的变化');
结果
1 19.001489
2 24.109594
1 19.001489
2 
3 21.879575
4 5 6 7 8 9 10 
11 22.349041
12 
13 24.854483
14 15 16 17 18 19 20 21 
22 24.854826
23 24 
25 24.854826
>> x
x =
     7.8585    24.8548
即经过25次迭代后，得到寻优结果
x=7.8585时
f(x)=24.8548