求 神经网络大神帮助

子任

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


a=max(p'); %#ok<UDIM>
for i=1:14
for j=1:3
ptest(i,j)=p(i,j)/a(i);
) \% P9 e) z0 Z8 z* gend
9 y/ c4 ~' j! j; Mend


/ o. w( M# Y2 g& H# jp1=ones(1,14);
p2=0.8.*ones(1,14);
p3=0.6.*ones(1,14);
1 m9 p2 l3 n, o! B' j8 fp4=0.4.*ones(1,14);
) B  P/ n/ I# ?: Pp5=0.2.*ones(1,14);
% U% `% f' S' n1 T/ V* L1 l5 `3 Lp6=0.*ones(1,14);
ptrain=[pl;p2;p3;p4;p5;p6];
t=[100,80,60,40,20,0];
% A+ l3 L) l" X( a/ Z3 z$ L! r6 \
4 N* A  g: d4 a! m, R
* p) F7 w# F* u/ l% s9 a9 [net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
net=init(net);
net.trainparam.epochs=100000;
) i! J, @2 R' G% V9 I/ Unet.trainparam.goal=1e-10;
[net,tr]=train(net,ptrain',t);
" d% h( n0 I# w+ T" f
& Q! d: C3 ^% W1 H
/ A9 S& p' s# c0 Q- `+ Z, {; y% S4 afor i=1:3
a=ptes(:,i);
score(i)=sim(net,a); %#ok<SAGROW>
& R' R  `7 n3 l& t9 ~end


' ]% H8 r9 c$ W& S- `4 ^
/ d9 t4 Q6 L% k1 F8 ~) }

数学中国—罂粟

你是要问什么啊？
( }  h- v- x0 a0 X! H

子任

大神能帮我运行一下代码吗
我不懂这个
能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
我这里有一个小模板，你可以看看，希望对你有用：
5 P4 P; f! Y- b1 uclc,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 ...
& z7 s! W; y8 a3 ?' w5 H% q  g    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
%人数 单位：万人
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...
. c3 d, s: X+ ?) s7 P1 E( N$ p. R    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
%机动车数 单位：万辆
  A' H  |2 k% o7 y  }, Osqglmj = [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];
4 t0 b" \1 b: i8 f: k7 \%公路面积 单位：万平方千米
. U  `1 h7 C& c" P# rglkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
5 P, \. w6 V0 Q8 O    19490 20433 22598 25107 33442 36836 40548 42927 43462];
%公路客运量 单位：万人
4 H' ^8 \3 t' |% V& \7 M. r6 lglhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
" x* {1 v" _* R1 {0 a9 C+ H: i' @    11115 13320 16762 18673 20724 20803 21804];
3 `3 j5 S- L# B0 F%公路货运量 单位：万吨
( [, N) B6 N, M. l6 Op = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
/ O5 j/ T6 h: jt = [glkyl;glhyl];          %目标矩阵

%%归一化处理
3 ~3 k4 f0 f" f: `" m: Y9 Q: V[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
dx = [-1,1;-1 1;-1,1];                       %
7 ^- ~2 C4 |0 j' a7 B# ]% Z
& F& P1 Q; Y* K( H/ r% t- y%BP网络训练
net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
* d6 U% r) P& b5 o! {
%%建立模型，并用梯度下降法训练
net.trainParam.show = 1000;          %轮回显示一次
net.trainParam.Lr = 0.05;            %学习速率
net.trainParam.epochs = 50000;       %最大训练轮回
  h, ~: [- c: s0 @+ }0 l  pnet.trainParam.goal = 0.65*10^(-3);  %均方误差
net = train(net,pn,tn);              %开始训练

* X! T/ T" Y% S: X; K% @%%利用原始数据对BP网络仿真
an = sim(net,pn);                    %用训练好的模型进行仿真
a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级

; _; j9 C; h4 [- L%%对比测试
x = 1990:2009;
! m1 ~1 c, x# N: `' bnewk = a(1,:);
" J; w# ~5 u1 e3 E! znewh = a(2,:);
% I4 t: A) j) T* m8 F& Q0 xfigure(2);
subplot(2,1,1);
' j# l- K, l  o6 y& W; wplot(x,newk,'r-o',x,glkyl,'b--+');
; w' D; K: ~' ]3 t) Q5 D2 N# s, Xlegend('预测输出客运量','实际客运量');
xlabel('年份');
ylabel('客运量/万人');
9 F7 |2 D2 o# q% _title('运用数据箱客运量学习和测试对比图');
! `$ _4 ^5 z, F; vsubplot(2,1,2);
plot(x,newh,'r-o',x,glhyl,'b--+');
, i( r7 w% D6 nlegend('预测输出货运量','实际货运量');
xlabel('年份');
/ F) u1 ]4 H: t. }ylabel('货运量/万吨');
title('运用数据箱货运量学习和测试对比图');
! M/ e8 }) v" N0 Y7 w5 I0 ~; G%利用训练好的网络进行预测
  E4 W: D4 ], j! N; W1 P6 ?' qpnew = [73.39 75.55
    3.9635 4.0975
' h- K( Z5 K/ w4 w1 y/ r    0.9880 1.0268];
pnewn = tramnmx(pnew,minp,maxp);
%利用原始输入数据的归一化参数对新参数进行归一化
anewn = sim(net,pnewn);
7 z& v- y5 }' g& @5 J6 J%利用归一化后的数据进行仿真
4 C  `! a9 S6 }/ t+ j* X3 M3 F* Banew = postmnmx(anewn,mint,maxt)
+ T0 j: V- w1 n0 W3 f* v%把仿真得到的数据还原原始的数量级


欢迎交流！
( ^5 \/ k7 R  \9 R5 ]8 q

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
/ T& X9 R: x' g# x  H& J' K非常感谢你
& n- Y8 v  `; |+ d- o% d0 M
: f. O8 \6 S8 G8 m# T3 }% t