第3章时域离散信号与系统幻灯片.pptVIP

5)线性卷积的MATLAB实现 MATLAB中可用conv函数求序列x(n)与h(n)的卷积y(n)，其调用格式为y=conv(x,h)。 例3-4 用MATLAB计算例3-3中两序列的卷积。x=[1,2,3]; h=[3,4,5,6,7]; %设置指标参数 y=conv (x,h) %计算线性卷积 程序运行结果与例3-3的卷积结果一致。需注意conv函数仅在序列起点为零点时适用。 一个差分方程对应一个滤波系统，因此可用MATLAB中的滤波函数filter求解差分方程（适用于零状态解），用函数impz求解系统的单位脉冲响应。 例3-7 某LTI系统差分方程为y(n)=0.5y(n-1)+1.5x(n)，且x(n)=δ(n)+δ(n-1)，初始条件为n0时，y(n)=0，用MATLAB求解系统的输出y(n)以及单位脉冲响应h(n)。 b=1.5; a=[1,-0.5]; %差分方程的系数组 x=[1,1,zeros(1,18)]; %输入序列(共输入20个点) y=filter(b,a,x); %使用函数filter求输出序列y(n) h=impz(b,a,20); %使用函数impz求单位采样响应h(n) 设系统的输入为 x(n)，先将x(n)表示成单位采样序列的移位加权和形式， 则系统输出为y=T[x(n)],所以 根据线性叠加原理， 因为T[δ(n)]=h(n),根据时不变性，T[δ(n-m)]=h(n-m), (3-39) 所以 T[a1·x1(n)+a2·x2(n)]= a1T[x1(n)]+a2T[x2(n)] 2.线性卷积计算 1)图解法 (1)将x(n),h(n)用x(m),h(m)表示，并将h(m)进行翻转，形成h(-m) ； (2)将h(-m)移位n，得到h(n-m) ； (3)将x(m)和h(n-m)相同m的序列值对应相乘后再相加。 注意！若两个序列长度为N、M，卷积后的序列长度为N＋M－1。 例3-2设x(n)=R4(n)，h(n)=R4(n)，求y(n)=x(n)*h(n)。 解 按照(3-39)式， m h ( 0 - m ) 0 1 2 3 1 m h ( 2 - m ) 0 1 2 3 1 m h ( 3 - m ) 0 1 2 3 1 m 0 1 2 3 1 h ( 1 - m ) 0 1 2 3 1 m x ( m ) 0 1 2 3 1 m h ( m ) m h ( 4 - m ) 0 1 2 3 m h ( 5 - m ) 0 1 2 3 m h ( 6 - m ) 0 1 2 3 m h ( 7 - m ) 0 1 2 3 1 1 1 1 4 5 4 4 5 6 7 4 5 6 0 1 2 3 1 m x ( m ) 两序列长度分别为N和M，则线性卷积的长度为N+M-1 n x ( n )* h ( n ) 0 1 2 3 4 5 6 7 1 2 3 4 5 图3-19 例3-2线性卷积 作图法 m -3 -2 -1 0 1 2 3 4 5 6 x(m) 1 1 1 1 h(m) 1 1 1 1 h(-m) 1 1 1 1 h(1-m) 1 1 1 1 h(2-m) 1 1 1 1 h(3-m) 1 1 1 1 h(4-m) 1 1 1 1 h(5-m) 1 1 1 1 h(6-m) 1 1 1 1 2)列表法 上式中矩形序列长度为4，求解上式主要是根据矩形序列的非零值区间确定求和的上、下限，R4(m)的非零值区间为：0≤m≤3， R4(n-m)的非零值区间为：0≤n-m≤3，其乘积值的非零区间，要求m同时满足下面两个不等式： 0≤m≤3且 n-3≤m≤n 因此， 3)解析法 例3-2 4)矩阵法 对于y(n)=x(n)*h(n)(设两序列长度分别为N和M)，可通过矩阵乘法Y=XH来求得，其中左乘矩阵X为由序列x(n)构成