DSP实验报告 经容-.doc

南京邮电大学
实验报告
实验名称熟悉MATLAB环境
快速傅里叶变换(FFT)及其应用
IIR数字滤波器的设计
FIR数字滤波器的设计
课程名称数字信号处理A
班级学号______12006903___________________
姓名_______经容__________________
开课时间 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。
(2) 用MATLAB实现下列序列:
a)
b)
c)
(4) 绘出下列时间函数的图形,对x轴、y轴以及图形上方均须加上适当的标注:
a)
b)
(6) 给定一因果系统,求出并绘制H(z)的幅频响应和相频响应。
(7) 计算序列{8 -2 -1 2 3}和序列{2 3 -1 -3}的离散卷积,并作图表示卷积结果。
(8) 求以下差分方程所描述系统的单位脉冲响应h(n),
每一小题均给出实验过程与结果(含实验程序、运行的数据结果和图形);
实验程序:
%(1)提示
clear
n = 0:1:3;
A=[1 2 3 4];
subplot(3,3,1)
stem(n,A)
xlabel('n')
ylabel('A')
B=[3,4,5,6];
subplot(3,3,2)
stem(n,B)
xlabel('n')
ylabel('B')
C=A+B;
subplot(3,3,3)
stem(n,C)
xlabel('n')
ylabel('C')
D=A-B;
subplot(3,3,4)
stem(n,D)
xlabel('n')
ylabel('D')
E=A.*B;
subplot(3,3,5)
stem(n,E)
xlabel('n')
ylabel('E')
F=A./B;
subplot(3,3,6)
stem(n,F)
xlabel('n')
ylabel('F')
G=A.^B;
subplot(3,3,7)
stem(n,G)
xlabel('n')
ylabel('G')
%(2)提示
%(a)
n=0:1:15;
x1=.^n;
stem(n,x1)
xlabel('n')
ylabel('x(n)')
title('2(a)')
n=0:1:15;
i=sqrt(-1);
a = +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)虚部')
title('2(b)')
n=0:1:15;
x3=3*cos(*pi*n+*pi) + 2*sin(*pi*n+*pi);
stem(n,x3)
xlabel('n')
ylabel('x(n)')
title('2(c)')
(4)提示
% (a)
t=0::10;
x=sin(2*pi*t);
plot(t,x,'r-')
xlabel('t'),ylabel('x(t)'),title('sin(2\pit)')
% (a)
t=0::10;
x=sin(pi*t)* cos(100*pi*t);
plot(t,x,'r-')
xlabel('t'),ylabel('x(t)'),title('cos(100\pit)')
(6)
%a
k=512;
num=[ - ];
den=[1 - -];
w=0:pi/k:pi;
h=freqz(num, den,w);
subplot(2,2,1);
plot(w/pi,real(h));grid
title('实部')
xlabel('\omega/ \pi');ylabel('幅度')
subplot(2,2,2);
pl