- 余辉
-
>>edit svmtrain
>>edit svmclassify
>>edit svmpredict
function [svm_struct, svIndex] = svmtrain(training, groupnames, varargin)
%SVMTRAIN trains a support vector machine classifier
%
% SVMStruct = SVMTRAIN(TRAINING,GROUP) trains a support vector machine
% classifier using data TRAINING taken from two groups given by GROUP.
% SVMStruct contains information about the trained classifier that is
% used by SVMCLASSIFY for classification. GROUP is a column vector of
% values of the same length as TRAINING that defines two groups. Each
% element of GROUP specifies the group the corresponding row of TRAINING
% belongs to. GROUP can be a numeric vector, a string array, or a cell
% array of strings. SVMTRAIN treats NaNs or empty strings in GROUP as
% missing values and ignores the corresponding rows of TRAINING.
%
% SVMTRAIN(...,"KERNEL_FUNCTION",KFUN) allows you to specify the kernel
% function KFUN used to map the training data into kernel space. The
% default kernel function is the dot product. KFUN can be one of the
% following strings or a function handle:
%
% "linear" Linear kernel or dot product
% "quadratic" Quadratic kernel
% "polynomial" Polynomial kernel (default order 3)
% "rbf" Gaussian Radial Basis Function kernel
% "mlp" Multilayer Perceptron kernel (default scale 1)
% function A kernel function specified using @,
% for example @KFUN, or an anonymous function
%
% A kernel function must be of the form
%
% function K = KFUN(U, V)
%
% The returned value, K, is a matrix of size M-by-N, where U and V have M
% and N rows respectively. If KFUN is parameterized, you can use
% anonymous functions to capture the problem-dependent parameters. For
% example, suppose that your kernel function is
%
% function k = kfun(u,v,p1,p2)
% k = tanh(p1*(u*v")+p2);
%
% You can set values for p1 and p2 and then use an anonymous function:
% @(u,v) kfun(u,v,p1,p2).
%
% SVMTRAIN(...,"POLYORDER",ORDER) allows you to specify the order of a
% polynomial kernel. The default order is 3.
%
% SVMTRAIN(...,"MLP_PARAMS",[P1 P2]) allows you to specify the
% parameters of the Multilayer Perceptron (mlp) kernel. The mlp kernel
% requires two parameters, P1 and P2, where K = tanh(P1*U*V" + P2) and P1
% > 0 and P2 < 0. Default values are P1 = 1 and P2 = -1.
%
% SVMTRAIN(...,"METHOD",METHOD) allows you to specify the method used
% to find the separating hyperplane. Options are
%
% "QP" Use quadratic programming (requires the Optimization Toolbox)
% "LS" Use least-squares method
%
% If you have the Optimization Toolbox, then the QP method is the default
% method. If not, the only available method is LS.
%
% SVMTRAIN(...,"QUADPROG_OPTS",OPTIONS) allows you to pass an OPTIONS
% structure created using OPTIMSET to the QUADPROG function when using
% the "QP" method. See help optimset for more details.
%
% SVMTRAIN(...,"SHOWPLOT",true), when used with two-dimensional data,
% creates a plot of the grouped data and plots the separating line for
% the classifier.
%
% Example:
% % Load the data and select features for classification
% load fisheriris
% data = [meas(:,1), meas(:,2)];
% % Extract the Setosa class
% groups = ismember(species,"setosa");
% % Randomly select training and test sets
% [train, test] = crossvalind("holdOut",groups);
% cp = classperf(groups);
% % Use a linear support vector machine classifier
% svmStruct = svmtrain(data(train,:),groups(train),"showplot",true);
% classes = svmclassify(svmStruct,data(test,:),"showplot",true);
% % See how well the classifier performed
% classperf(cp,classes,test);
% cp.CorrectRate
%
% See also CLASSIFY, KNNCLASSIFY, QUADPROG, SVMCLASSIFY.
% Copyright 2004 The MathWorks, Inc.
% $Revision: 1.1.12.1 $ $Date: 2004/12/24 20:43:35 $
% References:
% [1] Kecman, V, Learning and Soft Computing,
% MIT Press, Cambridge, MA. 2001.
% [2] Suykens, J.A.K., Van Gestel, T., De Brabanter, J., De Moor, B.,
% Vandewalle, J., Least Squares Support Vector Machines,
% World Scientific, Singapore, 2002.
% [3] Scholkopf, B., Smola, A.J., Learning with Kernels,
% MIT Press, Cambridge, MA. 2002.
%
% SVMTRAIN(...,"KFUNARGS",ARGS) allows you to pass additional
% arguments to kernel functions.
% set defaults
plotflag = false;
qp_opts = [];
kfunargs = {};
setPoly = false; usePoly = false;
setMLP = false; useMLP = false;
if ~isempty(which("quadprog"))
useQuadprog = true;
else
useQuadprog = false;
end
% set default kernel function
kfun = @linear_kernel;
% check inputs
if nargin < 2
error(nargchk(2,Inf,nargin))
end
numoptargs = nargin -2;
optargs = varargin;
% grp2idx sorts a numeric grouping var ascending, and a string grouping
% var by order of first occurrence
[g,groupString] = grp2idx(groupnames);
% check group is a vector -- though char input is special...
if ~isvector(groupnames) && ~ischar(groupnames)
error("Bioinfo:svmtrain:GroupNotVector",...
"Group must be a vector.");
end
% make sure that the data is correctly oriented.
if size(groupnames,1) == 1
groupnames = groupnames";
end
% make sure data is the right size
n = length(groupnames);
if size(training,1) ~= n
if size(training,2) == n
training = training";
else
error("Bioinfo:svmtrain:DataGroupSizeMismatch",...
"GROUP and TRAINING must have the same number of rows.")
end
end
% NaNs are treated as unknown classes and are removed from the training
% data
nans = find(isnan(g));
if length(nans) > 0
training(nans,:) = [];
g(nans) = [];
end
ngroups = length(groupString);
if ngroups > 2
error("Bioinfo:svmtrain:TooManyGroups",...
"SVMTRAIN only supports classification into two groups. GROUP contains %d different groups.",ngroups)
end
% convert to 1, -1.
g = 1 - (2* (g-1));
% handle optional arguments
if numoptargs >= 1
if rem(numoptargs,2)== 1
error("Bioinfo:svmtrain:IncorrectNumberOfArguments",...
"Incorrect number of arguments to %s.",mfilename);
end
okargs = {"kernel_function","method","showplot","kfunargs","quadprog_opts","polyorder","mlp_params"};
for j=1:2:numoptargs
pname = optargs{j};
pval = optargs{j+1};
k = strmatch(lower(pname), okargs);%#ok
if isempty(k)
error("Bioinfo:svmtrain:UnknownParameterName",...
"Unknown parameter name: %s.",pname);
elseif length(k)>1
error("Bioinfo:svmtrain:AmbiguousParameterName",...
"Ambiguous parameter name: %s.",pname);
else
switch(k)
case 1 % kernel_function
if ischar(pval)
okfuns = {"linear","quadratic",...
"radial","rbf","polynomial","mlp"};
funNum = strmatch(lower(pval), okfuns);%#ok
if isempty(funNum)
funNum = 0;
end
switch funNum %maybe make this less strict in the future
case 1
kfun = @linear_kernel;
case 2
kfun = @quadratic_kernel;
case {3,4}
kfun = @rbf_kernel;
case 5
kfun = @poly_kernel;
usePoly = true;
case 6
kfun = @mlp_kernel;
useMLP = true;
otherwise
error("Bioinfo:svmtrain:UnknownKernelFunction",...
"Unknown Kernel Function %s.",kfun);
end
elseif isa (pval, "function_handle")
kfun = pval;
else
error("Bioinfo:svmtrain:BadKernelFunction",...
"The kernel function input does not appear to be a function handle or valid function name.")
end
case 2 % method
if strncmpi(pval,"qp",2)
useQuadprog = true;
if isempty(which("quadprog"))
warning("Bioinfo:svmtrain:NoOptim",...
"The Optimization Toolbox is required to use the quadratic programming method.")
useQuadprog = false;
end
elseif strncmpi(pval,"ls",2)
useQuadprog = false;
else
error("Bioinfo:svmtrain:UnknownMethod",...
"Unknown method option %s. Valid methods are ""QP"" and ""LS""",pval);
end
case 3 % display
if pval ~= 0
if size(training,2) == 2
plotflag = true;
else
warning("Bioinfo:svmtrain:OnlyPlot2D",...
"The display option can only plot 2D training data.")
end
end
case 4 % kfunargs
if iscell(pval)
kfunargs = pval;
else
kfunargs = {pval};
end
case 5 % quadprog_opts
if isstruct(pval)
qp_opts = pval;
elseif iscell(pval)
qp_opts = optimset(pval{:});
else
error("Bioinfo:svmtrain:BadQuadprogOpts",...
"QUADPROG_OPTS must be an opts structure.");
end
case 6 % polyorder
if ~isscalar(pval) || ~isnumeric(pval)
error("Bioinfo:svmtrain:BadPolyOrder",...
"POLYORDER must be a scalar value.");
end
if pval ~=floor(pval) || pval < 1
error("Bioinfo:svmtrain:PolyOrderNotInt",...
"The order of the polynomial kernel must be a positive integer.")
end
kfunargs = {pval};
setPoly = true;
case 7 % mlpparams
if numel(pval)~=2
error("Bioinfo:svmtrain:BadMLPParams",...
"MLP_PARAMS must be a two element array.");
end
if ~isscalar(pval(1)) || ~isscalar(pval(2))
error("Bioinfo:svmtrain:MLPParamsNotScalar",...
"The parameters of the multi-layer perceptron kernel must be scalar.");
end
kfunargs = {pval(1),pval(2)};
setMLP = true;
end
end
end
end
if setPoly && ~usePoly
warning("Bioinfo:svmtrain:PolyOrderNotPolyKernel",...
"You specified a polynomial order but not a polynomial kernel");
end
if setMLP && ~useMLP
warning("Bioinfo:svmtrain:MLPParamNotMLPKernel",...
"You specified MLP parameters but not an MLP kernel");
end
% plot the data if requested
if plotflag
[hAxis,hLines] = svmplotdata(training,g);
legend(hLines,cellstr(groupString));
end
% calculate kernel function
try
kx = feval(kfun,training,training,kfunargs{:});
% ensure function is symmetric
kx = (kx+kx")/2;
catch
error("Bioinfo:svmtrain:UnknownKernelFunction",...
"Error calculating the kernel function: %s ", lasterr);
end
% create Hessian
% add small constant eye to force stability
H =((g*g").*kx) + sqrt(eps(class(training)))*eye(n);
if useQuadprog
% The large scale solver cannot handle this type of problem, so turn it
% off.
qp_opts = optimset(qp_opts,"LargeScale","Off");
% X=QUADPROG(H,f,A,b,Aeq,beq,LB,UB,X0,opts)
alpha = quadprog(H,-ones(n,1),[],[],...
g",0,zeros(n,1),inf *ones(n,1),zeros(n,1),qp_opts);
% The support vectors are the non-zeros of alpha
svIndex = find(alpha > sqrt(eps));
sv = training(svIndex,:);
% calculate the para- 苏州马小云
-
你的emd是%
g.
rilling,的那个吗?
len
=
size(imf,1);
for
k
=
1:len
len1
=
length(imf(k,:));
b(k)
=
sum(imf(k,:).*imf(k,:))/len1;%
时域均方值,能量
amp(k,:)
=
abs(imf(k,:));
b(k)
=
sqrt(b(k));
th
=
angle(hilbert(imf(k,:)));%hilbert变换的相位
d(k,:)
=
diff(th)/ts/(2*pi);%求导,得到频率:f
=
(1/2*pi)*d(th)/dt
end
你的频率公式用得有点问题,求出来不应是归一化频率
- 瑞瑞爱吃桃
-
能否把代码也发我一份,我想把这个方法用在地球化学方面
我邮箱258275248@q.com
- coco
-
求楼主解答 邮箱362382447@q.com 感谢万分
- 阿啵呲嘚
-
哥们 我现在在做这课程设计 能不能把你的给我看看一下了 727674595@q.com 谢谢了
- LocCloud
-
是啊,我也正在找EMD的代码。找到了,一定会告诉你的。