求 神经网络大神帮助

子任

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];
( C6 B& N, b  ~
- ]; t% k; _! M3 F/ e
a=max(p'); %#ok<UDIM>
for i=1:14
+ L  V) `( q- e% D- p6 f7 Qfor j=1:3
ptest(i,j)=p(i,j)/a(i);
end
end
# V- o; W" C6 D* B3 B9 A5 q
! u6 K3 S, u5 @" g% V& E( L
p1=ones(1,14);
: X) p( X* C0 yp2=0.8.*ones(1,14);
/ A9 H; k+ p9 M# E3 F- o, i) fp3=0.6.*ones(1,14);
p4=0.4.*ones(1,14);
p5=0.2.*ones(1,14);
! O+ N1 t+ b$ i1 c7 n2 Np6=0.*ones(1,14);
" y5 Q5 I/ u/ o0 Q7 i2 i: Fptrain=[pl;p2;p3;p4;p5;p6];
) I5 S' N3 j' l" Ht=[100,80,60,40,20,0];


# V7 Q' \/ T' G; z% [8 l+ cnet=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
net=init(net);
0 e$ }5 c" i9 Lnet.trainparam.epochs=100000;
net.trainparam.goal=1e-10;
3 c! t7 I7 e9 H) r4 J[net,tr]=train(net,ptrain',t);

  b: q  E- Z  @
for i=1:3
a=ptes(:,i);
$ G3 c; R) T, m# o$ rscore(i)=sim(net,a); %#ok<SAGROW>
end
* d# ]8 j0 Z+ l


, V+ x9 o( C* @% K8 q3 `

数学中国—罂粟

你是要问什么啊？
) {+ ]7 Z4 i% O& Y1 U; a% q

子任

大神能帮我运行一下代码吗
3 ?# S' f& ~4 N/ P+ r我不懂这个
能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
我这里有一个小模板，你可以看看，希望对你有用：
( a8 W5 P: u9 Q  kclc,clear
5 ~8 S+ I* p* N& g! p, r$ h1 ?4 o( jsqrs = [20.55 22.44 25.37 27.13 29.45 30.1 30.96 34.06 36.42 38.09 39.13 39.99 ...
; L- ?3 q8 Z. ^3 y; @' {+ K    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
. U9 O- f, S6 }# K- c" \0 f$ y%人数 单位：万人
/ n- S! D& Q; D, J# ^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];
+ n& k9 h2 K2 S1 U, K" @( G%机动车数 单位：万辆
5 F+ p3 S) F; bsqglmj = [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];
' H( ?; G3 \5 V) C%公路面积 单位：万平方千米
glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
, K6 N' `4 m2 t5 n4 G1 z9 x    19490 20433 22598 25107 33442 36836 40548 42927 43462];
  f% }; \) x" s; x2 D%公路客运量 单位：万人
glhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
- ?; [8 r- C6 S5 t+ \    11115 13320 16762 18673 20724 20803 21804];
! M5 E& e/ e6 N; ~" C# Q) |. Z%公路货运量 单位：万吨
p = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
t = [glkyl;glhyl];          %目标矩阵
3 ?2 j  U  ^' ~
# H1 r$ z4 z9 z6 c2 a%%归一化处理
6 g/ m& y/ K! D7 `# k; }& i( M[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
dx = [-1,1;-1 1;-1,1];                       %
  b1 |  E2 q6 w1 p7 H$ ^7 y( U/ X
%BP网络训练
6 g/ _" w: @2 x$ U: }, A8 G# Knet = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');

/ x3 Z, J! d9 z' i% f9 a4 D% x0 J%%建立模型，并用梯度下降法训练
: e. e) i+ B% z' Mnet.trainParam.show = 1000;          %轮回显示一次
$ N1 \* {/ X8 c* o8 A* I, |3 Onet.trainParam.Lr = 0.05;            %学习速率
net.trainParam.epochs = 50000;       %最大训练轮回
net.trainParam.goal = 0.65*10^(-3);  %均方误差
5 O1 N5 r: R+ C- h: e/ Dnet = train(net,pn,tn);              %开始训练

% I% K7 Y$ ~; z* ^: t%%利用原始数据对BP网络仿真
( D9 [2 ^/ c9 _! e9 K2 jan = sim(net,pn);                    %用训练好的模型进行仿真
/ u. W% z9 p- |/ }/ U1 t+ Va = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
  N& |6 O- r9 \
%%对比测试
x = 1990:2009;
4 X6 X; u% d" k* b% snewk = a(1,:);
2 k" I8 ?6 a! R/ W! W, b9 N3 \* knewh = a(2,:);
figure(2);
subplot(2,1,1);
plot(x,newk,'r-o',x,glkyl,'b--+');
& O. L2 }3 S! [legend('预测输出客运量','实际客运量');
* P4 B" g' l% kxlabel('年份');
ylabel('客运量/万人');
title('运用数据箱客运量学习和测试对比图');
subplot(2,1,2);
plot(x,newh,'r-o',x,glhyl,'b--+');
7 b" F" W' m, O1 D3 e2 T$ Wlegend('预测输出货运量','实际货运量');
xlabel('年份');
ylabel('货运量/万吨');
title('运用数据箱货运量学习和测试对比图');
2 ]6 I1 g1 ~% E" F3 h! ]: c: {%利用训练好的网络进行预测
pnew = [73.39 75.55
    3.9635 4.0975
( A7 y9 E, L9 v. D$ {0 g6 p6 @9 P    0.9880 1.0268];
pnewn = tramnmx(pnew,minp,maxp);
%利用原始输入数据的归一化参数对新参数进行归一化
anewn = sim(net,pnewn);
%利用归一化后的数据进行仿真
anew = postmnmx(anewn,mint,maxt)
%把仿真得到的数据还原原始的数量级
  O( m/ ~3 h  V7 k: x( C
6 `" X. _% ~; t$ f- W8 K$ x& c! h
欢迎交流！

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
非常感谢你
$ X' W2 {5 h2 {2 U" P
; e% Z9 B8 }4 w$ N* Q3 q0 M8 g7 X