QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3579|回复: 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];7 D6 C+ B9 g2 P4 X) M

    . ]9 e$ c. s3 Q* x! p  y0 }) h6 @, E* p5 a
    a=max(p'); %#ok<UDIM># b. k. z, B& x
    for i=1:14
    ' l  g. \- E0 Y7 }, Ffor j=1:3
    8 o3 Z' G8 m) W' Y% |' iptest(i,j)=p(i,j)/a(i);
    # w, |( ?* ~% Z; j7 Mend
    " ~& |8 v- |# c9 Y; [# e. Bend% n0 C/ e; R, T& k6 Z

    . u3 D) ?$ n7 L: Q; N3 K& F  h' L! L' R; m8 j9 K& d. X; p
    p1=ones(1,14);
    $ y( ^9 s/ \5 hp2=0.8.*ones(1,14);
    . T5 {2 q: U7 A& D: _( op3=0.6.*ones(1,14);0 x) `5 l' F0 m! ?8 i: P
    p4=0.4.*ones(1,14);
    8 t9 _1 b! l& c: vp5=0.2.*ones(1,14);. h( o: O( I; e* c+ a- o2 v1 S9 t3 Y8 x; q
    p6=0.*ones(1,14);
    / L. @3 G7 j# N5 e( i- v) Iptrain=[pl;p2;p3;p4;p5;p6];
    / l$ o3 [. b/ _, o  yt=[100,80,60,40,20,0];
    , H8 \7 S+ F+ \
      H+ g: g+ Z3 O7 C- t& F8 S0 ]+ t
    / |) I+ Q5 \+ W; C4 I& ]8 e5 znet=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
    7 k/ O/ i. _3 N+ tnet=init(net);
    3 o5 F- M  U/ V* \& P' k: Snet.trainparam.epochs=100000;
      H  |: p0 u& S0 ?+ L/ onet.trainparam.goal=1e-10;1 a# v9 \% ^5 _; V
    [net,tr]=train(net,ptrain',t);
    8 ~/ w/ J& u& w. R8 m9 t, a; |3 D# h# v% z- Y! V9 e; v

    6 X* a' L8 Q: g0 G& G# gfor i=1:3/ k3 i7 Z" N0 z7 Q1 }* R& d
    a=ptes(:,i);
    1 S: O7 J9 j$ z( j6 U8 |score(i)=sim(net,a); %#ok<SAGROW>
    ; U# ]/ @) m. {! {end
    & a" L& j; t2 k) F5 S2 |. D+ r7 ]& I! \8 v! W9 _

    ; o+ ]4 ^$ V% V. o3 H" t) D/ L# |+ B/ u: T$ z" H# T

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

    377

    主题

    148

    听众

    5980

    积分

    升级  19.6%

  • 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

    自我介绍
    比较喜爱数学
    大神能帮我运行一下代码吗
    4 ^# c9 l# k- V$ P) [- A( J( X我不懂这个3 ?. X$ M( y* p% B& ]0 R
    能留下联系方式吗
    : Z" e* j7 G" l( |
    回复

    使用道具 举报

    QYY1994        

    40

    主题

    17

    听众

    1233

    积分

    升级  23.3%

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

    [LV.7]常住居民III

    这程序是你自己编的吧!你能说说你想做的是什么吗?8 g7 y4 t0 n0 G8 N9 s9 u
    我这里有一个小模板,你可以看看,希望对你有用:
    4 }, t* X8 f' T1 z: cclc,clear( Z& Y* z8 L. H( W1 }
    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 ...
    5 i/ J/ o1 s( {- c" [+ R# s9 Y    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];7 a; f" S" e, W* O" B  f' y- i) O
    %人数 单位:万人( |% O+ z! b. |+ m; k
    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...* |( z- [4 i- E$ H9 _
        2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
    # Q2 A9 `* \, T, C7 C* c%机动车数 单位:万辆
    7 _+ b# s) b* V8 |6 ~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..." @  ?+ u( T9 a9 N
        0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];# Q; Q/ x$ v9 U" D, W$ L
    %公路面积 单位:万平方千米% H! }) r) ~" A8 |' i7 `$ Q
    glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....0 @3 p0 @$ E5 V+ K& G
        19490 20433 22598 25107 33442 36836 40548 42927 43462];
    5 ?7 X+ s$ W% `' i% d' n9 g$ s%公路客运量 单位:万人
    8 w  n; Z! k9 N, L- [. rglhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    : {8 [* J" I2 o* U2 Z1 P% i/ F* Y    11115 13320 16762 18673 20724 20803 21804];8 ?+ e% L- c+ z* _6 f
    %公路货运量 单位:万吨
    # y) `& N- E% X5 V, ]p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵3 r* e9 M  W& a/ P
    t = [glkyl;glhyl];          %目标矩阵
    ( \7 U* q1 A  J; E
      P! E* E, g' ]4 C%%归一化处理* S: k7 u; r& v8 t1 u" p9 A
    [pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理  z6 I) q; x5 G0 [9 i7 F& G
    dx = [-1,1;-1 1;-1,1];                       %/ ]2 @* ]& P, d1 ~$ l

    ) P$ B" }" m! V! V+ w, J%BP网络训练$ ~9 V# ?) D0 a
    net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');8 E: W. c& g5 Y& r; e

    $ a5 E  H9 |  y$ ^' f  w%%建立模型,并用梯度下降法训练
    ) w; I/ ?9 Z! }net.trainParam.show = 1000;          %轮回显示一次- |, K$ V7 D# i8 c+ Y% w# P
    net.trainParam.Lr = 0.05;            %学习速率
    4 {$ P" {/ m, W+ }. W$ }2 @net.trainParam.epochs = 50000;       %最大训练轮回3 V1 {+ g% a" X& x, M9 t
    net.trainParam.goal = 0.65*10^(-3);  %均方误差
    1 N- ]) Q. H# e* D4 }0 Nnet = train(net,pn,tn);              %开始训练- W1 t5 ]% R$ W" M/ u4 A4 l, c

    # w1 `# L2 @  t%%利用原始数据对BP网络仿真4 f5 S! `/ B) W+ W
    an = sim(net,pn);                    %用训练好的模型进行仿真+ I) l" \1 }) M( |  k& ~* @
    a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级( r) {+ A+ f6 C7 \4 a2 k) q

    ) O/ a( S1 O1 C' w# B%%对比测试. u7 Z4 u( X7 t6 u
    x = 1990:2009;
      W; A. W! ?5 K1 b9 y" n" ^- ynewk = a(1,:);
    3 Q4 v" f# b/ C3 `  W' b% \5 Xnewh = a(2,:);
    ( H: T2 y9 S5 J# g$ dfigure(2);
      D9 R1 n. R2 w. A: Tsubplot(2,1,1);  `0 _$ U2 l0 V, z/ d+ Q; K) O
    plot(x,newk,'r-o',x,glkyl,'b--+');
    : [5 b) u1 s5 Mlegend('预测输出客运量','实际客运量');
    ; D. z2 k0 f0 v; O9 l& d/ lxlabel('年份');8 h# [: @- r0 x" j0 e
    ylabel('客运量/万人');2 m7 r3 ]( P) f2 p" f  Z* S
    title('运用数据箱客运量学习和测试对比图');+ m1 T: x: C2 o6 ^  r  Y+ {& g
    subplot(2,1,2);8 t0 L$ L3 A  H. b$ {3 i
    plot(x,newh,'r-o',x,glhyl,'b--+');
    8 m7 s; f6 a6 ^+ Qlegend('预测输出货运量','实际货运量');
    & G) D5 m" Q1 `, S- J2 k7 Txlabel('年份');# O9 Z5 s# p" U5 U2 I8 T
    ylabel('货运量/万吨');
    # V; W! j9 w$ v  ~title('运用数据箱货运量学习和测试对比图');; I4 W5 Q" O4 A: X5 ], P2 j
    %利用训练好的网络进行预测
    # G% x+ r7 C6 D; `pnew = [73.39 75.550 c1 l: [: U; y2 u9 ]
        3.9635 4.0975) N2 ~" l9 K) `; E
        0.9880 1.0268];
    - f: g6 ~. y% p5 T& T. H1 ]pnewn = tramnmx(pnew,minp,maxp);
    ! a" Z; m* C" R" S%利用原始输入数据的归一化参数对新参数进行归一化  @* o! A1 h6 C1 K
    anewn = sim(net,pnewn);
    % x5 ?" G8 s+ i. F9 V%利用归一化后的数据进行仿真6 F$ V; ^  G3 Z' [# O% T/ s
    anew = postmnmx(anewn,mint,maxt)
    , V0 i7 m0 T" R%把仿真得到的数据还原原始的数量级  m/ j' o7 E4 c# j8 |: G, C

    * Q4 ~) s" t' H/ P" g
    , a  H/ m$ _* X欢迎交流!, l) t5 V) r" k. T- k
      P3 J% d, G6 M/ }) `' t# v# J
    回复

    使用道具 举报

    子任        

    4

    主题

    12

    听众

    360

    积分

    升级  20%

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

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    恩恩,我想做的是一个评估!就是对几个对象做一个评估,得出一个分数+ j  E; O$ h. w; w
    非常感谢你9 Q" y: R: m7 ~5 n& k; c
    " Y! ^9 ]3 b; l
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-4-27 10:46 , Processed in 0.410641 second(s), 77 queries .

    回顶部