数字信号实验报告.doc

《数字信号实验报告.doc》由会员分享，可在线阅读，更多相关《数字信号实验报告.doc（10页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、数字信号实验报告科目：数字信号处理姓名：殷超宇班级：14060142 学号：1406014226实验题目：Z Z 变换及离散时间系统分析p 指导教师：张志杰分数：一实验题目：Z 变换及离散时间系统分析p 二实验目的：1、通过本实验熟悉 Z 变换在离散时间系统分析p 中的地位和作用。2、掌握并熟练使用有关离散系统分析p 的 MATLAB 调用函数及格式，以深入理解离散时间系统的频率特性。三实验内容：给定系统 ) 8 .0 /( 2 .0 ) (2+ - = z z H ，编程并绘出系统的单位阶跃响应 y(n),频率响应 ) e (jwH ，并给出实验数据与代码。四参考代码：详见数字信号处理上机实

2、验指导(班群里有)五实验代码（代码从 B MATLAB ）软件复制粘贴于此处，教师检查重点）: :clear;_=ones(100); _(n)=1,n=1100；t=1:100; t 用于后面的绘图；b=0,0,-0.2; 形成向量 b；a=1,0,0.8; 形成向量 a；y=filter(b,a,_); 求所给系统的阶跃响应；plot(t,y,k-); grid on;ylabel( y(n)_label(n)六实验数据（图像或表格复制粘贴于此处，教师检查重点）：七实验心得与收获（可手写）：a,b 两个向量转化成符合的格式。这是一个震荡衰减的信号，信号越来越弱，通过实验，对 matlab

3、的使用，有了进一步的了解。明白不故作难懂。南京邮电大学实 验 报 告实验名称_熟悉MATLAB环境 _快速傅里叶变换及其应用 _IIR数字滤波器的设计_FIR数字滤波器的设计课程名称 数字信号处理A班级学号_9002111_姓 名 王都超 开课时间 2021/20_学年， 第 二 学期实验一熟悉MATLAB环境一、实验目的（1）熟悉MATLAB的主要操作命令。 （2）学会简单的矩阵输入和数据读写。 （3）掌握简单的绘图命令。（4）用MATLAB编程并学会创建函数。 （5）观察离散系统的频率响应。 二、实验内容(1) 数组的加、减、乘、除和乘方运算。输入A=1 2 3 4，B=3,4,5,6，求

4、C=A+B， D=A-B，E=A._B,F=A./B,G=A.B 。并用stem语句画出A、B、C、D、E、F、G。 D =-2-2-2-2 E =24F =0.33330.50000.60000.6667 G =2434096 (2) 用MATLAB实现下列序列： a) _(n)=0.8n 0n15n=0:1:15; _1=0.8.n; a=(0.2+3_i)_n; stem(_1) b) _(n)=e(0.2+3j)n 0n15n=0:1:15; _2=e_p(a); a=(0.2+3_i)_n; stem(_2)c) _(n)=3cos(0.125pn+0.2p)+2sin(0.25pn

5、+0.1p)0n15(4) 绘出下列时间函数的图形，对_轴、y轴以及图形上方均须加上适当的标注： a) _(t)=sin(2pt) 0t10sb) _(t)=cos(100pt)sin(pt) 0t4s t=0:0.01:4; _=cos(100_pi_t)._sin(pi_t); plot(t,_, r-); _label(t),ylabel(_(t),title(cos) (6) 给定一因果系统H(z)=(1+频响应和相频响应。2z-1+z-2)/(1-0.67z-1+0.9z-2)，求出并绘制H(z)的幅(7) 计算序列8 -2 -1 2 3和序列2 3 -1 -3的离散卷积，并作图表示

6、卷积结果。(8) 求以下差分方程所描述系统的单位脉冲响应h(n), 0n50y(n)+0.1y(n-1)-0.06y(n-2)=_(n)-2_(n-1)实验过程与结果（含实验程序、运行的数据结果和图形）； clear all; N=50; a=1 -2; b=1 0.1 -0.06; _1=1 zeros(1,N-1); n=0:1:N-1; h=filter(a,b,_1); stem(n,h) a_is(-1 53 -2.5 1.2)实验二快速傅里叶变换及其应用一、实验目的(1) 在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉MATLAB中的有关函数。 (2) 应用FFT对典型信

7、号进行频谱分析p 。(3) 了解应用FFT进行信号频谱分析p 过程中可能出现的问题，以便在实际中正确应用FFT。 (4) 应用FFT实现序列的线性卷积和相关。 二、实验内容实验中用到的信号序列 a) 高斯序列-(n-p)q_a(n)=e020n15 其他b) 衰减正弦序列e-ansin(2pfn)_b(n)=00n15其他c) 三角波序列 n_c(n)=8-n00n34n7 其他d) 反三角波序列4-n_d(n)=n-400n34n7 其他(1) 观察高斯序列的时域和幅频特性，固定信号_a(n)中参数p=8，改变q的值，使q分别等于2，4，8，观察它们的时域和幅频特性，了解当q取不同值时，对信

