matlab实验-卷积.pptVIP

实验二 用matlab进行信号卷积仿真 离散时间信号的卷积 我们试举一例来看conv的功能，已知序列f1(k)和f2(k)如下所示： f1(k)=1,(0≤k≤2) f2(k)=k,(0≤k≤3) 则调用conv( )函数求上述两序列的卷积和的MATLAB命令为： f1=ones(1,3); f2=0:3; f=conv(f1,f2) 运行结果为：f＝0 1 3 6 5 3 由这个例子可以看出，函数conv( )不需要给定序列f1(k)和f2(k)非零样值点的时间序号，也不返回序列f(k)=f1(k)*f2(k)的非零样值点的时间序号。因此，要正确地标识出函数conv( )的计算结果向量f，我们还必须构造序列f1(k)、f2(k)及f (k)的对应序号向量。我们可以看出f(k)的序号向量k是由序列f1(k)和f2(k)的非零样值点的起始序号及他们的时域宽度决定的 在此，我们自己可以构造一个实用一点的函数dconv( )，该函数实现的功能为： 可以在计算出卷积f (k)的同时，还可以绘出序列f1(k)、f2(k)及f (k)在时域内的波形图。并返回f (k)的非零样值点的对应向量 function [f,k]=dconv(f1,f2,k1,k2) %the function of compute f=f1*f2 %f: 卷积序列f(k)对应的非零样值向量 %k: 序列