求 神经网络大神帮助

子任

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];
2 `$ }, h3 k( B! a! _' n$ v
$ M% T2 j, l& r( D. ]+ P
+ R/ i3 M- u" M: R" C8 g% {a=max(p'); %#ok<UDIM>
9 c' X, C! H8 L# tfor i=1:14
! A6 V3 a( ]  G4 z" ^: Xfor j=1:3
  ^  _% M$ K. t/ vptest(i,j)=p(i,j)/a(i);
end
end

5 @& m* G9 l4 @0 z; Z9 l
  F2 c. N5 w) gp1=ones(1,14);
7 K0 X; }+ J7 V" c/ [+ R/ hp2=0.8.*ones(1,14);
p3=0.6.*ones(1,14);
+ h, _' Y: W2 {8 W% yp4=0.4.*ones(1,14);
p5=0.2.*ones(1,14);
p6=0.*ones(1,14);
ptrain=[pl;p2;p3;p4;p5;p6];
t=[100,80,60,40,20,0];
5 n9 ^; ?2 U# m! U6 P: x
& h  K% o% P* g- c, @
net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
net=init(net);
net.trainparam.epochs=100000;
net.trainparam.goal=1e-10;
[net,tr]=train(net,ptrain',t);
' V; q7 h  U, i- e9 L, A9 f
2 |7 p5 f, Z' j1 j6 Y5 e: d* p
/ V. [. n: p. F& l. N& P/ Hfor i=1:3
- |3 Z0 Q9 y7 a7 ^4 {9 |a=ptes(:,i);
! E/ c2 F5 X: c# mscore(i)=sim(net,a); %#ok<SAGROW>
+ V# G5 g8 n( R1 C( U+ f* T6 mend


- n( h+ P8 C4 ^) w

数学中国—罂粟

你是要问什么啊？

子任

大神能帮我运行一下代码吗
* g# D( {$ ^5 |9 t% M9 r4 s+ w我不懂这个
能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
9 _, D8 i7 p0 N$ }& A0 q我这里有一个小模板，你可以看看，希望对你有用：
8 K" x( f3 m2 A2 |: B1 ~4 Rclc,clear
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 ...
    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
( n/ B7 F3 l# d* ~0 ]: M" o7 Z%人数 单位：万人
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...
3 |5 @0 L4 x( P( h4 h% q$ v  S& H    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
%机动车数 单位：万辆
6 ]6 n. Z) M+ q' u( b% A) \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];
%公路面积 单位：万平方千米
3 a) V# ~' }% s  Xglkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    19490 20433 22598 25107 33442 36836 40548 42927 43462];
9 }( E7 h9 @# s1 l' D( {5 l( O5 F0 j%公路客运量 单位：万人
glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    11115 13320 16762 18673 20724 20803 21804];
1 s6 |# T* R  L/ t; n  v5 {%公路货运量 单位：万吨
p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
: W) d0 k+ V$ ft = [glkyl;glhyl];          %目标矩阵
& n2 {! p1 Y: c) L. e) ]8 B
%%归一化处理
[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
dx = [-1,1;-1 1;-1,1];                       %
$ Y% u/ {6 x7 l7 A5 z+ o+ Q
, I7 t; I/ ~4 o%BP网络训练
+ M0 ~$ K+ J) m. X" u% _. Gnet = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
5 X* C' B' z7 \# U# E4 |
%%建立模型，并用梯度下降法训练
net.trainParam.show = 1000;          %轮回显示一次
: Z8 f2 V2 F/ Z) p; ]net.trainParam.Lr = 0.05;            %学习速率
net.trainParam.epochs = 50000;       %最大训练轮回
4 t% z- Y; v! P- tnet.trainParam.goal = 0.65*10^(-3);  %均方误差
7 E  Y& y+ @- _0 {3 c' hnet = train(net,pn,tn);              %开始训练
' Z' O  [. J5 X( T( ^+ I2 t
$ A5 P2 A/ p9 b: e; i%%利用原始数据对BP网络仿真
6 X' _  s& s& K" Y( r: k% O0 G6 e$ @an = sim(net,pn);                    %用训练好的模型进行仿真
$ ^, Y% F: K5 D! \. _7 w& ja = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
- k* G$ s+ ?0 b" C+ z8 |
& F9 [: w, O1 a# i# C8 g%%对比测试
/ f9 g: T( _$ p6 }8 G: [x = 1990:2009;
2 C2 n9 v& e! G& K( ^/ Bnewk = a(1,:);
newh = a(2,:);
figure(2);
subplot(2,1,1);
6 K. l" x9 ^1 K& T0 P, y2 @9 ]plot(x,newk,'r-o',x,glkyl,'b--+');
+ ]  R; |0 ~" Z' mlegend('预测输出客运量','实际客运量');
xlabel('年份');
' ^" ~! X3 M& I5 aylabel('客运量/万人');
2 l, l# M, T, D$ L+ Ltitle('运用数据箱客运量学习和测试对比图');
subplot(2,1,2);
plot(x,newh,'r-o',x,glhyl,'b--+');
legend('预测输出货运量','实际货运量');
xlabel('年份');
ylabel('货运量/万吨');
5 D& V8 U) z# ^8 u" Ftitle('运用数据箱货运量学习和测试对比图');
0 C" q. \3 \# f0 F%利用训练好的网络进行预测
2 r: p2 ~4 f0 o' a$ H6 cpnew = [73.39 75.55
    3.9635 4.0975
    0.9880 1.0268];
$ x) R9 m8 _. o0 Q8 f! lpnewn = tramnmx(pnew,minp,maxp);
%利用原始输入数据的归一化参数对新参数进行归一化
8 A  n; J' z' H9 ]anewn = sim(net,pnewn);
%利用归一化后的数据进行仿真
0 b0 K8 B5 }% ]) u, N1 x  fanew = postmnmx(anewn,mint,maxt)
  M5 V. h  F+ i9 L5 E, U0 e%把仿真得到的数据还原原始的数量级
" U) A$ Z; J. V: n/ c2 F9 b
5 e# g+ L9 R; C0 g
欢迎交流！

$ z- w: E! I; t5 x$ R+ f

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
非常感谢你
  X4 {3 E* U( T) i2 B( U0 k8 \