QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3831|回复: 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];
    ' d! Z7 W" M0 E. m4 ~. ^3 {# @/ ]/ Q2 {6 n1 D8 K& x
    5 J9 E5 z# _" F/ N, _  U0 C( `
    a=max(p'); %#ok<UDIM>0 b2 z/ ^) ~, B8 b( K( A
    for i=1:14
    % I' }1 y& z, V* y/ E$ Rfor j=1:39 |( h! G+ u. j
    ptest(i,j)=p(i,j)/a(i);   Q. q# n4 D5 J7 T( H! M
    end
    2 w  w# e" D/ O2 Y, g0 Send
    $ J. T% l5 e+ f- ?3 ?8 p' ?; V' y9 X# B& W  Q, B4 ]
      r# |' S. |' q- A
    p1=ones(1,14);
    5 A/ U& ]4 d) s; Kp2=0.8.*ones(1,14);
    0 G0 V2 Q+ H" W' y) f, C6 A% Np3=0.6.*ones(1,14);
    1 s$ A4 g  D) ^6 \9 D; q0 Cp4=0.4.*ones(1,14);* o7 {0 \5 Q- p( I/ A6 W- i
    p5=0.2.*ones(1,14);/ H8 a! t" ]+ e/ m' h6 y9 ^, r5 k
    p6=0.*ones(1,14);
    7 T5 W8 w! R( B8 N* ~8 n( Jptrain=[pl;p2;p3;p4;p5;p6];: O4 [9 v2 w( u4 i: c( y. q
    t=[100,80,60,40,20,0];1 q0 T) O! h; @% D0 h6 Z
    . O, [0 y# p$ c' W% J' s% g/ q
    : J  G& n- u" D# p0 Q) S( T- B  c
    net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
    ; e0 d% l1 x- g$ @net=init(net);
    4 Z0 {+ a. i/ B" a* Znet.trainparam.epochs=100000;) s0 D' y; u/ ]/ R
    net.trainparam.goal=1e-10;6 z1 N, }8 Q0 \* X* a- ^# l$ F; A
    [net,tr]=train(net,ptrain',t);
    ; k1 _) v) G# v2 H& X. |
    7 c3 S# d+ {8 P% X: A' v* o  B
    8 }( k2 T9 z8 v5 h! X8 Pfor i=1:3& M% W0 u. t- q, X+ K3 `8 y
    a=ptes(:,i);+ t+ S& N) g. A; G% s3 z0 [
    score(i)=sim(net,a); %#ok<SAGROW>" P2 P% c7 r/ d% d- m6 o* h0 e
    end6 f( h9 {( {# z3 h( v

    2 l4 v% B" G0 q9 Z, J# z9 _0 i( m2 r% ~6 C

    , @' n1 D% r) a9 h4 r$ z# 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

    自我介绍
    比较喜爱数学
    大神能帮我运行一下代码吗9 c: W# b' P+ X
    我不懂这个7 \! Q2 I5 G# P9 n
    能留下联系方式吗- e" U- ^6 f! J/ J7 n  c; F
    回复

    使用道具 举报

    QYY1994        

    40

    主题

    17

    听众

    1233

    积分

    升级  23.3%

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

    [LV.7]常住居民III

    这程序是你自己编的吧!你能说说你想做的是什么吗?
    3 I8 Z1 M% T! U5 q我这里有一个小模板,你可以看看,希望对你有用:+ m% w" }! y, T' s, i; q1 M
    clc,clear: D- ^: w! \9 `5 j
    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 ...
    ' _% S" u, F' H$ K    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];9 c0 W7 y# @! F, |
    %人数 单位:万人8 F" R5 ?2 Y0 z$ G( s
    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..." E/ f- V* `; K, O
        2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];/ Y: N& B$ W9 m1 i2 I- l' Q: N9 G7 O
    %机动车数 单位:万辆
    0 P& h& O5 `0 R9 Gsqglmj = [0.09 0.11 0.11 0.14 0.20 0.23 0.23 0.32 0.32 0.34 0.36 0.36...
    2 Z" g/ E) Y2 e5 l" L6 P    0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];" V; x8 ]- Y0 p/ E( q
    %公路面积 单位:万平方千米
    $ J& `% f. W1 N4 e$ H! jglkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    ; ]+ t5 R" _$ h7 E! U) `    19490 20433 22598 25107 33442 36836 40548 42927 43462];
    3 q9 ?. d  ^% y- |! A: T$ g%公路客运量 单位:万人  ]" o2 \5 Q  H! a
    glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....2 Y- F5 Q$ J0 Z/ u0 L0 `/ O
        11115 13320 16762 18673 20724 20803 21804];
    . q0 y2 c. W9 a0 X% W; H%公路货运量 单位:万吨) R: n& G  ]0 ]1 G: Q
    p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵/ @* g$ i. f- l
    t = [glkyl;glhyl];          %目标矩阵
    7 L- m1 |& n7 V* g, [: ~! K8 u1 O% F9 F
    # U- F! |  E% c0 R%%归一化处理5 A; i. Y. u, O; }3 {& j8 u% {
    [pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理6 n# z+ i/ w3 d: Q( ^  U5 {, s
    dx = [-1,1;-1 1;-1,1];                       %5 d6 x# I  W3 r+ F# X6 L4 S+ |

    3 c/ U& O: p8 {  r) c& P; c%BP网络训练
    % L$ U! z, W$ c* C& Onet = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');1 I) f+ w: t# D. N0 v+ T
    8 t/ q' M$ d$ ?1 ]" Q( i' B% }, f" o
    %%建立模型,并用梯度下降法训练) N! O" N' [* h
    net.trainParam.show = 1000;          %轮回显示一次
    ; v) {0 R" i% q2 T9 K$ _net.trainParam.Lr = 0.05;            %学习速率8 K* l' M: a, Y; `
    net.trainParam.epochs = 50000;       %最大训练轮回7 Q6 b& ]. ?# s
    net.trainParam.goal = 0.65*10^(-3);  %均方误差
    . c9 D2 Q. Q7 p. |  hnet = train(net,pn,tn);              %开始训练5 _9 c8 G' R6 ^" a& A3 ~3 j% O+ }8 [
      [# g# Y4 V/ r* n& U/ W) O
    %%利用原始数据对BP网络仿真" p* Q) E' \; {
    an = sim(net,pn);                    %用训练好的模型进行仿真
    0 X5 P* t- r. ~# A8 v' v4 da = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
    : `% ^5 Z1 U2 f$ s$ n% l
    ' X" W, b) J" Q" t%%对比测试( ]( I  @3 @6 \' A& Y
    x = 1990:2009;8 A. ^, o! H8 `; S+ h" U8 E/ i8 c
    newk = a(1,:);3 R6 m7 R7 R* D) n# ?
    newh = a(2,:);
    " v& ~8 n8 }6 ?# \7 ]" C( ~% r7 k& Vfigure(2);
    ; f& \! V2 L8 L) n& Psubplot(2,1,1);8 @, r. v& @: |; b7 `& W4 x
    plot(x,newk,'r-o',x,glkyl,'b--+');' k3 q/ `; y1 c0 @: p3 A' M) z
    legend('预测输出客运量','实际客运量');; _  ~( K5 n* F* g8 l
    xlabel('年份');
    0 [; i* P/ p. `9 Dylabel('客运量/万人');6 n. V' ?" c( m' a  [8 |+ a
    title('运用数据箱客运量学习和测试对比图');
    , j" L! |4 z' d2 ^+ o6 xsubplot(2,1,2);
    5 b: e8 H2 _# Q) Rplot(x,newh,'r-o',x,glhyl,'b--+');: f* b, G9 {/ J5 c4 q
    legend('预测输出货运量','实际货运量');
      l- H) T0 r0 W1 Q1 K% ^" T# wxlabel('年份');
    * i- s8 ~3 K: d5 q( kylabel('货运量/万吨');
    & Q, K& }  S6 y( otitle('运用数据箱货运量学习和测试对比图');
    : T  q! f( j9 |0 E- i- Y9 j( X* m& ?%利用训练好的网络进行预测
    3 E, Q. j$ b! Mpnew = [73.39 75.55# {( D) }- ]2 m% K. X& I( C
        3.9635 4.09754 W/ D: ?5 V- }. b* M
        0.9880 1.0268];* L3 T/ ~8 b2 d) e* c
    pnewn = tramnmx(pnew,minp,maxp);
    ) x3 S8 O' d: x% b$ h6 W; ~% b%利用原始输入数据的归一化参数对新参数进行归一化* ~- ?, ?* K( `4 k' N- i0 n
    anewn = sim(net,pnewn);  }' F3 E; b3 s4 U2 P8 {% C
    %利用归一化后的数据进行仿真
    9 Q3 N( ]. {* t: I: }anew = postmnmx(anewn,mint,maxt)
    ( \* d' ]! ~7 _7 c8 D%把仿真得到的数据还原原始的数量级0 o* l* d: A1 k5 E2 d
    - B4 {, V% ~3 u' L9 C# `2 n

    0 l( G3 }0 f  Q& K. P" y欢迎交流!8 q5 B' f! r9 y3 N" g

    5 ]5 T2 ^4 U# P/ {! b! I6 o8 ?& A, L
    回复

    使用道具 举报

    子任        

    4

    主题

    12

    听众

    360

    积分

    升级  20%

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

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    恩恩,我想做的是一个评估!就是对几个对象做一个评估,得出一个分数/ O6 _# j1 {; N- U! a9 h( `
    非常感谢你
    6 K. T8 G  c( o- o* e0 i
    6 I0 p8 n' P# ^) H
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-9-10 18:11 , Processed in 0.623595 second(s), 77 queries .

    回顶部