心电信号系统的设计与仿真分析报告.doc

信号与线性系统课程设计报告 课题二 　　　　－－心电信号分析系统的设计与仿真 班级:　　通信c122班 姓名：　　张　怡 学号：　　125949 成绩： 指导教师：刘翠响 日期： 摘要：美国麻省理工学院提供的MIT-BIH数据库是一个权威性的国际　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　心电图检测标准库， 　　　近年来应用广泛，为我国的医学工程界所重视。MIT-BIH数据库共有48个病例，每个病例数据长30min，总计约有116000个心拍，包含有正常的心拍和异常的心拍内容丰富完整， 本报告意在对103.txt文件里的心电信号作为我们的原心电信号数据。利用matlab提供的文件textread或textscan函数，读取103.txt数据文件中的信号，并且还原实际波形。 因为在信号获取时可能会有其他声音的录入，所以本报告采取滤波器来获取比较纯正的心电信号，并对此分析。借助matlab和sinmlink来分析。最终通过GUIJIE界面的形式来展现。 关键词：MATLAB，心电信号，滤波器 １．课程设计的目的和意义 本设计课题主要研究数字心电信号的初步分析及滤波器的应用。通过完成本课题的设计，拟主要达到以下几个目的： 　　(1)了解MATLAB软件的特点和使用方法，熟悉基于Simulink的动态建模和仿真的步骤和过程； 　　(2)了解LabView虚拟仪器软件的特点和使用方法，熟悉采用LabView进行信号分析，系统设计　　　方法即仿真过程。 　　(3)了解人体心电信号的时域特征和频谱特征； 　　(4)进一步了解数字信号的分析方法； 　　(5)通过应用具体的滤波器进一步加深对滤波器理解； (6)通过本课题的设计，培养学生运用所学知识分析和解决实际问题的能力。 　　设计一个简单的心电信号分析系统。对输入的原始心电信号，进行一定的数字信号处理，进行频谱分析。采用Matlab语言设计，要求分别采用两种方式进行仿真，即直接采用Matlab语言编程的静态仿真方式、采用Simulink进行动态建模和仿真的方式。根据具体设计要求完成系统的程序编写、调试及功能测试。 3设计方案及论证 (1)心电信号的读取 txt格式的数据文件内容及格式： (2)利用读取和插值函数对信号分析 　 利用老师提供的函数对103.txt数据分析 [t,Xn]=duquexinhao1('103.txt'); baocun1(t,Xn); [t2,Xn2]=chazhi1(t,Xn); baocun2(t2,Xn2); subplot(2,1,1),plot(t,Xn),title('原信号时域波形') subplot(2,1,2),plot(t2,Xn2),title('插值后的时域波形') ％生成两个数据文件：Ｘn.txt(原信号有用的数据） 　　　　　　　　　Ｘn2,txt（插值后的数据） (3)对信号进行采样分析 　注：为了后面能够快速方面调用Xn2，t2数据，首先要生成Xn2.m,t2.m文件 　 程序如下： 　　　　　　save Xn2,Xn2; save t2,t2; 　（一）对Xn2点的采样 　　　　load Xn2 　　　　N=9998; 　　　　n=0:N-1 　　　　x=Xn2 　　　　plot(n,abs(x));title('信号图'); （二）对Xn2点的频谱分析 　　　　load Xn2 　　　　N=9998; 　　　　n=0:N-1 　　　　x=Xn2 　　　　y=fft(x,N) 　　　　plot(n,abs(y));title('信号的频谱'); 根据频谱图分析，设计高通滤波器和高通滤波器的参数 　 　低通滤波器的数字指标：wp=0.375*pi;ws=0.625*pi;Rp=1;As=15;T=1; 　高通滤波器的数字指标： wp=0.8*2*pi;ws=0.25*2*pi;Rp=0.1;As=40; 　 3.1利用巴特沃斯响应不变法设计低通滤波器 wp=0.375*pi;ws=0.625*pi;Rp=1;As=15;T=1;%数字指标 OmegaP=wp/T;OmegaS=ws/T;%模拟指标 [cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As)%butterworth低通滤波器原型设计函数; [b,a]=imp_invr(cs,ds,T)%利用脉冲响应不变法从模拟到数字滤波器变换函数 [db,mag,pha,w]=freqz_m(b,a); subplot(2,1,1);plot(w/pi,mag); title('响应幅度')