1、實驗1模糊聚類與模糊識別實驗目的：掌握模糊相似矩陣的建立方法與最大隸屬原則，會求傳遞閉包矩陣；會使用數學軟件MATLAB進行模糊矩陣的有關運算.實驗學時：2學時實驗內容：根據已知數據進行數據標準化.根據已知數據建立模糊相似矩陣，并求出其傳遞閉包矩陣.根據模糊等價矩陣繪制動態聚類圖.根據原始數據或標準化后的數據和的結果確定最佳分類.根據已知數據用最大隸屬原則進行模糊識別.實驗日期：2010年11月3日地點：104操作步驟：題目*5設有4種產品，它們的指標如下：=(3738JZ16.13.12).兀=(69,73,742Z64,17),=(7X86.49.27.6839).=(57.58.64.8

2、4.6X28).試用相關系數法建走相似知1陣，并丿刖專遞閉包法進行模糊聚類。解答：其特性指標矩陣為：373812161312、X*=697374226417738649276839A=373812161312;697374226417;738649276839;575864846328;n=4;m=6;a1=0;a2=0;fori=1:ma1=A(1,i);a2=A(1，i);forj=1:nif(a1A(j,i)a1=A(j,i);endif(a2c=0;s=0;b1=0;b2=0;b3=0;b4=0;fori=1:nforj=1:ns=0;c=0;fork=1:mb1=b1+A(i,k);

3、b2=b2+A(j,k);b3=b3+(A(i,k)-b1/m)A2;b4=b4+(A(j,k)-b2/m)A2;s=s+abs(A(i,k)-b1/m)*abs(A(j,k)-b2/m);c=sqrt(b3)*sqrt(b4);endif(j=i)R(ij)=s/c;elseR(ij)=1;endendendR1.0000NaNNaNNaN081811.00000.47870.47530.44400.35761.00000.31410.30420.24830.22971.0000傳遞閉包:tR(n,n)=l;s=0;whilesn*nfori=1:nforj=1:ntR(i,j)=0；for

4、k=1:nif(R(i,k)vR(k,j)s=R(i,k);%先取小elses=R(k,j);endif(tR(i,j)tRtR=1.00001.00001.00001.00000.81811.00000.81810.81810.44400.44401.00000.44400.30420.30420.30421.0000飛設有$個企業，記為1=輕衛2衛),對每個企業用以下m個指標來刻畫企業的技術密集水平：生產工人勞動生產率、每百萬元固定資產所容納的職工人數和技術管理人員在職工中的比例，令U.=(兀,工，再J,兀1，X，兀3依次表示上述m個指標坷=(1-&9黑0-15),込=(2丄99,0.21

5、),遇=(3.2,101,0.18),碼=(2.2,103,0.17),硝=(2.5,98,0.16),%=(2&102,0.20),嗎=(1.9,120,0.09),遍=(2.0,130,0.11).其中第1個指標的單位是萬元/(年人)，試用絕對值指標法進行聚類分析口注:56題為實驗室作業題。解答：每一步實現原理和上題大致同理，其中使用的不同的聚類公式見下。A=18950.15990.211010.181030.172.5980.162.81020.201.91200.092.01300.11;n=8;m=3;a1=0;a2=0;fori=1:ma1=A(1,i)；a2=A(1，i);for

6、j=1:nif(a1A(j,i)a1=A(j,i);endif(a2R(n,n)=1;fori=1:nforj=1:ns=0;fork=1:ms=s+abs(A(i,k)-A(j,k);endif(j=i)R(ij)=exp(-s);elseR(ij)=1;endendendR1.00000.43670.24140.50610.51220.26420.27650.22850.43671.00000.3353059510.48150.51220.17500.16690.24140.33531.00000.42540.47120.61820.10850.1034050610.59510.42541

7、.00000.64370.49310.25500.24310.51220.48150.47120.64371.00000.51590.19390.18490.26420.5122061820.49310.51591.00000.12570.11990.27650.17500.10850.25500.19390.12571.00000.59230.22850.16690.10340.24311.0000tR(n,n)=1;s=0;whilesG=1645186.4325.495124;1605381.47238448624446.72;A=84.3422.1147568

8、4.342211475684.3422.114756;B=abs(G-A)/10B=0.22800.20900.33800.36800.17000.28700.17300.10600.40600.31800.067000840C=B.C=0.22800.17000.40600.20900.28700.31800.33800.17300.06700.36800.10600.0840sum(C)/4ans=0.28580.18400.2188ans.ans=0.28580.18400.2188%藍色字段為兩次轉置，因為MATLAB中的求和與取大都是默認按照列進行計算的%所以為了需要方便，我把行變換