|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
clc
clear
close all
A=load('data.txt') %读取所有数据
%---------------------------------------------------------
%选取1-400个数据
x=1:1:400
y=A(1:1:400)
%选取1-400,间隔为2的200个数据作为样本数据
x_train=1:2:400
y_train=A(1:2:400)
%选取2-400,间隔为2的200个数据作为测试数据
x_test=2:2:400
y_test=A(2:2:400)
%定义SVR 参数
global p1 ;
p1=3
ker='erbf'; % 核函数 k = exp(-(u-v)*(u-v)'/(2*p1^2))
C=10000000;
e=0.1;
loss='einsensitive';
%训练支持向量机
[nsv beta bias] = svr(x_train',y_train,ker,C,loss,e); % 训练样本数据。
y_test1 = svroutput(x_train',x_test',ker,beta,bias); % 测试样本数据
plot(x_train,y_train,'g:',x_test,y_test,'b-',x_test,y_test1 ,'r-')
title('绿色——样本数据,蓝色——实际数据,红色——测试数据')
y_test1
grid on
[ 本帖最后由 lzh19840809 于 2009-5-10 19:30 编辑 ] |
|