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