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 [alpha,N,U]=youxianchafen2(r1,r2,up,under,num,deta)
%[alpha,N,U]=youxianchafen2(a,r1,r2,up,under,num,deta)
%該函數用有限差分法求解有兩種介質的正方形區域的二維拉普拉斯方程的數值解
%函數返回迭代因子、迭代次數以及迭代完成后所求區域內網格節點處的值
%a為正方形求解區域的邊長
%r1,r2分別表示兩種介質的電導率
%up,under分別為上下邊界值
%num表示將區域每邊的網格剖分個數
%deta為迭代過程中所允許的相對誤差限
n=num+1; %每邊節點數
U(n,n)=0; %節點處數值矩陣
N=0; %迭代次數初值
alpha=2/(1+sin(pi/num));%超松弛迭代因子
k=r1/r2; %兩介質電導率之比
U(1,1:n)=up; %求解區域上邊界第一類邊界條件
U(n,1:n)=under; %求解區域下邊界第一類邊界條件
U(2:num,1)=0;U(2:num,n)=0;
for i=2:num
U(i,2:num)=up-(up-under)/num*(i-1);%采用線性賦值對上下邊界之間的節點賦迭代初值
end
G=1;
while G>0 %迭代條件：不滿足相對誤差限要求的節點數目G不為零
Un=U; %完成第n次迭代后所有節點處的值
G=0; %每完成一次迭代將不滿足相對誤差限要求的節點數目歸零
for j=1:n
for i=2:num
U1=U(i,j); %第n次迭代時網格節點處的值
if j==1 %第n+1次迭代左邊界第二類邊界條件
U(i,j)=1/4*(2*U(i,j+1)+U(i-1,j)+U(i+1,j));
end
if (j>1)&&(j U2=1/4*(U(i,j+1)+ U(i-1,j)+ U(i,j-1)+ U(i+1,j));
U(i,j)=U1+alpha*(U2-U1); %引入超松弛迭代因子后的網格節點處的值
end
if i==n+1-j %第n+1次迭代兩介質分界面（與網格對角線重合）第二類邊界條件
U(i,j)=1/4*(2/(1+k)*(U(i,j+1)+U(i+1,j))+2*k/(1+k)*(U(i-1,j)+U(i,j-1)));
end
if j==n %第n+1次迭代右邊界第二類邊界條件
U(i,n)=1/4*(2*U(i,j-1)+U(i-1,j)+U(i+1,j));
end
end
end
N=N+1 %顯示迭代次數
Un1=U; %完成第n+1次迭代后所有節點處的值
err=abs((Un1-Un)./Un1);%第n+1次迭代與第n次迭代所有節點值的相對誤差
err(1,1:n)=0; %上邊界節點相對誤差置零
err(n,1:n)=0; %下邊界節點相對誤差置零
G=sum(sum(err>deta))%顯示每次迭代后不滿足相對誤差限要求的節點數目G
end
標簽：
有限差分
上傳時間：
2018-07-13
上傳用戶：Kemin
VIP專區-嵌入式/單片機編程源碼精選合集系列（120）資源包含以下內容：1. ARM的一個GPIO口的應用 它使能一個口 驅動燈亮.2. bit led2=P2^5 // led2對應接在P2.5腳
sbit led3=P2^6 // led3對應接在P2.6腳
sbit led4=P2^7 // led4對應接在P2.7腳.3. EA = 1 //開總中斷
ET0 = 1 //允許定時器0中斷
TMOD = 1 //定時器工作方式選擇
TL0 = 0x06
TH0 = 0xf8 //定時器賦予初值
.4. 嵌入式實驗代碼非常適合全面學習ARMS3C2440,有20多個練習代碼.5. 基于C8051F340單片機的串行flash芯片SST25VF016B的應用,配合一般的串口調試程序即可完成對板上器件的讀寫操作..6. 飛思卡爾MC9S12DG128的PORTP7中斷試驗.7. 微細間距QFP器件手工焊接指南,希望對大家有用.8. 金鵬OCMJ8*10D的觸摸屏驅動程序.9. CODE for embedded C ,hand coding version.10. konqueror3 embedded版本, KDE環境下的當家瀏覽器的嵌入式版本源碼包..11. 三菱PLC編寫用與點焊機上的PLC程序,另包括完整電氣原理圖..12. 24CXXX存儲數碼管
24C02讀寫(c)
24C02讀寫(匯編).13. 電機控制類
步進電機正反轉控制
步進電機調速控制.14. 個人設計的基于VHDL的數字電子日歷
在MAX+PLUSH軟件平臺上編譯、仿真.15. 超聲波倒車雷達的設計程序
利用SPCE061A單片機實現超聲波倒車雷達的測量計算方法。程序中采用軟件校正.16. PATTERNS FOR TIME-TRIGGERED EMBEDDED SYSTEMS by Michael J. Pont
This code is copyright (c) 200.17. EP9315 開發板手冊
cirrus logic(思睿邏輯)公司在2005年推出的一款ARM920T.18. 譯本嵌入式數據庫sqlite的電子書,值得大家參考..19. 1302是用來輸出北京時間的,18b20用來檢測溫度,輸出在ds1302上顯示,有仿真文件.20. 可將其需要的文件轉換成數組,主要應用在嵌入式或單片機編程中將要分析的文件生成一個數組參加程序編譯.21. 這是關于TDOA_AOA得一篇混合定位算法,有助于研究算發的人學習和使用..22. 51串口通訊程序.23. 一個開源的Modbus協議棧.24. 能夠詳細測量正負溫度的且小數點后四位的測溫系統.25. 一本關于嵌入式系統實時概念的著作.26. C語言寫的Ibutton的讀寫文件,硬件資料
AT89s51+at24c02,編譯通過,可以使用.27. arm的原理圖設計.28. 本書深入淺出的介紹了可重用的理念與實現.29. 德州儀器的關于cc2430評估模塊的原理圖及PCB資料.30. CC2430DB Reference Design TI的cc2430開發板設計資料.31. 嵌入式MCU可靠性接口及在車載定位系統中的應用.32. au1200下的boot代碼.33. ARM嵌入式系統開發：軟件設計與優化.34. 友善之臂－－－QQ2440V3原理圖.35. zigbee協議棧.36. 國外自制頻率計,很全,有原理圖,有源代碼,有說明..37. ECOS2.0的源碼.38. trf7960開發板原理圖,是TI公司的.39. c8051f的多機通信程序
在c語言環境下.40. 祥細介紹液晶顯示器原理祥細介紹液晶顯示器原理.
標簽：
通信
上傳時間：
2013-04-15
上傳用戶：eeworm
