求 神经网络大神帮助

子任

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];
+ x# Y- X0 E3 J0 O

, _" b2 C. t7 K7 Sa=max(p'); %#ok<UDIM>
4 j2 v+ P# e  {+ ?  Y* \4 Afor i=1:14
for j=1:3
4 C% X0 O7 @7 Q# q0 q( T4 j1 [ptest(i,j)=p(i,j)/a(i);
1 @4 f! c" c0 g& Z0 I. {( ?end
end
. t4 _! M6 b, I9 I# {0 F* h
% I! e1 c6 A4 j7 s
' ~  C% G$ ^) r" W  \% }p1=ones(1,14);
9 M, V: M; O+ L+ T4 W  [p2=0.8.*ones(1,14);
4 b6 n8 _; g% ]4 |( A) z* Np3=0.6.*ones(1,14);
p4=0.4.*ones(1,14);
p5=0.2.*ones(1,14);
; c- n! F3 D# q" }8 ?p6=0.*ones(1,14);
7 D1 [" `( K3 u* [) E2 xptrain=[pl;p2;p3;p4;p5;p6];
t=[100,80,60,40,20,0];
: `2 m- Z6 A, G$ c' [4 C
! w5 w) t3 d9 i; t) ]1 @7 ^
+ @% ], q6 h$ E, q9 c7 O& Ynet=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
+ G0 D& }  m, m/ F+ Z; c4 I. ~net=init(net);
net.trainparam.epochs=100000;
net.trainparam.goal=1e-10;
3 V/ V/ h( W* R$ W[net,tr]=train(net,ptrain',t);
: T0 \, ]6 x: u4 i( f- A' ]6 U8 z5 h
9 j# \  I0 u+ g7 p
% {  V5 I; [, C* }+ yfor i=1:3
6 f" q- B/ U) \- X* I7 oa=ptes(:,i);
score(i)=sim(net,a); %#ok<SAGROW>
end

数学中国—罂粟

你是要问什么啊？

子任

大神能帮我运行一下代码吗
8 T- n# j3 r1 v5 |+ u我不懂这个
, u8 n3 }  l9 P/ t% s* x能留下联系方式吗
1 ^$ M% _) O8 S! q

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
我这里有一个小模板，你可以看看，希望对你有用：
) E+ n. |) @6 V  Hclc,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 ...
* A  M2 [5 L& {" @* F5 L    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
/ c* h8 g% K! t% q( l0 Q& f%人数 单位：万人
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];
%机动车数 单位：万辆
2 N( _& z$ h+ V* Fsqglmj = [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];
%公路面积 单位：万平方千米
* v0 ^4 u% e& L; @7 C% Pglkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    19490 20433 22598 25107 33442 36836 40548 42927 43462];
%公路客运量 单位：万人
2 }* C) R% a/ P1 pglhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
) _* }8 T# e1 ~& _7 j3 O8 R    11115 13320 16762 18673 20724 20803 21804];
%公路货运量 单位：万吨
p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
t = [glkyl;glhyl];          %目标矩阵

%%归一化处理
1 g3 m+ _% [9 e- X[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
dx = [-1,1;-1 1;-1,1];                       %

9 U- F" u" W+ E%BP网络训练
1 I3 h7 I1 _4 knet = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
! H; I+ s/ q. R( `, F
5 L2 `  M$ B1 m7 ]9 C8 S%%建立模型，并用梯度下降法训练
net.trainParam.show = 1000;          %轮回显示一次
net.trainParam.Lr = 0.05;            %学习速率
net.trainParam.epochs = 50000;       %最大训练轮回
1 M: g# x4 ?& ~2 \7 ?net.trainParam.goal = 0.65*10^(-3);  %均方误差
: T  C$ c  w. A( s" F7 v% F* k4 anet = train(net,pn,tn);              %开始训练

%%利用原始数据对BP网络仿真
9 {  {, w3 n( ]3 }2 Q( ran = sim(net,pn);                    %用训练好的模型进行仿真
a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
) ^" V5 R7 H/ r( q2 Z
%%对比测试
x = 1990:2009;
, o9 U( a9 B- y" t& R- mnewk = a(1,:);
newh = a(2,:);
0 w6 D* I( ?+ F: g2 w! X9 c- ufigure(2);
) F  G6 H! f- w2 G# f, Qsubplot(2,1,1);
9 [) C- J' [0 d" rplot(x,newk,'r-o',x,glkyl,'b--+');
legend('预测输出客运量','实际客运量');
xlabel('年份');
6 E- l& J! y2 r" }& {) L: Jylabel('客运量/万人');
title('运用数据箱客运量学习和测试对比图');
& {( g( n) x# J4 e! ?$ Osubplot(2,1,2);
% Q/ _3 L" U8 Oplot(x,newh,'r-o',x,glhyl,'b--+');
- l( B; G5 w4 Ilegend('预测输出货运量','实际货运量');
9 Y6 X0 E$ B' L9 [: h# Mxlabel('年份');
5 a" N2 X7 {1 h' Y* nylabel('货运量/万吨');
' ^" b+ N) y7 S( K4 Ltitle('运用数据箱货运量学习和测试对比图');
/ f/ e- ^2 @  P, F$ y& K4 P%利用训练好的网络进行预测
pnew = [73.39 75.55
    3.9635 4.0975
, L( T8 n2 f7 e    0.9880 1.0268];
# P5 M$ J9 ~+ K# Y& zpnewn = tramnmx(pnew,minp,maxp);
%利用原始输入数据的归一化参数对新参数进行归一化
anewn = sim(net,pnewn);
%利用归一化后的数据进行仿真
$ H" ~$ q$ \+ Janew = postmnmx(anewn,mint,maxt)
%把仿真得到的数据还原原始的数量级
' G* ~" w* t. z
8 j; c. C/ f  E9 F6 |+ Z, i
欢迎交流！

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
4 C( \& L  ?! Y) I非常感谢你