function [U,center,result,w,obj_fcn]= fenlei(data) [data_n,in_n] = size(data) m= 2 ％ Exponent for U max_iter = 100 ％ Max. iteration min_impro =1e-5 ％ Min. improvement c=3 [center, U, obj_fcn] = fcm(data, c) for i=1:max_iter if F(U)>0.98 break else w_new=eye(in_n,in_n) center1=sum(center)/c a=center1(1)./center1 deta=center-center1(ones(c,1),:) w=sqrt(sum(deta.^2)).*a for j=1:in_n w_new(j,j)=w(j) end data1=data*w_new [center, U, obj_fcn] = fcm(data1, c) center=center./w(ones(c,1),:) obj_fcn=obj_fcn/sum(w.^2) end end display(i) result=zeros(1,data_n) U_=max(U) for i=1:data_n for j=1:c if U(j,i)==U_(i) result(i)=j continue end end end
標簽： data function Exponent obj_fcn
上傳時間： 2013-12-18
上傳用戶：ynzfm
function [U,V,num_it]=fcm(U0,X) ％ MATLAB (Version 4.1) Source Code (Routine fcm was written by Richard J. ％ Hathaway on June 21, 1994.) The fuzzification constant ％ m = 2, and the stopping criterion for successive partitions is epsilon =??????. ％*******Modified 9/15/04 to have epsilon = 0.00001 and fix univariate bug******** ％ Purpose:The function fcm attempts to find a useful clustering of the ％ objects represented by the object data in X using the initial partition in U0.
標簽： fcm function Version Routine
上傳時間： 2014-11-30
上傳用戶：二驅蚊器
兩臺處理機A 和B處理n個作業。設第i個作業交給機器 A 處理時需要時間ai，若由機器B 來處理，則需要時間bi。由于各作 業的特點和機器的性能關系，很可能對于某些i，有ai >=bi，而對于 某些j，j!=i，有aj<bj。既不能將一個作業分開由兩臺機器處理，也沒 有一臺機器能同時處理2 個作業。設計一個動態規劃算法，使得這兩 臺機器處理完成這n 個作業的時間最短（從任何一臺機器開工到最后 一臺機器停工的總時間）。研究一個實例：（a1,a2,a3,a4,a5,a6）= （2,5,7,10,5,2）；（b1,b2,b3,b4,b5,b6）=（3,8,4,11,3,4）
上傳時間： 2014-01-14
上傳用戶：獨孤求源
Euler函數: m = p1^r1 * p2^r2 * …… * pn^rn ai >= 1 , 1 <= i <= n Euler函數： 定義：phi(m) 表示小于等于m并且與m互質的正整數的個數。 phi(m) = p1^(r1-1)*(p1-1) * p2^(r2-1)*(p2-1) * …… * pn^(rn-1)*(pn-1) = m*(1 - 1/p1)*(1 - 1/p2)*……*(1 - 1/pn) = p1^(r1-1)*p2^(r2-1)* …… * pn^(rn-1)*phi(p1*p2*……*pn) 定理：若(a , m) = 1 則有 a^phi(m) = 1 (mod m) 即a^phi(m) - 1 整出m 在實際代碼中可以用類似素數篩法求出 for (i = 1 i < MAXN i++) phi[i] = i for (i = 2 i < MAXN i++) if (phi[i] == i) { for (j = i j < MAXN j += i) { phi[j] /= i phi[j] *= i - 1 } } 容斥原理：定義phi(p) 為比p小的與p互素的數的個數 設n的素因子有p1, p2, p3, … pk 包含p1, p2…的個數為n/p1, n/p2… 包含p1*p2, p2*p3…的個數為n/(p1*p2)… phi(n) = n - sigm_[i = 1](n/pi) + sigm_[i!=j](n/(pi*pj)) - …… +- n/(p1*p2……pk) = n*(1 - 1/p1)*(1 - 1/p2)*……*(1 - 1/pk)
上傳時間： 2014-01-10
上傳用戶：wkchong
//Euler 函數前n項和 /* phi(n) 為n的Euler原函數 if( (n/p) ％ i == 0 ) phi(n)=phi(n/p)*i else phi(n)=phi(n/p)*(i-1) 對于約數：divnum 如果i|pr[j] 那么 divnum[i*pr[j]]=divsum[i]/(e[i]+1)*(e[i]+2) //最小素因子次數加1 否則 divnum[i*pr[j]]=divnum[i]*divnum[pr[j]] //滿足積性函數條件 對于素因子的冪次 e[i] 如果i|pr[j] e[i*pr[j]]=e[i]+1 //最小素因子次數加1 否則 e[i*pr[j]]=1 //pr[j]為1次 對于本題： 1. 篩素數的時候首先會判斷i是否是素數。 根據定義，當 x 是素數時 phi[x] = x-1 因此這里我們可以直接寫上 phi[i] = i-1 2. 接著我們會看prime[j]是否是i的約數 如果是，那么根據上述推導，我們有：phi[ i * prime[j] ] = phi[i] * prime[j] 否則 phi[ i * prime[j] ] = phi[i] * (prime[j]-1) (其實這里prime[j]-1就是phi[prime[j]]，利用了歐拉函數的積性) 經過以上改良，在篩完素數后，我們就計算出了phi[]的所有值。 我們求出phi[]的前綴和 */
上傳時間： 2016-12-31
上傳用戶：gyq
全面詳細介紹了VHDL，英文版，作者Peter.J.Ashenden
上傳時間： 2017-01-02
上傳用戶：zaizaibang
/* * 五子棋.java * * Created on 2006年10月24日, 上午8:12 * * To change this template, choose Tools | Template Manager * and open the template in the editor. */ /** * * @author 05040045 */
標簽： template Created change choose
上傳時間： 2014-09-08
上傳用戶：jing911003
這是一本經典的設計模式書籍，我的寶貝哦 Copyright © 2002 by James W. Cooper IBM T J Watson Research Center February 1, 2002
上傳時間： 2014-08-16
上傳用戶：lmeeworm
粒子濾波程序，仿真實現自由度機器人對目標的跟蹤，使用kalman濾波估計總雅可比矩陣J，噪聲為非高斯噪聲
上傳時間： 2014-01-25
上傳用戶：頂得柱
Visual Studio 6.0是Microsoft公司推出的一個可視化應用程序集成開發環境（IDE Integrate Develop Environment）。 Visual Studio IDE不僅支持Visual C++，還支持Visual Basic、Visual J++、Visual InterDev等Microsoft系列開發工具。 本簡介介紹了通過 ActiveX Automation 接口顯示 AutoCAD 對象以及使用 Visual Basic for Applications 編程環境對這些對象進行編程的相關概念。
標簽： Visual Studio Environment Microsoft
上傳時間： 2013-12-16
上傳用戶：sjyy1001
