求 神经网络大神帮助

子任

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


" ?! o  B) [3 j1 ?6 R0 B% `a=max(p'); %#ok<UDIM>
for i=1:14
for j=1:3
ptest(i,j)=p(i,j)/a(i);
) p5 b& L9 m0 R; [/ Qend
end

- p4 P" w4 k7 f7 l  j
p1=ones(1,14);
p2=0.8.*ones(1,14);
p3=0.6.*ones(1,14);
8 T7 z8 R; Y5 F4 Y) ^; f7 T/ [p4=0.4.*ones(1,14);
p5=0.2.*ones(1,14);
: J% I, l' b% R+ m0 i$ w" Ap6=0.*ones(1,14);
! [& {" H( T7 Gptrain=[pl;p2;p3;p4;p5;p6];
4 L0 g9 S' v* ?3 M' @2 tt=[100,80,60,40,20,0];
  V2 y% L/ ^1 s4 l: _

net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
net=init(net);
net.trainparam.epochs=100000;
. }+ r6 N6 z6 W1 k; mnet.trainparam.goal=1e-10;
[net,tr]=train(net,ptrain',t);


for i=1:3
5 E0 Y  \+ @- @3 da=ptes(:,i);
score(i)=sim(net,a); %#ok<SAGROW>
end
+ }4 {! ?8 P" r, ?

0 U% t" D" n; i  [% _4 V5 K! i

数学中国—罂粟

你是要问什么啊？
& o5 y$ X1 Q) Q$ B7 |

子任

大神能帮我运行一下代码吗
$ z" d& h, ~) q7 A1 `& U我不懂这个
/ d& U0 z: j. I能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
我这里有一个小模板，你可以看看，希望对你有用：
clc,clear
* d0 t# R  F# @9 l7 Z1 xsqrs = [20.55 22.44 25.37 27.13 29.45 30.1 30.96 34.06 36.42 38.09 39.13 39.99 ...
$ t  N8 N) l: N: {. _" k    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...
    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
& V/ v- d" u5 R* X%机动车数 单位：万辆
- M/ Y0 t1 t6 K! t5 G! Qsqglmj = [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];
# ]7 A9 r7 t& X* B%公路面积 单位：万平方千米
- `2 G  P4 m+ G2 yglkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
7 Y' j) n# b( v, E. o. N2 f    19490 20433 22598 25107 33442 36836 40548 42927 43462];
* V- q& m8 ]/ s8 f2 X9 C% A%公路客运量 单位：万人
; H. _$ o$ Z4 Uglhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
* m! R  K, T% s3 Q    11115 13320 16762 18673 20724 20803 21804];
%公路货运量 单位：万吨
p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
t = [glkyl;glhyl];          %目标矩阵

: X0 e0 j! L! Z  d$ N) q%%归一化处理
7 b3 z. _! `% t6 a7 Y/ p5 v[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
! G) D+ G) |3 U; vdx = [-1,1;-1 1;-1,1];                       %

%BP网络训练
net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
3 J7 F/ r7 v+ w+ }
7 t& F! b6 E0 {4 H" O7 i%%建立模型，并用梯度下降法训练
0 h5 p. p2 ?. Q- p# Nnet.trainParam.show = 1000;          %轮回显示一次
6 p7 n1 r) G4 C, mnet.trainParam.Lr = 0.05;            %学习速率
net.trainParam.epochs = 50000;       %最大训练轮回
' q7 g! n4 G+ ]2 W0 O/ u3 ]3 vnet.trainParam.goal = 0.65*10^(-3);  %均方误差
net = train(net,pn,tn);              %开始训练
3 B" b' l, D6 M3 ]7 I
! [  G1 C1 L0 p( K%%利用原始数据对BP网络仿真
an = sim(net,pn);                    %用训练好的模型进行仿真
. x# `3 O: S" [& V0 ?, u; pa = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级

0 a7 q4 e1 e  [. b* _" I" z%%对比测试
x = 1990:2009;
newk = a(1,:);
newh = a(2,:);
figure(2);
subplot(2,1,1);
, t& k/ P" ^8 R! K! Pplot(x,newk,'r-o',x,glkyl,'b--+');
! G1 n$ e# u6 j9 @" d5 y- d+ Dlegend('预测输出客运量','实际客运量');
9 V& w$ [; t9 E, nxlabel('年份');
ylabel('客运量/万人');
title('运用数据箱客运量学习和测试对比图');
* a2 R# C* Q6 X( bsubplot(2,1,2);
plot(x,newh,'r-o',x,glhyl,'b--+');
legend('预测输出货运量','实际货运量');
xlabel('年份');
7 d# n+ Q  ^# d/ t- vylabel('货运量/万吨');
. e$ ~5 Y4 S$ M! btitle('运用数据箱货运量学习和测试对比图');
' B. b% C8 M+ h6 n, C4 S4 m%利用训练好的网络进行预测
pnew = [73.39 75.55
    3.9635 4.0975
    0.9880 1.0268];
pnewn = tramnmx(pnew,minp,maxp);
/ F: s7 F1 L, D/ M0 F. |0 f%利用原始输入数据的归一化参数对新参数进行归一化
. l& |4 f4 q+ j9 hanewn = sim(net,pnewn);
' V' B+ ]8 |* ^# }%利用归一化后的数据进行仿真
anew = postmnmx(anewn,mint,maxt)
& e! k$ o# A9 E3 E8 O%把仿真得到的数据还原原始的数量级
1 C# c0 r+ E' h% r6 j/ e6 z
" I! A% L! X. v: w
" J( d5 y+ E  Q3 U欢迎交流！

$ s4 W0 X, B1 n- j4 Y9 W

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
$ {1 e! v1 {# `# S& G3 F非常感谢你

. j6 [/ r- o5 W* P, g8 h* y