QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3722|回复: 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];1 ?3 W5 T- H0 B5 e

    7 q' ?- Q+ ]. o" g. \* j4 }- ]0 ?7 V# G7 b9 g
    a=max(p'); %#ok<UDIM>3 v, |0 a& ~4 L/ F  ~3 Z: s' m
    for i=1:14
    # {9 P4 p9 x/ j0 r, [, }- jfor j=1:33 A3 M6 T1 |& S) B
    ptest(i,j)=p(i,j)/a(i);
    $ ?9 c$ z1 R) h% _3 `end
    0 h* M1 a0 @' F' r  ~end# H' o8 c5 ?9 p. S

    # h; D! w+ T# Z! ^: G5 g, A  _5 f: \
    p1=ones(1,14);! R  }! R, ~) k1 P
    p2=0.8.*ones(1,14);
    * l9 x5 T0 W6 [% ~% ~' Ip3=0.6.*ones(1,14);
    4 u( z2 U/ p, Ap4=0.4.*ones(1,14);9 H5 s. G* p% {& k& t
    p5=0.2.*ones(1,14);: A0 g% }$ L" ]$ s1 Z
    p6=0.*ones(1,14);2 D5 S) m6 C2 {0 T+ m4 y3 d5 a6 w
    ptrain=[pl;p2;p3;p4;p5;p6];2 V/ N, x9 `# a+ i
    t=[100,80,60,40,20,0];
    ! \$ f2 V; m6 d8 N" t& m8 [6 |: g% o- Y7 x3 {

    : ^- I& W4 d* U7 t( H) z# Tnet=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');; A4 k% m  [  A& U7 p$ z* b8 f- ^
    net=init(net);
    " j9 i/ E: y2 `3 S; v% ?4 Bnet.trainparam.epochs=100000;* c# L! w' [1 z/ M4 C8 v1 Y
    net.trainparam.goal=1e-10;
    7 z$ B, [8 F! k; f[net,tr]=train(net,ptrain',t);0 Z4 N& y3 q+ Y5 o1 y. \
    / P) g- A  ~' K1 U; ]
    6 N; |$ }6 x; e
    for i=1:3* m# f7 @8 v, Z; C6 c' N
    a=ptes(:,i);3 T  R& q/ K2 E  P7 C
    score(i)=sim(net,a); %#ok<SAGROW>
    4 m2 `2 L. a3 e! B) F. Jend0 Z! B7 ?& w' W6 a" R

    9 n+ y0 H% U' \2 n3 y2 [
    3 q5 f% u- M7 ^& Q
    1 y( \2 x! Y  Z

    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

    自我介绍
    比较喜爱数学
    大神能帮我运行一下代码吗
    # r9 q% u: N+ e/ ]我不懂这个
    * B/ o& F& u# h% q' U. ?+ G& G, r能留下联系方式吗
    $ y/ Y7 z, K3 ~/ U+ s
    回复

    使用道具 举报

    QYY1994        

    40

    主题

    17

    听众

    1233

    积分

    升级  23.3%

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

    [LV.7]常住居民III

    这程序是你自己编的吧!你能说说你想做的是什么吗?7 z0 [$ _! q1 j
    我这里有一个小模板,你可以看看,希望对你有用:
    % G, X( ~( }9 Qclc,clear
      s5 o) i2 c6 c6 Ysqrs = [20.55 22.44 25.37 27.13 29.45 30.1 30.96 34.06 36.42 38.09 39.13 39.99 ...& ~+ m3 O6 Y: E3 h& ?  j
        41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
    % P( m9 L* r7 j2 y  L%人数 单位:万人4 J6 A. N% s! P4 d$ y% V
    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...5 H( s! W. n( e1 A# X  N4 B2 G4 m
        2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];$ Y) v3 G& {# j& l
    %机动车数 单位:万辆
    & r7 w5 n. C- @* d* Hsqglmj = [0.09 0.11 0.11 0.14 0.20 0.23 0.23 0.32 0.32 0.34 0.36 0.36...
    4 F6 v& C$ I6 ?8 k) K1 |    0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];
    % P$ @% J' r/ S: v%公路面积 单位:万平方千米
      i  g3 Z8 ~+ X# x! W$ C0 a: Zglkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....) f  K8 e. [5 Y1 S8 g& M, J# t3 p8 G) x
        19490 20433 22598 25107 33442 36836 40548 42927 43462];
    / u1 o/ p# a* x' _%公路客运量 单位:万人
      ?1 b; S7 ~  l; {% Q1 Dglhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....! ?. x8 m+ M  W+ ]
        11115 13320 16762 18673 20724 20803 21804];3 z2 \) C1 _, Y
    %公路货运量 单位:万吨: c) `( a  V) ]
    p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
    7 L% d# q$ F/ ^. y! i7 _: Rt = [glkyl;glhyl];          %目标矩阵
    5 h6 n0 ~$ }0 S  g( x" T9 @
    . E- S5 Q/ l1 h. @, H%%归一化处理# O. y& U; |' W/ V* ~3 `  a& x
    [pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理4 h. |- ~8 q6 X3 X9 i# T
    dx = [-1,1;-1 1;-1,1];                       %
    ) X  c1 |; |' k% G! [
    1 `" S8 X, x# A2 o0 @5 B%BP网络训练
    - o% j; c; b3 ]* ^8 {net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
    , Y5 S. q6 s/ ~, q$ P- j1 O, v; ^+ ~$ }3 A6 F, E+ A
    %%建立模型,并用梯度下降法训练
      g6 ^; E$ y3 m3 \$ Xnet.trainParam.show = 1000;          %轮回显示一次* c6 E  s0 v* G: _+ F6 X$ I' n3 a
    net.trainParam.Lr = 0.05;            %学习速率' Z" X8 ?' S- ~7 \/ E. X. E3 O+ D* I0 _
    net.trainParam.epochs = 50000;       %最大训练轮回
    $ R$ u% t0 e0 m8 q" m- H! s: Xnet.trainParam.goal = 0.65*10^(-3);  %均方误差
    5 O9 q7 @- C, @  w5 w: qnet = train(net,pn,tn);              %开始训练
    8 g" I( s# B3 }0 w
    # S4 l# y$ }! }, E0 j5 g9 K/ F%%利用原始数据对BP网络仿真" L9 H! a& T8 E$ t
    an = sim(net,pn);                    %用训练好的模型进行仿真2 S) f, }& H+ _- b+ P) N* n4 \; p* Z  U
    a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级7 x7 ^( t8 W5 _' q2 {+ @

    # q6 t7 F( K, u%%对比测试4 Q* o- s: Y8 B9 s3 `, Q9 _8 _
    x = 1990:2009;1 R' w* A  c1 [7 K+ a! p
    newk = a(1,:);
    & b3 ~8 ^6 \; S/ Q6 Jnewh = a(2,:);& ?, n+ p, C, a: I, I, X
    figure(2);' e6 t& x# p' j* L+ a$ F( {! j. X
    subplot(2,1,1);
    $ _7 e  {& I0 r. y0 D1 L( L; n: V' {- xplot(x,newk,'r-o',x,glkyl,'b--+');
    ! p+ n: D0 J* u8 k/ m2 K* k- J5 Q! @legend('预测输出客运量','实际客运量');
    3 J: g, q) w3 t1 xxlabel('年份');
    " c2 y+ g$ \- d& D$ b' Uylabel('客运量/万人');
    5 q( {# Q1 p# N- {0 Z) h4 W4 l5 |title('运用数据箱客运量学习和测试对比图');
    ) b% T/ W+ T: U' |% r6 u( |subplot(2,1,2);6 M8 O( x/ P3 |. a' M8 `  u$ A4 s
    plot(x,newh,'r-o',x,glhyl,'b--+');3 F. W* E! K, O/ D
    legend('预测输出货运量','实际货运量');5 ~- T8 s6 }9 X, }4 h7 u
    xlabel('年份');; m9 n# _5 X" S) e$ x/ G) W7 B
    ylabel('货运量/万吨');3 ]& u) Y' a) v
    title('运用数据箱货运量学习和测试对比图');. }' B; U7 ~( T. u5 e
    %利用训练好的网络进行预测
    ; K; W; t7 A6 Z' J6 wpnew = [73.39 75.552 F: N1 b3 P) A+ T* b
        3.9635 4.09758 |: O% z& |% e0 O/ e* H
        0.9880 1.0268];0 H  e. |! U4 v3 x3 n! \# ]+ P1 V1 [
    pnewn = tramnmx(pnew,minp,maxp);
    , Y( Q' v2 a, Y* m" O0 }%利用原始输入数据的归一化参数对新参数进行归一化! Q2 h: f8 D2 L
    anewn = sim(net,pnewn);
    7 Z  T0 s: h6 w0 w. Y%利用归一化后的数据进行仿真4 V; m: o: J6 p" C$ n  [
    anew = postmnmx(anewn,mint,maxt)' a% a$ A7 R1 d5 ]9 L4 r+ i
    %把仿真得到的数据还原原始的数量级  T" T! v' Q1 \, ]1 Q: k
    - |+ h4 T/ X* A1 I! x
    " A/ g2 N' o+ K: C$ w1 G/ E8 x- F
    欢迎交流!
    , f1 e6 q  L- s  ~- s& m8 Y& q4 K$ _
    回复

    使用道具 举报

    子任        

    4

    主题

    12

    听众

    360

    积分

    升级  20%

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

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    恩恩,我想做的是一个评估!就是对几个对象做一个评估,得出一个分数
    1 d6 h; }* [' F1 X非常感谢你
    ! j/ v- J3 K( |6 a6 M* n: g7 B5 `
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-7-27 02:54 , Processed in 0.289176 second(s), 77 queries .

    回顶部