randdata.m

本程序是基于船舶电力回归的SVM程序

clear all;clc;
%p=[3,4,5,6,7,8,9,1,2,4];
%a=[20,15.2,5,45,56,78,9,69,63,98,78,99];
load('lesspower.mat');load('lessprobability.mat');
p=p/10;
m=100,v=size(a,2);
fprintf('load rand data  ....\n')
fprintf('deal with data  ...\n');
    %st = cputime;
x=rand(m,v);e=0.001;ee=0.0001;
i=1:v;
t=max(x(:,i));
d=ones(m,v);
for n=1:v
    data=t(n);
d(:,n)=data.*d(:,n);
end
for sn=1:v
    r=find((d(:,sn)-x(:,sn))<e);%re=find((d-x)>e)
    rr=size(r,1);
while rr<p(sn)-1
    e=e+ee;
    r=find((d(:,sn)-x(:,sn))<e);
    rr=size(r,1);
end
s=size(r,1)-p(sn)+1;
if s>1
    e=e/10;
   ee=ee/10;
r=find((d(:,sn)-x(:,sn))<e);
rr=size(x,1);
while rr<p(sn)-1
    e=e+ee;
    r=find((d(:,sn)-x(:,sn))<e);
    rr=size(r,1);
end
end
%if size(rr,1)<size(r,1)
    %x(rr,sn)=1;
    %else 
    x(r,sn)=1;
    %end
%clear rr;clear r;
end
re=find(x<1);
x(re)=0;X=[];
for n=1:v
    X(:,n)=x(:,n).*a(n);
    %Y=Y+X(:,n);
end
Y=zeros(size(X,1),1);
for n=1:v
    Y=Y+X(:,n);
end
%fprintf('Execution time : %4.1f seconds\n',cputime - st);
k=1:m;
plot (k,Y(k),'r')
fprintf('Support vector machine regress  .... \n')
[nsv,beta,bias,error2,trnY]=powerregress(X,Y) ;