求 神经网络大神帮助

子任

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];
% [: s! b% V' D* U% c4 S% e3 b0 P
: S4 }4 P  U5 E: ^& r
1 h; d- Z, ?/ y* `9 g9 _1 \a=max(p'); %#ok<UDIM>
$ J- H( w8 W# K5 q# Hfor i=1:14
for j=1:3
ptest(i,j)=p(i,j)/a(i);
! k1 ?; f. d8 P* Rend
end
* M# m* _! b: _" m& i' I) n! @' _

; o8 q% G$ q$ C& k6 ]5 ]8 F3 V' Ep1=ones(1,14);
p2=0.8.*ones(1,14);
p3=0.6.*ones(1,14);
9 x# d+ P9 _6 ^" `p4=0.4.*ones(1,14);
8 r6 k) a7 t) g5 m( U( ep5=0.2.*ones(1,14);
p6=0.*ones(1,14);
3 p6 W) p$ s) j* Y5 G1 tptrain=[pl;p2;p3;p4;p5;p6];
, I( r; W/ e, H% ?' J; ~' Rt=[100,80,60,40,20,0];
9 S% {/ M- M; K% T

, D; w' O' i# N5 e1 a$ V5 Unet=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
5 \5 q3 B& S, O9 Vnet=init(net);
% X+ J% k, u9 j7 hnet.trainparam.epochs=100000;
6 P) k2 t5 L; N5 k/ l# @$ p3 anet.trainparam.goal=1e-10;
[net,tr]=train(net,ptrain',t);
; x5 ^0 X8 B# X
( C- u) C% m. \+ G$ x
for i=1:3
1 T% H9 ~  K+ Q8 Ya=ptes(:,i);
score(i)=sim(net,a); %#ok<SAGROW>
& \$ D2 D$ c+ Gend


' I& A% K7 C! U2 D0 o/ ]/ H" {; Q# Y

数学中国—罂粟

你是要问什么啊？
7 @% Y! Q% ]8 l) {9 d

子任

大神能帮我运行一下代码吗
% T- l, F! l8 m2 O+ a( X我不懂这个
  W( u, e& b, {' `* A) ?能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
我这里有一个小模板，你可以看看，希望对你有用：
clc,clear
) X  H; N3 J( c1 Y- {2 Ksqrs = [20.55 22.44 25.37 27.13 29.45 30.1 30.96 34.06 36.42 38.09 39.13 39.99 ...
. \; G$ S) Z5 _, Q( D( i- Y    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
6 d! i- @  G* D" o%人数 单位：万人
4 X& K3 U4 K% u1 W+ ]" M+ c: |5 g) A) osqjdcs = [0.6 0.75 0.85 0.90 1.05 1.35 1.45 1.60 1.70 1.85 2.15 2.20...
6 `8 U) }1 v/ O$ P2 e1 b3 v  i    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
' @. q$ B' P1 A5 W7 I% V& {%机动车数 单位：万辆
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...
    0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];
/ ^" F0 s3 ^* L+ X# d4 Y%公路面积 单位：万平方千米
9 w! x0 K5 g) J  u* B5 L/ M/ n2 Jglkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    19490 20433 22598 25107 33442 36836 40548 42927 43462];
%公路客运量 单位：万人
7 x' |' _1 S$ G9 p& dglhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    11115 13320 16762 18673 20724 20803 21804];
%公路货运量 单位：万吨
+ |, K: q8 v  ]0 U4 mp = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
t = [glkyl;glhyl];          %目标矩阵

1 r: g( v' L3 D. }/ F%%归一化处理
) t  u3 V% g7 `[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
& k3 {4 h" [; Fdx = [-1,1;-1 1;-1,1];                       %
  t; E" ~! `1 t! f
( J7 G% [) W. _8 Q%BP网络训练
net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');

%%建立模型，并用梯度下降法训练
; W; [, N; v( z# `net.trainParam.show = 1000;          %轮回显示一次
net.trainParam.Lr = 0.05;            %学习速率
# {  k1 W& T5 s) Pnet.trainParam.epochs = 50000;       %最大训练轮回
net.trainParam.goal = 0.65*10^(-3);  %均方误差
net = train(net,pn,tn);              %开始训练

- u2 m; Y% q4 r; c( [1 G%%利用原始数据对BP网络仿真
an = sim(net,pn);                    %用训练好的模型进行仿真
" H5 `7 V$ l4 [$ j7 ?! Y5 H' ~a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
# @/ U2 Y3 Q9 o; r! e' b2 E6 N
%%对比测试
x = 1990:2009;
: k2 I) _6 e% w3 X# `- J2 e7 Rnewk = a(1,:);
/ _1 R, r2 ]" Qnewh = a(2,:);
figure(2);
subplot(2,1,1);
plot(x,newk,'r-o',x,glkyl,'b--+');
legend('预测输出客运量','实际客运量');
( A6 [2 O6 Y9 C6 r/ A! v2 J3 \( yxlabel('年份');
ylabel('客运量/万人');
title('运用数据箱客运量学习和测试对比图');
, @$ M/ i/ S8 m7 A# a$ e$ Tsubplot(2,1,2);
plot(x,newh,'r-o',x,glhyl,'b--+');
legend('预测输出货运量','实际货运量');
% K, E$ W; `+ L, B2 I+ d$ c1 kxlabel('年份');
ylabel('货运量/万吨');
( ]. f, h) V1 \/ P# @1 Qtitle('运用数据箱货运量学习和测试对比图');
%利用训练好的网络进行预测
) K8 }8 a; y( W0 _pnew = [73.39 75.55
    3.9635 4.0975
) k$ Q6 c* G5 K7 C. L) x3 z3 S    0.9880 1.0268];
( P% i+ f# O2 g! H- S6 B7 i) Dpnewn = tramnmx(pnew,minp,maxp);
$ W1 ]7 w# c9 ~$ |5 z%利用原始输入数据的归一化参数对新参数进行归一化
anewn = sim(net,pnewn);
  q. P. S3 h% s) t" V) l%利用归一化后的数据进行仿真
9 b! X$ F9 w- [anew = postmnmx(anewn,mint,maxt)
3 A7 r/ ^2 c$ \8 Q$ b, r& a0 J- l- S%把仿真得到的数据还原原始的数量级
: E: _( w2 @& R6 L6 W

" e4 s" y" @9 e: `' ^欢迎交流！

1 y; K' g6 _. J) Q4 Z) ?% q, s( m9 X

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
非常感谢你
* N% E' V  g5 d5 ]0 O5 U# F; u