使用MATLAB實(shí)現(xiàn)關(guān)聯(lián)維算法及報(bào)告

1、NUPTJ地理與生物信息學(xué)院2011/2012 學(xué)年第二學(xué)期實(shí)驗(yàn)報(bào)告課程名稱：醫(yī)學(xué)信號(hào)處理實(shí)驗(yàn)名稱：關(guān)聯(lián)維算法實(shí)現(xiàn)班級(jí)學(xué)號(hào)：111111111學(xué)生姓名：AAAAAAA指導(dǎo)教師：BB期:2012 年 6 月一、實(shí)驗(yàn)?zāi)康暮鸵髮?shí)現(xiàn)關(guān)聯(lián)維算法并能夠在 Matlab 平臺(tái)上編寫(xiě)出程序， 用關(guān)聯(lián)維算法分析醫(yī)學(xué)信號(hào)并得出結(jié)果。二、實(shí)驗(yàn)準(zhǔn)備可以熟練使用 Matlab 進(jìn)行編程。三、實(shí)驗(yàn)內(nèi)容1、在 matlab 平臺(tái)上實(shí)現(xiàn)關(guān)聯(lián)維算法步驟：(1)取 m=48、10、20、25,(2)取T=3,應(yīng)用時(shí)間序列重構(gòu)一個(gè) m 隹相空間，(3)取 r=1:1:5,即：從 1 變化到 5,步長(zhǎng)為 1,(4)計(jì)算 C21NC

2、(r,m)=2*仇區(qū)Xjll)。*j)Ni,y(5)重復(fù)步驟 2)到 4),(6)當(dāng) m 曾大到 25 時(shí)計(jì)算時(shí)，在1nC2rm)-lnr圖里表現(xiàn)為曲線的直線部分的斜率不再隨 m 勺變化而變化，此時(shí) D2 即是所求的吸引子的關(guān)聯(lián)維。2、用關(guān)聯(lián)維算法分析醫(yī)學(xué)信號(hào)并得出結(jié)論。四、實(shí)驗(yàn)結(jié)果實(shí)現(xiàn)代碼：fid=fopen(C:vfdb-418.txt);f=fscanf(fid,%f,31000);start=1;t=3;m=48102025;n=100;Cr=ones(55);forc=1:5f_d=ones(m(c)n);fori=1:m(c)forj=1:nf_d(i,j)=f(2,start+(

3、i-1)*t+j-1)%重構(gòu)一個(gè) m 維相空間endendy=pdist(f_d);%距離(14950 欣巨陣B=squareform(y);%變成(100 x100)sum=0;forr=1:5fori=1:100forj=1:100if(r-B(i,j)0)sum=sum+1;endendendCr(c,r)=1/nA2*(sum-100)endendsave(Cr.mat,Cr)A=importdata(Cr.mat);figure,plot(log(1:5),log(A(1,:),r*,log(1:5),log(A(2,:),m+,log(1:5),log(A(3,:),b*,log(1:5),log(A(4,:),g.,log(1:5),log(A(5,:),ro)gridon,legend(m=4,m=8,m=10,m=20,m=25,2)xlabel(10g(r),ylabel(10g(Cr)源數(shù)據(jù)曲線：測(cè)試結(jié)果:0上+m=4+*日一!a11:I10.5苫0W-0.5.-1-1.5n+m=10L JU5mF25IQ卜oE-)I)0.204 口 60.Slog(r)1.21.41613五、實(shí)驗(yàn)小結(jié)本次實(shí)驗(yàn)中理解了什么是重構(gòu)相空間算法， 進(jìn)一步知道容量維數(shù)的計(jì)算，最后完成了本實(shí)驗(yàn)要求的關(guān)聯(lián)維算法，如何實(shí)現(xiàn)關(guān)聯(lián)維算法。通過(guò)對(duì)關(guān)聯(lián)維算法的理論學(xué)