求 神经网络大神帮助

子任

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];
* B' J5 V* I1 j/ B# h

. [# P: u& y! Y' n7 Wa=max(p'); %#ok<UDIM>
for i=1:14
, Y  T% W+ m' [, G/ _! Sfor j=1:3
ptest(i,j)=p(i,j)/a(i);
7 g2 ?% c/ `  U( b$ k% l7 V% @% M" [* [end
end
& I1 ]- \5 q* W2 V

* q8 g" s# v) S" D2 b( X1 ~* }p1=ones(1,14);
p2=0.8.*ones(1,14);
" E2 ~! f% }+ q( np3=0.6.*ones(1,14);
p4=0.4.*ones(1,14);
p5=0.2.*ones(1,14);
( \+ f% c9 h) S5 Rp6=0.*ones(1,14);
ptrain=[pl;p2;p3;p4;p5;p6];
9 N2 c$ ~( Q5 B/ P7 ]t=[100,80,60,40,20,0];


- J1 E( d+ `- Y1 E! {, anet=newff(minmax(ptrain'),[59,1],{tansig','purelin'},traingd');
net=init(net);
9 \0 s: p9 J" P* inet.trainparam.epochs=100000;
$ Q- e+ {* M( D9 B1 z  W$ y( g- ~net.trainparam.goal=1e-10;
[net,tr]=train(net,ptrain',t);

# n0 J) A/ o. Z! Z5 c9 \# t" Z7 B; V' N
for i=1:3
a=ptes(:,i);
9 }# v9 @, M& `9 t  p3 lscore(i)=sim(net,a); %#ok<SAGROW>
2 c, f' L% {, U/ Oend

9 L8 @2 C' X% b1 [7 S; H

数学中国—罂粟

你是要问什么啊？

子任

大神能帮我运行一下代码吗
我不懂这个
. t5 r& \3 Y, D- f: Q能留下联系方式吗

QYY1994

这程序是你自己编的吧！你能说说你想做的是什么吗？
/ j/ U" z7 A3 U4 S0 \. H* ]我这里有一个小模板，你可以看看，希望对你有用：
clc,clear
" }, W3 \' ~! k+ N4 Bsqrs = [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];
# F+ G" {; F+ Z4 U' z%人数 单位：万人
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...
& W# H2 ~; m" b* k) W8 ^    2.25 2.35 2.50 2.60 2.70 2.85 2.95 3.10];
+ ]/ y/ e  S" u7 c0 D, T%机动车数 单位：万辆
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...
    0.38 0.49 0.56 0.59 0.59 0.67 0.69 0.79];
%公路面积 单位：万平方千米
$ T3 a2 F$ ^% w% n* Eglkyl = [5126 6217 7730 9145 10460 11387 12353 15750 18304 19836 21024....
3 j5 ?% E4 o# G# w    19490 20433 22598 25107 33442 36836 40548 42927 43462];
. u) U2 l7 w1 a$ B3 K6 p1 |%公路客运量 单位：万人
8 L0 }  ]8 P: x: Kglhyl = [1237 1379 1385 1399 1663 1714 1834 4322 8132 8936 11099 11203 10524....
2 P* K5 M: J8 j% q. h    11115 13320 16762 18673 20724 20803 21804];
9 [) }2 u" t) H9 ^8 `- M0 w: M%公路货运量 单位：万吨
$ T% O- C4 }: c/ wp = [sqrs;sqjdcs;sqglmj];   %输入数据矩阵
t = [glkyl;glhyl];          %目标矩阵

" p+ ~* Y" ^1 Y/ v0 j%%归一化处理
3 R) ?  I& U6 N" {[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);  %对pt矩阵进行归一化处理
# h6 \% F% s( ^! D1 k; S9 O# c: Adx = [-1,1;-1 1;-1,1];                       %
# D' |% v5 H- T! f4 l; Q
%BP网络训练
net = newff(dx,[3,7,2],{'tansig','tansig','purelin'},'traingdx');
# E6 h0 p9 ^8 y! g7 a8 |5 W
%%建立模型，并用梯度下降法训练
net.trainParam.show = 1000;          %轮回显示一次
; L4 _4 {  v, S# E" P0 [  D- ]0 b$ Pnet.trainParam.Lr = 0.05;            %学习速率
net.trainParam.epochs = 50000;       %最大训练轮回
net.trainParam.goal = 0.65*10^(-3);  %均方误差
net = train(net,pn,tn);              %开始训练

0 |1 k6 b' {+ X! g%%利用原始数据对BP网络仿真
/ h# s# i6 |8 @an = sim(net,pn);                    %用训练好的模型进行仿真
a = postmnmx(an,mint,maxt);          %把仿真到的数据还原到原始的数量级

" h; S2 _- r! l2 L1 U+ ?%%对比测试
x = 1990:2009;
newk = a(1,:);
newh = a(2,:);
figure(2);
subplot(2,1,1);
2 g  |3 a2 X( xplot(x,newk,'r-o',x,glkyl,'b--+');
2 x2 H4 E3 c6 \. F4 H0 tlegend('预测输出客运量','实际客运量');
xlabel('年份');
  i5 x/ M% B2 z7 `3 Z, X8 q4 \, {: a  fylabel('客运量/万人');
title('运用数据箱客运量学习和测试对比图');
% E+ H0 e/ ]7 t4 x; F6 Vsubplot(2,1,2);
plot(x,newh,'r-o',x,glhyl,'b--+');
legend('预测输出货运量','实际货运量');
1 A5 T; b1 O9 O8 \! `, F) d; Cxlabel('年份');
ylabel('货运量/万吨');
title('运用数据箱货运量学习和测试对比图');
%利用训练好的网络进行预测
pnew = [73.39 75.55
: d; d% A2 ]' m8 A* s    3.9635 4.0975
6 e/ {6 T5 S  C3 v' U" v    0.9880 1.0268];
pnewn = tramnmx(pnew,minp,maxp);
%利用原始输入数据的归一化参数对新参数进行归一化
anewn = sim(net,pnewn);
%利用归一化后的数据进行仿真
anew = postmnmx(anewn,mint,maxt)
6 W7 {5 b5 ^/ F, P- I! q%把仿真得到的数据还原原始的数量级


- |6 q$ v# g. x1 o' B: J$ ~; W9 C欢迎交流！
1 d5 x9 w4 b: P8 I# c$ `+ H& w5 n
6 }+ J) W# }. I& t6 G0 \: i. ~

子任

恩恩，我想做的是一个评估！就是对几个对象做一个评估，得出一个分数
( q7 X! L5 {( B" w% {$ ]非常感谢你
1 B0 c1 ^/ f8 Q& D. K6 \9 o
- ~& G5 }# _1 N& q