南京邮电大学
实 验 报 告
实验名称  实验一熟悉MATLAB环境
实验二快速傅里叶变换(FFT)及其应用
实验三IIR数字滤波器的设计
实验四FIR数字滤波器的设计
课程名称      数字信号处理A     
班级学号__________________
  姓    名______ _______________
开课时间 2010 /2011 学年，  第 2学期
                    实验一熟悉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。
输入代码如下：
A=[1 2 3 4];
B=[3,4,5,6];
C=A+B;
stem(C)
C运行图片：
输入代码如下：
A=[1 2 3 4];
B=[3,4,5,6];
D=A-B
stem(D)
D运行图片：
输入代码如下：
A=[1 2 3 4];
B=[3,4,5,6];
E=A.*B
stem(E)
E运行图片：
输入代码如下：
A=[1 2 3 4];
B=[3,4,5,6];
F=A./B
stem(F)
F运行图片：
输入代码如下：
A=[1 2 3 4];
B=[3,4,5,6];
G=A.^B
stem(G)
G运行图片：
(2) 用MATLAB实现下列序列：
a)
输入代码如下：
n=0:1:15;
x1=0.8.^n
stem(x1)
运行图片：
b)
输入代码如下：
a=(0.2+3*i)*n;
matlab求傅里叶变换
x2=exp(a)
stem(x2)
运行图片：
c) 
输入代码如下：
x3=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi)
stem(x3)
运行图片：
d) 将c)中的x(n)扩展为以16为周期的函数，绘出四个周期。
输入代码如下：
x3=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi)
x4=[x3 x3 x3 x3];
stem(x4)
运行图片：
(4) 绘出下列时间函数的图形，对x轴、y轴以及图形上方均须加上适当的标注：
a)
输入代码如下：
t1=0:0.001:10;
x1=sin(2*pi*t1);
subplot(2,1,1)
plot(t1,x1,'r-')
xlabel('t'),ylabel('x(t)'),title('sin')
b)
输入代码如下：
t2=0:0.001:4;
x2=cos(100*pi*t2).*sin(pi*t2);
subplot(2,1,2)
plot(t2,x2)
xlabel('t'),ylabel('x(t)'),title('cossin')
运行图片：
(6) 给定一因果系统，求出并绘制H(z)的幅频响应和相频响应。
输入代码如下：
k=256;
num=[1 1.414 1];
den=[1 -0.67 0.9];
w=0:pi/k:pi;
h=freqz(num,den,w);
subplot(2,2,1);
plot(w/pi,abs(h));grid;
title('幅度谱');
xlabel('\omega/ \pi');ylabel('幅度')
subplot(2,2,2);
plot(w/pi,angle(h));grid
title('相位谱');
xlabel('\omega/ \pi');ylabel('幅度');
运行图片：
(7) 计算序列{8 -2 -1 2 3}和序列{2 3 -1 -3}的离散卷积，并作图表示卷积结果。

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。