1、模式識別上機作業隊別：研究生二隊姓名：孫祥威學號：11208213作業一：=(0,0),(0,l),O=(10),(l,l)。用感知器固定增量法求判別函數，設12w=(1,1,1)p=1。寫程序上機運行，寫出判別函數，打出圖表。lk解答：1、程序代碼如下：clc,clearw=001;011;-10-1;-1-1-1;W=111;rowk=1;flag=1;flagS=zeros(1,size(w,1);k=0;whileflagfori=1:size(w,1)ifisempty(find(flagS=0)flag=0;break;endk=k+1;pb=w(i,:)*W;ifpb=0flag

2、S(i)=0;W=W+rowk*w(i,:);elseflagS(i)=1;endendendW,kwp1=00;01;wp2=10;11;plot(wp1(:,1),wp1(:,2),o)holdonplot(wp2(:,1),wp2(:,2),*)holdony=-0.2:1/100:1.2;plot(1/3*ones(1,size(y),y,r-)axis(-0.251.25-0.251.25)2、判別函數。計算得到增廣權矢量為w*=(-3,0,1)t,故判別函數表達式為:3x+1=013、分類示意圖：1.20.80.60.40.2-0.2-0.200.20.40.60.811.2圖1感

3、知器算法分類結果圖作業二：在下列條件下，求待定樣本x=(2,0)T的類別，畫出分界線，編程上機。1、二類協方差相等；2、二類協方差不等。訓練樣本號k123123特征x1112-1-1-2特征x210-110-1類別12解答：經計算，兩類的協方差矩陣不相等。設P)二P(3),計算時相關項直接略12去。1、計算時，參考書中P96頁式(4-1-31)等。程序代碼如下：clc,clear,closeallD1=1,1,2;1,0,-1;D2=-1,-1,-2;1,0,-1;u1=mean(D1,2);u2=mean(D2,2);c1=zeros(size(D1,1),size(D1,1);fori=1

4、:size(D1,2)c1=c1+D1(:,i)*D1(:,i);endc1=c1/size(D1,2)-u1*u1;c2=zeros(size(D2,1),size(D2,1);fori=1:size(D2,2)c2=c2+D2(:,i)*D2(:,i);endc2=c2/size(D2,2)-u2*u2;I=eye(size(c1,1),size(c1,1);ic1=c1I;ic2=c2I;W1=-0.5*ic1;W2=-0.5*ic2;w1=ic1*u1;w2=ic2*u2;w10=-0.5*log(det(c1)-0.5*u1*ic1*u1;w20=-0.5*log(det(c2)-0

5、.5*u2*ic2*u2;symsx1x2;x=x1;x2;fprintfC決策界面方程為：)D=x*(W1-W2)*x+(w1-w2)*x+(w10-w20);pretty(D)fprintfC(2,0)代入決策面方程的值為：)value=subs(D,x1,x2,20)figureezplot(D)holdonplot(D1(1,:),D1(2,:),bo)plot(D2(1,:),D2(2,:),ks)plot(2,0,rp)運行結果顯示，決策面方程為：48x-18xx二0。1122、計算結果。從程序運行結果可以看出，x二(2,0)t帶入決策界面方程值為96,可見屬于第一類。見下圖：14

6、8x1-18x.1電=0642-4-6卜丨IIIII|-6-420246K1圖2兩類協方差不等分類結果圖圖中，黑色方塊為第二類，藍色圓圈為第一類，紅色五角星即為要判斷類21別的x二(2,0)t。綠色交叉的兩條直線即為決策面曲線。左上和右下為第一類，右上和左下為第二類。作業三：有訓練集資料矩陣如下表所示，現已知，N二9、N=N=N=3、n二2、123M二3，試問，X二(-2,2)t應屬于哪一類?訓練樣本號k123123123特征x1021-1-2-2001特征x201010-1-2-1-2類別123要求：用兩種解法：a、三類協方差不等；b、三類協方差相等。編程上機，畫出三類的分界線。解答：經計算

7、，三類的協方差矩陣不相等。設P()二P)二P()，計算時相關123項直接略去。1、計算時，參考書中P96頁式(4-1-31)等。程序代碼如下：clc,clear,closeallD1=021;010;D2=-1-2-2;10-1;D3=001;-2-1-2;X=-22;u1=mean(D1,2);u2=mean(D2,2);u3=mean(D3,2);c1=zeros(size(D1,1),size(D1,1);fori=1:size(D1,2)c1=c1+D1(:,i)*D1(:,i);endc1=c1/size(D1,2)-u1*u1;c2=zeros(size(D2,1),size(D2

8、,1);fori=1:size(D2,2)c2=c2+D2(:,i)*D2(:,i);endc2=c2/size(D2,2)-u2*u2;c3=zeros(size(D3,1),size(D3,1);fori=1:size(D3,2)c3=c3+D3(:,i)*D3(:,i);endc3=c3/size(D3,2)-u3*u3;I=eye(size(c1,1),size(c1,1);ic1=c1I;ic2=c2I;ic3=c3I;W1=-0.5*ic1;W2=-0.5*ic2;W3=-0.5*ic3;w1=ic1*u1;w2=ic2*u2;w3=ic3*u3;w10=-0.5*log(det(

9、c1)-0.5*u1*ic1*u1;w20=-0.5*log(det(c2)-0.5*u2*ic2*u2;w30=-0.5*log(det(c3)-0.5*u3*ic3*u3;symsx1x2real;x=x1;x2;D12=x*(W1-W2)*x+(w1-w2)*x+(w10-w20);D13=x*(W1-W3)*x+(w1-w3)*x+(w10-w30);D23=x*(W2-W3)*x+(w2-w3)*x+(w20-w30);fprintf(、二類決策界面方程為：)pretty(D12)fprintf(、三類決策界面方程為：)pretty(D13)fprintf(二、三類決策界面方程為：)

10、pretty(D23)DD1=x*W1*x+w1*x+w10;DD2=x*W2*x+w2*x+w20;DD3=x*W3*x+w3*x+w30;value=;fprintf(num2str(X(l),num2str(X(2)代入1-2決策面方程的值為：)v1=subs(DD1,x1,x2,X)fprintf(num2str(X(1),num2str(X(2)代入1-3決策面方程的值為：)v2=subs(DD2,x1,x2,X)fprintf(num2str(X(1),num2str(X(2)代入2-3決策面方程的值為：)v3=subs(DD3,x1,x2,X)V=v1v2v3;fprintf(可

11、見(num2str(X(1),num2str(X(2)屬于第num2str(find(V=max(V)類。n)figureholdonh=ezplot(D12);set(h,Color,r)h=ezplot(D13);set(h,Color,b)h=ezplot(D23);set(h,Color,y)plot(D1(1,:),D1(2,:),ko)plot(D2(1,:),D2(2,:),ks)plot(D3(1,:),D3(2,:),kp)plot(X(1),X(2),r*)title(分界線及樣本點分布)1-2類決策面方程為：9xx6x2一i2+33x一6x2一18x+24=0114073

12、74883553281221-3類決策面方程為：6x+6x+12xx6x2+6=0121222-3類決策面方程為：25x212xx6x227x2+24x18=01211225179981368524822、計算結果。從程序運行結果可以看出，x=（2,0）t帶入三個判決函數（略去P）相關項）計算的值分別為-95.3521、-17.3521、-29.3521,可見屬于第二類iW。分類界面參見下圖:2作業四：樣本12345678x102156567x202133445用對分法編程上機，分成兩類畫出圖形。解答:1、程序代碼如下：clc,clear,closeallExampleData=0,0,2,2

13、,4,4,5,6,6,7,-4,-2,-3,-3,-5,1,0,0,-1,-1,-3;6,5,5,3,4,3,1,2,1,0,3,2,2,0,2,1,-1,-2,-1,-3,-5;Data=0,2,1,5,6,5,6,7;0,2,1,3,3,4,4,5;N=size(Data,2);Data=1:N;Data;G1=Data;G2=;N1=size(G1,2);N2=size(G2,2);X1=mean(G1,2);X1=X1(2:end);X2=0;0;E=N1*N2/N*(X1-X2)*(X1-X2);RE=;ZY=;LE=E;whileLE=ELE=E;n1=size(G1,2);e1=

14、;fori=1:n1g1=G1;g2=G2;g2=g2g1(:,i);g1(:,i)=;nn1=size(g1,2);nn2=size(g2,2);x1=mean(g1,2);x1=x1(2:end);x2=mean(g2,2);x2=x2(2:end);ee1=nn1*nn2/N*(x1-x2)*(x1-x2);e1=e1;ee1;endE=max(e1);RE=RE;E;I=find(e1=E);nzy=G1(1,I);ZY=ZY;nzy;ifLEE3斗-*口-3-口-21-OO-1111O123456圖5K次平均算法聚類結果圖4、ISODATA算法實現程序為從網上下載的代碼，文件較長，參見附上的源代碼。調用計算作業數據代碼如下：clc,clear,closeallX=0124556111;Y=0113345456;ON=1;OC=4;OS=1;k=4;L=1;I=10;NO=1;min=50;centro,Xcluster,Ycluster,A,clustering=isodata(X,Y,k,L,I,ON,OC,OS,NO,min);clc;fprintf(聚成類的數目：dn