该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
clc;
clear;
p=[2.259;2.257;2.256;2.254;2.252;2.248;2.247;2.245;2.244;2.243;2.239;2.238;2.236;2.235;2.234;2.231;2.229;2.228;2.226;2.225;2.221;2.220;
2.219;2.217;2.216;2.211;2.209;2.208;2.207;2.206;2.202;2.201;2.199;
2.198;2.197;2.194;2.192;2.191;2.189;2.187]';
t=[298.24 ;298.01 ;297.79 ;297.55 ;297.32 ;296.42 ;296.19;295.97 ;295.73 ;295.48 ;294.59 ;294.36 ;294.12 ;293.89 ;293.66 ;292.74 ;292.52;292.31 ;
292.04 ;291.81;290.92 ;290.68 ;290.45 ;290.22 ;289.98 ;289.03 ;288.78 ;
288.51 ;288.28;288.02 ;287.13 ;286.87 ;286.59 ;286.36 ;286.14 ;285.21 ;
284.98 ;284.75;284.54 ;284.31];
P_test=[2.183;2.182;2.181;2.179;2.178;2.175;2.174;2.173;2.171;2.169]';
t1=[283.46 ;283.19 ;282.91 ;282.68 ;282.44 ;281.53 ;281.29;281.02 ;280.78 ;
280.55];
[pn,minp,maxp,tn,mint,maxt] = premnmx(p,t);
p2= tramnmx(P_test,minp,maxp);
net=newff(minmax(pn),[1,3 ,1],{'tansig','tansig','purelin'},'traingdm');
net.trainparam.show=50;
net.trainparam.mc=0.9;
net.trainparam.lr=0.05;
net.trainparam.epochs=2000;
net.trainparam.goal=0.001;
net=init(net);
[net,tr]=train(net,pn,tn);
PN=sim(net,p2);
[t2]= postmnmx(PN,mint,maxt);
E = t1 - t2;
MSE=mse(E);
%作图表示实测值和仿真值
figure(1);
X=1:5;
plot(X,t2,'r',X,t1,'b');
title(' b为真实值,r为预测值');
大佬 我这个该怎么修改呀? 我用的BP神经网络,求大佬回复一下 QAQ