连续时间信号卷积运算的MATLAB实现

1、连续时间信号卷积运算的MATLAB实现-、实验目的(1) 理解掌握卷积的概念及物理意义。(2) 理解单位冲激响应的概念及物理意义。二、实验原理根据前述知识，连续信号卷积元素按定义为f t =匚 t ” f2 t j； =J . f2 t 一 . d .-oO卷积计算可以通过信号分段求和来实现，即QO'、fi f2 tk :f t =匸 t “ f2 t 二匚.-oO如果只求当t(n为整数)时f t的值f n ，则由上式可得OOQOf nf1 k 二：f2 n. ： - k. ： - . ：f1 kL : f2 |jn - k - |kk =00cOZ人店心)十2( n心一 k心)ff式

2、中的k“：实际上就是连续信号fl t和f2 t经等时间间隔 厶均匀抽样的离散序列fl k，和f2 k人的卷积和。当厶足够小时，f n'：就是卷积积分的结 果一一连续时间信号f t的较好的数值近似。三、实验内容1、实验参考程序以下是MATLAB实现连续信号卷积的通用函数sconv()。fun ctio nf,k=sco nv(f1,f2,k1,k2,p)f=con v(f1,f2);f=f*p;k0=k1(1)+k2(1);k3=le ngth(f1)+le ngth(f2)-2;k=0:p:k3*p;subplot(2,2,1)plot(k1,f1)title( 'f1(t)&

3、#39;)xlabel( 't')ylabel( 'f1(t)')subplot(2,2,2)plot(k2,f2)title( 'f2(t)')xlabel( 't')ylabel( 'f2(t)')subplot(2,2,3)plot(k,f);h=get(gca, 'positi on' )h(3)=2.5*h(3);set(gca,'positi on',h)title( 'f(t)=f1(t)*f2(t)')xlabel( 't')ylabe

4、l( 'f(t)')end已知两连续时间信号如图所示，使用MATLAB求f t = f1 t " f2 t，并绘出f t的时域波形图。实现上述过程的MATLAB命令如下：p=0.5;k1=0:p:2;f1=0.5*k1;k2=k1;f2=f1;f,k=sc on v(f1,f2,k1,k2,p)运行得出结果。上述命令绘制的波形图在图中示出：图1.1以上给出了时间间隔为 p=0.5时的处理效果。而图1.2给出了抽样时间间隔 p=0.05时的处理 结果。可见，当抽样时间p足够小时，函数sconv()的计算结果就是连续时间卷积 f t "1 t “笃t的较好数值近

5、似。图1.22、习题练习(P34-1.11(a),(b)已知两连续时间信号如图所示， 使用MATLAB求f t = fl t " f2 t ,并绘出f t的 时域波形图。(a)输入以下 MATLAB命令： p=0.001;k1=0:p:2;f仁heaviside(k1)-heaviside(k1-1); k2=0:p:2;f2=0.5*k2;f,k=sco nv(f1,f2,k1,k2,p运行得出结果:*1I 1|HLJ_11_1IO.SJt1 1.6 :-o.$1 15 2 t(b)输入以下 MATLAB命令: p=0.001;k1=-1:p:1;f仁-abs(k1)+1;k2=0:p:2;f2=heaviside(k1)-heaviside(k1-1);f,k=sc on v(f1,f2,k1,k2,p)运行得出结果：