QQ登录

只需要一步,快速开始

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

[问题求助] 求 神经网络大神帮助

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

4

主题

12

听众

360

积分

升级  20%

  • TA的每日心情
    开心
    2016-10-11 16:51
  • 签到天数: 138 天

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    跳转到指定楼层
    1#
    发表于 2015-3-20 12:08 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    2体力
    p=[1 2 3;3 2 1;5 6 1;7 3 1;3 5 2;9 3 5;7 4 2;4 5 6;0 3 4;6 0 1;4 2 3;9 2 7;8 0 1;3 6 2];% R+ A8 d* L' L+ e/ N9 w

    * \: i; q) {+ Y, b2 @$ C! L! W* m* C( p- c0 D) @$ m
    a=max(p'); %#ok<UDIM>
    0 Z0 U' d' b/ p7 g4 b# E* bfor i=1:14* {: u3 F  ^5 o5 y, t! |
    for j=1:3: G) U# q- h: X4 @
    ptest(i,j)=p(i,j)/a(i);
    4 P8 ~9 l9 ^! N7 J4 {; W+ O3 Hend
    & s% X  \  [$ S) y9 {% _  t  [end
    9 u# X) M4 g9 C7 w, R
    4 a) }* m1 n2 F' S5 V
    + B; b! X# I1 J1 h1 s+ @# C$ Ap1=ones(1,14);
    # \# L- a9 o4 I; k& M' C3 xp2=0.8.*ones(1,14);
    # w5 t( b& o7 b! t: S& ^9 B( lp3=0.6.*ones(1,14);- w4 p' W% e2 Z- b3 {
    p4=0.4.*ones(1,14);
    " i9 @3 ]9 o% Q  xp5=0.2.*ones(1,14);' V! W9 _/ [  d8 m* B. _2 Y" D( [; W
    p6=0.*ones(1,14);+ B; V2 |) G/ G2 v% V: l
    ptrain=[pl;p2;p3;p4;p5;p6];) W* C: f3 J7 X& a8 |  T& J& R. b
    t=[100,80,60,40,20,0];
    ! B) b9 Z+ ~9 x0 j) _( x9 p$ z/ p! f% G! |+ Y/ b- L
    ! X& w* r( x2 V4 i6 m
    net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
    ( m6 u, j: S6 y, v6 R/ Z% Enet=init(net);
    ; k9 V  U+ x  C- C, r0 [0 b* T  J. dnet.trainparam.epochs=100000;( b6 g) y8 z) e6 }9 [) T6 o
    net.trainparam.goal=1e-10;
    0 M- c1 D/ \6 `[net,tr]=train(net,ptrain',t);8 h) P, N- }# R9 D$ D
    3 |: x; C3 m1 ^! e& {7 F5 ^+ @* ]

    ; ?: z$ r' ~* d$ \for i=1:3, _5 Y1 u; b8 @$ D' B0 C* x: B; |
    a=ptes(:,i);
    9 [* D3 l% c' k- R+ Q/ iscore(i)=sim(net,a); %#ok<SAGROW>
    5 s  {# Q$ j0 [; G2 ]1 ]- gend
    - h' P) B1 Z* x  ], V' R% o' K
    6 Z& y3 K# H- T4 J2 h1 C4 i5 w' z0 Y7 m
    . |! W7 M7 g' D# j+ o( ^4 s

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

    377

    主题

    148

    听众

    5992

    积分

    升级  19.84%

  • TA的每日心情
    奋斗
    2016-7-14 13:54
  • 签到天数: 219 天

    [LV.7]常住居民III

    社区QQ达人 元老勋章 发帖功臣 新人进步奖 优秀斑竹奖 金点子奖 原创写作奖 最具活力勋章 助人为乐奖 风雨历程奖

    群组2015年美赛冲刺

    群组高等数学串讲

    群组2014第三期英语写作

    群组SAS数据分析大赛冲刺

    群组2015美赛优秀论文解析

    回复

    使用道具 举报

    子任        

    4

    主题

    12

    听众

    360

    积分

    升级  20%

  • TA的每日心情
    开心
    2016-10-11 16:51
  • 签到天数: 138 天

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    大神能帮我运行一下代码吗
    " S" N/ q1 H, O, e9 ^( t我不懂这个/ ^* ~; g9 _8 \1 N' J
    能留下联系方式吗
    8 L. }* A  a+ p  k; v8 G4 {
    回复

    使用道具 举报

    QYY1994        

    40

    主题

    17

    听众

    1233

    积分

    升级  23.3%

  • TA的每日心情
    奋斗
    2015-4-4 20:23
  • 签到天数: 126 天

    [LV.7]常住居民III

    这程序是你自己编的吧!你能说说你想做的是什么吗?5 \7 i/ v' Z; ^; P7 P0 s( T" @
    我这里有一个小模板,你可以看看,希望对你有用:
    3 a4 o& O* i! \* _. c! [clc,clear. P8 |3 x+ [# E* n( q& n
    sqrs = [20.55 22.44 25.37 27.13 29.45 30.1 30.96 34.06 36.42 38.09 39.13 39.99 ...2 T( a* H( X# x$ L0 H  Z3 ^
        41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];& d1 n$ @0 w* y$ O4 o
    %人数 单位:万人0 N0 k1 Z, z1 M: c7 [
    sqjdcs = [0.6 0.75 0.85 0.90 1.05 1.35 1.45 1.60 1.70 1.85 2.15 2.20...
    6 e/ H6 N4 n& |, v. n) r: m5 }; {/ U    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];& Z8 `0 Z% B6 U! o8 V2 g) h
    %机动车数 单位:万辆7 {: Q& H* O  T$ P, h) R& A- j4 \
    sqglmj = [0.09 0.11 0.11 0.14 0.20 0.23 0.23 0.32 0.32 0.34 0.36 0.36...
    7 E! e$ d# \: m  @/ s" @( [    0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];
    4 ]( i" O; F) I1 ~' d%公路面积 单位:万平方千米* Q3 l3 Y7 E3 c! {4 d: r2 q
    glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....8 m& W. k  H, v0 [
        19490 20433 22598 25107 33442 36836 40548 42927 43462];
    5 [' f2 U5 K9 z# Z9 k7 ~! N%公路客运量 单位:万人
    8 w3 w: W5 H" o  W5 c4 ?& |& bglhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....0 d' n1 p+ f6 b6 }+ q3 a* [
        11115 13320 16762 18673 20724 20803 21804];5 X, ^9 l( ~1 s2 z8 ]. m3 }
    %公路货运量 单位:万吨% o! G& X) }% M. {
    p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
    " x! L- G# [4 {" T% n" ?( at = [glkyl;glhyl];          %目标矩阵7 x: Q- U! s( J3 {
    1 {) s  d/ G" R5 j/ m# r
    %%归一化处理( t) [3 q" u6 e5 U5 c
    [pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理" Q. Z8 L* Z% s/ K; n8 s
    dx = [-1,1;-1 1;-1,1];                       %5 R1 P/ e8 h( h' k

    + N2 G+ s( z: q+ s& W%BP网络训练
    5 J, t% |2 m( v: _" a" enet = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
    # x3 H9 u8 Z* O6 b* F" J
    $ M( Z6 E( I5 i# L9 b/ \%%建立模型,并用梯度下降法训练
    ; L& u( a3 _1 t1 `5 x% G/ n* x; Inet.trainParam.show = 1000;          %轮回显示一次
    3 i6 X) Q& |* ]% M, hnet.trainParam.Lr = 0.05;            %学习速率. o& x( c- I* H4 M7 G
    net.trainParam.epochs = 50000;       %最大训练轮回
    " f. q& i4 s+ y2 q! b9 [3 w1 Gnet.trainParam.goal = 0.65*10^(-3);  %均方误差
    4 @& |4 O, M3 G0 J* q# h# R  lnet = train(net,pn,tn);              %开始训练
    / x6 u3 p' Z3 `$ y1 _' B3 `, V3 p
    # ~5 n5 \4 W7 B+ O%%利用原始数据对BP网络仿真
    3 c4 e- A8 I% n, V" o( Xan = sim(net,pn);                    %用训练好的模型进行仿真
    & P9 |* G/ \* Y1 v2 _9 f, Ra = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
    - z6 j8 q9 |" q4 F/ J1 s% H$ P9 O6 s4 D* V. k% h5 H
    %%对比测试! a0 l$ Y6 b6 R* Z& x% J3 U
    x = 1990:2009;
    8 d# z% ?# @  O' Anewk = a(1,:);
      p6 k2 U- }1 l, {. dnewh = a(2,:);
    1 `$ i" A& ?+ [& m8 t. x: b7 Vfigure(2);/ @) m1 @$ V3 G3 `8 Y" i. L
    subplot(2,1,1);9 }# H5 k& f  H  S  r: H8 c0 Z
    plot(x,newk,'r-o',x,glkyl,'b--+');3 p  n* H& Y7 k# x  r, O3 A2 S
    legend('预测输出客运量','实际客运量');
    % a  U& x. g- V) [7 G4 N  s9 Y4 Txlabel('年份');* @& v: \- j- m% ]
    ylabel('客运量/万人');
    ) ~" N  t8 o+ K4 otitle('运用数据箱客运量学习和测试对比图');
    & O( }3 c+ R* ]1 `& p" k" Tsubplot(2,1,2);8 O, J. X* B# Y6 Q* G( m
    plot(x,newh,'r-o',x,glhyl,'b--+');( T* n- w6 ~$ {: J% `) d( W9 j0 D; G
    legend('预测输出货运量','实际货运量');
    : E& y2 E, K( M" kxlabel('年份');* F8 b' I; C' r, ]- A" I& K! r
    ylabel('货运量/万吨');
    : R, r' s# i1 ~+ Q; A3 h: C; Ktitle('运用数据箱货运量学习和测试对比图');
    * ]$ G9 Q% \: i5 n, k- e%利用训练好的网络进行预测
    % b. e; ]5 h9 L. e( epnew = [73.39 75.55. H' B! I4 U% q
        3.9635 4.0975
    + ]3 Q" D7 m: e7 O' z    0.9880 1.0268];
    : q; d1 o2 T0 q6 }, p2 [pnewn = tramnmx(pnew,minp,maxp);
    / \: A2 B9 L& p6 x%利用原始输入数据的归一化参数对新参数进行归一化- t# g# N# J4 w) x: e" b
    anewn = sim(net,pnewn);
    ! f8 ]2 S& S' k+ \& z%利用归一化后的数据进行仿真: a+ G7 A8 ~  @+ t# Q6 r1 `8 }
    anew = postmnmx(anewn,mint,maxt)
    ; ~5 N( `, N3 f$ d2 ?; v$ ~%把仿真得到的数据还原原始的数量级, |8 k4 [1 s7 R9 M3 [: }; q

    7 |( {  i. Y# J: m; Y& y& ?- R7 @
    & D, G' m  ^0 B# f" `* l$ I" {+ Y: ~欢迎交流!) A2 \' h- z* b+ f' t9 j
    " W3 L8 d( v+ F- e( N* N
    回复

    使用道具 举报

    子任        

    4

    主题

    12

    听众

    360

    积分

    升级  20%

  • TA的每日心情
    开心
    2016-10-11 16:51
  • 签到天数: 138 天

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    恩恩,我想做的是一个评估!就是对几个对象做一个评估,得出一个分数
    " ~+ Q' b0 n* @0 D8 M& d. D3 K- w非常感谢你* M. I% {/ S. l

    + E9 ~* w) R: O5 B% U
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-6-11 17:45 , Processed in 0.340491 second(s), 77 queries .

    回顶部