24小时热门版块排行榜    

查看: 539  |  回复: 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

银虫 (著名写手)


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的回帖
查看全部 4 个回答

yujunhui

银虫 (著名写手)

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

谢谢!万分感谢!
安静
2楼2009-10-26 14:49: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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上项目filecode邪修 +4 西山十月 2026-08-09 6/300 2026-08-10 00:16 by liufeng0619
[基金申请] 2026国自然放榜时间 +4 布布和一二 2026-08-08 4/200 2026-08-09 22:24 by otj2008
[基金申请] 关于代码变化问题,想知道的进来 +15 且听虎啸 2026-08-07 22/1100 2026-08-09 22:13 by 布布和一二
[基金申请] 这样的filecode谁见过 +10 布布和一二 2026-08-08 21/1050 2026-08-09 17:21 by 天神眷顾
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见