数字信号处课程设计二心电信号分析系统设计-matlab.docVIP

数字信号处理 课程设计报告 课题名称： 心电信号分析系统的设计 –matlab 班 级: 通信103 姓 名： 滕冲 学 号： 102189 成 绩： 指导教师： 王宝珠 日 期： 2011年12月28日 目录 摘要 ………………………………………………………………………3 课程设计的目的、意义…………………………………………………3 设计任务及技术指标 ……………………………………………………3 三、 设计方案论证……………………………………………………………3 相关函数及程……………………………………………………………3 读取txt形式的心电信号M文件………………………………………3 (2) 数据保存：分别用来保存插值前和插值后的数据………………………4 线性插值：………………………………………………………………4 (4)对以上步骤进行画图………………………………………………………4 (5) 两个滤波器………………………………………………………………5 （6）主程序……………………………………………………………………5 五、 实验结果与分析………………………………………………………10 六、 结论（只对设计课题内容进行总结并给出结论）……………………14 七、课程设计心得与体会……………………………………………………15 八、参考文献…………………………………………………………………16 摘要：设计一个心电信号分析系统。读取原始心电信号txt文件并将其前两行内容用程序的方法去除，对其做一次线性插值并保存在新的txt文件里面。由心电信号特点得知需要设计带通滤波器还有一个带阻滤波器，设置初步滤波技术指标画出初步滤波图形后看出的结果好不好，出图是否与设计指标相符，若计数指标不好应该怎样调节，对处理完的时域和频域波形进行分析，得出结论，做过程记录。该系统做完后，启用simulink功能，建立模型，进行连线，设置参数，最后启动仿真。观察示波器模块的输入和输出波形是否真确。 关键词：数字心电信号，线性插值，滤波器，功能仿真 课程设计的目的、意义 通过本课题的设计学会用matlab软件实现对心电信号进行读取，设计滤波器并进行滤波，分析心电信号的特点并反馈给滤波器的设计指标，同时强化对matlab语句应用的训练更好地掌握该软件；提高对实际问题的分析解决能力，锻炼综合能力；通过实际的动手过程发现问题解决问题，加深对所学知识的理解和掌握，重在过程。 二、 设计任务及技术指标 (1)对原始数字心电信号进行读取，由数字信号数据绘制出其时域波形。 (2)对数字信号数据做一次线性插值，使其成为均匀数字信号，以便后面的信号分析。 (3)根据心电信号的频域特征（自己查阅相关资料），设计相应的低通和高通滤波器。 (4)编程绘制实现信号处理前后的频谱，做频谱分析，得出相关结论。 (5)对系统进行综合测试，整理数据，撰写设计报告。 三、 设计方案论证 （1）采用线性插值函数interp1完成插值处理 （2）设计一个带通滤波器，经画图分析，需加上一个带阻滤波器，滤除60Hz处的干扰信号 （3）初始技术指标： 带通：wp=[1,90]*2*pi（滤波器通带截止频率）; ws=[0.25,99]*2*pi（阻带截止频率）; Rp=1db（通带衰减）； As=25db（阻带衰减）； 带阻：wp1=[58,62]*2*pi;（滤波器通带截止频率）; ws1=[59.9,60.1]*2*pi;（阻带截止频率）; Rp=1db（通带衰减）； As=25db（阻带衰减）； 四、 相关函数及程序 读取txt形式的心电信号M文件 %读取心电信号并转化成数组形式 function [t,Xn]=duquxinhao(w) fid=fopen(w); A=textscan(fid,%8c %f %*f,headerlines,2); %将其前两行内容去除 fclose(fid); a=A{1}; b=A{2}; k=length(a) for i=1:k c(i)=strread(a(i,:),%*s %f,delimiter,:); end c=c d=[c,b]; %为了减少分析数据的工作量截取2.5s(三个周期左右)的心电信号 for i=1:k if c(i)=2.5 e(i,:)=d(i,:); else break; end end t=e(:,1); %时间 Xn=e(:,2); %幅度 数据保存：分别用来保存插值前和插值后的数据 function sh