求 神经网络大神帮助

子任

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];
0 [0 y. W) ?9 E8 q7 k5 s' i" N
+ v0 |8 `, m4 S* I
! Y! j, B: K! p+ ^3 d& K, Ka=max(p'); %#ok<UDIM>
for i=1:14
% s9 t$ \3 \) f( H' Z2 i* f  W/ Zfor j=1:3
$ e4 f$ I8 h: O2 X2 @ptest(i,j)=p(i,j)/a(i);
end
- C! ]* R8 [4 I1 A  }end

% v( r+ x, K% N! I% _7 \
p1=ones(1,14);
7 g# m: ?2 f9 u* zp2=0.8.*ones(1,14);
p3=0.6.*ones(1,14);
! G5 C1 i0 J3 B  \: U. _p4=0.4.*ones(1,14);
p5=0.2.*ones(1,14);
: u) m$ F  ^, Dp6=0.*ones(1,14);
ptrain=[pl;p2;p3;p4;p5;p6];
9 T+ j9 T' B2 V! rt=[100,80,60,40,20,0];
" w/ K2 f2 L8 |9 w# m

net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
& i5 G+ w1 J( A( x' jnet=init(net);
net.trainparam.epochs=100000;
& d+ n7 k2 v; `! s6 Ynet.trainparam.goal=1e-10;
[net,tr]=train(net,ptrain',t);


for i=1:3
a=ptes(:,i);
score(i)=sim(net,a); %#ok<SAGROW>
# L& n$ ?& y* m1 Jend
( \' b0 b; T0 ?. p, b# H

% L6 c, h# M; j8 P9 v' T3 O* j

数学中国—罂粟

你是要问什么啊？

子任

大神能帮我运行一下代码吗
' ^& |# B* a, Z3 C1 d+ q1 A我不懂这个
能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
. S- h) N! J( K1 c) E8 ~我这里有一个小模板，你可以看看，希望对你有用：
! ^' g9 ], Z+ N  U, d' pclc,clear
1 @' I$ V/ y, ?* M1 g& T/ K* c6 osqrs = [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];
: ~! \# j7 i7 i' D7 b%人数 单位：万人
) [9 ~* F" {# y/ i+ {5 s1 Lsqjdcs = [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];
& \( U8 y. A1 x) x0 `8 P1 [2 Z8 `%机动车数 单位：万辆
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];
%公路面积 单位：万平方千米
% A7 N6 n' P& H3 Vglkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    19490 20433 22598 25107 33442 36836 40548 42927 43462];
  _7 A* l$ @* @& A" X%公路客运量 单位：万人
2 d0 c7 U6 N+ f! f! Q1 \glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    11115 13320 16762 18673 20724 20803 21804];
%公路货运量 单位：万吨
6 W) G: |% R1 n6 k5 C5 ]7 sp = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
+ e5 H  [4 w" y/ |t = [glkyl;glhyl];          %目标矩阵
5 P' ]9 ^6 K4 Y- Q8 K! P9 Q" Y
%%归一化处理
# {; K  @  W+ ]( A+ I2 Z( M[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
dx = [-1,1;-1 1;-1,1];                       %
: q8 X3 y* a+ t  X, U" Y- v6 i3 T
%BP网络训练
& W* r; H' F4 E  [net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
: [) a: z' V7 I& r# d) h6 [  \
%%建立模型，并用梯度下降法训练
+ v3 @4 H" H; p( jnet.trainParam.show = 1000;          %轮回显示一次
net.trainParam.Lr = 0.05;            %学习速率
2 O4 E$ Z% O$ R. {5 V4 z2 K6 ]; v4 _net.trainParam.epochs = 50000;       %最大训练轮回
net.trainParam.goal = 0.65*10^(-3);  %均方误差
net = train(net,pn,tn);              %开始训练

%%利用原始数据对BP网络仿真
an = sim(net,pn);                    %用训练好的模型进行仿真
, r/ q2 S0 o: ~" Q% G. Ia = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
, p4 R; m. _" I# |; a1 w7 P3 {
%%对比测试
$ z; f7 L! a9 v$ Kx = 1990:2009;
7 }4 j# d2 m; h, enewk = a(1,:);
newh = a(2,:);
figure(2);
subplot(2,1,1);
plot(x,newk,'r-o',x,glkyl,'b--+');
. Y% s2 S" \5 d+ Z, llegend('预测输出客运量','实际客运量');
# g6 a, r+ y2 `+ X: ?+ N  Wxlabel('年份');
& S3 z  g. M) w, j) F: eylabel('客运量/万人');
title('运用数据箱客运量学习和测试对比图');
# K2 p- Z9 }1 q5 E9 c: bsubplot(2,1,2);
plot(x,newh,'r-o',x,glhyl,'b--+');
6 }/ [) N2 y0 y+ T9 {$ slegend('预测输出货运量','实际货运量');
% A( ~: E& o4 K! u. Qxlabel('年份');
ylabel('货运量/万吨');
title('运用数据箱货运量学习和测试对比图');
%利用训练好的网络进行预测
$ u! N) x9 S. u( z+ Opnew = [73.39 75.55
% c# \: A! c+ s    3.9635 4.0975
/ X. ^3 O' W9 \8 m8 B0 c( ~8 J    0.9880 1.0268];
$ c( _) e' d% F# mpnewn = tramnmx(pnew,minp,maxp);
%利用原始输入数据的归一化参数对新参数进行归一化
anewn = sim(net,pnewn);
%利用归一化后的数据进行仿真
anew = postmnmx(anewn,mint,maxt)
* o0 h; E- ]3 W( g3 z4 r0 P. v%把仿真得到的数据还原原始的数量级
! w, w) M* R2 P+ z9 c& ~

7 j& {7 i* i+ H6 P欢迎交流！
  c% ?+ v) r7 |. ~, ]

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
非常感谢你
( x7 r$ ~! a5 e0 v- ?