mainprog.m

基于GA的BP改进算法

% 求下列函数的最大值????????????????????????????????????????????????%

%?? 将 x 的值用一个10位的二值形式表示为二值问题?????????????????????%
%----------------------------------------------
% 编程
%----------------------------------------------
%  主程序
%遗传算法主程序
clear
clf
popsize=20; %设置初始参数，群体大小
chromlength=10; %字符串长度（个体长度），染色体长度
pc=0.6; %设置交叉概率，本例中交叉概率是定值，若想设置变化的交叉概率可用表达式表示，或从写一个
        %交叉概率函数，例如用神经网络训练得到的值作为交叉概率
pm=0.001; %设置变异概率，同理也可设置为变化的

pop=initpop(popsize,chromlength); %运行初始化函数，随机产生初始群体
for i=1:30 %20为迭代次数
[objvalue]=calobjvalue(pop);%计算目标函数
fitvalue=calfitvalue(objvalue); %计算群体中每个个体的适应度
[newpop]=selection(pop,fitvalue); %复制
[newpop]=Crossover(pop,pc); %交叉
[newpop]=mutation(pop,pc);%变异
[bestindividual,bestfit]=best(pop,fitvalue);%求出群体中适应值最大的个体及其适应值
y(i)=max(bestfit);
n(i)=i;
pop5=bestindividual;
x(i)=decodechrom(pop5,1,chromlength)*10/1023;
pop=newpop;
end
y(i)
fplot('10*cos(3*x)-4*x',[0 30])

grid on
hold on
plot(x,y,'g*')
hold off