问题描述

matlab 代码错在哪 求解决

想要用12个数据去预测1个数据
P=[2.97 3.02 3.29 3.35 5.6 5.41 4.82 3.61 2.92 2.03 1.95 1.94 3.01 3.18 3.45 3.51 5.82 5.58 4.95 3.95 3.01 2.21 2.13 2.13 3.2 3.31 3.56 3.72 6.01 5.73 5.2 4.11 3.21 2.33 2.3 2.3 3.4 3.42 3.8 3.88 6.21 5.88 5.41 4.31 3.35 2.48 2.4 2.41 3.45 3.51 3.62 3.91 6.31 6.02 5.52 4.37 3.39 2.52 2.48 2.45 3.84 3.84 3.95 4.23 6.59 6.41 5.83 4.76 3.61 2.61 2.6 2.55 4.15 4.17 4.25 4.41 7.05 6.93 6.31 5.2 3.81 2.95 2.85 2.84 4.25 4.3 4.33 4.43 7.23 7 6.65 5.48 3.9 3.61 3.5 3.47 4.5 4.53 4.7 4.99 7.68 7.45 7.3 5.73 4.55 3.86 3.65 3.65 4.75 4.78 4.81 4.91 7.97 7.81 7.8 6.14 4.79 4.01 3.9 3.9 4.6 4.63 4.8 5 8.16 8.01 7.92 6.37 5.4 4.76 4.27 4.25 4.6 4.65 4.83 5 8.47 8 7.9 6.47 5.52 4.93 4.6 4.6 5.01 5.11 5.1 5.29 8.67 8.21 8.06 6.69 5.73 5.18 4.95 4.95 5.03 5.33 5.63 6.1 11.07 10.93 9.1 8.02 7.52 5.33 5.19 5.19 5.22 5.22 5.34 6.17 11.59 11.85 13.06 10.58 8.62 7.5 6.57 6.57 6.49 6.39 6.4 6.55 10.25 13.76 12.75 10.92 8.95 7.86 7.01 7];

n=0;
for i=1:180
for j=(n+1):(n+12)
a(i,j-n)=P(j);
end
n=n+1;
end
t(1,:)=P(13:192);
p=a';
[pn,minp,maxp,tn,mint,maxt]=premnmx(p,t);
net=newff([-1 1;-1 1;-1 1],[12,1],{'tansig','purelin'},'traingd');
net.trainParam.Lr=0.01;
net.trainParam.epochs=50000;
net.trainParam.goal=0.001;
net=train(net,pn,tn);
%利用原始数据对BP网络仿真
an=sim(net,pn);
a=postmnmx(an,mint,maxt);
%数据对比
x=13:192;
new=a;
figure(2);
plot(x,new,'r-o',x,t,'b--+');
legend('BP模拟值','实测值');
xlabel('年份');
ylabel('地下水埋深/m');
title('数据对比图');
%评价指标mse
%e=t-an
%perf=mse(e)
%预测
pnew=[6.49 6.39 6.40 6.55 10.25 13.76 12.75 10.92 8.95 7.86 7.01 7.00]';
[pnewn,minpnew,maxnew]=premnmx(pnew);
anewn=sim(net,pnewn);
anew=postmnmx(anewn,mint,maxt);
anew

解决方案

楼主 这是我积累的matlab的一些资料 希望对你有帮助吧
http://download.csdn.net/album/detail/1229