北邮dsp软件matlab仿真实验报告

《北邮dsp软件matlab仿真实验报告》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、信号处理与分析课程设计训练报告----------ECG信号高频噪声消除学号：姓名：指导老师：张新峰完成日期：2016年6月5日绪论本实验是信号处理相关方法在心电图ECG（electocardiogram）信号处理中的一个具体应用。通过心电图进行疾病诊断是心血管疾病诊断的一种重要方法。该方法可靠简便，在临床中应用很广泛。近年来给予数字信号处理和模式识别的ECG信号自动分析成为一个研究这点。ECG信号在采集过程中收到工频信号、肌电信号以及基线漂移等干扰，由此带来的噪声对自动分析的结果影响很大。因此在对ECG信号分析之前，必须进行噪声消除研究。本实验内容给定ECG信号，合理设计滤波

2、器，消除高频噪声的影响。一、实验内容心电信号是微弱低频人体生理电信号，通常频率在0.05~100Hz，幅值不超过4mv，通过安装在皮肤表面的电极来获取。由于人体是一个复杂的生命系统，存在50H工频干扰及基线漂移等其他生理电信号的干扰。噪声可能会影响到医生的临床诊断，因此，需对心电信号进行滤波，即必须做好前端数据采集的软硬件设计以保证心电数据的可靠和准确。（本次实验的心电数据是老师提供，为txt文档）传统医疗设备分别采用50Hz带阻滤波器和RC高通滤波器滤除工频干扰和基线漂移。但带阻滤波器电路复杂，其特性对元器件的精度敏感，而基线漂移本质上是一种缓慢变化的低频信号，采用RC滤波器

3、很难将高通滤波器的过渡带做得十分陡峭，基线漂移补偿效果不理想。因此，模拟方法往往不太容易获得很好的特性。在本实验中，我们以C语言为基础，来实现高频噪声的消除。二、实验原理1.低通巴特沃兹模拟滤波器的实现1）概述巴特沃兹滤波器的特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。在振幅的对数对角频率的波得图上，从某一边界角频率开始，振幅随着角频率的增加而逐步减少，趋向负无穷大。一阶巴特沃兹滤波器的衰减率为每倍频6分贝，每十倍频20分贝。二阶巴特沃兹滤波器的衰减率为每倍频12分贝、三阶巴特沃兹滤波器的衰减率为每倍频18分贝、如此类推。巴特沃兹滤波器的振幅对

4、角频率单调下降，并且也是唯一的无论阶数，振幅对角频率曲线都保持同样的形状的滤波器。只不过滤波器阶数越高，在阻频带振幅衰减速度越快。其他滤波器高阶的振幅对角频率图和低级数的振幅对角频率有不同的形状。2）设计巴特沃兹模拟低通滤波器1确定参数:，，，2求滤波器阶次N：==N=4确定系统函数N为奇数时：N为偶数时：双线性Z法把模拟滤波器转换成数字滤波器3.FFT的算法按时间抽选（DIT）的基-2FFT算法。当输入序列点数为（L为整数），先将序列x(n)按n的奇偶分成以下两组：则可以将DFT化为：k=0,1,2…,其中分别为点DFT。又根据系数的周期性可得到X(k)的完整表达：三、实现方

5、法流程图1.读取心电信号数据的实现floatagTemp2[1024+2];//定义一个数组用来存放数据FILE*p;//定义指针p=fopen("119_heart1.txt","r");//打开心电图数据文件for(i=0;i<1024+2;i++)//通过循环放入数组中{fscanf(p,"%f",&agTemp2[i]);}fclose(p);2.FFT变换的实现a.倒位序处理，先将x[n]进行倒位序后，才能进行fft变换voidchange(){complextemp;unsignedshorti=0,j=0,k=0;doublet;for(i=0;i

6、){k=i;j=0;t=(log(Num)/log(2));while((t--)>0){j=j<<1;j

7、=(k&1);k=k>>1;}if(j>i){temp=x[i];x[i]=x[j];x[j]=temp;}}}b.进行蝶形运算，将保存在全局复数数组x[n]中的时域信号转为频域信号:voidfft(){inti=0,j=0,k=0,l=0;complexup,down,product;change();for(i=0;i

8、=0;k