求 神经网络大神帮助

子任

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];

: H6 h( m/ @0 X. O: u6 }
a=max(p'); %#ok<UDIM>
, Y' Y  |- _6 E& D) bfor i=1:14
/ t$ m* Q( T" E9 t; y* z; l. Yfor j=1:3
ptest(i,j)=p(i,j)/a(i);
) u) n  \, f0 Y8 Wend
6 V* Z& o9 G: E& l/ [, J6 Cend

3 A1 N, P4 l: }; x% Q1 e. m/ ~
9 v9 r2 n8 v2 B2 M2 T" k3 y* j* Yp1=ones(1,14);
( v2 W- ^) p7 G% Y( op2=0.8.*ones(1,14);
p3=0.6.*ones(1,14);
p4=0.4.*ones(1,14);
" @( v2 t7 O9 Q( H" O6 Sp5=0.2.*ones(1,14);
p6=0.*ones(1,14);
2 W4 }3 j: c- H1 K) e9 K  kptrain=[pl;p2;p3;p4;p5;p6];
t=[100,80,60,40,20,0];


net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
, r: e5 }  b' d: ?+ vnet=init(net);
net.trainparam.epochs=100000;
) l+ L0 |9 z; F5 K7 Nnet.trainparam.goal=1e-10;
" O0 j& z& `  Y4 t[net,tr]=train(net,ptrain',t);
2 v! X) o# _# V: O3 }
6 W$ L6 {4 h$ [+ W8 O
for i=1:3
a=ptes(:,i);
score(i)=sim(net,a); %#ok<SAGROW>
( M* Y$ n7 c; aend
& `7 j+ E* C2 f

  n: {6 o8 r: ]+ y: D8 Z! R
: v  l  s2 ~! z" f" c

数学中国—罂粟

你是要问什么啊？

子任

大神能帮我运行一下代码吗
' y( P. `  m6 J7 s0 t) x. t' T& [我不懂这个
能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
我这里有一个小模板，你可以看看，希望对你有用：
clc,clear
7 b! I% ?, q. V- T' V. s. F: wsqrs = [20.55 22.44 25.37 27.13 29.45 30.1 30.96 34.06 36.42 38.09 39.13 39.99 ...
: B, S1 m7 \; A& p6 j& L    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
$ N% D9 P2 |: S- B3 |& 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...
    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
3 {7 m+ X* S+ M%机动车数 单位：万辆
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];
%公路面积 单位：万平方千米
glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    19490 20433 22598 25107 33442 36836 40548 42927 43462];
%公路客运量 单位：万人
( q. R  B6 x4 T- D& @glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    11115 13320 16762 18673 20724 20803 21804];
; N- L" i! X$ `5 a# k* f( S9 ?' ~%公路货运量 单位：万吨
p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
t = [glkyl;glhyl];          %目标矩阵
, K, @  H" U% D1 A6 Z. z9 E; M  V
%%归一化处理
0 m( t" h( B% g/ k7 d: G[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
& B1 H+ S5 [' a" `( C, b* {dx = [-1,1;-1 1;-1,1];                       %

$ |, F: N3 ~6 a%BP网络训练
net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
0 K; h8 v1 Q. N, q
%%建立模型，并用梯度下降法训练
net.trainParam.show = 1000;          %轮回显示一次
net.trainParam.Lr = 0.05;            %学习速率
net.trainParam.epochs = 50000;       %最大训练轮回
net.trainParam.goal = 0.65*10^(-3);  %均方误差
2 ?# J, m$ _  V% J1 vnet = train(net,pn,tn);              %开始训练

( [' |4 N& e- }& u%%利用原始数据对BP网络仿真
an = sim(net,pn);                    %用训练好的模型进行仿真
a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
% G4 R. }1 j! A: y
%%对比测试
x = 1990:2009;
newk = a(1,:);
newh = a(2,:);
figure(2);
5 j  X7 m/ N- q+ c7 B7 f- X% }subplot(2,1,1);
. |( Z0 [/ I+ J# jplot(x,newk,'r-o',x,glkyl,'b--+');
legend('预测输出客运量','实际客运量');
0 g- i1 v" u5 Y( ~/ xxlabel('年份');
ylabel('客运量/万人');
title('运用数据箱客运量学习和测试对比图');
" s- W5 n3 S" S: ?: e" hsubplot(2,1,2);
% ^6 I+ m0 a* N5 tplot(x,newh,'r-o',x,glhyl,'b--+');
legend('预测输出货运量','实际货运量');
2 Y/ l: S& E2 K5 C) S# A8 N7 kxlabel('年份');
ylabel('货运量/万吨');
+ m- m7 P+ @: Y2 stitle('运用数据箱货运量学习和测试对比图');
%利用训练好的网络进行预测
pnew = [73.39 75.55
5 t: G% j4 g+ @4 g  @4 l* I, h8 e1 a    3.9635 4.0975
) P  e* O- t& m$ D( ]    0.9880 1.0268];
pnewn = tramnmx(pnew,minp,maxp);
) G- N6 }! V3 E$ i( R%利用原始输入数据的归一化参数对新参数进行归一化
' f% ^/ a/ a% d9 j; S/ x6 canewn = sim(net,pnewn);
0 ~8 L' z( [4 A1 G, L%利用归一化后的数据进行仿真
anew = postmnmx(anewn,mint,maxt)
%把仿真得到的数据还原原始的数量级
5 z. g6 z5 g6 L2 Q9 m) y4 `5 y
, T( H5 w. n2 t3 G6 B; q
$ q: L8 P. S7 c欢迎交流！

: C# Z% u( t1 @5 g1 M! @

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
6 K5 u& |. F6 c; @+ x, B非常感谢你