QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3727|回复: 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];
    , r' Z1 j. Y  n3 S( b; e4 k4 Y: e, x8 `; H

    . W. q" C2 O0 m# I* Aa=max(p'); %#ok<UDIM>
    , Q4 D/ V+ T. J" o! h+ Afor i=1:14. w& D, i0 k$ v; `
    for j=1:3/ u9 H+ y  W1 x- y1 v
    ptest(i,j)=p(i,j)/a(i);
    1 Z' ?) c$ J" ^' xend. ~7 I( g, v$ x0 T# d! f6 t; ]. f! v
    end
    - N6 B- w: N7 U- a
    - V8 R4 |. K# C& @
    ( i( L" g" K2 e% |3 }p1=ones(1,14);
    ; S( F- [' }% Jp2=0.8.*ones(1,14);
    9 K! Q  N6 ?# s7 E9 b, {p3=0.6.*ones(1,14);
    + [5 f: x( ?* N# n4 ep4=0.4.*ones(1,14);* {; g) K6 J* t4 P
    p5=0.2.*ones(1,14);' B5 E, `: B4 p; J2 Z
    p6=0.*ones(1,14);
    / B  }" j8 F- D, F) \ptrain=[pl;p2;p3;p4;p5;p6];# [  U8 E5 x. X) D8 C$ m/ V
    t=[100,80,60,40,20,0];
    8 {0 ~& e, I* |- l6 p/ ^6 @* [7 u# Z6 r+ v; [' H

    0 w  M; S7 R. M- V0 B+ C1 L+ Knet=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
    ' M, z' B7 ~% M" [4 q) Xnet=init(net);
    6 v+ _, T4 w5 nnet.trainparam.epochs=100000;
    ) H/ ~' S; t$ v6 e3 q# T7 _( z# Ynet.trainparam.goal=1e-10;
    1 _( Z3 v4 s$ r, S* J' }; K% A1 @[net,tr]=train(net,ptrain',t);
    4 U; M  f3 q+ j; N& a$ H  Z3 D' d) I

      b; `& T$ w7 L& Y4 d4 G$ jfor i=1:3
    & b$ q6 }! i2 C' o; {' v5 \. y8 Ia=ptes(:,i);
    $ Q, P/ }" P3 k' Pscore(i)=sim(net,a); %#ok<SAGROW>4 i) i' w, V/ ^1 p. X
    end
    ! E) X! L1 Q8 L6 a7 U. T' k. v* J, o0 S" M8 m( M& F  c# ?
    . P# ^8 L4 O  L4 s$ k# P5 L
    0 \1 j* K4 i4 r6 Q

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

    377

    主题

    148

    听众

    6002

    积分

    升级  20.04%

  • 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$ n; N0 f; Z: B9 N  N) D. q+ L我不懂这个
    % r6 j% v, H, T$ S. @. _# h2 p$ u能留下联系方式吗
    " R, @% V8 O7 X& c7 P  k
    回复

    使用道具 举报

    QYY1994        

    40

    主题

    17

    听众

    1233

    积分

    升级  23.3%

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

    [LV.7]常住居民III

    这程序是你自己编的吧!你能说说你想做的是什么吗?
    / _4 {: f# k0 n7 X我这里有一个小模板,你可以看看,希望对你有用:* f+ W  o& t& O1 Z+ E" m5 B5 O- ~
    clc,clear' Y+ m3 w' u9 Z) r1 Z+ `! Q3 z, A
    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 ...
    0 s1 v( Z- S7 d, V6 h# s% Y9 S    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
    6 T$ Z. G5 e1 s) c%人数 单位:万人
    3 |$ l4 V) l5 w- K" H1 ~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...
    ! w. n( Q( h) A# S2 Z    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
    . w6 i$ @: m1 N%机动车数 单位:万辆
    3 v- _2 R4 i6 l/ Ssqglmj = [0.09 0.11 0.11 0.14 0.20 0.23 0.23 0.32 0.32 0.34 0.36 0.36...
    * Q4 O- S. {  Q    0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];
    & L& N+ J% M0 y: @  L' G) a) k%公路面积 单位:万平方千米
    1 t7 G& F% s' E+ Z7 ~0 Y' S  ?2 Bglkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    . [% r9 p2 s1 ~  r5 Q3 v    19490 20433 22598 25107 33442 36836 40548 42927 43462];
    % Y5 X9 h# \$ M' S) c%公路客运量 单位:万人0 m6 F8 O( D/ K0 L) p9 J
    glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    5 U0 S, w8 F& x1 [- P9 T3 _6 a4 r    11115 13320 16762 18673 20724 20803 21804];( \8 N* q) c: a2 l; g  N
    %公路货运量 单位:万吨* G8 d3 F2 C" }2 o# ~# o
    p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
    5 f+ c/ s1 x; k& qt = [glkyl;glhyl];          %目标矩阵/ {0 {, T% Y. u4 l5 w; O; p) `
    $ G) G; ?  ~3 z5 V$ l
    %%归一化处理# S) B" {5 U; z, f3 e5 m& W) v
    [pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理, C, Y1 P% e3 }7 m
    dx = [-1,1;-1 1;-1,1];                       %
      m5 v# p6 J6 w8 f! V) A4 y6 T8 `. }  i9 L2 X  P
    %BP网络训练* u; J5 T0 v6 w) u- K5 `; E" q- Z
    net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
    9 E# U6 T2 n0 _4 u
    2 Z" U" t/ \8 R5 ^' A%%建立模型,并用梯度下降法训练8 W, B! a2 B& n, s' f$ z$ c
    net.trainParam.show = 1000;          %轮回显示一次
    $ W- ~4 e. \8 g: R0 E$ g0 Wnet.trainParam.Lr = 0.05;            %学习速率  a4 W+ _9 c  e6 C/ E
    net.trainParam.epochs = 50000;       %最大训练轮回, K- q2 P5 d& M
    net.trainParam.goal = 0.65*10^(-3);  %均方误差8 C; M& W' W/ i" \! \3 b, E7 M/ r
    net = train(net,pn,tn);              %开始训练
    & n  P& f- P) ]7 Z& [, f+ ^7 _/ ^# |( M( n/ }0 G
    %%利用原始数据对BP网络仿真
    $ _- R" c* B5 @an = sim(net,pn);                    %用训练好的模型进行仿真
    ) J. j1 X/ A$ T+ Ba = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级" u' u$ S% A) ]2 i3 F( O8 l
    $ j- L* S/ W% y; d7 L( @  ~) r
    %%对比测试
    6 k; E. |0 D1 ~& N. {" D$ G( yx = 1990:2009;& |, F3 q% F8 G7 O6 [: E* q/ _& q9 i+ B
    newk = a(1,:);
    1 H$ o  \2 a- H/ W% ]3 c! Z" Fnewh = a(2,:);7 U# y' A  v! P9 L4 Q3 z2 @0 a
    figure(2);
    1 V% q- ?( s2 O! d" i; rsubplot(2,1,1);
    ) P8 J$ V! x9 rplot(x,newk,'r-o',x,glkyl,'b--+');9 S3 ^' y  M  E  z
    legend('预测输出客运量','实际客运量');
    3 O) l7 P+ Y& I" n1 C- a" {xlabel('年份');
    + Z8 z! l! B( i2 n/ I- gylabel('客运量/万人');9 h$ ]. K( \* k# \/ V' A
    title('运用数据箱客运量学习和测试对比图');
    ' Z: A' K1 S6 dsubplot(2,1,2);
    1 j2 j( `" ^8 i. K! }2 ^plot(x,newh,'r-o',x,glhyl,'b--+');" I! h' |4 S7 j" _; s5 i! j; {  O
    legend('预测输出货运量','实际货运量');
    9 T4 C! \% B4 Mxlabel('年份');* ~: z! I* B$ v; N* d
    ylabel('货运量/万吨');  o% m/ F, G& g% B- g# w
    title('运用数据箱货运量学习和测试对比图');
    # f6 n. S( t3 z+ E3 _1 S%利用训练好的网络进行预测0 K8 t# R3 [! d% b. b6 S
    pnew = [73.39 75.55' q9 n7 T" I+ ^  S
        3.9635 4.0975& M4 G% d' D. d- K% `
        0.9880 1.0268];
    6 N( M3 N+ }4 b1 D6 Fpnewn = tramnmx(pnew,minp,maxp);  f2 P- W) M) E% T
    %利用原始输入数据的归一化参数对新参数进行归一化
    2 v5 E% [9 e) y' o- }anewn = sim(net,pnewn);
    + u$ Q7 }1 A$ G' M%利用归一化后的数据进行仿真
    ( o5 x! A# S3 N8 k! canew = postmnmx(anewn,mint,maxt)6 Y) Z2 U9 L2 `' M5 n
    %把仿真得到的数据还原原始的数量级
    4 r: N# E! S7 ~+ a. B/ I
    0 G  A1 W4 y5 O2 C$ w3 A
    8 K5 N) Q+ a5 Q) w欢迎交流!
    2 W# _& a" d) Q7 _; F; `; f  w' c' f* k
    回复

    使用道具 举报

    子任        

    4

    主题

    12

    听众

    360

    积分

    升级  20%

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

    [LV.7]常住居民III

    自我介绍
    比较喜爱数学
    恩恩,我想做的是一个评估!就是对几个对象做一个评估,得出一个分数8 y# K4 L0 Y1 W" O
    非常感谢你: D/ U2 {: y" T9 n' W6 F9 n
    # e4 V% U; V7 p7 ?
    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-7-27 09:43 , Processed in 1.067845 second(s), 77 queries .

    回顶部