www.gusucode.com > LSSVMlabv1_8_R2009b_R2011a源码 > LSSVMlabv1_8_R2009b_R2011a源码/LSSVMlabv1_8_R2009b_R2011a/democonfint.m
% Copyright (c) 2011, KULeuven-ESAT-SCD, License & help @ http://www.esat.kuleuven.be/sista/lssvmlab clc; disp('This is a simple demo, solving a simple regression task using'); disp('LS-SVMlab and constructing confidence intervals. A dataset is constructed in the right formatting. The'); disp('data are represented as matrices where each row contains one'); disp('datapoint: '); disp(' '); disp('press <ENTER> key'); pause disp(' '); disp('>> X = (-3:0.02:3)'';'); X = (-3:0.02:3)'; disp('>> Y = sinc(X)+0.1.*randn(length(X),1);'); eval('Y = sinc(X)+0.1.*randn(length(X),1);',... 'Y = sin(pi.*X+12345*eps)./(pi*X+12345*eps)+0.1.*randn(length(X),1);'); disp('>> X'); X disp('>> Y'); Y disp('In order to make an LS-SVM model, we need 2 extra parameters: gamma'); disp('(gam) is the regularization parameter, determining the trade-off'); disp('between the fitting error minimization and smoothness of the'); disp('estimated function. sigma^2 (sig2) is the kernel function'); disp('parameter of the RBF kernel. These can be found via cross-validation:'); disp(' '); model = initlssvm(X,Y,'f',[],[],'RBF_kernel','o'); disp('>> model = tunelssvm(model,''simplex'',''crossvalidatelssvm'',{10,''mse''});'); model = tunelssvm(model,'simplex','crossvalidatelssvm',{10,'mse'}); disp(' '); disp('press <ENTER> key'); pause disp(' '); disp('Training the model '); disp(' '); disp('>> model = trainlssvm(model)'); model = trainlssvm(model); disp(' '); disp('Computation of Confidence Intervals '); disp(' '); disp('press <ENTER> key'); pause disp('ci = cilssvm(model);') ci = cilssvm(model); disp('The LS-SVM result and confidence intervals can be displayed if the dimension of the input'); disp('data is 1 or 2. '); disp(' '); disp('>> plotlssvm(model);'); figure; plotlssvm(model); disp(' '); hold all fill([X;flipud(X)],[ci(:,1);flipud(ci(:,2))],'c','FaceAlpha',0.5,'EdgeAlpha',1,'EdgeColor','w') disp('All plotting is done with this simple command. It looks for the'); disp('best way of displaying the result.') disp(' '); disp(' This concludes the demo');