南京邮电大学
实 验 报 告
实验名称 熟悉MATLAB环境
快速傅里叶变换(FFT)及其应用
IIR数字滤波器的设计 FIR数字滤波器的设计
课程名称 数字信号处理A
班级学号_ 12006311____ 姓 名_______张文欣_____________
开课时间 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。 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')
4A200123n10C500123n40E2000123n5000G00123n
(2) 用MATLAB实现下列序列： a) x(n)?0.8n 0?n?15 n=0:1:15; x1=0.8.^n; stem(n,x1) xlabel('n') ylabel('x(n)')
10B500123n0D-1-20123n10.500123nF
10.90.80.70.6x(n)2(a)0.50.40.30.20.1005n1015title('2(a)')
b) x(n)?e(0.2?3j)n 0?n?15 n=0:1:15; i=sqrt(-1); a = 0.2+3*i; x2=exp(a*n); figure
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)虚部')
12108106420-2-10-4-6-805n1015-15-2005n1015520152(b)x(n)实部x(n)虚部0-5
c) x(n)?3cos(0.125?n?0.2?)?2sin(0.25?n?0.1?) 0?n?15 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)') title('2(c)')
4321x(n)2(c)0-1-2-3-4-505n1015
(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)')
x(t)?sin(2?t) 0?t?10s