QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3836|回复: 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];
    ; o8 t: F" r7 |$ c! L. k0 p
    & b/ P' h$ U; e$ p
    & U4 A; e" ]( q* I; F: Qa=max(p'); %#ok<UDIM>/ |- @) |. o( t9 u, G- t
    for i=1:14" R0 c  c( o) q1 B3 e# d
    for j=1:3
    ( f& A$ X' K8 |. q6 }2 M- zptest(i,j)=p(i,j)/a(i);
    0 {: U0 g5 J# Y2 f# Nend
    % `7 R0 G* a. K2 J9 J  cend& ^' o$ N& M  Y$ g0 e2 x

    ) V: \1 n9 A: L1 W! ?, s: @! x+ ?1 m2 W% e: M. B
    p1=ones(1,14);
    . r( F( y% {3 _9 y6 o9 h& y+ t. Zp2=0.8.*ones(1,14);
    2 I2 U- U/ K! ~7 Sp3=0.6.*ones(1,14);
      ^3 \9 f9 L7 G' W3 I: Fp4=0.4.*ones(1,14);
    ( a2 d* O0 M6 E6 n2 u( k, Up5=0.2.*ones(1,14);' C* ~# L1 l1 p2 e' w
    p6=0.*ones(1,14);
    : g' ^( K- N* ]8 _0 iptrain=[pl;p2;p3;p4;p5;p6];
    # X4 r; f% x. o8 w' U1 It=[100,80,60,40,20,0];
    $ H1 @" ?6 S! U8 `0 g
      x  s& d' T! f2 c. J* a
    , G0 D- H; \& o' P8 unet=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');" L$ l- p. \6 G
    net=init(net);
    " v% v+ t& {7 O$ p7 rnet.trainparam.epochs=100000;+ J! |* U9 {- @/ {
    net.trainparam.goal=1e-10;
    8 Q2 G0 x4 |1 h[net,tr]=train(net,ptrain',t);" g5 i+ n+ s. f3 J. |
    # D8 d7 d; s6 q
    " z5 @8 B, `8 d$ v2 d
    for i=1:3
    # B' N: Z6 G/ [. H* q4 Ka=ptes(:,i);
    $ A/ @! u& h; O0 Z$ vscore(i)=sim(net,a); %#ok<SAGROW>
    9 t$ O7 R2 Z7 x! Gend
    8 r6 O6 {1 G) g3 W5 K' n
    ; P/ q, H* e" V2 a. s  A0 o# F- y& v0 t8 L" {& a+ @* E" i
      }9 D( t5 o. j5 V% U$ o  C8 s

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

    377

    主题

    148

    听众

    6015

    积分

    升级  20.3%

  • 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

    自我介绍
    比较喜爱数学
    大神能帮我运行一下代码吗
    # i' r" g$ ]/ C我不懂这个# g+ K+ g/ @) M$ @4 m+ [: a
    能留下联系方式吗
    : K8 q% k. b; s( A
    回复

    使用道具 举报

    QYY1994        

    40

    主题

    17

    听众

    1233

    积分

    升级  23.3%

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

    [LV.7]常住居民III

    这程序是你自己编的吧!你能说说你想做的是什么吗?( o, h. G7 j- z/ q1 M: M0 Y
    我这里有一个小模板,你可以看看,希望对你有用:. l$ N" `) `) n: F. i& `$ |2 _5 T" ?$ X
    clc,clear0 F( _5 @9 g+ b0 P
    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 ...1 A8 W3 ?6 t. `) k- Z2 G
        41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
    2 |4 ~7 s5 x0 ~, r0 E%人数 单位:万人( \% B( Q9 A  J: 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...
    & J. H6 M" _# U* G$ F    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];* F8 p$ M* z9 `9 _0 o' N
    %机动车数 单位:万辆8 }- R' l8 [, Y6 D! e2 o
    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...
    4 C3 b9 `% ~1 p  F4 G( B& P    0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];+ i( G$ v7 v+ v8 r  Z  B% w
    %公路面积 单位:万平方千米/ l% `$ j* j6 l. }. r2 }
    glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    7 t: z( P( R1 b    19490 20433 22598 25107 33442 36836 40548 42927 43462];2 B% z4 g- N; Z5 e; Q
    %公路客运量 单位:万人9 \- e0 s' e# s: \  N" t6 C
    glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    2 S" x3 u6 S0 x$ l. C) _( ]" Q$ N    11115 13320 16762 18673 20724 20803 21804];4 g, {4 R! ?; m! j
    %公路货运量 单位:万吨1 G1 q+ N" V1 _% b' r. t0 k
    p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
    ; h- C) o" p! Q0 H, @t = [glkyl;glhyl];          %目标矩阵
    - [+ l4 v; o# v
    + r; H4 {4 B5 L! v" n5 B, N6 M%%归一化处理
    0 i# o% _+ X. @5 `, t+ G2 `! X[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
    2 h) g% S( f2 `" \dx = [-1,1;-1 1;-1,1];                       %
    + t! v% H% b3 z1 g3 s4 U
    " h4 y, Q9 |2 Z8 h  {& U$ @; W$ i+ W7 L%BP网络训练! Q' b/ e# `) |
    net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');& f2 _5 ?, l6 y! u, Q" a

    9 K4 W2 b0 ~" ~%%建立模型,并用梯度下降法训练
    6 W/ c/ a6 A) Z: ]& j7 Znet.trainParam.show = 1000;          %轮回显示一次" \1 @4 u, J& d/ r7 F1 ^
    net.trainParam.Lr = 0.05;            %学习速率
    ! `, p) U4 d+ q1 j% znet.trainParam.epochs = 50000;       %最大训练轮回2 C* E! r& i4 G" n
    net.trainParam.goal = 0.65*10^(-3);  %均方误差7 c' g/ @; t/ J
    net = train(net,pn,tn);              %开始训练8 X# v, m% g% i9 c* ~5 l

    ( i  o" P( ]$ U9 @* O! C2 a%%利用原始数据对BP网络仿真9 g8 ~. r& a4 @/ T$ @. Y: X# x
    an = sim(net,pn);                    %用训练好的模型进行仿真
    6 t" g' ]% ]* f& ^a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
    8 H$ |( _! {! ]$ `
    . M7 m/ R7 a- E/ R4 f( Q: H- H%%对比测试
    / o- j$ u; f) Y! u& }x = 1990:2009;
    6 Z0 i# U( [  u/ o6 O' Lnewk = a(1,:);
    ( U# h, [$ R( Q, n0 Inewh = a(2,:);0 ~- `/ V, ?  T% b6 ?6 \) q
    figure(2);. D1 B$ e; Z9 ^. O3 r
    subplot(2,1,1);- R# D! a  u- k8 d
    plot(x,newk,'r-o',x,glkyl,'b--+');7 ]) I" e0 U; F0 ?* l# P% i: n
    legend('预测输出客运量','实际客运量');" F7 G4 S- ]: T; N5 ?& A
    xlabel('年份');
    , w  t1 Q8 M- w0 ^/ hylabel('客运量/万人');
    , D9 j3 Y1 x& m# C3 V# P0 Stitle('运用数据箱客运量学习和测试对比图');
    $ p& s; N9 `# Csubplot(2,1,2);
    ! g/ W1 V' |) J" |plot(x,newh,'r-o',x,glhyl,'b--+');  J, A: i) p5 K/ j2 M
    legend('预测输出货运量','实际货运量');& P0 V3 U* t( o. K
    xlabel('年份');
    : B1 R. B3 `, G- I/ G. I+ tylabel('货运量/万吨');
    ) N# A4 j' l/ m+ i8 [7 _title('运用数据箱货运量学习和测试对比图');
    : z6 \2 }8 q; _0 P" r5 s%利用训练好的网络进行预测/ E+ ?0 e2 I# N. n: L
    pnew = [73.39 75.55
      Y! D- Z' ]' V; Y) P9 h    3.9635 4.0975* R8 Z; B  v- I8 k; [
        0.9880 1.0268];# B6 B7 |; G) K1 \
    pnewn = tramnmx(pnew,minp,maxp);
    8 P/ Q) X+ n5 e: F%利用原始输入数据的归一化参数对新参数进行归一化6 s6 h& Z( D  A$ @1 a6 S0 v
    anewn = sim(net,pnewn);7 v( A6 t- T8 u5 B6 M' A
    %利用归一化后的数据进行仿真
    3 G' S/ H9 @' [6 j- x' Q5 @2 Xanew = postmnmx(anewn,mint,maxt)
    & q5 z- G. I' n, H# V%把仿真得到的数据还原原始的数量级
    5 W3 \: P$ \& i6 G$ T
    * t+ k6 W, `: L, A  e9 u; @4 P
    7 u# E+ \  z' \欢迎交流!) T$ G: M% y5 r$ `

    # z6 o. }  g0 c. e# y0 B+ S0 c
    回复

    使用道具 举报

    子任        

    4

    主题

    12

    听众

    360

    积分

    升级  20%

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

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    恩恩,我想做的是一个评估!就是对几个对象做一个评估,得出一个分数; I* q. L. \# X9 u" {
    非常感谢你
    5 h( ?& j% k) h  P: Q5 O
    8 ]' B" K: K3 C! l& d
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-9-10 20:09 , Processed in 0.702669 second(s), 75 queries .

    回顶部