求 神经网络大神帮助

子任

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

; P$ f8 a  h( Q/ `' p# N! W3 w7 Z
a=max(p'); %#ok<UDIM>
/ u) n) J& A# ]for i=1:14
, G) Q& h) `/ o* c" nfor j=1:3
ptest(i,j)=p(i,j)/a(i);
$ X/ Z7 y/ w. F: s  O/ v! L0 ^6 mend
' R* f! [2 F3 Z$ N: A( s; iend

6 k: d7 m- G. D& a
4 t' s- o( S( }p1=ones(1,14);
4 T" N1 w: h- w& ^9 @. Hp2=0.8.*ones(1,14);
6 P: F: M. [* xp3=0.6.*ones(1,14);
p4=0.4.*ones(1,14);
- C+ }$ P  V# d1 ep5=0.2.*ones(1,14);
9 Q5 G+ U+ _/ M2 q, h- x! l% Hp6=0.*ones(1,14);
: B2 _8 ?) u% i; l) S; A; U- o- W- Mptrain=[pl;p2;p3;p4;p5;p6];
6 u+ L) F. ?' b7 u/ Bt=[100,80,60,40,20,0];

* T4 A' A: y4 c/ c0 v
net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
net=init(net);
net.trainparam.epochs=100000;
net.trainparam.goal=1e-10;
4 O* W1 W# c- g5 L% i0 r[net,tr]=train(net,ptrain',t);
# A: U# F& c3 [1 I7 [: F7 {

for i=1:3
& x* N* r. y7 K9 n- Pa=ptes(:,i);
score(i)=sim(net,a); %#ok<SAGROW>
end

" s" K0 w2 d% |+ ]- I9 T
* P$ U8 k9 q# B5 L0 V4 s4 A3 F

数学中国—罂粟

你是要问什么啊？

子任

大神能帮我运行一下代码吗
! f  D: o7 D  I我不懂这个
能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
我这里有一个小模板，你可以看看，希望对你有用：
. G% b5 k/ @1 z4 F( S  bclc,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 ...
    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
, b2 c8 ^' k. \$ u- E& E! F%人数 单位：万人
7 ^7 ]# O9 \0 \5 ^7 E2 f6 gsqjdcs = [0.6 0.75 0.85 0.90 1.05 1.35 1.45 1.60 1.70 1.85 2.15 2.20...
% u7 q' t' I% X* [8 i0 q    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
+ E) Q% l. |: r%机动车数 单位：万辆
8 m4 L/ R4 x% J2 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...
) H0 ]( f2 N$ m* ~8 u    0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];
1 Q5 C! n7 G5 g3 x; E& r# y2 `%公路面积 单位：万平方千米
glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    19490 20433 22598 25107 33442 36836 40548 42927 43462];
9 z" L) H8 T; i* m%公路客运量 单位：万人
glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
# c! s8 s$ F; w8 @  d0 ^2 ?    11115 13320 16762 18673 20724 20803 21804];
# S# v0 R! O4 D0 p%公路货运量 单位：万吨
3 z: S1 w: Q. G( \. r& Hp = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
t = [glkyl;glhyl];          %目标矩阵
. Q! T$ n0 p' i! n8 ~# G
  }; B) b% t  B' n%%归一化处理
- A! N& N) g2 d! ^5 M[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
6 y4 c  k4 N/ O; r1 cdx = [-1,1;-1 1;-1,1];                       %
. `0 i1 g1 K, g- W7 V& U
%BP网络训练
/ c$ r& Z& Z& X8 c( Wnet = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
* o; H' U1 c- e& `5 F
3 J. K! s. D4 n7 w, e) c* {, X+ `%%建立模型，并用梯度下降法训练
net.trainParam.show = 1000;          %轮回显示一次
net.trainParam.Lr = 0.05;            %学习速率
, D( F  A+ d& O8 s1 z  U/ N9 knet.trainParam.epochs = 50000;       %最大训练轮回
net.trainParam.goal = 0.65*10^(-3);  %均方误差
net = train(net,pn,tn);              %开始训练
1 G$ G) }; k- d% L" t: `/ n- E
, A2 _1 i! z6 V/ c! @1 ]+ u%%利用原始数据对BP网络仿真
an = sim(net,pn);                    %用训练好的模型进行仿真
2 |  B- O/ r# p# l# Oa = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级

& k0 M3 n1 A; V4 R& E# u8 ^%%对比测试
x = 1990:2009;
newk = a(1,:);
4 k6 f2 ^  g( T4 `newh = a(2,:);
figure(2);
, g2 }9 I  f% B: j: ~# P+ e: ]subplot(2,1,1);
( `& B; X  f# L+ B( X" l! b! c. _plot(x,newk,'r-o',x,glkyl,'b--+');
legend('预测输出客运量','实际客运量');
xlabel('年份');
7 l6 Q+ B( Q3 w, T: B5 C+ sylabel('客运量/万人');
) I1 [; }! U! ]+ G5 G# utitle('运用数据箱客运量学习和测试对比图');
subplot(2,1,2);
plot(x,newh,'r-o',x,glhyl,'b--+');
* e) K5 J: f" M8 qlegend('预测输出货运量','实际货运量');
1 T+ G8 j, R1 `- q' f9 L  _8 xxlabel('年份');
0 f8 }+ ^& |9 ^) o# g: rylabel('货运量/万吨');
title('运用数据箱货运量学习和测试对比图');
2 _( x8 T- z$ s/ M( F$ u; X2 `%利用训练好的网络进行预测
) C: O# D0 f& c8 npnew = [73.39 75.55
    3.9635 4.0975
    0.9880 1.0268];
pnewn = tramnmx(pnew,minp,maxp);
%利用原始输入数据的归一化参数对新参数进行归一化
anewn = sim(net,pnewn);
  S0 Q5 \* E( S1 p# _; h%利用归一化后的数据进行仿真
& C2 I) x6 b; O4 h9 ?anew = postmnmx(anewn,mint,maxt)
6 t  T6 W3 J) y! q3 m# C/ J%把仿真得到的数据还原原始的数量级


" i8 I- K" N6 ?! i3 }7 T欢迎交流！

) u# h  k/ h) A$ s

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
非常感谢你

6 C, l+ r& |1 W$ F& H7 Z