QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3833|回复: 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];
    ' J) c/ z. s& k, t& l# J3 E' f% T. n" Q" F& @+ l

    4 c5 r" `1 k  H% ka=max(p'); %#ok<UDIM>( l( l8 D( w& S: o6 @! e/ m' w
    for i=1:14
    & O5 ~, K# J0 Yfor j=1:3
    # R) a  E6 f  f" @) Q, Q1 Tptest(i,j)=p(i,j)/a(i); & W  h1 s4 W+ u
    end2 ^5 V# N7 A+ q
    end
    2 s" a1 Q  h/ I  i% D! {
    ; u1 Z; O7 w2 V/ G3 t/ Z1 {& E7 n% [" `& C; {2 o( ]
    p1=ones(1,14);  p. T$ ~/ O% R% \% ^* Q
    p2=0.8.*ones(1,14);
    4 ~0 f2 [( a: y/ kp3=0.6.*ones(1,14);
    ) `% j1 e9 j( x: fp4=0.4.*ones(1,14);
    % o/ E7 C6 e, W3 gp5=0.2.*ones(1,14);
    & i* [: V" [2 H" K1 e5 b( i! B: W- Qp6=0.*ones(1,14);4 w* T; q) T0 Q/ ^
    ptrain=[pl;p2;p3;p4;p5;p6];
    / u; R  G, S- D! Y# o' Rt=[100,80,60,40,20,0];1 x9 B: J" l! U2 Q" r# k

    3 H9 D9 P* `" e" V- v$ D# e6 ?+ C3 G  Z) X& l  m2 A
    net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');0 y6 o/ B, {( J; a
    net=init(net);
    2 G+ D: v0 }( i; a) D( Knet.trainparam.epochs=100000;
    5 c9 H- I. J8 P1 Y+ t3 Hnet.trainparam.goal=1e-10;% [3 |5 T* c! O/ c( G3 j! `8 L
    [net,tr]=train(net,ptrain',t);
    . ^" o4 ]: c+ }* C$ ~, Q2 V9 T, i& Z

    ) j9 j8 n5 c' i8 J: w! pfor i=1:3
    ; Y6 _* p& j( Ha=ptes(:,i);
    ( _/ E: c* V  ]score(i)=sim(net,a); %#ok<SAGROW>& Q0 v' M8 j. r' R2 [8 G9 K( k5 r
    end$ Z% x/ p, z5 R1 E* Y

      \$ j6 W# A; F( x
    : S+ n3 w) Y$ T9 S, U
    0 a5 z1 h$ Q+ c0 M* B1 q" |& A

    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

    自我介绍
    比较喜爱数学
    大神能帮我运行一下代码吗
    5 Q! |4 G: |: J: }4 g5 _我不懂这个
    ) s1 _7 a' {. c% H能留下联系方式吗8 ]/ g8 h% z; M! T
    回复

    使用道具 举报

    QYY1994        

    40

    主题

    17

    听众

    1233

    积分

    升级  23.3%

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

    [LV.7]常住居民III

    这程序是你自己编的吧!你能说说你想做的是什么吗?0 Y" }, i5 R) Q6 s1 [
    我这里有一个小模板,你可以看看,希望对你有用:1 a3 X. b6 h0 a1 o/ x
    clc,clear
    8 O( }4 c3 a( z! Nsqrs = [20.55 22.44 25.37 27.13 29.45 30.1 30.96 34.06 36.42 38.09 39.13 39.99 ...
    4 F5 V& n' _! y    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
    7 x, ]3 ]0 d+ M$ L0 F( ^) r2 H- |7 M1 P%人数 单位:万人
    : Y, ~6 B5 g2 H6 `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..." o" x0 X9 |$ \, \2 F
        2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];# B; X: v3 D6 ]
    %机动车数 单位:万辆
    $ t8 H- l' L6 P7 n+ Wsqglmj = [0.09 0.11 0.11 0.14 0.20 0.23 0.23 0.32 0.32 0.34 0.36 0.36.... C, R2 c0 s2 ^- m5 q4 ?
        0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];" L1 ?% o0 o5 }' K/ N
    %公路面积 单位:万平方千米/ y4 n' L) H: Y& t) a
    glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....9 n. @1 M; d( a) C* ?  B# o
        19490 20433 22598 25107 33442 36836 40548 42927 43462];
    1 o* _  T+ q7 H( k%公路客运量 单位:万人0 p) I( C8 p* A( z" F
    glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    ! h. s- o) u! z. @' H! X    11115 13320 16762 18673 20724 20803 21804];
    4 l( H) Y7 B- b% {' t* D. ^; D- m%公路货运量 单位:万吨
    ; q6 z/ _, \; e/ C5 K. I# [% n1 ?p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
    : ^0 B: F: n, W! Vt = [glkyl;glhyl];          %目标矩阵
    - {6 ^) L8 T: \/ i" `$ L& c( A/ ]$ @  {
    %%归一化处理$ @( I: \9 p; d3 v) Q( {* n
    [pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理0 j* V/ L0 I5 }! u1 ~: i5 k
    dx = [-1,1;-1 1;-1,1];                       %
    % P1 x* p5 G# }& F  o9 S
    + \6 {, z5 O! o1 f, |& l# y7 ?%BP网络训练; ]* L% `0 B6 U8 l0 k
    net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');: E5 V* W, m4 ~% f# ]+ X5 Z

    * i( ~" K9 ~4 H# U# K7 O%%建立模型,并用梯度下降法训练
    0 ]. }" R4 ~# ]  Jnet.trainParam.show = 1000;          %轮回显示一次# J( S" Y& t2 a1 T  u
    net.trainParam.Lr = 0.05;            %学习速率
    0 E. W( V$ a+ Enet.trainParam.epochs = 50000;       %最大训练轮回
    ! h" u* i- W! z, ?: z, o% Gnet.trainParam.goal = 0.65*10^(-3);  %均方误差
    / C" n! v9 k2 P& v8 i8 dnet = train(net,pn,tn);              %开始训练5 j1 ]0 R) e' B1 q$ l- a  P3 j8 W1 Q

    - ]  L( N8 ]. h3 q# D* z4 B%%利用原始数据对BP网络仿真* \6 v4 f2 l# z* w3 W! N$ F
    an = sim(net,pn);                    %用训练好的模型进行仿真
    $ h7 R' k' a$ M7 b( I3 xa = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
    " Z) O# ], u+ I4 I0 R! F9 v8 s, Y7 l; U% a
    %%对比测试$ y2 i/ O8 }, v, |) y7 y9 q) @
    x = 1990:2009;: X4 B' [$ m4 A' ?! H' C. ~5 f9 V
    newk = a(1,:);9 G4 O  q; d. y" f+ L
    newh = a(2,:);4 w$ _2 ~* i5 V* R
    figure(2);( B* n2 c0 m! X) M8 x
    subplot(2,1,1);  {# J5 G8 T+ |
    plot(x,newk,'r-o',x,glkyl,'b--+');, A( ]: e* I+ K+ K1 F2 f% j
    legend('预测输出客运量','实际客运量');( |, E: z4 T& H- W1 O" q% h
    xlabel('年份');1 l6 ]% c2 D6 P. ?( @5 p
    ylabel('客运量/万人');, K# g* P( b! Y" q
    title('运用数据箱客运量学习和测试对比图');
    & ~0 D. Z8 x% a7 N6 `subplot(2,1,2);
    / V2 Y2 W: ^* C' u' x& |plot(x,newh,'r-o',x,glhyl,'b--+');
    5 a; S, s# R: z' x( g9 D5 i( |legend('预测输出货运量','实际货运量');; a- Z- o3 h& ^* X. e
    xlabel('年份');/ p3 I, D- H  z- `* j
    ylabel('货运量/万吨');  a3 c, U' U* E, Q4 d. E& n4 Z" j- M+ p
    title('运用数据箱货运量学习和测试对比图');
    ; x5 Q% A! y' }%利用训练好的网络进行预测8 t$ _7 n, n  g+ h( I
    pnew = [73.39 75.554 h% l, Y# G, h, i3 \9 _# z
        3.9635 4.0975
    ) O+ y- }3 ]- B  v3 ^- u& j    0.9880 1.0268];
    % A% k" J8 u) Kpnewn = tramnmx(pnew,minp,maxp);2 d0 p; E: l9 }0 ~& S; c; d' k
    %利用原始输入数据的归一化参数对新参数进行归一化/ e9 R1 w, i' R$ w" W& v) F
    anewn = sim(net,pnewn);0 r& h& u9 a) {7 c9 j0 P
    %利用归一化后的数据进行仿真, U- k$ |" H* B1 t5 E4 T
    anew = postmnmx(anewn,mint,maxt)
    9 ], d! X; ]! l% A& }0 D" i; ^%把仿真得到的数据还原原始的数量级+ Q7 z8 C& ~( t% }# e

    ( a" H! W2 F" C9 V5 }
    ' C7 R( `3 M0 t& u" j" K欢迎交流!! W0 |6 `! ?' `
    4 R2 u7 o  \* f
    回复

    使用道具 举报

    子任        

    4

    主题

    12

    听众

    360

    积分

    升级  20%

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

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    恩恩,我想做的是一个评估!就是对几个对象做一个评估,得出一个分数
      [2 d, t) T* Z非常感谢你
    % L- c& ?2 g* w8 p  p  a" c
    $ _0 R+ i* i' K+ k: Y
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-9-10 19:30 , Processed in 0.404994 second(s), 76 queries .

    回顶部