求 神经网络大神帮助

子任

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

8 s, L" k- B) D" V6 L' m
a=max(p'); %#ok<UDIM>
for i=1:14
for j=1:3
1 H( j9 V' V# k) }ptest(i,j)=p(i,j)/a(i);
+ r$ L$ }8 ^8 D2 n0 E2 |end
* i7 C) h2 v/ g, q  rend
+ T* h  P9 D) o5 }
. d+ ?& [( C3 k8 b
* T: Y1 E9 m/ f2 |p1=ones(1,14);
& o7 S' I, P9 fp2=0.8.*ones(1,14);
p3=0.6.*ones(1,14);
p4=0.4.*ones(1,14);
2 N0 k7 c. k4 ~9 \0 yp5=0.2.*ones(1,14);
p6=0.*ones(1,14);
  e, K  e3 ?, c: m8 k5 F% _ptrain=[pl;p2;p3;p4;p5;p6];
t=[100,80,60,40,20,0];
& _9 G8 L. t. g' \7 P

net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
net=init(net);
/ _8 z3 L$ q" b0 O2 i$ Y$ J( dnet.trainparam.epochs=100000;
net.trainparam.goal=1e-10;
[net,tr]=train(net,ptrain',t);
" A6 a/ i" A) m1 P1 I: O

" }  d3 s0 n+ B% w2 I5 [3 xfor i=1:3
a=ptes(:,i);
score(i)=sim(net,a); %#ok<SAGROW>
end



3 E# K5 z2 x: [' z) B

数学中国—罂粟

你是要问什么啊？

子任

大神能帮我运行一下代码吗
; T) J* \/ K# t* e我不懂这个
能留下联系方式吗
" b0 M0 ^5 W" S

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
我这里有一个小模板，你可以看看，希望对你有用：
clc,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 ...
6 A, X4 A; [* ~  c3 g, D    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
! c( u! L! P' U7 |6 p%人数 单位：万人
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...
. K& G/ P1 b$ {- o    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
%机动车数 单位：万辆
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...
( L6 j4 H7 ~  l, ]    0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];
# q8 I  y  f, X6 I/ z" }# j%公路面积 单位：万平方千米
glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    19490 20433 22598 25107 33442 36836 40548 42927 43462];
1 u; I* h1 X& H5 V. W%公路客运量 单位：万人
glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
' d9 T! h! G9 \, n/ M2 n, i* d1 S7 k    11115 13320 16762 18673 20724 20803 21804];
%公路货运量 单位：万吨
5 [# f% q/ _! w& x9 |- [p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
& [* q0 x/ g4 ~1 _1 ?& o: ?t = [glkyl;glhyl];          %目标矩阵
' T3 n9 V" g5 O( ?0 z% c- n
7 q8 r5 w% i$ S5 k. v  S& s3 r%%归一化处理
[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
0 c0 V$ f9 C7 B( z! R$ Fdx = [-1,1;-1 1;-1,1];                       %
3 V& f0 I" c$ Y" Y
0 d% k( C. K# V) M% [9 Q%BP网络训练
net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');

7 ~* B% e' M* C1 Z%%建立模型，并用梯度下降法训练
net.trainParam.show = 1000;          %轮回显示一次
net.trainParam.Lr = 0.05;            %学习速率
net.trainParam.epochs = 50000;       %最大训练轮回
& h7 @8 ~6 H$ W; ?- Rnet.trainParam.goal = 0.65*10^(-3);  %均方误差
3 Y* ]! [' v/ {8 f/ P0 }net = train(net,pn,tn);              %开始训练

: `$ M) O  d+ ^* U6 t4 U%%利用原始数据对BP网络仿真
an = sim(net,pn);                    %用训练好的模型进行仿真
a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
: p7 H9 }% M* D8 O1 {
%%对比测试
$ [4 N5 H1 {, o4 ^x = 1990:2009;
& O+ Q+ n6 ?# K5 a0 P0 Knewk = a(1,:);
# Y/ r6 g& U1 S/ p+ hnewh = a(2,:);
$ {" r* K6 y% K" H/ l2 F, wfigure(2);
5 e' F7 t+ \; ~# J& M( T+ n, Rsubplot(2,1,1);
plot(x,newk,'r-o',x,glkyl,'b--+');
% U) S& f( g$ c, Slegend('预测输出客运量','实际客运量');
xlabel('年份');
- \$ F  t+ ~( G6 x# ~ylabel('客运量/万人');
" T% M& C; H+ K1 ]2 ]' v8 Ptitle('运用数据箱客运量学习和测试对比图');
subplot(2,1,2);
7 G% P0 G2 o3 H4 Kplot(x,newh,'r-o',x,glhyl,'b--+');
( [0 b# Q4 G* F9 \0 ulegend('预测输出货运量','实际货运量');
xlabel('年份');
2 {1 p: E5 s+ W; m& ^( Mylabel('货运量/万吨');
6 C: L6 l; Y  [0 P+ L& xtitle('运用数据箱货运量学习和测试对比图');
4 K# e- A( h" B0 H5 a%利用训练好的网络进行预测
pnew = [73.39 75.55
    3.9635 4.0975
  y: h, j0 ~5 t; |    0.9880 1.0268];
pnewn = tramnmx(pnew,minp,maxp);
; x% c; h0 _& }%利用原始输入数据的归一化参数对新参数进行归一化
anewn = sim(net,pnewn);
$ x/ I- V, X( U$ q%利用归一化后的数据进行仿真
anew = postmnmx(anewn,mint,maxt)
%把仿真得到的数据还原原始的数量级

: N3 n+ l: A' H- w7 U1 ]& i
) Y  W: {, c/ |7 m+ h2 G; |$ i9 a6 r; K欢迎交流！

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
+ n$ m/ X: C! ^6 E; V/ V! g" ^# n非常感谢你
) k1 U$ T2 S- |8 A8 u' b
! f1 F" |- a5 [! v$ j% q2 a; F9 i