信号系统实验三.docVIP

信号与系统实验 实验三：信号的卷积 （第三次实验） 【实验目的】 1. 理解卷积的物理意义； 2. 掌握运用计算机进行卷积运算的原理和方法；3. 熟悉卷积运算函数conv的应用； 实验原理说明 1.卷积的定义 2.卷积的计算 由于计算机技术的发展，通过编程的方法来计算卷积积分和卷积和已经不再是冗繁的工作，并可以获得足够的精度，因此信号的时域卷积分析法在系统分析中得到了广泛的应用。 卷积积分的数值运算可以应用信号的分段求和来实现，即： 数值运算只求当时的信号值，则由上式可以得到： 上式中实际上就是连续信号等间隔均匀抽样的离散序列的卷积和，当足够小的时候就是信号卷积积分的数值近似。因此，在利用计算机计算两信号卷积积分时，实质上是先将其转化为离散序列，再利用离散卷积和计算原理来计算。 3.卷积的应用 1 求解系统响应 卷积是信号与系统时域分析的基本手段，主要应用于求解系统响应，已知一LTI系统的单位冲激响应和系统激励信号则系统响应为激励与单位冲激响应的卷积。 需要注意的是利用卷积分析方法求得的系统响应为零状态响应。 2 相关性分析 相关函数是描述两个信号相似程度的量。两信号之间的相关函数一般称之为互相关函数或者互关函数，定义如下： 若是同一信号，此时相关函数称为自相关函数或者自关函数： 对于相关函数与卷积运算有着密切的联系，由卷积公式与相关函数比较得： 可见，由第二个信号反转再与第一个信号卷积即得到两信号的相关函数。 4.涉及的Matlab函数 4.1 conv函数 格式w = conv(u,v)，实现两有限序列u，v的卷积运算。长度为长度和减一。 给定如下因果线性时不变系统： y[n]+0.71y[n-1]-0.46y[n-2]-0.62y[n-3=0.9x[n]-0.45x[n-1]+0.35x[n-2]+0.002x[n-3] (1)不用impz函数，使用filter命令，求出以上系统的单位冲激响应h[n]的前20个样本; 图像如下： (2)得到h[n]后，给定x[n],计算卷积输出y[n]；并用滤波器h[n]对输入x[n]滤波，求得y1[n]; 图像如下： y[n]和)y1[n]有差别吗？为什么要对x[n]进行补零得到的x1[n]来作为输入来产生y1[n]？ (4)思考：设计实验，证明下列结论 单位冲激信号卷积：? ? 图像如下： 代码如下： n = [0: 30]; t0 = 5; d = (n - t0 == 0); f = sin(n); f1 = conv(d, f); subplot(3, 1, 1); f1 = f1(1:31); stem(n, f1); title((t-t0)*f(t)); grid; subplot(3, 1, 2); f = [zeros(1, t0) f]; f = f(1:31); stem(n, f); title(f(t-t0)); grid; subplot(3, 1, 3); stem(n, f - f1); title((t-t0)*f(t) - f(t-t0)); grid; 图像如下： 卷积交换律 代码如下： n = 1:50; f1 = (-1).^n; f2 = cos(n); f3 = sin(n); y1 = conv(f1, (f2 + f3)); y1 = y1(1:50); y2 = conv(f1, f2) + conv(f1, f3); y2 = y2(1:50); subplot(3, 1, 1); stem(n, y1); title(f1(t)*[f2(t)+f3(t)]); grid; subplot(3, 1, 2); stem(n, y2); title(f1(t)*f2(t)+f1(t)*f3(t)); grid; subplot(3, 1, 3); y3 = (y1 - y2) 10 ^ -14; stem(n, y3); title(f1(t)*[f2(t)+f3(t)]-[f1(t)*f2(t)+f1(t)*f3(t)]); grid; 图像如下： ③卷积分配律 代码如下： n = 0:30; f1 = sin(n); f2 = cos(n); y1 = conv(f1, f2); y1 = y1(1:31); y2 = conv(f2, f1); y2 = y2(1:31); subplot(3, 1, 1); stem(n, y1); title(f1(t)*f2(t)); grid; subplot(3, 1, 2); stem(n, y2); title(f2(t)*f1(t)); grid; subplot(3, 1, 3); y3 = (y1 - y2) 10 ^ -14; stem(n