QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3736|回复: 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];
    ( C6 B& N, b  ~
    - ]; t% k; _! M3 F/ e5 b) m' c7 o. t. |8 r; c2 C) q
    a=max(p'); %#ok<UDIM>! T8 ~, M7 d* I' p
    for i=1:14
    + L  V) `( q- e% D- p6 f7 Qfor j=1:3& h* {7 g4 S- j
    ptest(i,j)=p(i,j)/a(i); 9 l# }1 b/ K3 a( t
    end" U8 n4 i& V, w/ H$ k
    end
    # V- o; W" C6 D* B3 B9 A5 q
    ! u6 K3 S, u5 @" g% V& E( L% g! G; j2 k! Q, \
    p1=ones(1,14);
    : X) p( X* C0 yp2=0.8.*ones(1,14);
    / A9 H; k+ p9 M# E3 F- o, i) fp3=0.6.*ones(1,14);4 M- v8 R6 s& B, h( i3 H
    p4=0.4.*ones(1,14);# \/ a, h9 I3 t1 ?" P, `
    p5=0.2.*ones(1,14);
    ! O+ N1 t+ b$ i1 c7 n2 Np6=0.*ones(1,14);
    " y5 Q5 I/ u/ o0 Q7 i2 i: Fptrain=[pl;p2;p3;p4;p5;p6];
    ) I5 S' N3 j' l" Ht=[100,80,60,40,20,0];8 f! q6 x$ T$ H  u
    ) f0 F4 a4 {+ ~- H; Q4 ~

    # V7 Q' \/ T' G; z% [8 l+ cnet=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');/ ]- R9 Q$ D0 M0 K* a* q
    net=init(net);
    0 e$ }5 c" i9 Lnet.trainparam.epochs=100000;( l8 J; b6 X9 U. k. o3 A
    net.trainparam.goal=1e-10;
    3 c! t7 I7 e9 H) r4 J[net,tr]=train(net,ptrain',t);1 ?' M$ v. `4 I

      b: q  E- Z  @# u3 f; v# M' P7 W6 X, ]
    for i=1:38 ~2 U* h' d. P& ~6 K# c# Z' g
    a=ptes(:,i);
    $ G3 c; R) T, m# o$ rscore(i)=sim(net,a); %#ok<SAGROW>% `, d; t- a# Y8 I6 ?
    end
    * d# ]8 j0 Z+ l  k) ~/ g# H8 C$ d& f  |
    " C4 W0 e# L* B0 Y

    , V+ x9 o( C* @% K8 q3 `

    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

    自我介绍
    比较喜爱数学
    大神能帮我运行一下代码吗
    3 ?# S' f& ~4 N/ P+ r我不懂这个% [% p  n6 d7 l5 E! K
    能留下联系方式吗/ c3 V5 K7 x9 A
    回复

    使用道具 举报

    QYY1994        

    40

    主题

    17

    听众

    1233

    积分

    升级  23.3%

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

    [LV.7]常住居民III

    这程序是你自己编的吧!你能说说你想做的是什么吗?  c7 u  y& w  m# w" }7 \4 K& h8 Z
    我这里有一个小模板,你可以看看,希望对你有用:
    ( a8 W5 P: u9 Q  kclc,clear
    5 ~8 S+ I* p* N& g! p, r$ h1 ?4 o( jsqrs = [20.55 22.44 25.37 27.13 29.45 30.1 30.96 34.06 36.42 38.09 39.13 39.99 ...
    ; L- ?3 q8 Z. ^3 y; @' {+ K    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
    . U9 O- f, S6 }# K- c" \0 f$ y%人数 单位:万人
    / n- S! D& Q; D, J# ^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...) R* h7 l' F& D/ @7 G/ w* X
        2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
    + n& k9 h2 K2 S1 U, K" @( G%机动车数 单位:万辆
    5 F+ p3 S) F; bsqglmj = [0.09 0.11 0.11 0.14 0.20 0.23 0.23 0.32 0.32 0.34 0.36 0.36...* h' f( Y. `4 u- N
        0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];
    ' H( ?; G3 \5 V) C%公路面积 单位:万平方千米$ g- S2 I7 J- x# x3 B. n9 P
    glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    , K6 N' `4 m2 t5 n4 G1 z9 x    19490 20433 22598 25107 33442 36836 40548 42927 43462];
      f% }; \) x" s; x2 D%公路客运量 单位:万人3 Z/ n2 H, ^0 r& E
    glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    - ?; [8 r- C6 S5 t+ \    11115 13320 16762 18673 20724 20803 21804];
    ! M5 E& e/ e6 N; ~" C# Q) |. Z%公路货运量 单位:万吨! k) t( }+ |# \8 o% Z" u  D
    p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵8 }, C# i; ^; ^- G
    t = [glkyl;glhyl];          %目标矩阵
    3 ?2 j  U  ^' ~
    # H1 r$ z4 z9 z6 c2 a%%归一化处理
    6 g/ m& y/ K! D7 `# k; }& i( M[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理1 }, X- O: d# l
    dx = [-1,1;-1 1;-1,1];                       %
      b1 |  E2 q6 w1 p7 H$ ^7 y( U/ X3 b6 Y6 s* }+ C& z5 s
    %BP网络训练
    6 g/ _" w: @2 x$ U: }, A8 G# Knet = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');; d5 e- S8 h2 F* K, \

    / x3 Z, J! d9 z' i% f9 a4 D% x0 J%%建立模型,并用梯度下降法训练
    : e. e) i+ B% z' Mnet.trainParam.show = 1000;          %轮回显示一次
    $ N1 \* {/ X8 c* o8 A* I, |3 Onet.trainParam.Lr = 0.05;            %学习速率( K0 s' A. z0 i  W
    net.trainParam.epochs = 50000;       %最大训练轮回( b# c# X4 b( p: s: @
    net.trainParam.goal = 0.65*10^(-3);  %均方误差
    5 O1 N5 r: R+ C- h: e/ Dnet = train(net,pn,tn);              %开始训练; z; O6 `4 q# E- T; ^, v8 a# b& m

    % I% K7 Y$ ~; z* ^: t%%利用原始数据对BP网络仿真
    ( D9 [2 ^/ c9 _! e9 K2 jan = sim(net,pn);                    %用训练好的模型进行仿真
    / u. W% z9 p- |/ }/ U1 t+ Va = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
      N& |6 O- r9 \; u* c7 s6 c+ a9 B; L5 u
    %%对比测试; c- U  G6 r: n
    x = 1990:2009;
    4 X6 X; u% d" k* b% snewk = a(1,:);
    2 k" I8 ?6 a! R/ W! W, b9 N3 \* knewh = a(2,:);- d' r7 z  f- x+ e& e
    figure(2);" c" q$ i+ g: x# k6 p/ I  t
    subplot(2,1,1);9 |7 N9 D# U: G7 O% C" @
    plot(x,newk,'r-o',x,glkyl,'b--+');
    & O. L2 }3 S! [legend('预测输出客运量','实际客运量');
    * P4 B" g' l% kxlabel('年份');! f  W; ^/ d1 n- B, }
    ylabel('客运量/万人');; u; W9 G* ]! C. P( J( ^
    title('运用数据箱客运量学习和测试对比图');7 E. P8 P' y% S6 N. V
    subplot(2,1,2);3 _, l4 p4 u( n8 i0 _4 d5 o7 a
    plot(x,newh,'r-o',x,glhyl,'b--+');
    7 b" F" W' m, O1 D3 e2 T$ Wlegend('预测输出货运量','实际货运量');8 h, f  w4 m/ n  i  W- n
    xlabel('年份');/ r* E  i' a* {: A3 R& R
    ylabel('货运量/万吨');' h+ j% H7 {' m* R
    title('运用数据箱货运量学习和测试对比图');
    2 ]6 I1 g1 ~% E" F3 h! ]: c: {%利用训练好的网络进行预测$ l! Y; s2 a- }$ ?' P: Y4 O
    pnew = [73.39 75.55, N6 x* {( T7 W, ~9 S
        3.9635 4.0975
    ( A7 y9 E, L9 v. D$ {0 g6 p6 @9 P    0.9880 1.0268];) F7 W/ K/ R% b- F2 Q
    pnewn = tramnmx(pnew,minp,maxp);, [$ B' @6 B3 B6 |; F
    %利用原始输入数据的归一化参数对新参数进行归一化$ k/ Y- S# F) `$ }* r% e( k
    anewn = sim(net,pnewn);* [7 F; P" P+ U6 w  j
    %利用归一化后的数据进行仿真( `- }' A: U5 J* |4 C' J
    anew = postmnmx(anewn,mint,maxt)4 t0 O( ]1 n2 k$ O  h% S
    %把仿真得到的数据还原原始的数量级
      O( m/ ~3 h  V7 k: x( C
    6 `" X. _% ~; t$ f- W8 K$ x& c! h+ S+ x: t; k0 P0 ]
    欢迎交流!+ r) U0 `( T4 e# y
    ( V' v2 u! ]9 r$ k* f
    回复

    使用道具 举报

    子任        

    4

    主题

    12

    听众

    360

    积分

    升级  20%

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

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    恩恩,我想做的是一个评估!就是对几个对象做一个评估,得出一个分数7 b; I& s: r8 W% R$ J# g) g
    非常感谢你
    $ X' W2 {5 h2 {2 U" P
    ; e% Z9 B8 }4 w$ N* Q3 q0 M8 g7 X
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-7-28 03:24 , Processed in 0.360630 second(s), 76 queries .

    回顶部