实验四 线性卷积验实验.doc

实验四 线性卷积实验 [实验目的] 熟悉并验证卷积的性质 利用卷积生成新的波形，建立波形间的联系 验证卷积定理 [实验原理] 信号的卷积是针对时域信号处理的一种分析方法。信号的卷积一般用于求取信号通过某系统后的响应。在信号与系统中，我们通常求取某系统的单位冲激响应，所求得的h(k)可作为系统的时域表征。任意系统的系统响应可用卷积的方法求得： 两个序列的线性卷积可经过下列步骤： 将x(n)和h(n)用x(m)和h(m)表示，并将h(m)进行翻转，形成h(-m)； ②将h(-m)移位n，得到h(n-m)。当n0时，序列右移；n0时，序列左移； ③将x(m)和h(n-m)相同m的序列值对应相乘后 求和，将以上所有对应点的乘积累加起来，可得到卷积结果y(n)。 [实验内容] MATLAB提供了一个内部函数conv来计算两个有限长序列的卷积。conv函数假定两个序列都从n=0开始。给出序列x=[3,11,7,0,-1,4,2]； h=[2,3,0,-5,2,1]；求两者的卷积y 。 解： x=[3,11,7,0,-1,4,2] x = 3 11 7 0 -1 4 2 h=[2,3,0,-5,2,1] h = 2 3 0 -5 2 1 y=conv(x,h) y = 6 31 47 6 -51 -5 41 18 -22 -3 8 2 2．将函数conv稍加扩展为函数conv_m，它可以对任意基底的序列求卷积。格式如下： function [y,ny]=conv_m(x1,x2) % 信号处理的改进卷积程序 % [y,ny]=conv_m(x1,x2) % [y,ny]=卷积结果 % [x,nx1]=第一个信号 % [h,nx2]=第二个信号 解： function [y,ny]=conv_m(x1,x2) x1=input(x1=); x2=input(x2=); N1=length(x1); M=length(x2); L=N1+M-1; ny=0:L-1 for(n=1:L) y(n)=0; for(m=1:M) k=n-m+1; if(k=1k=N1) y(n)=y(n)+x2(m)*x1(k); end end end y y1=conv(x1,x2); nx1=0:N1-1; nx2=0:M-1; subplot(221); stem(nx1,x1,.k);ylabel(x1(n));xlabel(n);grid on; title(x1); subplot(222); stem(nx2,x2,.k);ylabel(x2(n));xlabel(n);grid on; title(x2); subplot(223); stem(ny,y,.k);ylabel(y(n));xlabel(n);grid on; title(y); subplot(224); stem(y1);ylabel(y1);xlabel(n);grid on; title(y1); 根据线性卷积分析，参考程序流程为 3.对下面三个序列，用conv_m函数，验证卷积特性（交换律、结合律、分配 律、同一律） 交换律 结合律 分配律 同一律 其中：x1(n)=n[u(n+10)-u(n-20)] x2(n)= [u(n)-u(n-30)] 　　 x3(n)=(1.2)n[u(n+5)-u(n-10)]