clear
X1=[2.5 3.5 4.5 5.5 6.5 7.5 8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5
3.5 4.5 5.5 6.5 7.5 8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5
4.5 5.5 6.5 7.5 8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5
5.5 6.5 7.5 8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5
6.5 7.5 8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5
7.5 8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5 24.5
8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5 24.5 25.5
9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5 24.5 25.5 26.5
10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5 24.5 25.5 26.5 27.5
11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5 24.5 25.5 26.5 27.5 28.5
]';
Y1 = [12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5 24.5 25.5 26.5 27.5 28.5 29.5
]';
X2=[20.5 21.5 22.5 23.5 24.5 25.5 26.5 27.5 28.5 29.5
21.5 22.5 23.5 24.5 25.5 26.5 27.5 28.5 29.5 30.5
22.5 23.5 24.5 25.5 26.5 27.5 28.5 29.5 30.5 31.5
23.5 24.5 25.5 26.5 27.5 28.5 29.5 30.5 31.5 32.5
24.5 25.5 26.5 27.5 28.5 29.5 30.5 31.5 32.5 33.5
25.5 26.5 27.5 28.5 29.5 30.5 31.5 32.5 33.5 34.5
26.5 27.5 28.5 29.5 30.5 31.5 32.5 33.5 34.5 35.5
27.5 28.5 29.5 30.5 31.5 32.5 33.5 34.5 35.5 36.5
28.5 29.5 30.5 31.5 32.5 33.5 34.5 35.5 36.5 37.5
29.5 30.5 31.5 32.5 33.5 34.5 35.5 36.5 37.5 38.5
]';
X3=[2.5 3.5 4.5 5.5 6.5 7.5 8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5
3.5 4.5 5.5 6.5 7.5 8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5
4.5 5.5 6.5 7.5 8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5
5.5 6.5 7.5 8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5
6.5 7.5 8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5
7.5 8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5 24.5
8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5 24.5 25.5
9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5 24.5 25.5 26.5
10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5 24.5 25.5 26.5 27.5
11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5 24.5 25.5 26.5 27.5 28.5
]';
ystar= [30.5 31.5 32.5 33.5 34.5 35.5 36.5 37.5 38.5 39.5
]';
gam = 5000;
sig2 = 10000;
type = 'function estimation';
[alpha,b] = trainlssvm({X1,Y1,type,gam,sig2});
Y2 = simlssvm({X1,Y1,type,gam,sig2,'RBF_kernel','preprocess'},X2)
Y3 = simlssvm({X1,Y1,type,gam,sig2,'RBF_kernel','preprocess'},X3)
Xh=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28];
figure(1);
Y4=[Y3',Y2'];
Y5=[Y1',ystar'];
plot(Xh,Y4,'ro');
hold on
plot(Xh,Y4,'b-')
plot(Xh,Y5,'bp')
plot(Xh,Y5,'b--')
legend('回归值', '实线','实测值')%图例
xlabel('时间(月)')
ylabel('边坡变形累计值(mm)')
title('支持向量机实测与回归曲线');
figure(2);
[m,b,r]=postreg(Y2,ystar) %相关系数
sgs=1
评论0