QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3731|回复: 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];
    + x# Y- X0 E3 J0 O- v: e8 R. @- D* w7 X. E

    , _" b2 C. t7 K7 Sa=max(p'); %#ok<UDIM>
    4 j2 v+ P# e  {+ ?  Y* \4 Afor i=1:14) p% ]4 l# S# B) q9 Q: A# h
    for j=1:3
    4 C% X0 O7 @7 Q# q0 q( T4 j1 [ptest(i,j)=p(i,j)/a(i);
    1 @4 f! c" c0 g& Z0 I. {( ?end: G3 |9 W2 }' A$ S- I. ?
    end
    . t4 _! M6 b, I9 I# {0 F* h
    % I! e1 c6 A4 j7 s
    ' ~  C% G$ ^) r" W  \% }p1=ones(1,14);
    9 M, V: M; O+ L+ T4 W  [p2=0.8.*ones(1,14);
    4 b6 n8 _; g% ]4 |( A) z* Np3=0.6.*ones(1,14);& V3 O$ a/ S, f6 `0 p  C+ j6 v
    p4=0.4.*ones(1,14);6 K1 ~# j& x9 h9 h/ C" K+ }/ f
    p5=0.2.*ones(1,14);
    ; c- n! F3 D# q" }8 ?p6=0.*ones(1,14);
    7 D1 [" `( K3 u* [) E2 xptrain=[pl;p2;p3;p4;p5;p6];  H, ]  {8 y1 q# a  s, \! f
    t=[100,80,60,40,20,0];
    : `2 m- Z6 A, G$ c' [4 C
    ! w5 w) t3 d9 i; t) ]1 @7 ^
    + @% ], q6 h$ E, q9 c7 O& Ynet=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
    + G0 D& }  m, m/ F+ Z; c4 I. ~net=init(net);8 f1 ]) G6 n5 X" b6 i
    net.trainparam.epochs=100000;. P. E! @* J( C* C2 m0 @
    net.trainparam.goal=1e-10;
    3 V/ V/ h( W* R$ W[net,tr]=train(net,ptrain',t);
    : T0 \, ]6 x: u4 i( f- A' ]6 U8 z5 h
    9 j# \  I0 u+ g7 p
    % {  V5 I; [, C* }+ yfor i=1:3
    6 f" q- B/ U) \- X* I7 oa=ptes(:,i);* f6 @  W- U4 w8 I& Q8 E, P
    score(i)=sim(net,a); %#ok<SAGROW>9 |; c" k  }4 Z5 @/ a* J3 F; W0 W
    end2 |9 n" F$ o6 f% [- t& x( n
    # f1 |+ g5 g1 Q  O
    / z. M) O  O$ |0 T
    : `' x, e3 T0 }

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

    377

    主题

    148

    听众

    6002

    积分

    升级  20.04%

  • 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

    自我介绍
    比较喜爱数学
    大神能帮我运行一下代码吗
    8 T- n# j3 r1 v5 |+ u我不懂这个
    , u8 n3 }  l9 P/ t% s* x能留下联系方式吗
    1 ^$ M% _) O8 S! q
    回复

    使用道具 举报

    QYY1994        

    40

    主题

    17

    听众

    1233

    积分

    升级  23.3%

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

    [LV.7]常住居民III

    这程序是你自己编的吧!你能说说你想做的是什么吗?! x, D, D& ?1 K( Q9 S3 `
    我这里有一个小模板,你可以看看,希望对你有用:
    ) E+ n. |) @6 V  Hclc,clear% b; p0 ^( L; m* v$ Q8 o
    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 ...
    * A  M2 [5 L& {" @* F5 L    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
    / c* h8 g% K! t% q( l0 Q& f%人数 单位:万人# ]- t; a5 J4 u. P3 x8 c8 b6 c0 \+ n
    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..." H  b5 [4 P5 ?/ Q5 _
        2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];! _3 j) t0 U6 f( v. s- g
    %机动车数 单位:万辆
    2 N( _& z$ h+ V* Fsqglmj = [0.09 0.11 0.11 0.14 0.20 0.23 0.23 0.32 0.32 0.34 0.36 0.36...* }/ N$ ]2 ^  @
        0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];5 N  p: [3 d8 ]% l
    %公路面积 单位:万平方千米
    * v0 ^4 u% e& L; @7 C% Pglkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024...., z; x4 j! U. ]+ S2 P  X
        19490 20433 22598 25107 33442 36836 40548 42927 43462];. p  H4 r6 a( z
    %公路客运量 单位:万人
    2 }* C) R% a/ P1 pglhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    ) _* }8 T# e1 ~& _7 j3 O8 R    11115 13320 16762 18673 20724 20803 21804];" I0 }1 f# s/ r  y2 _
    %公路货运量 单位:万吨& [# C1 u1 v" v" O/ C/ r" W
    p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵9 G* F! ^! D2 ]5 a' v. p6 e! s
    t = [glkyl;glhyl];          %目标矩阵5 L4 [! I* {. N+ m4 V6 G: p
    # T- ~* y5 U5 _
    %%归一化处理
    1 g3 m+ _% [9 e- X[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理( W8 \. O; z& x% u: a2 b- x" d( d
    dx = [-1,1;-1 1;-1,1];                       %2 W# w7 l& |2 a( l; M! [. W

    9 U- F" u" W+ E%BP网络训练
    1 I3 h7 I1 _4 knet = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
    ! H; I+ s/ q. R( `, F
    5 L2 `  M$ B1 m7 ]9 C8 S%%建立模型,并用梯度下降法训练9 l: q2 a2 v; i) b3 b) }
    net.trainParam.show = 1000;          %轮回显示一次& k7 v1 C( ^( n
    net.trainParam.Lr = 0.05;            %学习速率, R; t6 [  d# d
    net.trainParam.epochs = 50000;       %最大训练轮回
    1 M: g# x4 ?& ~2 \7 ?net.trainParam.goal = 0.65*10^(-3);  %均方误差
    : T  C$ c  w. A( s" F7 v% F* k4 anet = train(net,pn,tn);              %开始训练6 D( o" W4 M1 \( n, f
    6 F6 Q8 _; G) U9 H, C
    %%利用原始数据对BP网络仿真
    9 {  {, w3 n( ]3 }2 Q( ran = sim(net,pn);                    %用训练好的模型进行仿真" ~" C, ~; O3 r$ m& N9 j
    a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
    ) ^" V5 R7 H/ r( q2 Z/ o7 p! ~5 p! Z& m2 P
    %%对比测试& X+ w: N) Z' S' _& j8 d% t) M
    x = 1990:2009;
    , o9 U( a9 B- y" t& R- mnewk = a(1,:);; u/ a. i5 W9 C& \
    newh = a(2,:);
    0 w6 D* I( ?+ F: g2 w! X9 c- ufigure(2);
    ) F  G6 H! f- w2 G# f, Qsubplot(2,1,1);
    9 [) C- J' [0 d" rplot(x,newk,'r-o',x,glkyl,'b--+');" z! M" I" Z3 B2 Q7 E+ k  T* L  ?
    legend('预测输出客运量','实际客运量');0 Y8 K, w* F' b# c( a4 c
    xlabel('年份');
    6 E- l& J! y2 r" }& {) L: Jylabel('客运量/万人');$ w# |  _8 t( m  ^
    title('运用数据箱客运量学习和测试对比图');
    & {( g( n) x# J4 e! ?$ Osubplot(2,1,2);
    % Q/ _3 L" U8 Oplot(x,newh,'r-o',x,glhyl,'b--+');
    - l( B; G5 w4 Ilegend('预测输出货运量','实际货运量');
    9 Y6 X0 E$ B' L9 [: h# Mxlabel('年份');
    5 a" N2 X7 {1 h' Y* nylabel('货运量/万吨');
    ' ^" b+ N) y7 S( K4 Ltitle('运用数据箱货运量学习和测试对比图');
    / f/ e- ^2 @  P, F$ y& K4 P%利用训练好的网络进行预测* O! A! q# m0 x$ a
    pnew = [73.39 75.551 N1 E3 n8 y4 ^  }6 \5 b
        3.9635 4.0975
    , L( T8 n2 f7 e    0.9880 1.0268];
    # P5 M$ J9 ~+ K# Y& zpnewn = tramnmx(pnew,minp,maxp);# p+ @; K- j' f! n. D) I7 ?/ o* w
    %利用原始输入数据的归一化参数对新参数进行归一化* s7 e( r" y; I. I
    anewn = sim(net,pnewn);1 V  Q- c- Q  a; N4 q8 x
    %利用归一化后的数据进行仿真
    $ H" ~$ q$ \+ Janew = postmnmx(anewn,mint,maxt)* R0 a& ?! L0 o0 R0 \' t. t
    %把仿真得到的数据还原原始的数量级
    ' G* ~" w* t. z
    8 j; c. C/ f  E9 F6 |+ Z, i, L; u# r$ M2 f) e, l2 Y+ K
    欢迎交流!/ e! O0 |, F9 i
    ; x$ K4 [1 Y8 h6 ~7 f9 m
    回复

    使用道具 举报

    子任        

    4

    主题

    12

    听众

    360

    积分

    升级  20%

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

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    恩恩,我想做的是一个评估!就是对几个对象做一个评估,得出一个分数
    4 C( \& L  ?! Y) I非常感谢你- P, J7 p) r) F8 z, ]7 k
    3 {! E) A) t4 H) _
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-7-27 12:42 , Processed in 0.546833 second(s), 77 queries .

    回顶部