QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 10446|回复: 15
打印 上一主题 下一主题

在MATLAB中如何将神经网络和遗传算法联系起来

[复制链接]
字体大小: 正常 放大
trok_2005        

3

主题

2

听众

61

积分

升级  58.95%

该用户从未签到

新人进步奖

跳转到指定楼层
1#
发表于 2005-1-22 22:42 |只看该作者 |倒序浏览
|招呼Ta 关注Ta

我因科研需要刚开始接触MATLAB,发现确实如各位大侠所说是个好东西,目前我在处理一组实验数据,希望用神经网络对其进行训练,可以实现仿真,另外我还希望采用遗传算法寻优,如何将神经网络的输出转化成遗传算法需要适应度函数呢?请各位大侠多多指点,先这里谢过了。

zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
retin        

25

主题

5

听众

679

积分

升级  19.75%

该用户从未签到

新人进步奖

回复

使用道具 举报

ilikenba 实名认证       

1万

主题

49

听众

2万

积分

  • TA的每日心情
    奋斗
    2024-6-23 05:14
  • 签到天数: 1043 天

    [LV.10]以坛为家III

    社区QQ达人 新人进步奖 优秀斑竹奖 发帖功臣

    群组万里江山

    群组sas讨论小组

    群组长盛证券理财有限公司

    群组C 语言讨论组

    群组Matlab讨论组

    希望你能将如何用神经网络做的仿真说清楚,既然已经仿真成功了,那所说的寻优又是指的对什么寻优呀?神经网络又作了怎样的仿真?问题不清楚,这个很难回答!

    回复

    使用道具 举报

    trok_2005        

    3

    主题

    2

    听众

    61

    积分

    升级  58.95%

    该用户从未签到

    新人进步奖

    以下是引用ilikenba在2005-1-23 12:26:35的发言:/ Y- q7 E) s. k' M0 O% b$ n1 }1 M 1 u( S+ N* i/ o6 w' ^" O7 v

    希望你能将如何用神经网络做的仿真说清楚,既然已经仿真成功了,那所说的寻优又是指的对什么寻优呀?神经网络又作了怎样的仿真?问题不清楚,这个很难回答!

    1 H- S6 K1 B! g9 k- a

    谢谢您的关注,是这样,比如说我利用正交设计获得18组实验数据,我用BP神经网络建立了网络模型,我还希望能利用这个模型求解最优,即寻找合适的网络输入,以使网络输出值达到最大。但由于神经网络模型的建立是通过神经元间的连接权值和阈值来实现的,很难用简单的函数形式表达,因此我希望采用遗传算法来寻找最优的实验条件。但我不知道怎样将神经网络的输出转化为遗传算法所需要的适应度函数。 . W6 |9 j# ^" v g5 V, o4 A

    我是用NET=NEWFF(MINMAX(P),[20,1],{'TANSIG','PURELIN'},'TRAINLM');进行建立网络,用NET=TRAIN(NET,P,T)对网络进行训练,用ANEW=SIM(NET,P)对网络进行仿真。 , ~+ b4 D4 c- B/ c) t) ?: F. a

    不知我的问题说清楚了没有,希望您或其他大侠帮忙看一下,小妹这里谢过了!

    回复

    使用道具 举报

    ilikenba 实名认证       

    1万

    主题

    49

    听众

    2万

    积分

  • TA的每日心情
    奋斗
    2024-6-23 05:14
  • 签到天数: 1043 天

    [LV.10]以坛为家III

    社区QQ达人 新人进步奖 优秀斑竹奖 发帖功臣

    群组万里江山

    群组sas讨论小组

    群组长盛证券理财有限公司

    群组C 语言讨论组

    群组Matlab讨论组

    按照你的说法,你可以完全利用遗传算法寻优,并不一定要和神经网络结合的,这里只要找到一种合适的编码规则就行。对于遗传算法你完全可以随意地利用一种选择、交叉、变异规则,因为它是按照概率1收敛的,所以这些规则的选择好坏只影响收敛速度,我想你所说的适应度函数就是表示每个种群的概率密度函数吧!你完全可以按照网络输出值的数值大小计算每个种群所占的概率!
    回复

    使用道具 举报

    trok_2005        

    3

    主题

    2

    听众

    61

    积分

    升级  58.95%

    该用户从未签到

    新人进步奖

    谢谢 ilikenba 的回复,我是新手,不太明白,研究了半天,编了一点东西,可运行出了问题,帮忙看看是怎么回事,如何修改,谢谢!

    p=[ 40 21 2.5 6. 5 ; l- C0 f. W# l& T' D, F 40 25 3. 6.5 10 9 X% K" Y9 i0 Q3 F 40 30 3.5 7. 20 ) }' Z6 s2 p) a" i/ y 45 21 2.5 6.5 10 - I& O& [ r5 [% V( d7 g0 t# a0 j 45 25 3. 7. 20 4 L# h- T0 j" v" U/ U, o 45 30 3.5 6. 5 6 j: O- X) j7 s( `; A2 V/ c+ s 50 21 3. 6. 20 5 b, R4 }" M% ] 50 25 3.5 6.5 5 2 z* D; o6 ^/ N9 Y X' P 50 30 2.5 7. 10 * p: s/ L$ p6 I 40 21 3.5 7. 10 & s" E7 M% {, M! U; f7 ~# U" ] 40 25 2.5 6. 20 / s: C2 U; X/ w5 W0 i 40 30 3. 6.5 5 7 p7 n/ ` Y4 g6 J- o 45 21 3. 7. 5 7 T6 s3 K. H8 p. ? 45 25 3.5 6. 10 / t' q; f. X7 N! V 45 30 2.5 6.5 20 ( S# S+ ^8 o/ c6 C, X3 v) E. e 50 21 3.5 6.5 20 ( a7 W; V5 \& t! Y 50 25 2.5 7. 5 - L- K' c/ h! U4 K: G; o 50 30 3. 6. 10];

    t=[2.673;3.412;1.315;2.019;1.996;0.885;9.317;4.484;1.768;5.720;2.316;0.678;1.042;1.492;0.765;8.812;1.204;3.130]; ) O3 | e( Y2 l4 l( bT=t';' ^( E+ I0 a9 t0 @0 F% i0 V P=p';# w) _+ w8 d; G* U% c+ Z' T net=newff(minmax(P),[12,1],{'tansig','purelin'},'trainlm');+ t/ l0 \8 F! e$ J* C8 X %训练网络

    net.trainParam.show=10;( Q: ]. Y( ^! X+ @8 t %net.trainParam.lr=0.05;7 J. L# v" M; C# X% i9 } %net.trainParam.lr_inc=1.05;/ n v! H6 V- p6 z6 H net.trainParam.epochs=10000; + k+ W r- {# Z+ ^4 ~7 j# R0 @5 r5 M1 enet.trainParam.goal=1e-5; ! b3 g4 y, J h: {; Y % randn('seed',192736547); 6 F# n* q2 Z f0 Q, M' o% net=init(net);$ [- B+ d# J% P5 _0 h. W [net,tr]=train(net,P,T);; t( I* g* U) j. `9 Y minmax_var=minmax(P);9 D( _& @3 K% @ minmax_target=minmax(T);/ t, s; N5 x+ `& k% Y0 t save('result','net','minmax_var','minmax_target');

    % 将网络输 转换成5 I& }' h# v3 U, u9 X$ @% c% ]1 n) l* t + g8 `/ K4 Q, r0 L- l load('result'); 7 d& z4 ~( f; I4 \& b: ~ [Pnew,endPop,bestSols,trace]=ga(minmax(P),'fitness');5 }/ J; `& |. M7 r : X+ k e) Z8 l5 O( \9 I( S %性能跟踪6 D7 R' R! J [& u! F6 M plot(trace(:,1),trace(:,3),'y-');1 [; T" k5 G- h hold on / V, ~" d2 w. r& h plot(trace(:,1),trace(:,2),'r-'); - D4 E* h2 P2 h) V$ R! Q xlabel('Generation'); 6 b4 H' q+ W9 {7 H+ Y2 D# ^/ u ylabel('Fitness');" ^, S* p2 P( b9 y+ @7 J legend('change of solution','average change of population');: _ {8 }% y3 N1 m TRAINLM, Epoch 0/10000, MSE 12.2801/1e-005, Gradient 1739.63/1e-010 ; M. l( q# l6 g3 Z9 S# oTRAINLM, Epoch 10/10000, MSE 0.694955/1e-005, Gradient 93.6508/1e-0108 F/ {4 b) U! }6 k TRAINLM, Epoch 20/10000, MSE 0.0242391/1e-005, Gradient 2.89095/1e-010) E( b* N- J# X6 u$ ~+ ~2 O9 R8 } TRAINLM, Epoch 30/10000, MSE 0.0206875/1e-005, Gradient 4.2655/1e-010 : e2 u N. }3 s- H2 dTRAINLM, Epoch 40/10000, MSE 0.0185878/1e-005, Gradient 18.249/1e-0108 s/ x( m+ {3 r8 ^/ @ TRAINLM, Epoch 50/10000, MSE 0.00947447/1e-005, Gradient 55.0854/1e-010 5 ^) v* Y' a7 [7 LTRAINLM, Epoch 53/10000, MSE 1.24279e-006/1e-005, Gradient 0.504667/1e-0101 F( H% [" ]5 A& i0 B TRAINLM, Performance goal met.

    ??? Undefined function or variable 'minmax_target'.

    Error in ==> D:\MATLAB6p5p1\work\fitness.m * l% p( F* n$ O9 |( H5 oOn line 2 ==> min_target=minmax_target(1);

    Error in ==> D:\MATLAB6p5p1\work\initializega.m 7 g+ n: E( m5 m8 @6 ROn line 41 ==> eval(estr);

    Error in ==> D:\MATLAB6p5p1\work\ga.m 8 n, h+ D- q/ y( ^. a5 P6 B1 d r* n9 \* LOn line 148 ==> startPop=initializega(80,bounds,evalFN,evalOps,opts(1:2));' f6 z0 P) ?+ Q: F r5 D! Y- D 所使用的适应度函数是

    function [sol,eval]=fitness(P,options)3 m4 X }, i# r+ [, C min_target=minmax_target(1);. h) |8 O1 \0 _ max_target=minmax_target(2); ) [- Q9 V' W8 f8 {5 F- B5 U+ A eval=sim(net,P) ' j9 P: h5 q' _6 f; A9 f" d if isformax' s; Q$ Q/ Y' s! W, H! U J7 X eval=eval-min_target+(max_target-min_target); 7 E2 J6 a: [, g+ G else * e! {, j- {0 i' M- [- B, ^- p2 h# x eval=-eval+max_target+(max_target-min_target);, G0 j* v& W5 G$ h% j end

    回复

    使用道具 举报

    cdqxc        

    0

    主题

    0

    听众

    4

    积分

    升级  80%

    该用户从未签到

    回复

    使用道具 举报

    lzh0601        

    2

    主题

    3

    听众

    89

    积分

    升级  88.42%

    该用户从未签到

    回复

    使用道具 举报

    dsh1001        

    0

    主题

    1

    听众

    2

    积分

    升级  40%

    该用户从未签到

    我也是你这种问题,就是通过实验获得数据,然后建立模型进行训练、仿真和预测,这些都可以" X+ q( O: p0 q* d& f5 L4 x* s0 L
    就是最后如何实现参数的优化,就是给定了输出如何确定最优的输入
    . _' ?" O1 L4 b; R, Y! N好像是用模拟退火算法之类的可以实现,目前还没有解决,期待你解决以后多交流
    回复

    使用道具 举报

    0

    主题

    4

    听众

    5

    积分

    升级  0%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-8-25 02:52 , Processed in 0.642584 second(s), 106 queries .

    回顶部