求 神经网络大神帮助

子任

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

3 ~' R8 m3 u% _! I5 G5 ^6 u
a=max(p'); %#ok<UDIM>
- a9 n5 P. O( ]+ u! X7 g5 Kfor i=1:14
for j=1:3
! J3 Q$ t( o' b6 ~! N5 J+ B( Rptest(i,j)=p(i,j)/a(i);
end
$ V+ e& ?7 r, T- @end


p1=ones(1,14);
- r; M" N6 B+ W$ O) ^9 dp2=0.8.*ones(1,14);
2 q+ u8 w+ s/ e1 g; Zp3=0.6.*ones(1,14);
p4=0.4.*ones(1,14);
p5=0.2.*ones(1,14);
; h3 K0 e2 o) b# R- t4 x) G1 kp6=0.*ones(1,14);
ptrain=[pl;p2;p3;p4;p5;p6];
t=[100,80,60,40,20,0];

* G3 W/ F3 s1 ?1 J6 ~) p1 Z, f! U
net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
' j* l# E/ ~. L) Wnet=init(net);
6 i& n) V' X9 v# |net.trainparam.epochs=100000;
4 z/ y0 n( m  S) b# i) V  T" Unet.trainparam.goal=1e-10;
' k: P3 |  c! c[net,tr]=train(net,ptrain',t);
2 r$ c3 W! w0 Z/ w/ Q) K4 {* P  L

  v' N! W+ X  lfor i=1:3
a=ptes(:,i);
) p; m7 [4 X4 X; x) ^( fscore(i)=sim(net,a); %#ok<SAGROW>
end

9 A! c, }* G# y; q% L/ B$ j

数学中国—罂粟

你是要问什么啊？
! V! T% q; `9 M" \! G

子任

大神能帮我运行一下代码吗
* ?0 h4 Q$ L% r) H我不懂这个
能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
/ l0 }  I4 l* S3 c6 s% w7 V# X我这里有一个小模板，你可以看看，希望对你有用：
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 ...
  }" N7 P! `4 C2 S! D7 [' \    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
%人数 单位：万人
) [+ B4 a4 p& P! U* nsqjdcs = [0.6 0.75 0.85 0.90 1.05 1.35 1.45 1.60 1.70 1.85 2.15 2.20...
# V& S. r6 E# v" |4 q7 f    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
%机动车数 单位：万辆
( p7 _; I7 m5 H* A1 Jsqglmj = [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];
; h4 d; t% Q' j: c%公路面积 单位：万平方千米
7 f5 c6 `, A; Q$ l0 U5 Xglkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    19490 20433 22598 25107 33442 36836 40548 42927 43462];
%公路客运量 单位：万人
, A) E9 V+ A+ f: G: Eglhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    11115 13320 16762 18673 20724 20803 21804];
" D; `! K, x6 q%公路货运量 单位：万吨
9 F  m6 ~! D! u: `" Kp = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
t = [glkyl;glhyl];          %目标矩阵
' v8 Z5 w/ x- h4 E5 n* i, P
! ?& `9 x$ [( {. b%%归一化处理
[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
dx = [-1,1;-1 1;-1,1];                       %
& I3 c* c1 ]. W5 P) I
4 C( `  e' Q# d$ A+ i4 {6 M%BP网络训练
net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');

6 g/ C' \# B4 e9 u6 `" S%%建立模型，并用梯度下降法训练
net.trainParam.show = 1000;          %轮回显示一次
& N% g, g5 i  l9 Fnet.trainParam.Lr = 0.05;            %学习速率
net.trainParam.epochs = 50000;       %最大训练轮回
6 ^7 q- q! ^, U* r' Dnet.trainParam.goal = 0.65*10^(-3);  %均方误差
4 v( G4 B& ?" Y& r! t: Z$ @net = train(net,pn,tn);              %开始训练

%%利用原始数据对BP网络仿真
0 y* k6 `7 h# P# r" F& v, O: Xan = sim(net,pn);                    %用训练好的模型进行仿真
  H9 Q  f0 \. na = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级

* m4 Z. `3 l. ?0 J% P8 F( X0 Y%%对比测试
; s5 N( I2 o4 f1 S2 Vx = 1990:2009;
newk = a(1,:);
+ o! Z# m5 a- g9 y% m  Z1 Xnewh = a(2,:);
# H: D1 N$ G  ]  S- k# T3 Nfigure(2);
; |! l' ?1 N7 t$ `& @! usubplot(2,1,1);
" G' j- R8 p5 T2 t; A9 j5 yplot(x,newk,'r-o',x,glkyl,'b--+');
legend('预测输出客运量','实际客运量');
" ^& S- q( \& M8 Rxlabel('年份');
) U3 r- x5 `, V+ [  P# xylabel('客运量/万人');
title('运用数据箱客运量学习和测试对比图');
subplot(2,1,2);
plot(x,newh,'r-o',x,glhyl,'b--+');
! J4 j1 o$ e5 nlegend('预测输出货运量','实际货运量');
% @- m- y8 `& jxlabel('年份');
& h9 [! z/ M3 [0 U/ E; lylabel('货运量/万吨');
% D- }5 v0 H) R! E% qtitle('运用数据箱货运量学习和测试对比图');
%利用训练好的网络进行预测
- u; O. _2 R( ^: e+ F8 J# q0 bpnew = [73.39 75.55
    3.9635 4.0975
    0.9880 1.0268];
0 A4 V2 X  X; F- N: A1 |+ upnewn = tramnmx(pnew,minp,maxp);
%利用原始输入数据的归一化参数对新参数进行归一化
$ R1 o% I  f( tanewn = sim(net,pnewn);
%利用归一化后的数据进行仿真
; _5 }, u& G( _# h( [0 l) i( ?7 xanew = postmnmx(anewn,mint,maxt)
7 @3 V, I! |! R$ N%把仿真得到的数据还原原始的数量级

2 B0 Q) K  u  d2 x7 O
+ P* q7 E; P% l  h: M欢迎交流！
2 O! g$ m2 O: Z/ t
5 y* u2 z. y4 z" P2 O4 e

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
非常感谢你
; v0 g, f" Z! z$ o& q, A& a
* l0 B( I; O& Z: |