DSP实验报告_封面及要求.docVIP

南京邮电大学 实 验 报 告 实验名称 熟悉MATLAB环境 快速傅里叶变换(FFT)及其应用 IIR数字滤波器的设计 FIR数字滤波器的设计 课程名称 数字信号处理A 班级学号_________________________ 姓 名_________________________ 开课时间 2014/2015学年， 第 二 学期 实验一 熟悉MATLAB环境 实验目的 （1）熟悉MATLAB的主要操作命令。 （2）学会简单的矩阵输入和数据读写。 （3）掌握简单的绘图命令。 （4）用MATLAB编程并学会创建函数。 （5）观察离散系统的频率响应。 实验内容 (1) 数组的加、减、乘、除和乘方运算。输入A=[1 2 3 4]，B=[3,4,5,6]，求 C=A+B， D=A-B，E=A.*B,F=A./B,G=A.^B 。并用stem语句画出A、B、C、D、 E、F、G。 clear n = 0:1:3; A=[1 2 3 4]; subplot(4,2,1) stem(n,A) xlabel(n) ylabel(A) B=[3,4,5,6]; subplot(4,2,2) stem(n,B) xlabel(n) ylabel(B) C=A+B; subplot(4,2,3) stem(n,C) xlabel(n) ylabel(C) D=A-B subplot(4,2,4) stem(n,D) xlabel(n) ylabel(D) E=A.*B subplot(4,2,5) stem(n,E) xlabel(n) ylabel(E) F=A./B subplot(4,2,6) stem(n,F) xlabel(n) ylabel(F) G=A.^B subplot(4,2,7) stem(n,G) xlabel(n) ylabel(G) (2) 用MATLAB实现下列序列： a) n=0:1:15; x1=0.8.^n; stem(n,x1) xlabel(n) ylabel(x(n)) title(2(a)) b) n=0:1:15; x1=0.8.^n; stem(n,x1) xlabel(n) ylabel(x(n)) title(2(a)) figure n=0:1:15; i=sqrt(-1); a = 0.2+3*i; x2=exp(a*n); stem(n,x2) subplot(1,2,1) stem(n,real(x2)) xlabel(n) ylabel(x(n)实部) subplot(1,2,2) stem(n,imag(x2)) xlabel(n) ylabel(x(n)虚部) n=0:1:15; x3=3*cos(0.125*pi*n+0.2*pi) + 2*sin(0.25*pi*n+0.1*pi); stem(n,x3) xlabel(n) ylabel(x(n)) c) n=0:1:15; x3=3*cos(0.125*pi*n+0.2*pi) + 2*sin(0.25*pi*n+0.1*pi); stem(n,x3) xlabel(n) ylabel(x(n)) (4) 绘出下列时间函数的图形，对x轴、y轴以及图形上方均须加上适当的标注： a) t=0:0.001:10; x=sin(2*pi*t); plot(t,x,r-) xlabel(t),ylabel(x(t)),title(sin(2\pit)) b) t=0:0.001:10; x=cos(100*pi*t)*sin(pi*t) plot(t,x,r-) xlabel(t),ylabel(x(t)),title(sin(2\pit)) (6) 给定一因果系统，求出并绘制H(z)的幅频响应和相频响应。 (7) 计算序列{8 -2 -1 2 3}和序列{2 3 -1 -3}的离散卷积，并作图表示卷积结果。 x=[8 -2 -1 2 3]; kx=0:4; h=[2 3 -1 -3]; kh=0:3; y=conv(x,h); k=kx(2)+kh(1):kx(end)+kh(end); stem(k,y); xlabel(n);ylabel(y(n)) (8) 求以下差分方程所描述系统的单位脉冲响应h(n), 每一小题均给出实验过程与结果（含实验程序、运行的数据结果和图形）； 实验二 快速傅里叶变换(FFT)及其应用 实验目的 (1) 在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉MATLAB中的有关函数。 (2) 应用FFT对典型信号进行频谱分析。 (3) 了解应用FFT进行信号频谱分析过程中可能出现的问题，以便在实