数字信号处理综合性实验.docVIP

数字信号处理综合性、设计性实验报告 实验室： 实验日期： 年 月 日 院（系） 年级、专业、班 姓名 成绩 课程 名称 数字图像处理 实验项目 名 称 综合性试验 指导 教师 教师 评语 教师签名： 年 月 日 一 ，实验目的 1、掌握综合试验的能力 2、复习所学内容，综合所学知识 二，实验原理 利用MATLAB计算序列的离散卷积。 上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。 2. Matlab信号处理工具箱中提供了设计巴特沃思模拟滤波器的函数buttord、buttap和 butter，其中butter格式如下： 用于计算巴特沃思模拟滤波器系统函数中分子和分母多项式系数向量B和A，其中N和分别是滤波器的阶和3dB截止频率，返回向量B和A中的元素和分别是上面的表示式中的分母和分子系数。ftype缺省时，设计低通滤波器;ftype=high时，设计高通滤波器；ftype=stop时，设计带阻滤波器，此时为向量，且ftype缺省时，设计带通滤波器，带通的频率区间为。S缺省时，设计数字滤波器。 三，使用仪器，材料 微型计算机 MATLAB7.0 四，实验内容与步骤 １、用MATLAB计算序列{-2 0 1 –1 3}和序列{1 2 0 -1}的离散卷积。 ２、求下列直接型系统函数的零、极点，并将它转换成二阶节形式 ３、设采样周期T=250μs（采样频率fs =4kHz），用脉冲响应不变法和双线性变换法设计一个三阶巴特沃兹滤波器，其3dB边界频率为fc =1kHz。 五，实验过程原始记录（数据，图表，计算等） 1． a=[-2 0 1 -1 3]; b=[1 2 0 -1]; c=conv(a,b); M=length(c)-1; n=0:1:M; stem(n,c); xlabel(n); ylabel(幅度); 2． num=[1 -0.1 -0.3 -0.3 -0.2]; den=[1 0.1 0.2 0.2 0.5]; [z,p,k]=tf2zp(num,den); m=abs(p); disp(零点);disp(z); disp(极点);disp(p); disp(增益系数);disp(k); sos=zp2sos(z,p,k); disp(二阶节);disp(real(sos)); zplane(num,den) 0.9615 -0.5730 -0.1443 + 0.5850i -0.1443 - 0.5850i 极点 0.5276 + 0.6997i 0.5276 - 0.6997i -0.5776 + 0.5635i -0.5776 - 0.5635i 增益系数 1 二阶节 1.0000 -0.3885 -0.5509 1.0000 1.1552 0.6511 1.0000 0.2885 0.3630 1.0000 -1.0552 0.7679 3． [B,A]=butter(3,2*pi*1000,s); [num1,den1]=impinvar(B,A,4000); [h1,w]=freqz(num1,den1); [B,A]=butter(3,2/0.00025,s); [num2,den2]=bilinear(B,A,4000); [h2,w]=freqz(num2,den2); f=w/pi*2000; plot(f,abs(h1),-.,f,abs(h2),-); grid; xlabel(频率/Hz ) ylabel(幅值/dB) butter的边界频率2π×1000，为脉冲响应不变法原型低通滤波器的边界频率；第二个butter的边界频率2/T=2/0.00025，为双线性变换法原型低通滤波器的边界频率.图3.1给出了这两种设计方法所得到的频响，虚线为脉冲响应不变法的结果；实线为双线性变换法的结果。脉冲响应不变法由于混叠效应，使得过渡带和阻带的衰减特性变差，并且不存在传输零点。同时，也看到双线性变换法，在z=