求 神经网络大神帮助

子任

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];
( H) F( ^7 z0 [5 A$ s" y) m: l
4 n# l6 j6 d9 R- h' N
6 D4 ]! [7 H3 R: `! q  n) @6 Ia=max(p'); %#ok<UDIM>
; ?' b$ i: u5 O( g# {0 A7 H, }for i=1:14
for j=1:3
. m# w# ?" z: c4 d8 j& z  iptest(i,j)=p(i,j)/a(i);
5 u! H: S5 ~: U( L4 a* Aend
end


& M# S; B: v! Z7 b- ?/ \# j8 hp1=ones(1,14);
p2=0.8.*ones(1,14);
p3=0.6.*ones(1,14);
7 B3 B# ?' F3 k" g6 }' tp4=0.4.*ones(1,14);
p5=0.2.*ones(1,14);
p6=0.*ones(1,14);
ptrain=[pl;p2;p3;p4;p5;p6];
/ k+ }# _2 R7 |  @. c$ a3 ft=[100,80,60,40,20,0];
( b- v/ x( g: H2 [

net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
  f6 ]+ M( l% C6 M. b0 Pnet=init(net);
net.trainparam.epochs=100000;
net.trainparam.goal=1e-10;
[net,tr]=train(net,ptrain',t);


$ U. s. |$ l/ s0 h1 A. t& ffor i=1:3
2 L! B4 }8 J3 O$ ga=ptes(:,i);
1 ], K7 Z  E$ ?8 O  \score(i)=sim(net,a); %#ok<SAGROW>
end

6 j2 P% |, R% l# Q. d7 P
5 H2 {% C0 N$ P' F: K

数学中国—罂粟

你是要问什么啊？

子任

大神能帮我运行一下代码吗
: U, a6 D2 g( q9 t8 G1 [: q我不懂这个
能留下联系方式吗
* \0 A# f5 d' c6 {# O

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
: N! C# |& p7 f1 y  `' R/ ?+ s$ T! V我这里有一个小模板，你可以看看，希望对你有用：
clc,clear
/ {# W0 w8 u% L/ j$ j1 \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 ...
: o" ^% T* o' }" ~8 E! r! S6 d/ w    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
%人数 单位：万人
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...
    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
%机动车数 单位：万辆
" p# t4 N( s# r) isqglmj = [0.09 0.11 0.11 0.14 0.20 0.23 0.23 0.32 0.32 0.34 0.36 0.36...
9 ]( L  L0 Y5 u" m. E7 L& ?    0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];
%公路面积 单位：万平方千米
$ U7 v; ], m) u8 \glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    19490 20433 22598 25107 33442 36836 40548 42927 43462];
8 d9 o& S0 m: f, X( f%公路客运量 单位：万人
glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
- r) o2 z% S& M. x# w( x* P    11115 13320 16762 18673 20724 20803 21804];
! M( c: ]3 w0 f# ^%公路货运量 单位：万吨
p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
6 {9 H& Y6 G: m- O6 ~t = [glkyl;glhyl];          %目标矩阵

; {5 H( ]* |6 B( I9 p4 a- R%%归一化处理
[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
9 I# O, a& @& `5 b& ^9 s/ bdx = [-1,1;-1 1;-1,1];                       %
$ `# O# M/ C" \/ Q5 y+ A$ |6 S' D! q% ]
%BP网络训练
net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');

$ |8 Z1 r2 r' [. i%%建立模型，并用梯度下降法训练
& W4 J9 T! w7 P0 F7 \net.trainParam.show = 1000;          %轮回显示一次
5 u6 B4 v; D5 m3 rnet.trainParam.Lr = 0.05;            %学习速率
net.trainParam.epochs = 50000;       %最大训练轮回
4 b* Q! K% ^2 xnet.trainParam.goal = 0.65*10^(-3);  %均方误差
net = train(net,pn,tn);              %开始训练

9 U, z, A9 f) K% {+ {+ }/ C%%利用原始数据对BP网络仿真
an = sim(net,pn);                    %用训练好的模型进行仿真
( _( V9 P% O9 n; X3 w: j. ]a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级

%%对比测试
$ T! d+ u+ J  F; Y7 lx = 1990:2009;
newk = a(1,:);
newh = a(2,:);
figure(2);
7 x# o" G- @; N! U3 K4 A( Msubplot(2,1,1);
4 [# ^( S7 W$ s% l7 A, qplot(x,newk,'r-o',x,glkyl,'b--+');
legend('预测输出客运量','实际客运量');
xlabel('年份');
ylabel('客运量/万人');
; K4 U" G+ p, e* m) d% @  s8 U1 ititle('运用数据箱客运量学习和测试对比图');
subplot(2,1,2);
( A# f. w: I8 K% I8 m2 Vplot(x,newh,'r-o',x,glhyl,'b--+');
- i3 @  ?# k" L; ~  x; P* l; slegend('预测输出货运量','实际货运量');
xlabel('年份');
ylabel('货运量/万吨');
title('运用数据箱货运量学习和测试对比图');
6 |. Z. I7 T. y. B  J7 l. ~3 J%利用训练好的网络进行预测
pnew = [73.39 75.55
3 d; H+ O; X+ U0 C' j5 n    3.9635 4.0975
1 v+ k/ c, L$ g5 {    0.9880 1.0268];
pnewn = tramnmx(pnew,minp,maxp);
1 ]2 F3 y6 A8 x8 [%利用原始输入数据的归一化参数对新参数进行归一化
anewn = sim(net,pnewn);
0 Y% J0 `% p9 o0 S/ u' ~- ^) k4 K& _%利用归一化后的数据进行仿真
6 {- I5 F0 b/ k. o' ?anew = postmnmx(anewn,mint,maxt)
3 C4 _$ p  v) g# k" p" D( t%把仿真得到的数据还原原始的数量级

' X' y0 ^' d' S% W( u+ m
  l: X; v; I+ M欢迎交流！
0 [& A0 A) j9 W9 J0 A9 n; ?  y
: d7 Y* Z$ u! l2 H$ ?9 E# }0 ^9 @

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
+ w- c' K, _' j( _: }非常感谢你

( ?" M: }6 ^: {7 k/ l; V/ j1 d4 d