求 神经网络大神帮助

子任

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 K) m) V& x9 D9 t3 G) Q

& m) \1 e. _  o+ ka=max(p'); %#ok<UDIM>
0 k2 l4 ^* H9 k) T) i6 I0 V9 I' @; \& Dfor i=1:14
for j=1:3
. \2 x' V- B' o* E( Tptest(i,j)=p(i,j)/a(i);
end
end

5 N# Z# ^# I' D
$ K5 p0 [5 W8 }p1=ones(1,14);
p2=0.8.*ones(1,14);
p3=0.6.*ones(1,14);
* [) g1 I, u/ ?/ z& W) up4=0.4.*ones(1,14);
1 x6 L+ O. ~" D. u7 Rp5=0.2.*ones(1,14);
. s' K5 ?! \! b) B0 R  F/ j! _p6=0.*ones(1,14);
ptrain=[pl;p2;p3;p4;p5;p6];
# o" j8 ?2 _) y$ gt=[100,80,60,40,20,0];

" @) U# [: J4 q+ g$ f8 Y
net=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
net=init(net);
0 }! D" D: {) j/ \  gnet.trainparam.epochs=100000;
net.trainparam.goal=1e-10;
[net,tr]=train(net,ptrain',t);


: p& y. E% G" ~: @1 Sfor i=1:3
a=ptes(:,i);
score(i)=sim(net,a); %#ok<SAGROW>
end



0 M% j$ ]2 f, {+ X

数学中国—罂粟

你是要问什么啊？

子任

大神能帮我运行一下代码吗
0 U) T8 c# K, b5 d7 G) ]. u我不懂这个
能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
我这里有一个小模板，你可以看看，希望对你有用：
clc,clear
1 d$ J% \% F: h1 j7 o8 d: }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 ...
! q$ l0 g6 z5 c    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
; Q, y& h# i, C( N%人数 单位：万人
, v! ]2 @. ~* W) t: i; m" X7 usqjdcs = [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];
5 j$ ^* R3 l% V; N" S# ^# P%机动车数 单位：万辆
sqglmj = [0.09 0.11 0.11 0.14 0.20 0.23 0.23 0.32 0.32 0.34 0.36 0.36...
5 S; U* G9 s5 K0 t    0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];
8 J0 H, X7 v6 u/ t, b: V% c: d%公路面积 单位：万平方千米
glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
5 @; u/ o% M. V) z    19490 20433 22598 25107 33442 36836 40548 42927 43462];
%公路客运量 单位：万人
5 p9 {  x0 }+ pglhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
    11115 13320 16762 18673 20724 20803 21804];
0 u7 N% F, y6 r& C. O- i$ j8 G%公路货运量 单位：万吨
! S$ M* z+ O! Pp = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
' [+ \- x$ t+ ot = [glkyl;glhyl];          %目标矩阵
1 p/ M& U0 ?5 k6 |
%%归一化处理
+ r* l1 }3 _$ }) H  E' B[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
dx = [-1,1;-1 1;-1,1];                       %

+ Z/ L) N7 C' Q" J/ [7 v" v%BP网络训练
net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');

: g! D  t1 [) B- H%%建立模型，并用梯度下降法训练
" c9 ^4 ]2 L+ a. L8 I# @. qnet.trainParam.show = 1000;          %轮回显示一次
net.trainParam.Lr = 0.05;            %学习速率
net.trainParam.epochs = 50000;       %最大训练轮回
2 z9 M$ r  Z% W3 V7 e$ Nnet.trainParam.goal = 0.65*10^(-3);  %均方误差
' Q0 G1 `5 A. ^net = train(net,pn,tn);              %开始训练
" i( a" {% d/ @8 G- l9 S
%%利用原始数据对BP网络仿真
2 t' g( j. q( `8 j# Oan = sim(net,pn);                    %用训练好的模型进行仿真
6 n2 L) @! @1 @5 ]; z$ Za = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级

1 j+ I, I+ ^8 R3 A% S9 ?- I* H6 W%%对比测试
4 R" \$ J7 R- \3 ^# ux = 1990:2009;
newk = a(1,:);
* d1 n$ r& {' ~! G9 s" |* `newh = a(2,:);
1 }& y: b( ]) p9 q8 l; Bfigure(2);
subplot(2,1,1);
9 Q' ~7 a5 U; mplot(x,newk,'r-o',x,glkyl,'b--+');
6 m/ ]1 c( S' P  C& Llegend('预测输出客运量','实际客运量');
* I  V" v& U0 V+ r! z, h: {xlabel('年份');
+ |8 r# t  R; I/ N, Dylabel('客运量/万人');
title('运用数据箱客运量学习和测试对比图');
% H; S. S( d4 ]" H' _& ksubplot(2,1,2);
* e) Q) H; R% B( F% s. ^3 ~; pplot(x,newh,'r-o',x,glhyl,'b--+');
legend('预测输出货运量','实际货运量');
xlabel('年份');
ylabel('货运量/万吨');
9 J& z; X1 @- W/ S( Utitle('运用数据箱货运量学习和测试对比图');
; P2 W5 [; c" F+ V. l; l1 D" L%利用训练好的网络进行预测
& k( k. K7 ]5 \4 q( kpnew = [73.39 75.55
6 t( a# k+ p' ^0 I6 U0 V- ~    3.9635 4.0975
* ?" Y8 v6 p) Z. @8 O1 }1 M) M    0.9880 1.0268];
pnewn = tramnmx(pnew,minp,maxp);
%利用原始输入数据的归一化参数对新参数进行归一化
anewn = sim(net,pnewn);
6 r- O6 B/ R  S8 |4 W) r5 j%利用归一化后的数据进行仿真
anew = postmnmx(anewn,mint,maxt)
$ n* S) e- |+ k: }5 C8 n0 h1 _2 T5 T%把仿真得到的数据还原原始的数量级


欢迎交流！

" J  z4 W9 A* d6 n

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
5 h9 D/ y- }% ~7 x非常感谢你

! s. E. `; F0 |% c7 ?