现代信号处理教案分析.docVIP

?现代信号处理 课程设计实验报告 实验课题：现代信号处理 专业班级： 学生姓名： 学生学号： 指导老师： 完成时间： 目录 前言-------------------------------------------------2 课程设计题目、课程设计题目 使用FFT对产生的测试信号进行频谱分析并展示其幅频特性与相频特性，指定需要滤除的频带，通过选择滤波器类型，确定对应的滤波器低通、高通技术指标。FIR 数字滤波器设计原理 FIR 数字滤波器可以用下面的差分方程来描述 其中x（n）是输入序列，y（n）是输出序列，h（n）是系统的单位脉冲响应。由于其采用卷积的方法得到输出，故可以采用快速算法FFT进行卷积计算，提高信号处理速度。FIR 数字滤波器满足线性相位的条件是其乘法器的系数关于中点对称，又因为其结构无反馈的回路，所以是无条件的稳定系统。 2 窗函数法设计及MATLAB 实现 FIR 滤波器目前常用的设计方法有窗函数法和频率采样法，窗函数法是从时域进行设计，而频率采样法是从频域进行设计。窗函数法由于简单、物理意义清晰，因而得到了较为广泛的应用。窗函数法设计的基本思想是：首先根据技术指标要求，选取合适的阶数N和窗函数的类型，使其幅频特性逼近理想滤波器幅频特性。其次，因为理想滤波器的是无限长的，所以需要对进行截断，数学上称这种方法为窗函数法。 其中是希望逼近的滤波器的单位脉冲响应，是希望逼近的滤波器的幅频特性，是所设计的滤波器的单位脉冲响应，是窗函数。在截断的过程中要保证FIR滤波器的稳定性与线性相位的特性。由于截断效应使所设计滤波器的幅频特性存在误差，称之为吉布斯效应。最后，验证所设计的滤波器是否满足要求，若不满足，重新设计。常用的窗函数有矩形窗、三角窗、汉宁窗和凯塞窗等。窗函数的选择一般要满足以下三个条件： （1）具有较低的旁瓣幅度，尤其是第一旁瓣幅度。 （2）旁瓣幅度下降速度要快，以利于增加阻带衰减。 （3）主瓣宽度要窄，以获得较陡的过渡带。IIR数字滤波器的设计 IIR数字滤波器冲激响应无限宽，与模拟滤波器相匹配。因此可采取在模拟滤波器设计的基础上进一步变换的方法设计IIR滤波器，即通过已知的模拟滤波器系统的系统函数H(s)来设计数字滤波器的系统函数H(z)，主要是通过脉冲响应不变法或双线性变换法完成s平面到Z平面的变换。通过典型的模拟滤波器（如巴特沃斯滤波器、切比雪夫滤波器等）可实现一定参数要求的数字滤波器。 IIR数字滤波器的设计具体步骤： 按照既定规则将数字滤波器技术指标转为模拟低通滤波器技术指标。以转换后的技术指标为依据设计模拟低通滤波器H(s)（传递函数) 。 再按照一定规则将H(s)转换成H(z)（数字滤波器的传递函数）。若设计的是高通、带通或带阻滤波器，还需进行以下步骤：将高通、带通或带阻数字滤波器的技术指标转换为低通模拟滤波器的技术指标, 然后设计出低通H(s), 再将H(s)转换为H(z)。 t=0:pi/256:2*pi; f1=7;f2=37;f3=39;f4=43; x=100*sin(2*pi*f1*t)+100*sin(2*pi*f2*t)+100*sin(2*pi*f3*t)+100*sin(2*pi*f4*t); plot(t,x,g);grid;title(原始输入信号); xlabel(时间/t);ylabel(振幅); X=fft(x,N); plot(X);grid; title(输入信号频谱); 设计输入信号并对其进行傅里叶变换。 输入信号如下图 （2）输入信号频谱分析 幅频特性 Fs=100; len=512; f=Fs*(0:len/2-1)/len; plot(f,abs(X(1:len/2)));grid title(原始信号幅频特性曲线); 相频特性 Fs=100; len=512; f=Fs*(0:len/2-1)/len; plot(f,angle(X(1:len/2)));grid; title(相频特性曲线); 幅频特性 相频特性 2、巴特沃斯高通滤波器 高通滤波器采用Butterworth滤波器, 因为它具有较平坦的幅频特性, 而且其元件值比较合乎实际、对精度要求稍宽。用[N,Wn]=buttord(Wp,Ws,Rp,Rs,s)语句实现Butterworth滤波器阶数和边缘频率的确定；[b,a]=butter(N,Wn,s)语句来实现Butterworth滤波器的设计。其中Wp为通带截止频率；Rp为通带截止频率处