8、号序列的时域幅频特性的影响；固定q=8，改变p，使p分别等于8，13，14，观察参数p变化对信号序列的时域及幅频特性的影响，观察p等于多少时，会发生明显的泄漏现象，混叠是否也随之出现？记录实验中观察到的现象，绘出相应的时域序列和幅频特性曲线。 (3) 观察三角波和反三角波序列的时域和幅频特性，用N=8点FFT分析p 信号序列_c(n)和观察两者的序列形状和频谱曲线有什么异同？绘出两序列及其幅频特性_d(n)的幅频特性，曲线。在_c(n)和_d(n)末尾补零，用N=32点FFT分析p 这两个信号的幅频特性，观察幅频特性发生了什么变化？两种情况的FFT频谱还有相同之处吗？这些变化说明了什么？ (5

9、) 用FFT分别实现_a(n)（p8，q2）和_b(n)（a0.1，f0.0625）的16点循环卷积和线性卷积。n=0:15; p=8; q=2;_a=e_p(-(n-p).2/q); subplot(2,3,1); stem(n,_a,.); title(_a波形);_a=fft(_a,16); subplot(2,3,4); stem(abs(_a),.);title(_a(k)=FFT_a(n)的波形 ); A=1; f=0.0625; a=0.1;_b=e_p(-a_n)._sin(2_pi_f_n); subplot(2,3,2); stem(n,_b,.); title(_b波形)

10、; _b=fft(_b,16); subplot(2,3,5); stem(abs(_b),.);title(_b(k)=FFT_b(n)的波形 );实验过程与结果（含实验程序、运行的数据结果和图形）； 实验三 IIR数字滤波器的设计一、实验目的(1) 掌握双线性变换法及脉冲响应不变法设计IIR数字滤波器的具体设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR数字滤波器的计算机编程。(2) 观察双线性变换及脉冲响应不变法设计的滤波器的频域特性，了解双线性变换法及脉冲响应不变法的特点。 (3) 熟悉巴特沃思滤波器、切比雪夫滤波器和椭圆滤波器的频率特性。 二、实验内容

