求 神经网络大神帮助

子任

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


: x! T/ k8 s( k! D; aa=max(p'); %#ok<UDIM>
/ o" M# \3 o2 ]( H8 Q1 P6 M# c, vfor i=1:14
$ [" {2 j$ g" O- H! v9 w% k6 bfor j=1:3
4 y/ B' b+ B( o! i: X1 H! m# Tptest(i,j)=p(i,j)/a(i);
end
end
8 E6 w/ {1 ?3 a9 {5 h6 e$ R. Y

p1=ones(1,14);
p2=0.8.*ones(1,14);
3 C- p/ G& r2 B1 R$ O1 u  ep3=0.6.*ones(1,14);
( W$ B: f. y. W+ wp4=0.4.*ones(1,14);
p5=0.2.*ones(1,14);
p6=0.*ones(1,14);
ptrain=[pl;p2;p3;p4;p5;p6];
. C% o0 i, {* V, wt=[100,80,60,40,20,0];
6 k: s5 ~- j+ A2 k: G
2 O" Q- ?8 Z8 }  n& C8 @
net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
net=init(net);
' Q+ j8 v, `: v9 e* p6 p/ Nnet.trainparam.epochs=100000;
. k; z# ~- e# u9 f3 qnet.trainparam.goal=1e-10;
8 x3 B, \+ }/ U$ @[net,tr]=train(net,ptrain',t);
% _* Z& c, |) E0 q+ Q
  g2 u/ Q* S5 Z. w# y/ T
/ j. J9 i* C! P; U) yfor i=1:3
a=ptes(:,i);
, ?8 F( f6 Y; r1 ^) D2 `6 \7 jscore(i)=sim(net,a); %#ok<SAGROW>
end


% @. Y3 u1 B$ H/ E9 C! P

数学中国—罂粟

你是要问什么啊？

子任

大神能帮我运行一下代码吗
. J9 U7 {9 [4 G4 A我不懂这个
, H8 y% P1 W7 i# U! z3 n能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
我这里有一个小模板，你可以看看，希望对你有用：
2 e! [7 u+ w; Gclc,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 ...
3 k: q9 c1 @9 D' m, \" d1 Z    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
%人数 单位：万人
, u9 [  b; U' C5 X, @- Bsqjdcs = [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];
%机动车数 单位：万辆
* g% @* {$ j; X" Z) I: w+ `% ~$ Rsqglmj = [0.09 0.11 0.11 0.14 0.20 0.23 0.23 0.32 0.32 0.34 0.36 0.36...
) U( H) w* E, ?6 a4 P, v    0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];
7 c7 ^- H  H1 g, g2 v- N%公路面积 单位：万平方千米
glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
' E* o- T& g& }' j5 Y8 ]    19490 20433 22598 25107 33442 36836 40548 42927 43462];
%公路客运量 单位：万人
glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    11115 13320 16762 18673 20724 20803 21804];
%公路货运量 单位：万吨
p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
t = [glkyl;glhyl];          %目标矩阵

%%归一化处理
[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
- d5 I: d8 a4 V- mdx = [-1,1;-1 1;-1,1];                       %

%BP网络训练
net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
2 O1 O+ |5 p  ]6 G4 O  l
- z* H% U& }8 ^% q%%建立模型，并用梯度下降法训练
1 p4 c' }8 F5 tnet.trainParam.show = 1000;          %轮回显示一次
net.trainParam.Lr = 0.05;            %学习速率
' a2 i7 J9 O% C  Q2 M* Dnet.trainParam.epochs = 50000;       %最大训练轮回
+ h* v! V- U0 ?' T! knet.trainParam.goal = 0.65*10^(-3);  %均方误差
net = train(net,pn,tn);              %开始训练

5 X' U$ P' B4 Z6 k$ N- q& V%%利用原始数据对BP网络仿真
an = sim(net,pn);                    %用训练好的模型进行仿真
a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
; F, y, n- v- Q. d' X0 Z
%%对比测试
x = 1990:2009;
/ D1 |& K5 H2 x, J! G1 @/ U2 ?newk = a(1,:);
3 B- [! w7 }* M8 _# Rnewh = a(2,:);
figure(2);
4 t7 [6 y) |% ]- g1 Isubplot(2,1,1);
& v0 d% H; p8 z* jplot(x,newk,'r-o',x,glkyl,'b--+');
legend('预测输出客运量','实际客运量');
" x5 ~  F0 H6 I7 o% Lxlabel('年份');
3 a* o+ x& C, H" zylabel('客运量/万人');
9 B6 f, [! v4 j2 I( B; Ztitle('运用数据箱客运量学习和测试对比图');
subplot(2,1,2);
2 B0 g8 i/ e4 w' B5 r0 zplot(x,newh,'r-o',x,glhyl,'b--+');
5 j5 B$ s8 @5 B% K2 tlegend('预测输出货运量','实际货运量');
& N9 _" e0 I( e7 Ixlabel('年份');
" {$ I# s4 b, J6 C5 aylabel('货运量/万吨');
: Z# F5 f+ R# H6 n, @' O2 S& dtitle('运用数据箱货运量学习和测试对比图');
%利用训练好的网络进行预测
pnew = [73.39 75.55
, W& k8 s  u7 C' a. |& l  `    3.9635 4.0975
1 F1 o6 N1 H2 h4 G& C6 w    0.9880 1.0268];
7 x" g, l9 i% gpnewn = tramnmx(pnew,minp,maxp);
%利用原始输入数据的归一化参数对新参数进行归一化
. E& }8 E  @/ h0 L* G, `, oanewn = sim(net,pnewn);
%利用归一化后的数据进行仿真
- D  F  |5 c8 Aanew = postmnmx(anewn,mint,maxt)
%把仿真得到的数据还原原始的数量级


0 k' w, h2 m9 I& V( O0 g" v4 ]欢迎交流！

子任

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

& ?0 g! ~' S+ p8 J