| 查看: 3278 | 回复: 6 | |||
[求助]
遗传算法优化BP神经网络权值阈值
|
|
我从网上下了一个遗传算法优化BP神经网络的程序，可是运行后总出现错误，不知道哪里错了，请大侠们帮忙给我检查一下；源代码如下： 程序一：GA训练BP权值的主函数 function net=GABPNET(XX,YY) %-------------------------------------------------------------------------- % GABPNET.m % 使用遗传算法对BP网络权值阈值进行优化，再用BP算法训练网络 %-------------------------------------------------------------------------- %数据归一化预处理 nntwarn off XX=premnmx(XX); YY=premnmx(YY); %创建网络 net=newff(minmax(XX),[19,25,1],{'tansig','tansig','purelin'},'trainlm'); %下面使用遗传算法对网络进行优化 P=XX; T=YY; R=size(P,1); S2=size(T,1); S1=25;%隐含层节点数 S=R*S1+S1*S2+S1+S2;%遗传算法编码长度 aa=ones(S,1)*[-1,1]; popu=50;%种群规模 initPpp=initializega(popu,aa,'gabpEval');%初始化种群 gen=100;%遗传代数 %下面调用gaot工具箱，其中目标函数定义为gabpEval [x,endPop,bPop,trace]=ga(aa,'gabpEval',[],initPpp,[1e-6 1 1],'maxGenTerm',gen,... 'normGeomSelect',[0.09],['arithXover'],[2],'nonUnifMutation',[2 gen 3]); %绘收敛曲线图 figure(1) plot(trace(:,1),1./trace(:,3),'r-'); hold on plot(trace(:,1),1./trace(:,2),'b-'); xlabel('Generation'); ylabel('Sum-Squared Error'); figure(2) plot(trace(:,1),trace(:,3),'r-'); hold on plot(trace(:,1),trace(:,2),'b-'); xlabel('Generation'); ylabel('Fittness'); %下面将初步得到的权值矩阵赋给尚未开始训练的BP网络 [W1,B1,W2,B2,P,T,A1,A2,SE,val]=gadecod(x); net.LW{2,1}=W1; net.LW{3,2}=W2; net.b{2,1}=B1; net.b{3,1}=B2; XX=P; YY=T; %设置训练参数 net.trainParam.show=1; net.trainParam.lr=1; net.trainParam.epochs=50; net.trainParam.goal=0.001; %训练网络 net=train(net,XX,YY); 程序二：适应值函数 function [sol, val] = gabpEval(sol,options) % val - the fittness of this individual % sol - the individual, returned to allow for **rckian evolution % options - [current_generation] load data2; nntwarn off XX=premnmx(XX); YY=premnmx(YY); P=XX; T=YY; R=size(P,1); S2=size(T,1); S1=25;%隐含层节点数 S=R*S1+S1*S2+S1+S2;%遗传算法编码长度 for i=1:S, x(i)=sol(i); end; [W1, B1, W2, B2, P, T, A1, A2, SE, val]=gadecod(x); 程序三：编解码函数 function [W1, B1, W2, B2, P, T, A1, A2, SE, val]=gadecod(x) load data2; nntwarn off; XX=premnmx(XX); YY=premnmx(YY); P=XX; T=YY; R=size(P,1); S2=size(T,1); S1=25;%隐含层节点数 S=R*S1+S1*S2+S1+S2;%遗传算法编码长度 % 前R*S1个编码为W1 for i=1:S1, for k=1:R, W1(i,k)=x(R*(i-1)+k); end end % 接着的S1*S2个编码（即第R*S1个后的编码）为W2 for i=1:S2, for k=1:S1, W2(i,k)=x(S1*(i-1)+k+R*S1); end end % 接着的S1个编码（即第R*S1+S1*S2个后的编码）为B1 for i=1:S1, B1(i,1)=x((R*S1+S1*S2)+i); end % 接着的S2个编码（即第R*S1+S1*S2+S1个后的编码）为B2 for i=1:S2, B2(i,1)=x((R*S1+S1*S2+S1)+i); end % 计算S1与S2层的输出 A1=tansig(W1*P,B1); A2=purelin(W2*A1,B2); % 计算误差平方和 SE=sumsqr(T-A2); val=1/SE; % 遗传算法的适应值 这些命令流直接复制粘贴到Matlab里面就可以用吗？ |
回复此楼» 收录本帖的淘帖专辑推荐
 资源代码 |
» 猜你喜欢
东南大学364求调剂
已经有5人回复
-
pva水凝胶电池
已经有0人回复
-
金属材料论文润色/翻译怎么收费?
已经有151人回复
-
分析位错线
已经有5人回复
-
上海高校教授长年招收硕士和博士研究生-动力工程与工程热物理
已经有5人回复
-
天然气锅炉
已经有0人回复
-
讲师身份申请面上项目中标率有影响吗
已经有32人回复
-
希望面上好运
已经有75人回复
-
本科211能动类专业 ，初试302分，六级已过，求调剂！
已经有3人回复
-
论文投稿
已经有0人回复
-
303分 0807学硕求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab中MSE 过大
已经有8人回复
- 利用PSO算法训练BP神经网络时，如何确定粒子位置的取值范围？
已经有3人回复
- 如何人工神经网络来预测下一个数值？？？
已经有8人回复
- matlab中的BP神经网络归一化问题
已经有3人回复
- BP网络初始值与权值问题
已经有3人回复
- 关于spss和神经网络结合的问题
已经有3人回复
- 人工神经网络BP算法源代码与演示程序
已经有222人回复
- DE值的具体算法 真的是不知道了
已经有9人回复
- matlab 的 bp神经网络 预测 问题..小作业..
已经有9人回复
2楼2012-08-20 10:07:38
3楼2012-08-20 10:09:17
xesgue
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1123.9
- 散金: 810
- 红花: 1
- 帖子: 233
- 在线: 298.4小时
- 虫号: 2126184
- 注册: 2012-11-11
- 性别: GG
- 专业: 食品科学基础
4楼2012-11-16 14:07:41
1007695994
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 434.5
- 散金: 6
- 红花: 1
- 帖子: 80
- 在线: 18.3小时
- 虫号: 2331967
- 注册: 2013-03-09
- 性别: GG
- 专业: 电力系统
5楼2013-04-07 15:16:24
6楼2013-10-27 11:07:12
1007695994
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 434.5
- 散金: 6
- 红花: 1
- 帖子: 80
- 在线: 18.3小时
- 虫号: 2331967
- 注册: 2013-03-09
- 性别: GG
- 专业: 电力系统
7楼2013-10-28 00:55:16
