基于小波变换的ecg信号的压缩代码

《基于小波变换的ecg信号的压缩代码》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、MIT-BIHECG信号的数据读取方法和Matlab程序收藏最近在写一篇基于小波变换的ECG信号压缩算法的论文，遇到了怎样获取ECG信号测试数据的问题，在百度和专业论坛里搜索了一番，发现也有很多朋友为此发愁。现在论文写好了，投稿中，顺便也把怎样获取和处理ECG信号数据的方法写出来,供有需要的朋友参考，省却在百度和论坛里苦苦求索的麻烦，呵呵二一、首先，如果是对ECG心电信号进行观察、分析和诊断使用的话，有两个方法：(1)从MIT-BIH数据库下载请参考我前些天发布的文章《MIT-BIHECG心电数据的下载和读取图解》,里面有详尽介绍。http://bl

2、og.csdn.net/chenyusiyuan/archive/2008/01/06/2027887.aspx(2)用专门的Matlab心电数据读取程序我10FI在浩惠电子论坛(http://www.hhee匸com/bbs/)的“医疗器械”版块找到了读取ECG心电数据的Matlab程序(rddata.m),如获至宝啊！这个程序是由外国人写的，能够读取MIT-BIH数据库.atr>.dat>.hea三种文件的数据，根据这些数据计算出实际的心电信号值，并绘制出信号波形。程序不大，注释也算齐全,不过是英文的，需要这个程序的朋友请按以下链接下载。匿名提取

3、文件连接htB：〃pickuD.mofile.com/3497080791233097或登录Mofile,使用提取码3497080791233097提取文件PS:关于rddata.m的下载，在打开http://pickup.mo行le.com/3497080791233097后，下载链接是在“文件标签"和“文件介绍"之后的地方，夹在两块广告图片之间，共有3个链接：“推荐快车(flashget)高速下载文件下载文件(IE浏览器)下载文件(非IE浏览器广，非常隐蔽，一般点击中间的那个链接“下载文件(IE浏览器广就会弹出保存对话框了。二、如果是要对ECG信

4、号进行压缩、编码等信号处理操作上面程序获得的数据就不便于使用了，因为那是转换为具有实际意义的心电数据，信号数据值一般在・2〜2之间，单位是mV。那么，要找新的ECG读取程序来获取数据吗？不用！实际上，程序rddata.m中本身就是把MIT.dat文件中存储的二值数据转换为十进制数据，然后再进一步处理转换成具有实际意义的心电信号值。我们进行信号处理时，需要用到的就是从二值数据转换来的初始十进制数据，由于.dat文件中是三个字节存储2个数，即每个数12bits,转换后得到的十进制数范围应该是0-2048o我所理解的数据存储方式图示如下，不知是否正确，仅供

5、参考：由于rddata.m程序中的注释是英文的，且有些地方也说明不清楚，我从程序中截取出将二值数据转换为十进制数据的部分代码，将注释转换为中文，并根据自己的理解作一些补充说明，希望对大家有所帮助！具体的程序代码如下：%%程序Fun.ReadECGData用于读取ECG信号数据，将原始的二值数据转换为十进制数%输入参数及其示例：%PATH=,D:MATLABR2007bworkECGData1;%指定数据的储存路径%HEADERFILE=417.hea*;%.hea格式，头文件，可用记事本打开%DATAFILE='117.dat';%.dat格

6、式，ECG数据%SAMPLES2READ=204&%指定需耍读入的样本数%%若.dat文件中存储有两个通道的信号：%%则读入2*SAMPLES2READ个数据%输岀参数：M——一个SAMPLES2READ行2列的数据矩阵，每列数据代表一个通道的信号值%functionM=Fun_ReadECGData(PATH,HEADERFILE,DATAFILE,SAMPLES2READ)%——LOADHEADERDATA%——读入头文件数据%%示例：用记事木打开的117.hea文件的数据%%1172360650000%117.dat21220011102483

7、9311700MLII%117.dat212200111024930280830V2%#69M950654x2%#None%%%%【注】函数fprintf的功能将格式化的数据写入到指定文件中。%表达工弋：count=fprintf(fid,format,A,...)%在字符串Tormaf的控制下，将矩阵A的实数数据进行格式化，并写入到文件对象fid中。该函数返回所写入数据的字节数counto%fid是通过函数fopen获得的整型文件标识符。fid=l,表示标准输出(即输岀到屏幕显示)；fid=2,表示标准偏差。%fprintf(1；$>WORK

8、INGON%s...HEADERFILE);%在Matlab命令行窗口提示当前工作状态signalh=