《模式识别》课程实验线性分类器设计实验

1、模式识别课程实验线性分类器设计实验一、实验目的：1掌握线性分类器设计方法；2、掌握感知准则函数分类器设计方法。二、实验内容：1对下列两种情况，求采用判决准则时的投影向量和分类界面，并做图。1(2,0)T,(2,2)T,(2,4)T,(3,3)T(0,3)t,(B2,2)t,(,)t,(1,(3,)t21(1,1)T,(2,0)T,(2,1)T,(0,2)T,(1,3)T(,2)t,(0,0)t,(，0),(,)t,(0,亶)T22、对下面的两类分类问题，采用感知准则函数，利用迭代修正求权向量的方法求两类的线性判决函数及线性识别界面，并画出识别界面将训练样本区分的结果图。1(1,1)T,(2,0

2、)T,(2,1)T,(0,2)T,(1,3)T(,2)t,(0,0)t,(，0),(,)t,(0,J)T2三、实验原理：()判决准则投影方向：w*()w12(2)感知准则函数：J(v)(ikz)pz!Zk当Zk为空时，即J(v)0，v*即为所求p四、解题思路:i线性判决器：用函数求两类样本的均值求两类样本的均值的类内离散矩阵利用类内离散矩阵求总类内离散矩阵求最佳投影方向定义阈值，并求得分界面2感知准则函数分类器：获得增广样本向量和初始增广权向量对样本进行规范化处理获得解区，并用权向量迭代修正错分样本集，得到最终解区五、实验结果：1线性判决分类器：条件：取阈值系数为第一种情况第.二种情况Fish

3、er性强类-器西喘示第1类表示第2类2、感知准则函数判决条件：取步长row为1惑気口彳隹贝II竽IJ决-1判决结果：六、结果分析：线性判决器中，调整阈值系数时，分界面会随之平行上下移动，通过调整阈值系数的大小，就能比较合理的得到分界面。而投影方向不会随阈值系数的变化而变化。线性判决器中，由于题目中样本数目少，所以其先验概率统一取为。而改变其先验概率对结果的影响不是很大。3感知准则函数分类判决中，调整步长的大小可使分界面的方向和位置发生变动，但是变动范围不大，仅限于小范围的调整。七、源程序functionvarargout=fisherganzhi(varargin)gui_Singleton=

4、1;gui_State=struct(gui_Name,gui_Singleton,gui_OpeningFcngui_OutputFcn,gui_LayoutFcn,gui_Callback,mfilename,.gui_Singleton,.,fisherganzhi_OpeningFcn,.fisherganzhi_OutputFcn,.,.);ifnargin&ischar(varargin1)gui_State.gui_Callback=str2func(varargin1);endifnargoutvarargout1:nargout=gui_mainfcn(gui_State,v

5、arargin:);elseendgui_mainfcn(gui_State,varargin:);%-Executesjustbeforefisherganzhiismadevisible.functionfisherganzhi_OpeningFcn(hObject,eventdata,handles,varargin)handles.output=hObject;%Updatehandlesstructureguidata(hObject,handles);%-Outputsfromthisfunctionarereturnedtothecommandline.functionvarar

6、gout=fisherganzhi_OutputFcn(hObject,eventdata,handles)varargout1=handles.output;%-Executesonbuttonpressinpushbutton2.functionpushbutton2_Callback(hObject,eventdata,handles)pw1=str2num(get(handles.edit1,String);pw2=str2num(get(handles.edit2,String);k=str2num(get(handles.edit3,String);w1=11;20;21;02;1

7、3;w2=-12;00;-10;-1-1;0-2;u1=mean(w1);u2=mean(w2);s1=zeros(2,2);s2=zeros(2,2);forn=1:length(w1)s1=s1+(w1(n,:)-u1)*(w1(n,:)-u1);end;form=1:length(w2)s2=s2+(w2(m,:)-u2)*(w2(m,:)-u2);end;sw=pw1*s1+pw2*s2;wo=inv(sw)*(u1-u2);y0=k*(wo*u1+wo*u2);x1=-4:0.3:6;x2=(y0-wo(1)*x1)/wo(2);x3=wo(2)*x1/wo(1);cla(handl

8、es.axes1);axes(handles.axes1);holdon;plot(x1,x2,b,x1,x3,k);plot(w1(:,2),w1(:,1),*r);plot(w2(:,2),w2(:,1),og);title(Fisher线性分类器(2);12);handles)text(3.5,13,*表示第1类,FontSize,10,Color,rtext(3.5,11,o表示第2类,FontSize,10,Color,gtext(wo(1),wo(2),leftarrow投影方向,FontSizeholdoff;%-Executesonbuttonpressinpushbutton

9、1.functionpushbutton1_Callback(hObject,eventdata,pw1=str2num(get(handles.edit1,String);pw2=str2num(get(handles.edit2,String);k=str2num(get(handles.edit3,String);w1=20;22;24;33;w2=03;-22;-1-1;1-2;3-1;u1=mean(w1);u2=mean(w2);s1=zeros(2,2);s2=zeros(2,2);forn=1:length(w1)s1=s1+(w1(n,:)-u1)*(w1(n,:)-u1);

10、end;form=1:length(w2)s2=s2+(w2(m,:)-u2)*(w2(m,:)-u2);end;sw=pw1*s1+pw2*s2;wo=inv(sw)*(u1-u2);y0=k*(wo*u1+wo*u2);x1=-4:0.3:6;x2=(y0-wo(1)*x1)/wo(2);x3=wo(2)*x1/wo(1);cla(handles.axes1);axes(handles.axes1);holdon;plot(x1,x2,b,x1,x3,k);plot(w1(:,2),w1(:,1),*r);plot(w2(:,2),w2(:,1),og);title(Fisher线性分类器

11、(2);text(3.5,9,*表示第1类,FontSize,10,Color,r);text(3.5,8,o表示第2类,FontSize,10,Color,g);text(wo(1),wo(2),leftarrow投影方向,FontSize,12);holdoff;functionedit3_Callback(hObject,eventdata,handles)%-Executesduringobjectcreation,aftersettingallproperties.functionedit3_CreateFcn(hObject,eventdata,handles)ifispc&ise

12、qual(get(hObject,BackgroundColor),get(0,defaultUicontrolBackgroundColor)set(hObject,BackgroundColor,white);endfunctionedit2_Callback(hObject,eventdata,handles)%-Executesduringobjectcreation,aftersettingallproperties.functionedit2_CreateFcn(hObject,eventdata,handles)ifispc&isequal(get(hObject,Backgro

13、undColor),get(0,defaultUicontrolBackgroundColor)set(hObject,BackgroundColor,white);endfunctionedit1_Callback(hObject,eventdata,handles)functionedit1_CreateFcn(hObject,eventdata,handles)ifispc&isequal(get(hObject,BackgroundColor),get(0,defaultUicontrolBackgroundColor)set(hObject,BackgroundColor,white

14、);end%-Executesonbuttonpressinpushbutton3.functionpushbutton3_Callback(hObject,eventdata,handles)row=str2num(get(handles.edit4,String);w1=11;20;21;02;13;w2=-12;00;-10;-1-1;0-2;y1=ones(1,5);w1;y2=-ones(1,5);w2;w=y1y2;v=y1(:,1);sumz=zeros(3,1);z=ones(1,10);whilesum(z)0v=v+row*sumz;h=v*w;z=(abs(h)h);Z=

15、w.*z;z;z;sumz=sum(Z);endx1=-2:0.5:2.5;x2=(-v(1)-v(2)*x1)/v(3);cla(handles.axes1);axes(handles.axes1);holdon;plot(x1,x2);plot(w1(1,:),w1(2,:),*r)plot(w2(1,:),w2(2,:),og)title(感知准则函数分类器);text(1.5,5,*表示第1类,FontSize,10,Color,r);text(1.5,4,o表示第2类,FontSize,10,Color,g)holdoff;functionedit4_Callback(hObject

16、,eventdata,handles)%hObjecthandletoedit4(seeGCBO)%eventdatareservedtobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%Hints:get(hObject,String)returnscontentsofedit4astext%str2double(get(hObject,String)returnscontentsofedit4asadouble%-Executesduringobjectcreation,aftersettingallproperties.functionedit4_CreateFcn(hObject,eventdata,handles)%hObjecthandletoedit4(seeGCBO)%eventdatareservedtobedefinedinafutureversionofM