用MATLAB实现DFA算法及报告_第1页
用MATLAB实现DFA算法及报告_第2页
用MATLAB实现DFA算法及报告_第3页
用MATLAB实现DFA算法及报告_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

.PAGE.地理与生物信息学院2011/2012学年第二学期实验报告课程名称：医学信号处理实验名称：DFA算法实现班级学号：********** 学生姓名：****指导教师：****日期：2012年6月实验题目：Matlab软件处理生物医学数字信号二、实验要求：采用DFA算法实现实验内容：能使用Matlab软件编写、实现DFA算法并处理生物医学信号实验过程：〔1针对某个心电时间序列<总长度为>,首先进行求和其中,是第个数据,而是所分析的心电时间序列的平均值。这种求和步骤可以映射原始时间序列到一个自相似的过程。接着,量度求和后的时间序列的垂直特征尺度。方法是：求和后的时间序列被分成等长为的许多个小片段。画出每一个长度为的小片段的最小二乘拟和直线〔它用来代表该片段中的趋势。直线段的y坐标标记为。〔3其次,对求和后的时间序列进行去趋势,即在每一个片段中把减去局域趋势。对一个给定的片段长度,这种求和并且去趋势后的时间序列的涨落的特征尺寸可以由下式计算：〔4在所有时间尺度〔片段大小上重复上面的计算,可以得出和片段大小的关系曲线。对的斜率决定尺度指数〔自相似参量。五、主要实验程序：Clearall;A=textread<'d:\nsrdb-16786.txt'>;N=1000;n=[48163264128];l=length<n>;r=fix<rand<1>*10000>;b<1:N>=A<r:<r-1+N>,2>*100+2100;plot<1:N,b<1:N>>;c=b-mean<b>;fori=1:Ny<i>=sum<c<1:i>>;endfork=1:l%k=1n=10k=2n=50...d=N/n<k>;forj=1:dp=polAit<<j-1>*n<k>+1:j*n<k>,y<<j-1>*n<k>+1:j*n<k>>,1>;y1<k,<j-1>*n<k>+1:j*n<k>>=polyval<p,<j-1>*n<k>+1:j*n<k>>;endfigure,plot<1:N,y>;holdon;plot<1:N,y1<k,1:N>,'r'>;set<gca,'XGrid','on'>;%set<gca,'XTick',[0:n<k>:N]>;title<['n=',num2str<n<k>>]>;f<k>=sqrt<<sum<<y<1:N>-y1<k,1:N>>.^2>>/N>;endfigure,plot<log10<n<1:k>>,log10<f<1:k>>,'o'>;holdon;a=polAit<log10<n<1:k>>,log10<f<1:k>>,1>;plot<log10<n<1:k>>,polyval<a,log10<n<1:k>>>>;title<['拟合直线斜率',num2str<a<1>>]>;原始数据曲线：n=4、8、16、32、64、128拟合直线：六、实验小结：DFA方法可以检测出一个非平稳时间序列的内在的自相似性,同时避免散杂地检测可明显看出

人人文库> 全部分类> 应用文书 > 合同范本

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论