中国地质大学（武汉）
数字信号处理上机实习
学生姓名： 班 级：071132 学 号：2013100
指导老师：王晓莉
题目一 离散卷积计算
一、实验题目
设线性时不变（LTI）系统的冲激响应为h(n)，输入序列为x(n) 1、h(n)=(0.8)n，0≤n≤4； x(n)=u(n)-u(n-4) 2、h(n)=(0.8)nu(n), x(n)=u(n)-u(n-4) 3、h(n)=(0.8)nu(n), x(n)=u(n)
求以上三种情况下系统的输出y(n)，显示输入和输出波形。
二、实验目的
1.理解和掌握离散卷积计算；
2.学习如何用Mtalab实现离散卷积计算。
三、算法设计
离散卷积定义为：
y(n)?k????x(k)h(n?k)
n1、h(n)=(0.8)n，0?n?4，x(n)=u(n)-u(n-4)， y(n)?x(n)?h(n)?m????x(m)h(n?m)
?(a) 当n?0 时，y(n)?0；
(b) 当0?n?3时,y(n)??(0.8)n ；
m?0n(c) 当4?n?7时，y(n)?m?n?3?4(0.8)n；
(d) 当n?7时，y(n)?0；
2、h(n)?(0.8)u(n)，x(n)?u(n)-u(n?4)，y(n)?x(n)?h(n)?(a) 当n?0 时，y(n)?0；
nm????x(m)h(n?m)
?(b) 当0?n?3时,y(n)??(0.8)n；
m?0n
(c) 当4?n?20时，y(n)?m?n?320?n(0.8)n；
(d) 当21?n?23时，y(n)?(e) 当n?23时，y(n)?0；
m?n?3?(0.8)n；
3、h(n)?(0.8)u(n)，x(n)?u(n)，y(n)?x(n)?h(n)?(a) 当n?0 时，y(n)?0；
(b) 当0?n?70时,y(n)??(0.8)n；
m?0nnm????x(m)h(n?m)
?
(c) 当71?n?140时，y(n)?(d) 当n?140时，y(n)?0；
m?n?70?70(0.8)n；
四、程序分析
所用到的函数：
（1）y=conv（x.,h）：卷积运算函数，计算y(n)?x(n)*h(n); （2）n1=0:4：n1取0~4；
（3）subplot(m,n,p):subplot()函数是将多个图画到一个平面上的工具。其中，m表示是图排成m行，n表示图排成n列，也就是整个figure中有n个图是排成一行的，一共m行，如果m=2就是表示2行图。p表示图所在的位置，p=1表示从左到右从上到下的第一个位置。
（4）title(‘content ’):title()函数的功能是为当前坐标系添加标题“content ”。
五、程序设计
n=0:4; h=0.8.^n;