求 神经网络大神帮助

子任

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];
) ]9 Q& K+ f4 u, f4 y

a=max(p'); %#ok<UDIM>
for i=1:14
. b$ k+ D: c, h1 z' K2 d' h* `; y! R4 ofor j=1:3
ptest(i,j)=p(i,j)/a(i);
7 g2 a% m, l# [8 j: a- Zend
2 i2 y7 M) S% Q0 n# mend

. x' h0 L  K" B! S
p1=ones(1,14);
p2=0.8.*ones(1,14);
p3=0.6.*ones(1,14);
( |' P% `5 ~" _2 N* kp4=0.4.*ones(1,14);
p5=0.2.*ones(1,14);
p6=0.*ones(1,14);
ptrain=[pl;p2;p3;p4;p5;p6];
$ i7 z2 ^6 {0 u% }; k1 }6 Kt=[100,80,60,40,20,0];
' R, X8 [6 ?" K% U5 Z; V; K

. r3 m& Q" X. n! H) u: r& Lnet=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
* ?2 S- B: {% j) {% w2 i1 F" Gnet=init(net);
net.trainparam.epochs=100000;
net.trainparam.goal=1e-10;
[net,tr]=train(net,ptrain',t);
' B: x$ B3 Q" u6 t% g
( ]6 e& t  }/ S/ e5 w$ Q
for i=1:3
a=ptes(:,i);
score(i)=sim(net,a); %#ok<SAGROW>
end
( S. q& r: `2 i7 ]: X

/ k3 C% {0 n- v& i# f, `
$ [/ @% l; I6 j* G

数学中国—罂粟

你是要问什么啊？
- u3 q2 B! \) C: z; U2 h

子任

大神能帮我运行一下代码吗
# }" _8 u9 y( L0 X我不懂这个
; _$ z! Y& |+ U& G2 A+ e" p能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
我这里有一个小模板，你可以看看，希望对你有用：
clc,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 ...
0 w/ A; ~2 C% S# I    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
2 [) s8 Y$ D0 g; c% W; Y3 Z%人数 单位：万人
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];
* F9 ^2 C8 p6 `% n%机动车数 单位：万辆
" K' I" E. ]6 \- x! A3 k! tsqglmj = [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];
3 k% E  Y4 K/ Q/ {%公路面积 单位：万平方千米
# l( Z: x  X+ I# C  oglkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
. q9 W: w3 `! s9 A) s) G, \8 O    19490 20433 22598 25107 33442 36836 40548 42927 43462];
$ v& N6 U3 k4 ]- F%公路客运量 单位：万人
* Y4 r8 B8 Y5 wglhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    11115 13320 16762 18673 20724 20803 21804];
4 T1 `+ @) v" v, E3 o3 C/ Z+ }" k%公路货运量 单位：万吨
p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
' r6 y- i" r: E: @. c. it = [glkyl;glhyl];          %目标矩阵

%%归一化处理
7 v3 G- Y& W3 P) @1 M" p7 f[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
, A: m  a5 M/ p/ }. S; fdx = [-1,1;-1 1;-1,1];                       %

%BP网络训练
4 A# U& h/ e3 e7 D& Z' J# f+ Knet = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');

%%建立模型，并用梯度下降法训练
net.trainParam.show = 1000;          %轮回显示一次
6 W2 K) a! j5 o3 ^net.trainParam.Lr = 0.05;            %学习速率
0 L) u* `! m3 I  hnet.trainParam.epochs = 50000;       %最大训练轮回
net.trainParam.goal = 0.65*10^(-3);  %均方误差
net = train(net,pn,tn);              %开始训练
- O: a- B0 |! `4 @) H
: H& U9 y1 Q0 [. J%%利用原始数据对BP网络仿真
an = sim(net,pn);                    %用训练好的模型进行仿真
a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级

1 n& y. k3 T, s9 I* C5 |) C/ l%%对比测试
) i, ]# g+ {/ k- C" l/ V/ jx = 1990:2009;
, X& R" e' w" O( H2 S. P) x" q# ?newk = a(1,:);
' a. q4 a, V* enewh = a(2,:);
figure(2);
; i' c) L+ V' ?/ V4 H. csubplot(2,1,1);
; N+ i" F% p. Y- wplot(x,newk,'r-o',x,glkyl,'b--+');
; r2 }* l& S; W! y3 |( U" qlegend('预测输出客运量','实际客运量');
5 i/ u3 e7 A( N4 A3 lxlabel('年份');
+ w9 I/ m( N6 X8 eylabel('客运量/万人');
title('运用数据箱客运量学习和测试对比图');
4 J7 f0 A3 \" z3 tsubplot(2,1,2);
plot(x,newh,'r-o',x,glhyl,'b--+');
- s8 I5 j5 N; Z/ o- M: E& qlegend('预测输出货运量','实际货运量');
xlabel('年份');
ylabel('货运量/万吨');
/ `; J. J% S* ]* stitle('运用数据箱货运量学习和测试对比图');
2 U2 z; V5 Q. f%利用训练好的网络进行预测
8 F1 q: _. b6 J. r& X7 ipnew = [73.39 75.55
" `: C9 V+ n2 g$ e) k- X4 Y    3.9635 4.0975
' x; _, q4 `/ [/ p: {6 n( J% j    0.9880 1.0268];
pnewn = tramnmx(pnew,minp,maxp);
/ r7 x  s4 L3 C9 j0 U" b1 k9 V3 ^%利用原始输入数据的归一化参数对新参数进行归一化
anewn = sim(net,pnewn);
%利用归一化后的数据进行仿真
anew = postmnmx(anewn,mint,maxt)
3 e8 S* ^! E! q+ Z' G/ |%把仿真得到的数据还原原始的数量级
( q; _1 q: F3 P/ F$ a& c! V

欢迎交流！

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
- \3 m+ p  o5 y0 }, o0 F+ u3 t. v非常感谢你
" P+ w$ l& W2 v# }3 j" E