模式識(shí)別第二次作業(yè)

模式識(shí)別第二次作業(yè)模式識(shí)別第二次作業(yè)模式識(shí)別第二次作業(yè)模式識(shí)別第二次作業(yè)編制僅供參考審核批準(zhǔn)生效日期地址：電話：傳真:郵編:模式識(shí)別第二次作業(yè)1、調(diào)用函數(shù)(matlab工具箱函數(shù))替代svm_matlab中的fmincon函數(shù)。seqminopt函數(shù)的調(diào)用方法如下形式：[ALPHAOFFSET]=seqminopt(TRAINING,Y,BOXC,KERNELFUNC,SMOOPTIONS)輸入輸出各變量含義如下：TRAINING待訓(xùn)練數(shù)據(jù)Y列向量，代表響應(yīng)數(shù)據(jù)，為±1BOXC松弛變量，列向量KERNELFUNC句柄SMOOPTIONS可選，用于設(shè)置迭代次數(shù)等ALPHA優(yōu)化得到的αOFFSET超平面偏移使用fmincon函數(shù)結(jié)果為：W_fmincon=用seqminopt替換fmincon后結(jié)果為：W_seqminopt=可以發(fā)現(xiàn)兩種形式下結(jié)果是一樣的。設(shè)置100個(gè)離散點(diǎn)，分類結(jié)果為：程序closeall;clearall;clc%生成訓(xùn)練數(shù)據(jù)=100;=1*[3;5];=randn(2,+kron,ones(1,);=50;=randn(2,+kron,ones(1,);=100;=-1*[3;5];=randn(2,+kron,ones(1,);=50;=randn(2,+kron,ones(1,);TestMatrix=zeros(2);cnt_c1=1;cnt_c2=1;tt=1:1000;tt=(tt-500)*;figure;plot(1,:),(2,:),'bo','linewidth',2,'markersize',5)holdon;plot(1,:),(2,:),'ko','linewidth',2,'markersize',5)%要訓(xùn)練的數(shù)據(jù)及其類別cnt=1;X=zeros+,3);foriii=1:X(cnt,1:2)=(:,iii);X(cnt,3)=1;y(cnt,1)=1;cnt=cnt+1;endforiii=1:X(cnt,1:2)=(:,iii);X(cnt,3)=1;y(cnt,1)=-1;cnt=cnt+1;end%替換為seqminopt后不用下面這個(gè)cnt=1;Xsvm=zeros+,3);foriii=1:Xsvm(cnt,1:2)=-1*(:,iii);Xsvm(cnt,3)=-1*1;ysvm(cnt,1)=-1;cnt=cnt+1;endforiii=1:Xsvm(cnt,1:2)=(:,iii);Xsvm(cnt,3)=1;ysvm(cnt,1)=-1;cnt=cnt+1;end%%%%%%%%%%%%%%%%%%%%%%%%%%下面就是替換為seqminopttraindata=X(:,1:2);%用X而不是Xsvmtrainlabel=y;%用y而非ysvmboxConstraints=ones(size(trainlabel));%松弛變量，我默認(rèn)為1kfun=@linear_kernel;%線性函數(shù)句柄[alphas,bias]=seqminopt(traindata,trainlabel,boxConstraints,kfun);svIndex=find(alphas>sqrt(eps));%尋找支持向量的索引sv=traindata(svIndex,:);%對(duì)應(yīng)x(i)alphaHat=trainlabel(svIndex).*alphas(svIndex);%對(duì)應(yīng)y(i)*α(i)weight=zeros(1,2);%W=sum(y(i)*α(i)*x(i))fori=1:length(svIndex)weight=weight+alphaHat(i)*sv(i,:);endW=zeros(3,1);W(1:2)=weight;W(3)=bias;W_seqminopt=W%用來(lái)做對(duì)比%%%%%%%%%%%%%%%%%%%%%%fmincon_options=optimset('Algorithm','interior-point',...'MaxIter',50);[W1,val,exflag,output]=fmincon(@(x)(norm(x(1:2),2)),[0;0;1],...Xsvm,ysvm,[],[],[],[],[],fmincon_options);W_fmincon=W1%做對(duì)比%%%%%%%%%%%%%%%%%%%%%%%%%forkkk=1:TMP_Tr(1:2,1)=(:,kkk);TMP_Tr(3,1)=1;yy=W.'*TMP_Tr;if(yy>0)Result_c1(:,cnt_c1)=(:,kkk);cnt_c1=cnt_c1+1;TestMatrix(1,1)=TestMatrix(1,1)+1;elseResult_c2(:,cnt_c2)=(:,kkk);cnt_c2=cnt_c2+1;TestMatrix(1,2)=TestMatrix(1,2)+1;endendforkkk=1:TMP_Tr(1:2,1)=(:,kkk);TMP_Tr(3,1)=1;yy=W.'*TMP_Tr;if(yy>0)Result_c1(:,cnt_c1)=(:,kkk);cnt_c1=cnt_c1+1;TestMatrix(2,1)=TestMatrix(2,1)+1;elseResult_c2(:,cnt_c2)=(:,kkk);cnt_c2=cnt_c2+1;TestMatrix(2,2)=TestMatrix(2,2)+1;endendTestMatrix(1,:)=TestMatrix(1,:)/;TestMatrix(2,:)=TestMatrix(2,:)/;TestMatrixfigure;tryplot(Result_c1(1,:),Result_c1(2,:),'bo','linewidth',2,'markersize',5)endholdontryplot(Result_c2(1,:),Result_c2(2,:),'ko','linewidth',2,'markersize',5)end%plot(1,:),(2,:),'r+','linewidth',2,'markersize',10)%%pl