24小时热门版块排行榜    

CyRhmU.jpeg
查看: 453  |  回复: 3
当前主题已经存档。

yujunhui

银虫 (著名写手)

[交流] 【求助 遗传算法优化神经网络时候 initializega 函数出错,请高手指点?【已完成】

%% 清除环境变量
clear all
clc
warning off
nntwarn off

%% 声明全局变量
global p14     % 训练集输入数据
global t14     % 训练集输出数据
global R     % 输入神经元个数
global S2    % 输出神经元个数
global S1    % 隐层神经元个数
global S     % 编码长度
S1=9;



% 410S+20R   
x1=[0.036    12.56     0         0.315      0.473   0.002  0.002   0.14        0     0          0.58      0.23    0.005     0.015 ];
%0Cr13Al+16MnR
x2=[ 0.08     13.5       0         1.00       1.00     0.03  0.035        0.170   0.230     0.15       1.540       0.26    0.007     0.022 ];
%0Cr13+20R         
x3=[0.080   13.100     0.6       0.800      0.600   0.020 0.030    0.14    0          0        0.58       0.23      0.005    0.015 ];
%0Cr13+16MnR  
x4=[ 0.080   13.100     0.6       0.800      0.600   0.020     0.030   0.170   0.230     0.15       1.540       0.26    0.007     0.022 ];
%  410S+16MnR       
x5=  [ 0.036    12.56     0         0.315      0.473   0.002     0.002 0.170     0.230     0.15      1.540      0.26    0.007     0.022];

        

p=[x1' x2' x3' x4' x5'];

t1=[ 20             40              192.0                 262             182.7                320                  440             30.5           10     ];
t2=[ 18             25            165.0                 222.0             149.8                390                  525             30.5           62   ];
t3=[ 20             40              141.0                 208.0              155.2              310                  455             29.5           91    ];
t4=[ 23             30              145.0                  201.0             139.5                385                  535              27           58    ];
t5=[ 20        40              147.0                 228.0             154.3                340                  525             28.5           62    ];                  


t=[t1' t2' t3' t4' t5'];


[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);
p14=pn(:,1:4);
t14=tn(:,1:4);
p5=pn(:,5);
tt5=tn(:,5);



% 网络创建
net=newff(minmax(p14),[S1,9],{'tansig','purelin'},'trainlm');
% 设置训练参数
net.trainParam.goal=0.01;%训练目标误差
net.trainParam.show=200;%训练时候每次展示的结果
net.trainParam.epochs=5000;% 总的迭代次数
net.trainParam.mc=0.95;%动量因子
% 网络训练
[net,tr]=train(net,p14,t14);
% 仿真测试
y1=sim(net,p5);    % BP神经网络的仿真结果
[a] = postmnmx(y1,mint,maxt);



%% GA-BP神经网络

R=size(p14,1);
S2=size(t14,1);
S=R*S1+S1*S2+S1+S2;
aa=ones(S,1)*[-1,1];   % 变量上下项矩阵
popu=50;  % 种群规模
initPpp=initializega(popu,aa,'gabpEval');  % 初始化种群


我在使用initializega 函数过程中,电脑提示错误,



initPpp=initializega(popu,aa,'gabpEval');  % 初始化种群



??? Error using ==> mtimes
Inner matrix dimensions must agree.

Error in ==> gadecod at 29
A1=tansig(W1*p,B1);

Error in ==> gabpEval at 6
[W1,B1,W2,B2,val]=gadecod(x);

Error in ==> initializega at 80
  eval(estr);


请高手指点?

谢谢了 !
最好把怎么修改给我说下 !

谢谢!万分感谢!

[ Last edited by nono2009 on 2009-10-26 at 16:24 ]
回复此楼
安静
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yujunhui

银虫 (著名写手)

最好把怎么修改给我说下 !

谢谢!万分感谢!
安静
2楼2009-10-26 14:49:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yujunhui

银虫 (著名写手)


wuguocheng(金币+0,VIP+0): 问题解决了就好 10-26 16:23
nono2009(金币+1,VIP+0):私下指点的吗?谢谢分享经验。 10-26 16:23
谢谢大家,经过   yuthreestone  书童 的指点 我找到原因了!

经过  高人 的指点 我现在找到原因了.%  
%% 声明全局变量
global p    % 训练集输入数据
global t     % 训练集输出数据
global R     % 输入神经元个数
global S2    % 输出神经元个数
global S1    % 隐层神经元个数
global S     % 编码长度
S1=9;

全局变量的设定 代表的参数 不能乱改!
如果修改了话,在以后的运算过程出现问题 !

程序分享,希望大家 共同进步!


在次感谢 高人!
谢谢!
安静
3楼2009-10-26 15:52:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yujunhui

银虫 (著名写手)


nono2009(金币+1,VIP+0):鼓励分享经验!欢迎常来交流. 10-26 17:54
恩 ,私下解决的!
  分享下经验,希望 大家一起 进步!
安静
4楼2009-10-26 17:39:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yujunhui 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见