数字信号处理实验报告孙正.docVIP

课程设计(综合实验)报告 ( 2013 -- 2014 年度第 2 学期) 名 称：数字信号处理综合实验 院 系：电子与通信工程系 班 级： 学 号： 学生姓名： 指导教师： 孙正 设计周数： 1 成 绩： 日期： 2014 年 12 月 日 综合实验的目的与要求 1．1 实验目的 掌握《数字信号处理基础》课程的基本理论； 掌握应用MATLAB进行数字信号处理的程序设计方法。 1．2 实验要求 设计题目及设计要求： 已知低通数字滤波器的性能指标如下： ，，， 要求： 选择合适的窗函数，设计满足上述指标的数字线性相位FIR低通滤波器。用一个图形窗口，包括四个子图，分析显示滤波器的单位冲激响应、相频响应、幅频响应和以dB为纵坐标的幅频响应曲线。 用双线性变换法，设计满足上述指标的数字Chebyshev I型低通滤波器。用一个图形窗口，包括三个子图，分析显示滤波器的幅频响应、以dB为纵坐标的幅频响应和相频响应。 已知模拟信号 其中，，，，取采样频率。 要求： 以对进行取样，得到。用一个图形窗口，包括两个子图，分别显示以及()的波形； 用FFT对进行谱分析，要求频率分辨率不超过5Hz。求出一个记录长度中的最少点数，并用一个图形窗口，包括两个子图，分别显示以及的幅值； 用要求1中设计的线性相位低通数字滤波器对进行滤波，求出滤波器的输出，并用FFT对进行谱分析，要求频率分辨率不超过5Hz。求出一个记录长度中的最少点数，并用一个图形窗口，包括四个子图，分别显示（）、、（）和的幅值； 用要求2中设计的Chebyshev低通数字滤波器对进行滤波，求出滤波器的输出，并用FFT对进行谱分析，要求频率分辨率不超过5Hz。求出一个记录长度中的最少点数，并用一个图形窗口，包括四个子图，分别显示（）、、（）和的幅值 设计原理及方法 设计思路： 1）窗函数设计数字线性相位FIR低通滤波器 1、根据设计给出的技术指标，计算允许的过渡带宽度及阻带衰减，初步选择窗函数和N值。 2、若选用理想低通逼近，则计算积分，求出理想低通的冲激响应hd(n);理想低通的截止频率选择为wc=(ws+wp)/2,hd(n)=sin[wc(n-a)]/pi*n-a); 3、将hd(n)与窗函数相乘的FIR数字滤波器的冲激响应h(n)：h(n)=hd(n)w(n) 4、计算FIR数字滤波器的频率响应，并验证是否达到所要求的指标。 2）双线性变换法设计的数字Chebyshev I型低通滤波器 1、根据设计给出的数字滤波器的技术指标，计算模拟原型低通的指标 2、模拟切比雪夫1型滤波器的计算，利用function[b,a]=afd_cheb1(wp,ws,Rp,As) 进行切比雪夫1型模拟滤波器的设计；利用function[b,a]=cheb1ap_o(N,Rp,omega) 非归一化的切比雪夫1型模拟原型低通滤波器的设计 3、双线性变换法将Ha(s)转换为H(z) 4、直接型转换为级联型，计算数字滤波器的频率响应，并验证是否达到所要求的指标。 3）模拟信号的波形产生和采样及FFT频谱分析和FIR滤波及IIR滤波 N：一个记录长度中的取样数 F：频率分辨率 fs：采样周期 最小记录长度的选择： tp=NT=1/F N=fs/F 为了提高fft和ifft的计算效率,应尽量使序列长度N为2的幂次方,或将序列补零使其长使其长度为2的幂次方。 filter(b,a,x):可采用此函数给出滤波器的输出 b=H(z)的分多项式系数(对FIR b=h) a=H(z)的分母多项式系数(对FIR a=[1]) 程序代码 数字线性相位FIR低通滤波器 1 %设计FIR部分 wp=0.26*pi; ws=0.41*pi; tr_width=ws-wp; N=35; %哈明窗窗宽度 n=0:1:N-1; wc=(ws+wp)/2;%理想低通的截止频率 hd=ideal_lp(wc,N);%理想低通的冲激响应 w_ham=(hamming(N)); %选用哈明窗 h=hd.*w_ham; %FIR滤波器的冲激响应 [db,mag,pha,grd,w]=freqz_m(h,[1]); delta_w=2*pi/1000; Rp=-(min(db(1:wp/delta_w+1))) %通带衰减 As=-round(max(db(ws/delta_w+1:501))) %阻带衰减 figure(1); subplot(