| 查看: 3279 | 回复: 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人回复
-
2026考研调剂-上海工程技术大学材料机械
已经有0人回复
-
金属材料论文润色/翻译怎么收费?
已经有139人回复
-
2026考研调剂上海工程技术大学材料学院
已经有0人回复
-
2026上海工程技术大学李老师化学材料机械计算机方向招硕士7-10名
已经有0人回复
-
上海工程技术大学材料方向招收硕士7-10名
已经有10人回复
-
296材料专硕求调剂
已经有22人回复
-
锂优电池测试仪如何实现电池性能检测的精准性与可靠性
已经有0人回复
-
西安石油大学新能源学院接收材料类、能源动力类、机械类、计算机类等专业专硕调剂生！
已经有11人回复
-
西安石油大学新能源学院接收材料类、能源动力类、机械类、计算机类等专业专硕调剂生！
已经有10人回复
-
JMPT期刊有字数限制吗
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
