| 查看: 3243 | 回复: 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里面就可以用吗？ |
回复此楼» 收录本帖的淘帖专辑推荐
 资源代码 |
» 猜你喜欢
钢厂耐磨辊子激光熔覆修复加工
已经有0人回复
-
轴上高速激光熔覆铜技术
已经有0人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有281人回复
-
大型设备核心耐磨件激光淬火修复
已经有0人回复
-
北京高校副校长团队招收机械类，环境类学硕和专硕
已经有10人回复
-
液压杆激光熔覆修复加工
已经有0人回复
-
抽油杆激光熔覆修复加工优势及应用实例
已经有0人回复
-
抽油杆激光熔覆修复加工流程详解
已经有0人回复
-
西安航空学院机械工程学院杨武成院长一行莅临西安国盛激光参观交流
已经有0人回复
-
西安交通大学张志芬老师一行莅临西安国盛激光参观交流
已经有0人回复
-
北京高校副校长团队招收机械类，环境类学硕和专硕
已经有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
5