基于DSP的数字滤波器的设计2

-.z---..--总结资料—ECG信号的FIR滤波组员：110921035王敏(组长〕110921045李小娟1109210蔡佳佳生物医学工程学院11级医学仪器目录TOC\o"1-2"\h\z\u摘要I第1章课程设计的目的和要求1第2章系统功能介绍及总体设计方案2第3章主要内容和步骤33.1滤波器原理33.2操作步骤3第4章实验过程54.1汇编语言实验步骤与内容5第5章结论与体会10：源程序清单11参考文献12ECG信号的FIR滤波摘要：随着我国人民生活条件的逐渐改善，心血管病逐渐成为威胁人类**和生命的主要疾病之一，而心电图现已成为诊断此类疾病的常规检查手段。**是人全面开展的根底，也是社会开展的必要条件。心血管疾病是威胁人类生命的主要疾病。因此，对心血管疾病的诊断和治疗一直是各国医学界所重视，准确的进展心电信号提取，为医生提供有效的辅助分析是重要而有意义的课题。在电子技术迅速开展的今天，医用电子监护系统己逐渐应用于临床诊断中。心电监护仪就是现代临床医学中不可或缺的辅诊仪器,通过对心血管病人和危重病人实施长时间不连续的监护,为疾病的诊断、分析、治疗和研究提供重要的依据。针对心电信号的特点进展心电信号的采集、数据转换模块的设计与开发。设计一种用于心电信号采集的电路，然后进展A/D转换，使得心电信号的频率到达采样要求。心电数据采集系统是心电信号从人体体表进入计算机的必要通道,将直接影响后续的信号处理。由于心电信号会受到各种电磁干扰以及电路本身噪声的影响,所以如何有效抑制各种干扰和噪声、提高信号的完整性,是心电数据采集系统研究的重点和难点。目前对心电信号的降噪有多种方法，本文主要从滤波的方面介绍将噪声从信号中别离。关键词：A/D转换采样噪声降噪滤波ECGfilterWiththegradualimprovementofpeople'slivingconditionsinourcountry,cardiovasculardiseasehasbeeoneofthemajordiseasesthreateninghumanhealthandlife,andhasbeearoutinemethodofECGindiagnosingthisdisease.Healthisthefoundationforpeople'sall-rounddevelopment,butalsoanecessaryconditionforsocialdevelopment.Cardiovasculardiseasesarethemaindiseasesthreateninghumanlife.Therefore,thediagnosisandtreatmentofcardiovasculardiseasehasbeenvaluedallmedicalcircle,accurateECGsignale*traction,provideeffectiveau*iliaryanalysisfordoctorsisanimportantandmeaningfultask.Withtherapiddevelopmentofelectronictechnologytoday,medicalelectronmonitoringsystemhasbeengraduallyappliedinclinicaldiagnosis.ECGmonitorisindispensableinmodernclinicalmedicineau*iliaryinstrument,thelongtimemonitoringimplementationforcardiovascularpatientsandcriticalpatientsuninterrupted,provideanimportantbasisfordiseasediagnosis,analysis,treatmentandresearch.Designanddevelopmentofdataacquisition,dataconversionmoduleofECGsignalbasedonthecharacteristicsofECGsignal.DesignofacircuitforECGsignalacquisition,thentheA/Dconversion,thefrequencyofECGsamplingrequirementstoachieve.ECGdataacquisitionsystemisanecessarychannelECGsignalsfromthehumanbodyintotheputer,willdirectlyaffectthesubsequentsignalprocessing.BecauseoftheinfluenceoftheECGsignalwillbedisturbedbyelectromagneticinterferenceandcircuitnoise,sohowtoeffectivelyrestrainallkindsofinterferenceandnoise,improvesignalintegrity,istheemphasisanddifficultyoftheECGdataacquisitionsystem.ThenoisereductionofECGsignalsinavarietyofways,thisarticlemainlyfromthefilterpapertoseparatenoisefromsignal.Keywords:A/Dconversionsamplingnoisenoisereductionfiltering第1章课程设计的目的和要求随着科学技术的开展，现代对信号的要求越来越高。而作为当代大学生，并且所学专业为生物医学工程，我们不能仅局限于理论知识的学习，还要掌握信号的提取，读取，转化，滤波以及信号的分析。所以此次我们选择未经处理的心电信号，让我们基于MATLAB和DSP对心电信号进展读取和滤波，并对滤波之后的信号进展分析，增强我们对信号的了解以及实际动手能力。并且要通过课程设计，加深对DSP芯片TMS320C54*的构造、工作原理的理解，获得DSP应用技术的实际训练，掌握设计较复杂DSP系统的根本方法。通过使用汇编语言编写具有完整功能的图形处理程序或信息系统，使学生加深对所学知识的理解，进一步稳固汇编语言讲法规则。学会编制构造清晰、风格良好、数据构造适当的汇编语言程序，从而具备解决综合性实际问题的能力。第2章系统功能介绍及总体设计方案2.1系统功能介绍一个实际的应用系统中，总存在各种干扰。数字滤波器在语音信号处理、信号频谱估计、信号去噪、无线通信中的数字变频以及图像信号等各种信号处理中都有广泛的应用，数字滤波器也是使用最为广泛的信号处理算法之一。心电信号中工频干扰比拟严重，所以我们设计一个滤波器滤除工频干扰。在本设计中，我们选择教师给的第一个下载./physiobank/database/mitdb/，下载的为213号心电数据，利用MATLAB读取数据、转化格式并且求解滤波器系数，然后利用CCS进展滤波。设计一个FIR低通滤波器，其参数为：滤波器名称：FIR低通滤波器采样频率：Fs=360HZ通带/阻带截止频率：55Hz～65Hz通带最大衰减：0.5dB阻带最少衰减：50dB滤波器级数：N=160滤波器系数：由MATLAB编程求得。2.2总体设计方案流程图图1总体设计方案第3章主要内容和步骤3.1滤波器原理对于一个FIR滤波器系统，它的冲击响应总是又限长的，其系统函数可记为：其中是FIR的滤波器的阶数，为延时结，为端口信号函数。最根本的FIR滤波器可用下式表示：其中输入采样序列，是滤波器系数，是滤波器的阶数表示滤波器的输出序列，也可以用卷积来表示输出序列与、的关系，如下：3.2操作步骤1、从MIT-BIH中获取经过AD转换后的数据，将其以t*t文档的形式保存〔此时是双通道数据〕，用MATLAB程序将其读取出来，并选择其中一组通道的3000个数据。2、将获得的的数据用MATLAB程序转换成DSP能读取的dat格式文件3、用MATLAB编程得到滤波器系数4、在CCS中用TMS320C54*来进展编程实现滤波功能第4章实验过程4.1.实验步骤：1.数据下载：我们选择教师给的第一个下载./physiobank/database/mitdb/，下载的为213号心电数据〔经ADC转换后的数据〕。数据转化：用MATLAB编程读取心电数据并将心电数据转化为DSP能识别的dat格式：通过查资料，我们知道了DSP能识别的dat格式，根据我们的滤波程序，要想对心电信号进展滤波，必须将t*t数据转化为DSP能识别的dat格式。我们应用MATLAB对其进展转化。首先转化通道一的数据，利用M=M(:,1)取通道一的数据,则可以将通道一的t*t数据转化为dat格式的数据。t*t数据：dat格式数据：通道一信号波形通道二信号波形3.系数计算：A.选用低通滤波器：采样频率为360HZ，截止频率100HZ，阶数160用MATLAB编程求出滤波器的系数B.选用带阻滤波器：采样频率360HZ，阶数160,通带/阻带截止频率：55Hz～65Hz用MATLAB编程求出滤波器的系数4滤波检测：〔选用DSP〕〔1〕.首先创立新的工程，并将我们写好的asm及cmd文件添加到工程，导入rts.lib文件；〔2〕.对参加的文件进展编译，无错，生成out文件，在file中选择loadprogram,并翻开这个out文件，在asm文件中设置断点；〔3〕.选择File->Data->Load…翻开之前Matlab生成的input.dat文件；〔4〕.翻开View->Graph->Time/Frequency，修改参数〔5〕.运行程序,查看图形；5.滤波结果：信号1低通滤涉及频谱信号1带阻滤涉及频谱6.可行性论证：在确定我们要做工频干扰滤波后，我们将主要精力放在查找滤波程序上，通过查找资料，我们找到了设计滤波器并能进展低通滤波以及带阻滤波的程序〔在DSP上运行〕，通过王敏的阅读和思考，她确定我们可以运用这段程序，只需将滤波器的系数予以修改，并且将心电数据转化为DSP可以识别的dat格式即可。所以我们利用MATLAB求出了我们所需的滤波器系数，并且将心电数据转化为了DSP可以识别的dat格式。最后，以DSP为平台，实现滤波。我们认为这样一步步转化是可行的。第五章结论与体会1.从这个实验中最先是实验数据的问题，我们开场的思路是但在其中的一步中出现问题，我们在网上下载的是三种格式的数据〔.hea,.dat,.art)，在我们用程序转换之后得到的数据是最原始的的电压值，数值非常小，而且是double类型的数据，这在转换成DSP所需dat格式文件时非常不方便，因为这种格式要求为头文件+数据的形式，其中头文件信息包含五个局部：magiumbei〔固定为1651〕，数据存取格式，数据存放的起始地址，页码以及数据长度。其中数据格式可以取1-4，分别代表十六进制，整数，长整数，浮点数。此时获得的数据非常小，而且是double类型的数据，转换成以上类型是都会存在误差，是的输入波形发生变化，这会影响到比拟输入和滤波输出后图像的比拟。而改用经AD转换后的数据很方便的采用整数格式，因此我们采用的头文件格式为：1651201000fff。2.在这次实验中我们学会了求滤波系数的两种方法：一是用MATLABGUI中的SPTOOL设计滤波器系数，二是用MATLAB编程来得到滤波器系数，在实验中我们采用的是第二种方法，直接采用编程得到的特定格式数据的系数文件，在粘贴在汇编语言相应的coffe系数取。3.在ccs软件中实现滤波功能是我们遇到的两个问题：一是DSP中TMS320C54*汇编语言的不熟悉，二是DSP芯片的存储器的构造不了解，这导致我们在编写asm文件和cmd文件时出现了很多问题，在不断的查资料和纠错过程中我们对这两个方面的知识掌握的更全面。4.我们的动手能力真的很差，并且遇到问题不会自己思考，所以我们不仅要掌握知识，加强动手能力，就读于理工科，更要学会思考，思考哪里出错了，思考这里为什么不行，不要依赖于同学或者教师。：源程序清单1.利用MATLAB读取心电数据并转化为DSP能识别的dat格式：clearall;temp=importdata('D:\matlab\work\213.t*t');%读取网上下载的213号心电数据a=temp.data;%从temp这个构造体中得到数据ecg=a(1:3000,2);%得到通道1中的3000个数据fid=fopen('213ecg.dat','w');%翻开文件,'w'是将此文件定义为可写的，fid是此文件的整数标示fprintf(fid,'1651201000fff\n');%输出文件头，文件头必须是dsp所能识别的，就如此句程序所设定的fprintf(fid,'0*%*\n',ecg);%输出ecg数组，并写到与fid标示符一样的文件fclose(fid);%关闭fid标示符的文件。滤波系数求解程序：n=160;fsample=360;fp=100;fs=65;mode=1;w=hamming(n+1);wp=fp/(fsample/2);ws=fs/(fsample/2);ifmode==1h=fir1(n,wp,'low',w);endifmode==2h=fir1(n,[wp,ws],'stop',w);endm=0:n;figure(1)subplot(2,1,1);plot(m,w);grid;a*is([0n01.1*ma*(w)]);subplot(2,1,2);plot(m,h);grid;a*is([0n1.1*min(h)1.18*ma*(h)]);figure(2)freqz(h)subplot(2,1,2)h=[h,0];b=round(h*(2^15-1));fid=fopen('firlwindow.t*t','w');fori=1:32;fprintf(fid,'\r\t%s\t','.word');forj=(i-1)*32+1:(i-1)*32+31fprintf(fid,'%6.0f',b(j));endfprintf(fid,'%6.0f',b((i-1)*32+32));endfclose(fid);DSP实现滤波：fir1window.asm.mmregsin_buffer.usect"in_buf",160input.usect"in_dat",