医学信号处理课程设计.pptVIP

医学信号处理课程设计 指导教师：侯海燕 符志鹏 一、课程设计基本要求 1 巩固和加深对数字信号处理的基本概念、理论、方法的理解，提高综合运用这些知识进行程序设计的能力。 2 能够根据需要查阅相关的手册和资料，提高独立分析问题和解决问题的能力。 3 掌握基于MTLAB的数字信号处理领域的计算机辅助分析和设计的基本方法。 4 能独立撰写设计报告，准确分析实验结果。 5 培养自学意识，严肃认真的治学态度和严谨求实的工作作风。 二、课程设计目的 1 熟悉和巩固数字滤波器的设计方法和原理 2 熟练掌握使用MATLAB程序设计方法 3 掌握序列离散傅里叶变换的MATLAB实现，并进行频谱分析 4 熟练掌握使用MATLAB设计IIR或FIR数字滤波器 三、设计原理 频谱分析（傅里叶变换） 滤波器设计 IIR（冲激响应不变法，双线性变换法） FIR（窗函数法，频率抽样法） 信号滤波（时域卷积，频域相乘） 双线性变换法IIR数字滤波器的设计 （1）利用双线性变换的频率预畸变公式，把所要求的数字频率指标w转换为相应的模拟频率指标Ω。 （2）将模拟滤波器的频率指标Ω归一化η，并转换成归一化原型低通滤波器的频率指标λ； （3）按照巴特沃思、切必雪夫低通滤波器的设计方法，得到模拟低通滤波器的归一化传输函数G(p) ； （4）进行模拟频率逆变换，将G(p)转换所要求的模拟滤波器的传输函数H(s)； （5）利用双线性变换法将得到的模拟滤波器传输函数H(s)转换成所要求的数字滤波器的传输函数H(z)。 冲激响应不变法IIR数字滤波器的设计 （1）把所要求的数字频率指标w转换为相应的模拟频率指标Ω=w/T。 （2）将模拟滤波器的频率指标Ω归一化η，并转换成归一化原型低通滤波器的频率指标λ； （3）按照巴特沃思、切必雪夫低通滤波器的设计方法，得到模拟低通滤波器的归一化传输函数G(p) ； （4）进行模拟频率逆变换，将G(p)转换所要求的模拟滤波器的传输函数H(s)； （5）利用冲激响应不变法法将得到的模拟滤波器传输函数H(s)转换成所要求的数字滤波器的传输函数H(z)。 窗函数法FIR数字滤波器的设计 (1)根据技术要求确定待求滤波器的单位取样响应hd(n)。如果给出待求滤波器的频响为Hd(ejω)，那么单位取样响应用下式求出 (2)根据阻带衰减选择窗函数w(n) 根据过渡带宽度Δω确定N值。 (3) 加窗，计算滤波器的单位取样响应h(n)， h(n)=hd(n)w(n) 验算技术指标是否满足要求，如不符合要求修改适当改变窗函数形状或长度，重复上述步骤直到满意为止。 四、选题 1-7 基于MATLAB的心电信号的分析与处理设计 8,16 使用FFT实现任意三个同频带信号的频分复用 9-15 基于MATLAB的声音信号的分析与处理设计 基于MATLAB的心电信号分析与处理设计（1） 1．找到一段心电信号，画出心电信号的时域波形和频谱图（幅频和相频）； 2．给该段心电信号加上随机噪声信号（噪声幅值为心电信号最大幅值的3%），画出加噪信号的时域波形和频谱图（幅频）。 3．根据心电信号频率范围及噪声的频率范围选择一种滤波器（高通，低通，带通，带阻中的一种）并设置滤波器性能指标，用冲激响应不变法设计一个IIR巴特沃斯滤波器以消除2中加噪心电信号中的噪声干扰，输出所设计的滤波器的阶次和分子、分母多项式的系数，画出滤波器的频率响应（幅频响应和相频响应）曲线； 4．用该滤波器对2中的加噪心电信号进行滤波，画出滤波以后心电信号的时域波形和频谱（幅频）；分析信号滤波前后时域和频域的变化； 5．运用GUI设计一个心电信号处理系统界面（选作）。 基于MATLAB的声音信号分析与处理设计（1） 1．录制2秒自己的声音，如发音“河南科技大学”等简短的语句，抽样率是7000Hz/s，共得14000个数据。 2．抽取7000个数据进行频谱分析，得到幅度谱。 3．取这14000个数据进行频谱分析，得到幅度谱。比较二者异同并分析原因。 4．分析人声的频率范围，根据所录制的声音信号的频谱特性，选择合适的滤波器参数，选择矩形窗用窗函数法设计一个FIR低通滤波器；要求输出所设计的滤波器的阶次，画出滤波器的单位抽样响应，分析其对称性；画出滤波器的频率响应（幅频响应和相频响应）曲线，根据相频曲线判断滤波器是否为线性相位滤波器； 5．用滤波器对信号进行滤波，画出滤波以后信号的时域波形和频谱，与原信号进行时域和频域的对比。 6．再把处理后的所有数据储存为声音文件，与原始声音进行比较。 使用FFT实现任意三个同频带信号的频分复用 1．利用MATLAB产生三路正弦信号，这三路信号的频率范围设置不同（低频，中频，高频），（注意采样频率，中频信号频率设置在