QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4017|回复: 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];
    / R6 O7 B) R% ]" l/ \9 t9 U3 E- ]# [! j: K/ R# X' |& w
      E4 w, X) s- C% M& \0 T
    a=max(p'); %#ok<UDIM>/ K1 \3 H) y2 I- f5 j
    for i=1:148 T3 e. n( h/ {: o- D$ I+ M
    for j=1:3
    # b9 q( q3 Z/ f* o4 Xptest(i,j)=p(i,j)/a(i); 0 k/ Q$ g, v, Y' B. ?4 o# w$ ?
    end; t. x! L( B$ t' ~- w+ Y. U
    end
    3 \% y: r  W9 T' E. g
    ( l) u+ ?: ]2 M& |$ l3 l1 R  J/ G4 J; t( k
    p1=ones(1,14);
    ' m5 c: @4 @3 w! S9 Hp2=0.8.*ones(1,14);
    . [' w/ `. H) G1 Ep3=0.6.*ones(1,14);. C0 Z$ J  w* l+ P5 U4 n
    p4=0.4.*ones(1,14);! q2 @9 }1 }6 W; X* W' N
    p5=0.2.*ones(1,14);# \* H6 Y( h: j( T( _% W5 S8 j. x
    p6=0.*ones(1,14);: Q1 Z' V+ t5 n3 l5 e. t
    ptrain=[pl;p2;p3;p4;p5;p6];) C1 t; T7 `/ [
    t=[100,80,60,40,20,0];
    4 p9 ]  D: K8 K* @6 Y4 Q2 E# {, w/ z; z
    ' f3 X0 C$ N( J8 D. c: P
    net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');2 Z% v7 c. x8 B6 I
    net=init(net);
    1 b' n; @9 M+ Unet.trainparam.epochs=100000;( `$ ^2 t& k4 N6 B( |$ p: |
    net.trainparam.goal=1e-10;
    2 \0 q5 n- N3 d5 ]0 \6 l0 v3 G) i[net,tr]=train(net,ptrain',t);1 N! B7 N9 W) g- ?  Z

    * i+ d3 A" J0 e/ I: ^9 {7 W+ t6 f% |3 J; ]# a* R3 ?  u5 o
    for i=1:3& F+ z% c; Y9 m4 U, t" Q
    a=ptes(:,i);
    + G3 g6 d% F8 k  a5 [% [score(i)=sim(net,a); %#ok<SAGROW>; t! i" `! `7 r' I
    end) [! w# a( m/ W) N2 I+ }
    & N3 `: ]' F5 h* E; q& m* _/ V

    " m* T* \) [( F$ `, p: o" z1 l4 \7 w

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

    377

    主题

    148

    听众

    6036

    积分

    升级  20.72%

  • 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

    自我介绍
    比较喜爱数学
    大神能帮我运行一下代码吗: u( z# a' O  G0 `9 q. H( D6 O, F
    我不懂这个  C: M$ A+ x0 K1 ?/ q( s
    能留下联系方式吗4 N2 h( O5 W/ a, P7 ^
    回复

    使用道具 举报

    QYY1994        

    40

    主题

    17

    听众

    1233

    积分

    升级  23.3%

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

    [LV.7]常住居民III

    这程序是你自己编的吧!你能说说你想做的是什么吗?( o7 ~- _' C% `1 ^5 t$ V# C# Z9 b
    我这里有一个小模板,你可以看看,希望对你有用:5 @- G2 ]2 q' B( ?, d3 N- M/ D
    clc,clear0 Q( Q! ]7 \3 i2 N, w
    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 ...
    $ V2 N# P. {! l; _5 B    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
    8 C' a1 I. m" M' X, f%人数 单位:万人3 {  Q. P! e' u, n. r4 W  v" b
    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...
    1 d2 F5 _1 s9 n( b* g    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
    & A* s) a* r; }%机动车数 单位:万辆
    : Z" ^9 A8 ]! t! M3 ^8 I+ t" d! @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..., N; \0 T  g9 |7 ]7 ~) e
        0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];
    # K1 [0 o8 E$ R9 @% z8 W%公路面积 单位:万平方千米7 i+ p  K  j! [( V' l
    glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    ' Y) n- o$ v4 p1 W% j    19490 20433 22598 25107 33442 36836 40548 42927 43462];+ F- k6 W: x5 Q$ |. c$ m
    %公路客运量 单位:万人) e# w, }- b8 O* e! z. W  X/ g2 w# x0 u
    glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524..../ D! @5 V% K' R, H- ~
        11115 13320 16762 18673 20724 20803 21804];
    * l8 z' i* W1 v7 f1 R+ I8 d%公路货运量 单位:万吨7 H% V! T1 ~& P& H5 s; J, @% Y2 ^
    p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
    , q+ t0 Y! F8 X) H! h; W9 w$ vt = [glkyl;glhyl];          %目标矩阵1 ?( t0 {3 \7 r% W4 N

    2 f: j" J! E" _) w7 }! A' }%%归一化处理1 \+ ~* [/ C$ g& s2 E8 R3 m1 n
    [pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
    " P# q* i. g) z' f8 ndx = [-1,1;-1 1;-1,1];                       %
    7 h, Q- L/ |6 T* n3 p$ n! ^: J! i" G1 ?, w7 W( B; u& r. O
    %BP网络训练
    / a1 s; X; A; w/ G. anet = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
    ' x' X" Q8 |8 R, y
    $ I0 Z) A0 Z+ p$ n%%建立模型,并用梯度下降法训练
    # [/ ?; A9 T1 }( v( f( N9 onet.trainParam.show = 1000;          %轮回显示一次
    " ?. [4 ~0 |8 q0 h' C; n# Wnet.trainParam.Lr = 0.05;            %学习速率4 V- p7 M# l! |( O. Y
    net.trainParam.epochs = 50000;       %最大训练轮回3 F" _0 g& I; I8 W% B
    net.trainParam.goal = 0.65*10^(-3);  %均方误差
    ! |$ V* J$ l* R. r. s7 inet = train(net,pn,tn);              %开始训练; a% `" N/ o2 I% e4 Z, Q

    , t. L8 s8 K$ j* _7 W%%利用原始数据对BP网络仿真
    3 b; u( c" n1 |) N0 a/ n$ Zan = sim(net,pn);                    %用训练好的模型进行仿真. V8 k( A$ H# V9 ^3 c$ T% c% k8 S
    a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
    - }+ e; t& {2 r  N  h) _/ Z4 {' f7 z
    %%对比测试& Y- L* Q2 t9 }  o3 E# n
    x = 1990:2009;
    8 n- O5 k, j2 d$ @newk = a(1,:);8 i/ o2 e9 H: x1 E/ d! {
    newh = a(2,:);/ P( t6 w% X. T2 z1 _
    figure(2);# }% i6 L( A( ~6 z9 _* c
    subplot(2,1,1);
    ; b- ]6 G, j) B9 D7 N6 n+ Nplot(x,newk,'r-o',x,glkyl,'b--+');0 U" Z, b/ N) `" y/ h% n2 p
    legend('预测输出客运量','实际客运量');
    8 d3 Z# R% Y: n& y6 j: x+ `xlabel('年份');) Z# g7 M9 f4 @7 ~" x/ w
    ylabel('客运量/万人');6 x% K1 B- j0 k& y5 V9 V
    title('运用数据箱客运量学习和测试对比图');
    1 p' y2 J! b! D+ i& Isubplot(2,1,2);$ W1 b5 Q$ Z9 i) c
    plot(x,newh,'r-o',x,glhyl,'b--+');
    * }/ Y9 G5 R6 E( w* D2 Xlegend('预测输出货运量','实际货运量');3 p1 u/ p1 N0 Q0 G
    xlabel('年份');
    6 W  s) x: W( u3 D" rylabel('货运量/万吨');
    # Z! `3 I) Z% G" N) ?  o9 wtitle('运用数据箱货运量学习和测试对比图');7 w8 I- T' D- @0 E8 p# X, ]5 _
    %利用训练好的网络进行预测
    * }8 W# _$ X; @$ N& Q; s2 Bpnew = [73.39 75.55
    . {0 C& H2 ^1 s( [1 k    3.9635 4.09759 j3 n; Q6 F& a/ J, ?+ p0 [# x
        0.9880 1.0268];7 m( w; ], p3 K2 [% C$ G& d: f1 b; m
    pnewn = tramnmx(pnew,minp,maxp);/ C7 w/ f. C5 Y5 w# C8 v
    %利用原始输入数据的归一化参数对新参数进行归一化" U1 G/ T% `. F3 e+ V( v) b
    anewn = sim(net,pnewn);
    0 _1 z: Z% O1 T; Q* n%利用归一化后的数据进行仿真
    - p9 e1 O1 A8 p% \) R3 E# P. eanew = postmnmx(anewn,mint,maxt)
    , R6 n5 b& N: }% r%把仿真得到的数据还原原始的数量级
    ' H8 u0 j/ L6 _- f2 v
    2 F& b2 |3 ^1 ]5 _# Q" P6 z0 `0 n* {- ^! G9 N5 G- i0 s
    欢迎交流!4 B1 p) w( B3 }: Q

    2 i2 }; B! S2 H: T
    回复

    使用道具 举报

    子任        

    4

    主题

    12

    听众

    360

    积分

    升级  20%

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

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    恩恩,我想做的是一个评估!就是对几个对象做一个评估,得出一个分数
    7 K4 U( S7 t6 ?# O7 ^! U非常感谢你6 y* Y+ U: ^4 u' K6 _
    % F2 Q( |9 {% s" X" A3 `7 Q' N
    回复

    使用道具 举报

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

    fastpost qq
    收缩
    • 电话咨询

    • 04714969085

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-10-31 14:25 , Processed in 0.768984 second(s), 76 queries .

    回顶部