求 神经网络大神帮助

子任

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


  L! z9 C9 P& }/ Ba=max(p'); %#ok<UDIM>
. k; X4 j" z2 _, S) r8 @for i=1:14
! D$ E4 N0 X; g) W7 z9 K1 pfor j=1:3
ptest(i,j)=p(i,j)/a(i);
end
end
! H% z. W, W- K' Q5 x

p1=ones(1,14);
p2=0.8.*ones(1,14);
p3=0.6.*ones(1,14);
p4=0.4.*ones(1,14);
* Y; B! w; e* @% d! y" dp5=0.2.*ones(1,14);
p6=0.*ones(1,14);
' q- }8 _8 a* A" c' hptrain=[pl;p2;p3;p4;p5;p6];
t=[100,80,60,40,20,0];
( `* Y' M" u3 y5 I

3 x# x; K% y  O# n- c% [9 M# w+ `: t  X( F9 Qnet=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
) F, t* Q: G6 y* x+ W, ?net=init(net);
net.trainparam.epochs=100000;
net.trainparam.goal=1e-10;
3 D8 B/ h( {, `! g[net,tr]=train(net,ptrain',t);
9 k$ ]: ~% W" W7 Y
/ h/ s; U0 k5 Z8 g  D$ W; `* U
; S. L1 E6 d5 D0 y# L. E5 X8 C8 _for i=1:3
" b- c$ y/ [" Z1 c) E' `a=ptes(:,i);
3 P1 [$ ]) B% E. yscore(i)=sim(net,a); %#ok<SAGROW>
, I: O3 @: u; [6 K; E6 gend

/ ^9 d) n- u! R/ x3 ]  ^6 o

  s- L* k+ R, m: j, z7 a7 X

数学中国—罂粟

你是要问什么啊？

子任

大神能帮我运行一下代码吗
& J0 e5 t2 \, z. _: U( K7 f8 m* ?: j; b$ m我不懂这个
8 Q6 @: P! z) o能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
" z4 ^: B/ \) a$ J$ Y6 l2 ]我这里有一个小模板，你可以看看，希望对你有用：
% ^7 V& x/ H8 ~( Eclc,clear
; u, j9 O; m- R- \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 ...
    41.93 44.59 47.30 52.89 55.73 56.76 59.17 60.63];
%人数 单位：万人
( J% j7 h" X( W! X0 T; a3 _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];
%机动车数 单位：万辆
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...
9 P: B& C% I; L% A    0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];
%公路面积 单位：万平方千米
glkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
    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];
$ G* V3 X, R& ?  y! c& G%公路货运量 单位：万吨
: Q) g; W" b+ y) U( Z1 q& J# Ep = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
t = [glkyl;glhyl];          %目标矩阵
* O- s" S" u2 j2 c5 w. @
%%归一化处理
' R7 l8 D% S) b! N[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
7 U& ~- H; T: Z( _dx = [-1,1;-1 1;-1,1];                       %

%BP网络训练
net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');

) q5 X: P" {2 m3 b, y%%建立模型，并用梯度下降法训练
+ ?( M4 I+ g2 e) h- A; N& Rnet.trainParam.show = 1000;          %轮回显示一次
& b. r8 y" r+ h% ^2 f- D4 w4 ?8 Xnet.trainParam.Lr = 0.05;            %学习速率
6 o1 Y" E6 G* C5 n8 Z6 onet.trainParam.epochs = 50000;       %最大训练轮回
/ p+ W( q% D" f# dnet.trainParam.goal = 0.65*10^(-3);  %均方误差
1 Y5 X( v8 H$ |. F' E- X2 |0 i5 ~net = train(net,pn,tn);              %开始训练
4 H- F2 Z0 U' l. V1 ~" Y
%%利用原始数据对BP网络仿真
* T2 d7 c7 @" t7 W  `! Oan = sim(net,pn);                    %用训练好的模型进行仿真
a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级
- x, F0 @) y/ T3 K
%%对比测试
( S& F% C- ^% b5 |& ?0 ?3 ox = 1990:2009;
newk = a(1,:);
newh = a(2,:);
figure(2);
subplot(2,1,1);
plot(x,newk,'r-o',x,glkyl,'b--+');
legend('预测输出客运量','实际客运量');
9 j$ b1 q3 D$ m5 Q) l8 N: }xlabel('年份');
ylabel('客运量/万人');
title('运用数据箱客运量学习和测试对比图');
& r( b/ e9 F- zsubplot(2,1,2);
( C' D1 C+ p) @. m- R  V# v& Rplot(x,newh,'r-o',x,glhyl,'b--+');
legend('预测输出货运量','实际货运量');
xlabel('年份');
$ o, B7 k/ g; _' K% G. ?! iylabel('货运量/万吨');
title('运用数据箱货运量学习和测试对比图');
%利用训练好的网络进行预测
* B! z) b- S0 W7 j% W" X% `: ]pnew = [73.39 75.55
    3.9635 4.0975
    0.9880 1.0268];
pnewn = tramnmx(pnew,minp,maxp);
1 }" A3 [; X  v2 G$ J%利用原始输入数据的归一化参数对新参数进行归一化
anewn = sim(net,pnewn);
%利用归一化后的数据进行仿真
anew = postmnmx(anewn,mint,maxt)
1 B1 v; ^$ V  P& B5 f%把仿真得到的数据还原原始的数量级
& \3 k, i9 H$ ~% t5 X# G5 @) t; z
" f  |9 v6 h) [$ p. w/ p
欢迎交流！

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
, Z/ ^7 t" f9 Z% I- x; }. h( Q) q非常感谢你
; y2 ^4 U9 \. L% U$ p1 d
2 S5 j# e: `6 G' H4 t