
MATLAB中做了一个RBF网络,想加一步优化,不知道该怎么下手,请高手们帮忙分析一下。最好是用PSO优化。 10
m_data=[112.6016.82.980.99922;3.3017.32.880.99924;0124.215.62.920.99937;077.619.22.87...
m_data=[112.6 0 16.8 2.98 0.99922;
3.3 0 17.3 2.88 0.99924;
0 124.2 15.6 2.92 0.99937;
0 77.6 19.2 2.87 0.99939;
4.6 64.7 16.5 2.96 0.99938;
3.3 103.5 18.2 3.02 0.99934;
3.3 129.4 20.5 3.08 0.99939;
0 64.7 20.7 3.11 0.99938;
0 0 20.1 3.06 0.99922;
0 62.9 19.8 3.01 0.9996;
0 0 20.3 2.88 0.9999;
0 129.6 18.7 2.83 0.9994;
0 142.3 16.6 2.94 0.9993;
0 0 18.2 2.82 0.99921;
3.9 0 17.9 2.78 0.99925;
2.6 0 17.3 2.74 0.99924;
31.8 0 16.9 2.79 0.99925;
128.5 0 17.2 2.82 0.99923;
6.6 0 18.6 2.86 0.99922;
4.6 0 19.4 2.93 0.9993;
0 51.8 19.6 2.99 0.9994;
2.6 0 19.1 3.03 0.99928;
0 0 20.3 3.07 0.9993;
0 51.7 21.2 3.13 0.99935;
1.3 0 20.9 3.11 0.99925;
1.3 77.6 18.5 3.06 0.99936;
1.3 25.6 17.9 2.96 0.99938;
33.1 0 18.3 2.91 0.99928;
4.6 0 18.2 2.93 0.99925;
2 0 19.4 2.84 0.9993;
78.1 77.6 19.7 2.82 0.99939;
29.8 51.7 20.3 2.95 0.99937;
1.3 25.9 19.5 2.93 0.99936;
20.5 126.8 19.1 2.97 0.99938;
0 0 18.7 3.02 0.99925;
101.9 0 18.4 3.06 0.99924;
65.6 0 18.2 3.04 0.99926;
13.2 0 17.9 3.01 0.99925;
13.2 0 18.1 3.02 0.99928];
X=m_data(:,1:4);T=m_data(:,5);T=T';
C=X;
delta=cov(X');
delta=sum(delta);
for i=1:1:39
for j=1:1:39
H(i,j)=((X(i,:)-C(j,:)))*((X(i,:)-C(j,:))');
H(i,j)=exp(-H(i,j)./delta(j));
end
end
p=H;
err_goal=0.001;
sc=3;
net=newrb(p,T,err_goal,sc,1000,1);
Y=sim(net,p);
E=T-Y;
SSE=sse(E);
MSE=mse(E);
figure;
plot(T);
hold on;
plot(Y,'r');
title('RBF拟合曲线');
ylabel('氯乙烯的纯度');
xlabel('输入数据样本点');
legend('离线测量值','软测量值') 展开
3.3 0 17.3 2.88 0.99924;
0 124.2 15.6 2.92 0.99937;
0 77.6 19.2 2.87 0.99939;
4.6 64.7 16.5 2.96 0.99938;
3.3 103.5 18.2 3.02 0.99934;
3.3 129.4 20.5 3.08 0.99939;
0 64.7 20.7 3.11 0.99938;
0 0 20.1 3.06 0.99922;
0 62.9 19.8 3.01 0.9996;
0 0 20.3 2.88 0.9999;
0 129.6 18.7 2.83 0.9994;
0 142.3 16.6 2.94 0.9993;
0 0 18.2 2.82 0.99921;
3.9 0 17.9 2.78 0.99925;
2.6 0 17.3 2.74 0.99924;
31.8 0 16.9 2.79 0.99925;
128.5 0 17.2 2.82 0.99923;
6.6 0 18.6 2.86 0.99922;
4.6 0 19.4 2.93 0.9993;
0 51.8 19.6 2.99 0.9994;
2.6 0 19.1 3.03 0.99928;
0 0 20.3 3.07 0.9993;
0 51.7 21.2 3.13 0.99935;
1.3 0 20.9 3.11 0.99925;
1.3 77.6 18.5 3.06 0.99936;
1.3 25.6 17.9 2.96 0.99938;
33.1 0 18.3 2.91 0.99928;
4.6 0 18.2 2.93 0.99925;
2 0 19.4 2.84 0.9993;
78.1 77.6 19.7 2.82 0.99939;
29.8 51.7 20.3 2.95 0.99937;
1.3 25.9 19.5 2.93 0.99936;
20.5 126.8 19.1 2.97 0.99938;
0 0 18.7 3.02 0.99925;
101.9 0 18.4 3.06 0.99924;
65.6 0 18.2 3.04 0.99926;
13.2 0 17.9 3.01 0.99925;
13.2 0 18.1 3.02 0.99928];
X=m_data(:,1:4);T=m_data(:,5);T=T';
C=X;
delta=cov(X');
delta=sum(delta);
for i=1:1:39
for j=1:1:39
H(i,j)=((X(i,:)-C(j,:)))*((X(i,:)-C(j,:))');
H(i,j)=exp(-H(i,j)./delta(j));
end
end
p=H;
err_goal=0.001;
sc=3;
net=newrb(p,T,err_goal,sc,1000,1);
Y=sim(net,p);
E=T-Y;
SSE=sse(E);
MSE=mse(E);
figure;
plot(T);
hold on;
plot(Y,'r');
title('RBF拟合曲线');
ylabel('氯乙烯的纯度');
xlabel('输入数据样本点');
legend('离线测量值','软测量值') 展开
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询