使用MATLAB实现关联维算法及报告.docVIP

地理与生物信息学院 2011 / 2012 学年第 二 学期 实 验 报 告 课程名称： 医 学 信 号 处 理 实验名称： 关 联 维 算 法 实 现 班 级 学 号 ： 111111111 学 生 姓 名 ： AAAAAAA 指 导 教 师 ： B B 日 期： 2012 年 6 月 一、实验目的和要求 实现关联维算法并能够在Matlab平台上编写出程序，用关联维算法分析医学信号并得出结果。 二、实验准备 可以熟练使用Matlab进行编程。 三、实验内容 1、在matlab平台上实现关联维算法 步骤： （1）取m=4、8、10、20、25， （2）取=3，应用时间序列重构一个m维相空间， （3）取r=1：1：5，即：从1变化到5，步长为1， （4）计算C2， （） （5）重复步骤2）到4）， （6）当m增大到25时计算时，在图里表现为曲线的直线部分的斜率不再随m的变化而变化，此时D2即是所求的吸引子的关联维。 2、用关联维算法分析医学信号并得出结论。 四、实验结果 实现代码： fid=fopen('C: \vfdb-418.txt'); f=fscanf(fid,'%f',[3 1000]); start=1; t=3; m=[4 8 10 20 25]; n=100; Cr=ones([5 5]); for c=1:5 f_d=ones([m(c) n]); for i=1:m(c) for j=1:n f_d(i,j)=f(2,start+(i-1)*t+j-1) %重构一个m维相空间 end end y=pdist(f_d'); %距离(1 4950)矩阵 B=squareform(y); %变成(100 x 100) sum=0; for r=1:5 for i=1:100 for j=1:100 if(r-B(i,j)>0) sum=sum+1; end end end Cr(c,r)=1/n^2*(sum-100) end end save('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') grid on,legend('m=4','m=8','m=10','m=20','m=25',2) xlabel('log(r)'),ylabel('log(Cr)') 源数据曲线： 测试结果： 实验小结 本次实验中理解了什么是重构相空间算法，进一步知道容量维数的计算，最后完成了本实验要求的关联维算法，如何实现关联维算法。通过对关联维算法的理论学习和老师的讲解，用matlab实现。由于本实验循环较多，除了减少实验数据以外，尽量优化代码，简化算法，利用技巧完成实验。