11、（1）P162 例4.4 设采样周期T=250ms(采样频率fs=4kHz)，分别用脉冲响应不变法和双线性变换法设计一个三阶巴特沃思低通滤波器，其3dB边界频率为fc=1kHz。脉冲响应不变法： fc=1000; fs=4000; OmegaC=2_pi_fc; B,A=butter(3, OmegaC,s); num1,den1=impinvar(B,A,fs); h1,w=freqz(num1,den1); f = w/pi_fs/2; plot(f,abs(h1); 双线性变换法： fc=1000; fs=4000;OmegaC=2_fs_tan(pi_fc/fs); B,A=butte

12、r(3, OmegaC,s); num2,den2=bilinear(B,A,fs); h2,w=freqz(num2,den2); f = w/pi_fs/2; plot(f,abs(h2); 同一图中画两条曲线： fc=1000; fs=4000; OmegaC=2_pi_fc; B,A=butter(3, OmegaC,s); num1,den1=impinvar(B,A,fs); h1,w=freqz(num1,den1); f = w/pi_fs/2; OmegaC=2_fs_tan(pi_fc/fs); B,A=butter(3, OmegaC,s); num2,den2=bili

13、near(B,A,fs); h2,w=freqz(num2,den2); f = w/pi_fs/2; plot(f,abs(h1),r-.); hold on; plot(f,abs(h2),g-);(选做) (2)fc=0.2kHz,d=1dB,fr=0.3kHz,At=25dB,T=1ms；分别用脉冲响应不变法及双线性变换法设计一巴特沃思数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。实验过程与结果（含实验程序、运行的数据结果和图形）； 实验四FIR数字滤波器的设计一、实验目的(1) 掌握用窗函数法，频率采样法及优化设计法设

14、计FIR滤波器的原理及方法，熟悉相应的计算机编程； (2) 熟悉线性相位FIR滤波器的幅频特性和相频特性；(3) 了解各种不同窗函数对滤波器性能的影响。 二、实验内容(1) 生成一个长度为20的矩形窗，画出其时域和幅频特性曲线。 n=0:1:19; N=20; (1:20)=1; H,w=freqz(,1); subplot(2,1,1); stem(n,) subplot(2,1,2); plot(w,abs(H); (2) 用矩形窗设计一个21阶的线性相位低通FIR数字滤波器，截止频率Wc0.25，求出滤波器系数，并绘出滤波器的幅频特性。修改程序，分别得到阶次为N41，61的滤波器，并显示

15、其各自的幅频曲线。a) 在上面所得的几幅图中，在截止频率两边可以观察到幅频响应的摆动行为。 请问波纹的数量与滤波器脉冲响应的长度之间有什么关系？ b) 最大波纹的高度与滤波器脉冲响应的长度之间有什么关系？ 实验过程与结果（含实验程序、运行的数据结果和图形）； 21阶的线性相位低通FIR数字滤波器： Wc=0.25_pi; N=21; M=(N-1)/2;位移量for n=0:(N-1)if (n= fi_(M)中间的点单独算hd(n+1)=Wc/pi;else hd(n+1)=sin(Wc_(n-M) /(pi_(n-M); end; end; =bo_car(N); 不同窗函数h=hd._;

16、 H,w=freqz(h,1); n=0:1:N-1; subplot(3,1,1); stem(n,h) subplot(3,1,2); plot(w,abs(H); subplot(3,1,3); plot(w,angle(H);41阶的线性相位低通FIR数字滤波器： Wc=0.25_pi; N=41; M=(N-1)/2;位移量for n=0:(N-1)if (n= fi_(M)中间的点单独算hd(n+1)=Wc/pi;else hd(n+1)=sin(Wc_(n-M) /(pi_(n-M); end; end; =bo_car(N); 不同窗函数h=hd._; H,w=freqz(h,

17、1); n=0:1:N-1; subplot(3,1,1); stem(n,h) subplot(3,1,2); plot(w,abs(H); subplot(3,1,3); plot(w,angle(H);61阶的线性相位低通FIR数字滤波器： Wc=0.25_pi; N=61; M=(N-1)/2;位移量for n=0:(N-1)if (n= fi_(M)中间的点单独算hd(n+1)=Wc/pi;else hd(n+1)=sin(Wc_(n-M) /(pi_(n-M); end; end; =bo_car(N); 不同窗函数h=hd._; H,w=freqz(h,1); n=0:1:N-1

18、; subplot(3,1,1); stem(n,h) subplot(3,1,2); plot(w,abs(H); subplot(3,1,3); plot(w,angle(H);数字信号处理实验小结及心得体会：通过这次实验，我对MATLAB语言有了一定的认识，虽然还不能完全用MATLAB独立编写程序，但对这种语言环境有了新的了解。我知道了一般的加减乘除在MATLAB中不同的意义。知道输入、输出语句怎么形成。通过快速傅里叶变换及其应用的实验，加深了我对FFT的理解，还有对各典型信号的频谱分析p ，改变参数后时域和幅频特性的变化。IIR数字滤波器的设计让我知道了巴特沃思滤波器和切比雪夫滤波器的频率特性，还有双线性变换及脉冲响应不变法设计的滤波器的频率特性。做这个实验的时候程序有点困难，很多细节问题不能考虑清楚，导致图形出不来。FIR数字滤波器的设计出来的是三种窗的图形，通过三种窗的比较，我了解了他们各自的特点，幅频和相频特性。我在这次实验中的收获很大，接触了很多新的知识，但在实验写程序时，我发现自己还有很多不足。很多程序写不完全。这是自己今后要加强的地方。作者继续努力。第 10 页 共 10 页