QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3837|回复: 4
打印 上一主题 下一主题

[问题求助] 求 神经网络大神帮助

[复制链接]
字体大小: 正常 放大
子任        

4

主题

12

听众

360

积分

升级  20%

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

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    跳转到指定楼层
    #
    发表于 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];' O9 z+ }+ x! D1 f9 K: Y

    . D5 O3 \7 I4 _' B, s4 j! d5 c$ L
    , S- Q# e5 o: \8 Za=max(p'); %#ok<UDIM>
    $ L- f% r8 T' L( H" U5 I- J' ofor i=1:14# v4 ~9 [; O9 O6 e; r* [8 `
    for j=1:3
    . C+ L, @$ C. {7 p$ s, bptest(i,j)=p(i,j)/a(i); . D/ {+ b/ w5 J, Z
    end
    - u/ s0 \- ]8 I4 ]- O0 u- d# tend6 k* T7 V. r  x. c  }
    / V0 h% w. \' M2 O" j7 k+ T9 S2 g) Y

    ) b# s7 s, h1 J* Q, X! Y" o1 Lp1=ones(1,14);/ |, X! k+ r7 i' z8 |
    p2=0.8.*ones(1,14);
    # z5 E+ b: X; J0 R3 K, O9 ?" |p3=0.6.*ones(1,14);
    4 V5 Z& [" W# P& C3 H  Rp4=0.4.*ones(1,14);% x6 Q" i# `1 ?* A: S7 [$ F, S
    p5=0.2.*ones(1,14);
    3 Y2 a  _" l, ~, J: q2 op6=0.*ones(1,14);# X; Y& r6 l" c; S- a9 E% j3 c
    ptrain=[pl;p2;p3;p4;p5;p6];3 j  F6 [3 X) e' s0 @
    t=[100,80,60,40,20,0];* i4 U4 ]! }' k: P
    9 m% d! g+ h# [0 O7 O$ g* f. n8 r2 d8 r( h

    # `' v6 D+ _$ l. ~net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');+ s$ }& H. W1 R
    net=init(net);
    4 x: ~5 i; `( y3 N3 P: e( ynet.trainparam.epochs=100000;
    : s' \/ p; A0 V# e% E; A/ |net.trainparam.goal=1e-10;
    9 j( a7 n  L8 P) f[net,tr]=train(net,ptrain',t);8 N9 B; j3 W! f5 W& ^1 ^5 U8 o- B
    8 Y% Q3 @) r( J/ \4 q
    ( J: @% A! f5 S% \" P
    for i=1:3
    * B$ h, ~8 N4 E' a( V! d# ca=ptes(:,i);! t3 o1 o5 G) i, q" K
    score(i)=sim(net,a); %#ok<SAGROW>
      @$ i+ K1 X4 Y- V/ rend7 ~# S6 w  t2 b
    . d0 {# R1 R+ k5 s7 {. x8 A
    : ]$ E/ r* h; d5 ]# R/ o  i

    2 r8 M! k' ?. G$ q. C2 M

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

    4

    主题

    12

    听众

    360

    积分

    升级  20%

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

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    恩恩,我想做的是一个评估!就是对几个对象做一个评估,得出一个分数
    : a7 w% b% P- \  a: f9 m0 }非常感谢你
    4 |2 `1 P( `0 Q
    3 c! ?2 z8 t7 U1 r& e& s1 r
    回复

    使用道具 举报

    QYY1994        

    40

    主题

    17

    听众

    1233

    积分

    升级  23.3%

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

    [LV.7]常住居民III

    这程序是你自己编的吧!你能说说你想做的是什么吗?0 U% _5 r: K0 t# s' G1 E& [5 F
    我这里有一个小模板,你可以看看,希望对你有用:$ Y/ T# X; U+ g: d, u7 ?9 }/ B
    clc,clear. h; ?8 E( @4 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 ...& X% M5 E" R0 X8 @$ }
        41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
    6 [1 O* K8 b& C3 A%人数 单位:万人7 ^* U, k( i' t; c! [# F
    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...
    . M1 m/ _) D6 T! |# u( O0 u    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];7 X% J* T( {% V: `) _( o
    %机动车数 单位:万辆/ O6 f& c9 A$ X
    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...* D" a/ p# d- k* N% X7 E
        0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];. W# e( J5 x  n) [( A* `! d: x' P
    %公路面积 单位:万平方千米
    ) ~1 Q9 N% J9 X7 [glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    * m; G9 z1 H: f3 O5 w' F: I    19490 20433 22598 25107 33442 36836 40548 42927 43462];; e3 k$ V- H7 Q7 n3 Y/ k( ]
    %公路客运量 单位:万人
    7 W6 S8 y/ g! W1 O( D! A( Nglhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    - S2 ^# a. d% }+ y1 [0 ?: q    11115 13320 16762 18673 20724 20803 21804];
    5 \% X5 ?) N3 W8 U% L/ H: O%公路货运量 单位:万吨7 @6 P8 m. {* P4 s
    p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵: i1 @, X, [2 J: p/ s2 x
    t = [glkyl;glhyl];          %目标矩阵
      n+ m+ q( e: F3 l
    4 A4 X% r, i4 }7 @* R% E! `%%归一化处理+ @! ]- z* d) D5 G5 U, |' i
    [pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理' A" V: f$ M+ \. X- E
    dx = [-1,1;-1 1;-1,1];                       %' ~! O( `& `5 N3 |* E0 j7 z. g

    6 n3 ~  O9 G: a! S%BP网络训练1 c: |) F- H8 Q# O
    net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');! R; S9 V3 x5 u
    " k% F" L- C9 @4 I% d
    %%建立模型,并用梯度下降法训练1 C( f' ?* J: U7 g% g6 V/ m: r
    net.trainParam.show = 1000;          %轮回显示一次% D0 K4 Y* X7 e, G6 j
    net.trainParam.Lr = 0.05;            %学习速率8 p! f. \  E" W0 x; E
    net.trainParam.epochs = 50000;       %最大训练轮回
    4 @0 v$ b$ x  D' C$ Xnet.trainParam.goal = 0.65*10^(-3);  %均方误差; G4 k  o1 t2 ^
    net = train(net,pn,tn);              %开始训练
      I4 `* v1 ]+ ~- r+ e( |" n" Z; q6 s4 C  m
    %%利用原始数据对BP网络仿真
    , y' c6 T' v) t, _/ z; Xan = sim(net,pn);                    %用训练好的模型进行仿真7 D8 D4 w; |" e% j2 G: I& a
    a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
    * l9 L' N" y. M8 f9 a; l1 B0 ~# Z; e  X) j& H4 S3 B1 ]5 t
    %%对比测试
    1 z* _! m0 I$ p. A# C& g! Qx = 1990:2009;  v. q0 w3 a  k0 F8 M
    newk = a(1,:);
    + ~! H" q/ d0 g" Snewh = a(2,:);! d# r- E1 u% q( Z" C6 e3 M
    figure(2);, Q3 n: a( N) \/ n) I$ w$ r
    subplot(2,1,1);
      \: J' c6 A2 Y3 M5 A- rplot(x,newk,'r-o',x,glkyl,'b--+');
    ; _! ?* u8 p# h3 y7 S9 T7 @3 ]legend('预测输出客运量','实际客运量');
      |  U  s5 }! Kxlabel('年份');7 q  |4 I- _( r$ c, O2 Z, n
    ylabel('客运量/万人');
    9 G- B1 L7 \+ C. utitle('运用数据箱客运量学习和测试对比图');7 m1 U8 x0 J) t! c" ?" U
    subplot(2,1,2);
    ' S; {& G- v; M" ]0 Iplot(x,newh,'r-o',x,glhyl,'b--+');$ g; j* t/ N& q: ~
    legend('预测输出货运量','实际货运量');
    / D. @8 a# [) a( v( Pxlabel('年份');' m$ _" P$ q# L; m2 [9 {7 G
    ylabel('货运量/万吨');
    6 H: y9 y$ y4 p6 h6 X) Q  vtitle('运用数据箱货运量学习和测试对比图');1 o' [. ?& ~% ^5 A1 L+ X
    %利用训练好的网络进行预测4 l( M4 n. A; [3 m3 K$ r
    pnew = [73.39 75.55
    . q/ h: V& t. b( [; b    3.9635 4.0975
    ) S( u/ ~, v0 i1 r' I    0.9880 1.0268];
    , B% E: l0 J( ~; U0 [8 [2 N" s6 f! cpnewn = tramnmx(pnew,minp,maxp);
      F3 f$ M1 C; w' O%利用原始输入数据的归一化参数对新参数进行归一化$ Z0 D7 P* r. B4 p8 r1 q0 a; ]
    anewn = sim(net,pnewn);/ E9 H$ M) T. _/ w
    %利用归一化后的数据进行仿真
    4 N5 |& X9 x6 \4 S3 Y0 }5 ], \anew = postmnmx(anewn,mint,maxt): @* U6 `$ h. {: b
    %把仿真得到的数据还原原始的数量级
    1 O1 b7 z/ z9 V) e+ N
    , F8 d  J$ V2 |: ?8 U1 P# \2 x( j7 _$ B! [! T! P  b
    欢迎交流!
    $ c4 A+ ]) d* E3 y! ?1 E" C2 V0 h8 G" o/ Q
    回复

    使用道具 举报

    子任        

    4

    主题

    12

    听众

    360

    积分

    升级  20%

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

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    大神能帮我运行一下代码吗
    5 `, `( j1 e/ T我不懂这个/ f5 o" k- N. w( i; m! F+ M) Q
    能留下联系方式吗
    $ O, `3 Y/ g5 o% N7 w- i
    回复

    使用道具 举报

    377

    主题

    148

    听众

    6015

    积分

    升级  20.3%

  • TA的每日心情
    奋斗
    2016-7-14 13:54
  • 签到天数: 219 天

    [LV.7]常住居民III

    社区QQ达人 元老勋章 发帖功臣 新人进步奖 优秀斑竹奖 金点子奖 原创写作奖 最具活力勋章 助人为乐奖 风雨历程奖

    群组2015年美赛冲刺

    群组高等数学串讲

    群组2014第三期英语写作

    群组SAS数据分析大赛冲刺

    群组2015美赛优秀论文解析

    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-9-10 20:10 , Processed in 0.492356 second(s), 78 queries .

    回顶部