資源描述:
《Matlab 源程序如下》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。
1、functionfisher_testN=10;%個(gè)類(lèi)的樣本數(shù)range=2;%范圍offset=1.5;%偏移:決定了有樣本重疊的可能性while(1)%生成屬于omega1類(lèi)的二維樣本,范圍:[02]x1=range*rand(2,N);%生成屬于omega2類(lèi)的二維樣本,范圍:[13]x2=offset+range*rand(2,N);%定義幾個(gè)參量m1=(sum(x1')/N)';m2=(sum(x2')/N)';%樣本均值S1=zeros(2,2);fori=1:size(x1,2)S1=S1+(x1(:,i)-m1)*(x1(:,i)
2、-m1)';endS2=zeros(2,2);fori=1:size(x2,2)S2=S2+(x2(:,i)-m2)*(x2(:,i)-m2)';endSw=S1+S2;%總的類(lèi)內(nèi)離散度矩陣if(inv(Sw)3、-deltabound],'YLim',[-deltabound]);plot(x1(1,:),x1(2,:),'Marker','x','LineStyle','none');plot(x2(1,:),x2(2,:),'Marker','o','LineStyle','none');legend([{'第一類(lèi)'},{'第二類(lèi)'}]);startX=(range+offset)/2;startY=startX;xlen=1;endX=startX-xlen;endY=startY+xlen*w_(1)/w_(2);%控制繪圖的幾個(gè)量annotat
4、ion('arrow',[startXendX]/3.9,[startYendY]/3.9);text(startX,startY,'最佳投影方向')holdoff;%于是在一維空間上有:y1=w_'*x1;y2=w_'*x2;mm1=sum(y1)/N;mm2=sum(y2)/N;%類(lèi)均值%假設(shè)沒(méi)有任何先驗(yàn)知識(shí):threshold=(mm1+mm2)/2;%%%%%%%%%%%%%%%%%%%%%%%%%下面開(kāi)始試驗(yàn)這個(gè)閾值%和投影方向到底怎樣%%%%%%%%%%%%%%%%%%%%%%%%correct1=0;correct2=0;for(i
5、=1:200)sa1=range*rand(2,1);sa2=offset+range*rand(2,1);if(w_'*sa1>threshold)correct1=correct1+1;endif(w_'*sa2