基于svm人脸识别代码

基于svm人脸识别代码

ID:31264424

大小:59.27 KB

页数:8页

时间:2019-01-07

基于svm人脸识别代码_第1页
基于svm人脸识别代码_第2页
基于svm人脸识别代码_第3页
基于svm人脸识别代码_第4页
基于svm人脸识别代码_第5页
资源描述:

《基于svm人脸识别代码》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、基于SVM人脸识别Matlab代码function[nsv,al3,bi,T]=svml68(x,y,ker,C,par,tol);%支持向量机在模式识别方面的应用.%[nsv,al3,bi,T]=svml68(x,y,ker,C,par,tol);%输入参数:%x:输入矩阵,每行代表一个样本,列数表示样本特征数。矩阵行数代表样本%Y:由-1或+1组成的列向量.每行表示相应样本的类别。只有两类。%ker:字符串'linear','poly_h','poly_i'and'rbf'。%用于传递内积函数的类型。%C:惩罚因子。%par:内积函数的参数。

2、如多项式内积表示多项式的次数。查看kernel函数。%tol:算法的迭代误差。默认值是10=5。%输出:%nsv:支持向量个数.%al3:alpha值.%bi:偏差值.%T:运行时间.%clear;clc;closeall:%dispC输入样本矩阵(每行表示一个测试数据):’);%x=[randn(l,20)-3randn(1,20)+3:randn(1,40)]J%dispC样本所属类别:');%y二[-ones(1,20)ones(1,20)]'%[nsv,alpha,bias,T]=svml68(x,y,'linear',10);%fpri

3、ntf('运行时间T=%g,,T);%fprintfC支持向量个数:%g,nsv);%svcplot(x,y,'linear,,alpha,bias);%W=x'*alpha;%disp(,权系数向量:’);%W%dispC测试数据:');%X_test=[randn(1,20)-3randn(1,20)+3;randn(1,40)]'%dispC测试数据类别:');%Y_test=[-ones(l,20)ones(1,20)]'%dispC测试数据类别判决:’);%Output二X_test*W+bias%dispC错误率:');%nu

4、mber_error=sum(abs(Y_test-sign(Output)))/2;%error=number_error/20%径向基(或多项式)内积函数:(testrbf.m)%clear;clc;closeall;%disp('输入样本矩阵(每行表示一个测试数据):’);%x=[randn(l,20)-3randn(1,20)+3;randn(1,40)]*%dispC样本所属类别:');%y=[-ones(1,20)ones(1,20)]'%[nsv,alpha,bias,T]=svml68(x,y,'rbf',10,4);%fprin

5、tf('运行时间T=%g,,T);%fprintf('支持向量个数:%g,,nsv);%svcplot(x,y,'linear,,alpha,bias);%dispC测试数据:');%X_test=[randn(1,20)-3randn(1,20)+3;randn(1,40)]'%dispC测试数据类别:');%Y_test=[-ones(1,20)ones仃,20)]'%dispC测试数据类别判决:');Output二kernel('rbf,,X_test,x(find(alpha),:),4)*alpha(find(alpha))+b

6、ias%dispC错误率:');%number_error=sum(abs(Y_test-sign(Output)))/2;%error=number_error/20T=clock;if(nargin<4

7、nargin>6)helpsvml68elseif(nargin==4&ker(l)='p')dispC将多项式内积的次数置为2。');par=2;・3f,,sqrt(size(x,2)));elseif(nargin==4&ker(1)==,x)fprintf(l,'径向基函数内积的参数(方差)置为par=sqrt(size(x,2))

8、;elseif(nargin==4)par=O;endif(nargin<6)tol=10'-5;endK=1(T6;N=size(x,1);ns=-l;hacer=l;bi=O;H=kernel(ker,x,x,par);ilp=l:2:N;iln=2:2:N;il=[ilpilnj,;a=zeros(N,1);a(il)=C;al3=zeros(N,1);bi=O;i2p=[];i2n=[]:Lp⑴二N*C;e=y;while(hacer)al3_a=al3;bi_a=bi;al3=zeros(N,1);al3(i2p)二C;al3(i2n)

9、=-C;if(1ength(il))Xi=inv([H(il,il)+diag(l・/(a(订)))ones(length(il),1)

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。