求 神经网络大神帮助

子任

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
! G  X* g) t( l+ efor j=1:3
ptest(i,j)=p(i,j)/a(i);
end
end
' q. Q" t" g" \9 q3 j7 Y1 P
4 r4 m0 x6 g" A0 S- U) @  \
0 w# p! u, q- u. xp1=ones(1,14);
p2=0.8.*ones(1,14);
# g' J. g2 r  ~: f: h+ ?, M" d, Np3=0.6.*ones(1,14);
p4=0.4.*ones(1,14);
p5=0.2.*ones(1,14);
p6=0.*ones(1,14);
ptrain=[pl;p2;p3;p4;p5;p6];
t=[100,80,60,40,20,0];


: [, U+ J; z' P- H# a. |net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
$ d7 s! v; W/ vnet=init(net);
net.trainparam.epochs=100000;
# }' n  a: d& Z4 m; c* _5 wnet.trainparam.goal=1e-10;
[net,tr]=train(net,ptrain',t);


, M. R. _7 p9 O, M( ?$ V- bfor i=1:3
a=ptes(:,i);
& ~7 g7 S4 u7 S1 v# }score(i)=sim(net,a); %#ok<SAGROW>
" x. {& `& `; Eend
% V% v" g6 @3 T5 P
0 Q* _' D, H9 {* f& G& Z7 B. h: ^1 ]

( f( r8 x$ _- u9 H

数学中国—罂粟

你是要问什么啊？
. z1 H. G4 B: e4 W

子任

大神能帮我运行一下代码吗
! y- i  J9 S5 `% e我不懂这个
. P3 C3 Q  _7 R) W: `, R& C能留下联系方式吗
2 `0 c  Q6 J- o% ?/ N7 V' x

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
我这里有一个小模板，你可以看看，希望对你有用：
  s4 h( |; i; L6 J7 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];
* {! r, b3 L; C6 U%人数 单位：万人
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...
% ~+ `" c6 E% _+ D* x    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
/ w" ~  j' G& v9 l, `! e$ i%机动车数 单位：万辆
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];
  Q# w5 q3 p) ^$ U%公路面积 单位：万平方千米
glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
. a+ Y( t+ C/ u& q6 o$ T+ g7 ]/ W    19490 20433 22598 25107 33442 36836 40548 42927 43462];
8 D% l4 `5 z8 h6 {  y7 W( G: G%公路客运量 单位：万人
glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
" \" a" e6 E2 m% B/ n) \( W8 X    11115 13320 16762 18673 20724 20803 21804];
/ s4 ~  P* ~, F) @8 `; G%公路货运量 单位：万吨
p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
t = [glkyl;glhyl];          %目标矩阵

%%归一化处理
[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
$ M: V, e7 G/ W0 S' Q" z& ldx = [-1,1;-1 1;-1,1];                       %
- n2 T5 n/ H0 @$ J! S7 L! Z. b) {. t
%BP网络训练
net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');

%%建立模型，并用梯度下降法训练
5 c2 S; r, q. ?9 k! i; jnet.trainParam.show = 1000;          %轮回显示一次
net.trainParam.Lr = 0.05;            %学习速率
net.trainParam.epochs = 50000;       %最大训练轮回
net.trainParam.goal = 0.65*10^(-3);  %均方误差
% e8 b7 b. [% A, |0 i" \net = train(net,pn,tn);              %开始训练
/ z8 \# @# M# O& a. c2 l1 A
% p3 i% ], ^+ Z8 S%%利用原始数据对BP网络仿真
" z0 H/ I0 d3 _& W) O# M% `an = sim(net,pn);                    %用训练好的模型进行仿真
a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级

%%对比测试
/ ?, [4 u/ @* M; j: X' _- F) r- I6 ox = 1990:2009;
newk = a(1,:);
newh = a(2,:);
figure(2);
2 E! z8 O$ |0 R, Q1 i- q# W$ Asubplot(2,1,1);
# n& R4 |  {* M9 ~& X' yplot(x,newk,'r-o',x,glkyl,'b--+');
legend('预测输出客运量','实际客运量');
' [4 H9 E8 U9 _/ f1 v& R  Q' Xxlabel('年份');
% O9 @3 ?+ h% u& ~" D/ \+ Sylabel('客运量/万人');
title('运用数据箱客运量学习和测试对比图');
6 r- K4 B4 M& H+ Q& Osubplot(2,1,2);
  F- L8 q( K  r: v0 D0 jplot(x,newh,'r-o',x,glhyl,'b--+');
legend('预测输出货运量','实际货运量');
xlabel('年份');
3 `3 o+ w0 y, d7 ~8 ?! W$ ]% d0 [: {1 yylabel('货运量/万吨');
# O5 Y. A4 x( P( dtitle('运用数据箱货运量学习和测试对比图');
. j: [8 I  V: y- Q4 |%利用训练好的网络进行预测
$ d" S" o, O! |6 ^# Ppnew = [73.39 75.55
    3.9635 4.0975
4 K. c, O3 H3 T. T( R$ m, W    0.9880 1.0268];
2 b. o' Y9 ?( P5 G( D/ p+ m  Upnewn = tramnmx(pnew,minp,maxp);
%利用原始输入数据的归一化参数对新参数进行归一化
8 L$ x( w+ _5 x' wanewn = sim(net,pnewn);
  ?  F% F$ T; J' E2 N# q" s%利用归一化后的数据进行仿真
% m% i, a* A! H4 sanew = postmnmx(anewn,mint,maxt)
%把仿真得到的数据还原原始的数量级

" ^* Y/ R" d- |- n/ q2 o5 [
欢迎交流！

7 i, I9 E) r8 ~. M+ F2 r: h3 |  `

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
# m. t5 X+ z% N, |非常感谢你
0 b- ~$ l0 |* Y& O& a7 {