基于SVM人臉識(shí)別代碼(共7頁(yè))

1、精選優(yōu)質(zhì)文檔-傾情為你奉上基于SVM人臉識(shí)別Matlab代碼function nsv,al3,bi,T=svm168(x,y,ker,C,par,tol);% 支持向量機(jī)在模式識(shí)別方面的應(yīng)用.% nsv,al3,bi,T=svm168(x,y,ker,C,par,tol);% 輸入?yún)?shù):% x: 輸入矩陣，每行代表一個(gè)樣本，列數(shù)表示樣本特征數(shù)。矩陣行數(shù)代表樣本數(shù) 。% y: 由 -1 或 +1組成的列向量. 每行表示相應(yīng)樣本的類別。只有兩類。% ker: 字符串 linear, poly_h, poly_i and rbf。% 用于傳遞內(nèi)積函數(shù)的類型。% C: 懲罰因子。% par: 內(nèi)積函數(shù)

2、的參數(shù)。如多項(xiàng)式內(nèi)積表示多項(xiàng)式的次數(shù)。查看kernel函數(shù)。% tol: 算法的迭代誤差。默認(rèn)值是 10-5。% 輸出:% nsv: 支持向量個(gè)數(shù).% al3: alpha值.% bi: 偏差值.% T: 運(yùn)行時(shí)間.% clear;clc;close all;% disp(輸入樣本矩陣(每行表示一個(gè)測(cè)試數(shù)據(jù)):);% x=randn(1,20)-3 randn(1,20)+3;randn(1,40)% disp(樣本所屬類別:);% y=-ones(1,20) ones(1,20)% nsv,alpha,bias,T=svm168(x,y,linear,10);% fprintf(運(yùn)行時(shí)間T=

3、%gn,T);% fprintf(支持向量個(gè)數(shù)：%gn,nsv);% svcplot(x,y,linear,alpha,bias);% W=x*alpha;% disp(權(quán)系數(shù)向量:);% W% disp(測(cè)試數(shù)據(jù):);% X_test=randn(1,20)-3 randn(1,20)+3;randn(1,40)% disp(測(cè)試數(shù)據(jù)類別:);% Y_test=-ones(1,20) ones(1,20)% disp(測(cè)試數(shù)據(jù)類別判決:);% Output=X_test*W+bias% disp(錯(cuò)誤率:);% number_error=sum(abs(Y_test-sign(Output)

4、/2;% error=number_error/20% 徑向基(或多項(xiàng)式)內(nèi)積函數(shù):(testrbf.m)% clear;clc;close all;% disp(輸入樣本矩陣(每行表示一個(gè)測(cè)試數(shù)據(jù)):);% x=randn(1,20)-3 randn(1,20)+3;randn(1,40)% disp(樣本所屬類別:);% y=-ones(1,20) ones(1,20)% nsv,alpha,bias,T=svm168(x,y,rbf,10,4);% fprintf(運(yùn)行時(shí)間T=%gn,T);% fprintf(支持向量個(gè)數(shù)：%gn,nsv);% svcplot(x,y,linear,al

5、pha,bias);% disp(測(cè)試數(shù)據(jù):);% X_test=randn(1,20)-3 randn(1,20)+3;randn(1,40)% disp(測(cè)試數(shù)據(jù)類別:);% Y_test=-ones(1,20) ones(1,20)% disp(測(cè)試數(shù)據(jù)類別判決:);Output=kernel(rbf,X_test,x(find(alpha),:),4)*alpha(find(alpha)+bias% disp(錯(cuò)誤率:);% number_error=sum(abs(Y_test-sign(Output)/2;% error=number_error/20-T=clock;if (na

6、rgin6)help svm168elseif(nargin=4 & ker(1)=p)disp(將多項(xiàng)式內(nèi)積的次數(shù)置為2。);par=2;elseif(nargin=4 & ker(1)=r)fprintf(1,徑向基函數(shù)內(nèi)積的參數(shù)（方差）置為 %1.3fn,sqrt(size(x,2);par=sqrt(size(x,2);elseif(nargin=4)par=0;endif(nargin0).*y(find(e.*y0);if(abs(Lp(k)-Lp(k-1)/Lp(k-1)(ns+1)nsv=length(i1)+length(i2n)+length(i2p);hacer=0;el

7、seif(Lp(k)-Lp(k-1)/Lp(k-1)tol)LL=Lp(k-1);HH=Lp(k);posL=0;posH=1;while(posL=0)jj=(posH+posL)/2;al3b=al3*jj+al3_a*(1-jj);bib=bi*jj+bi_a*(1-jj);I=find(al3b-al3_a);if(length(I)eb=e_a-H(:,I)*(al3b(I)-al3_a(I)-(bib-bi_a);endif(HHLL)HHa=HH;posHa=posH;HH=al3b*H*al3b/2+C*sum(eb(find(eb.*y0).*y(find(eb.*y0);p

8、osH=jj;elseLL=al3b*H*al3b/2+C*sum(eb(find(eb.*y0).*y(find(eb.*y0);posL=jj;endendjj=posH+posL;al3b=al3*jj+al3_a*(1-jj);bib=bi*jj+bi_a*(1-jj);I=find(al3b-al3_a);if(length(I)eb=e_a-H(:,I)*(al3b(I)-al3_a(I)-(bib-bi_a);endvalor=al3b*H*al3b/2+C*sum(eb(find(eb.*y0).*y(find(eb.*y0);posL=posH;LL=HH;posH=jj;H

9、H=valor;if(valorLp(k-1)while(abs(valor-Lp(k-1)/Lp(k-1)10-3)jj=(posH+posL)/2;al3b=al3*jj+al3_a*(1-jj);bib=bi*jj+bi_a*(1-jj);I=find(al3b-al3_a);if(length(I)eb=e_a-H(:,I)*(al3b(I)-al3_a(I)-(bib-bi_a);endvalor=al3b*H*al3b/2+C*sum(eb(find(eb.*y0).*y(find(eb.*y0);if(valorLp(k-1)HH=valor;posH=jj;elseLL=val

10、or;posL=jj;endendendpos=jj;al3=al3*pos+al3_a*(1-pos);bi=bi*pos+bi_a*(1-pos);I=find(al3-al3_a);if(length(I)e=e_a-H(:,I)*(al3(I)-al3_a(I)-(bi-bi_a);endLp(k)=al3*H*al3/2+C*sum(e(find(e.*y0).*y(find(e.*y0);if(pos=0 & y=1);i1n=find(eK)=K;i2p=find(al3=(C-C/100) & al3=(C-C/100) & al3_a0 & y=1 & 15*e14*e_a);if(length(setdiff(i1p,i2p)i1p=setdiff(i1p,i2p);elseval,pos=min(e(i2p);i2p=setdiff(i2p,i2p(pos);i1p=setdiff(i1p,i2p);endi2n=find(al3=(-C-C/100) & al3=(-C-C/100) & al3_a=(-C+C/100) & e0 & y=-1 & 15*e14*e_a);if