系统的卷积响应.docVIP

山东轻工业学院实验报告 课程名称：信号与系统教程 指导教师：贺爱玲 实验日期：2011.11.3___ 院（系）： 理学院 专业班级：信息09-2 实验地点：机电楼A座 学生姓名： 学号： 200903021076 同组人：无________ 实验项目名称： 实验二：系统的卷积响应 ____________________________ 一、实验目的： 已知有系统： 求该系统的冲激相应、阶跃响应。 并求当 时的零状态响应。 2、求以下卷积： 若需取时间间隔，全取为 0.001 二、实验器材： 安装有matlab软件的计算机一台。 三、实验内容： 1、求冲激相应，阶跃响应及零状态响应。 为求该系统的冲击响应与阶跃响应，现在m文件编辑窗口中，写入一下程序： 代码如下： b=[4,16]; a=[1,7,12]; subplot(1,2,1) impulse(b,a) title(冲激响应) subplot(1,2,2) step(b,a) title(阶跃响应) 通过编译连接，可得到冲激相应与阶跃响应图形如下： 同样通过matlab程序，可得到零状态响应图形如下，其源代码附于图形后。 零状态响应图形： 代码： yzs=dsolve(D2y+7*Dy+12*y=-8*exp(-t)+32*exp(-t),y(0)=0,Dy(0)=0); ezplot(yzs,[0 8]); 2、求以下各卷积： （1）ε(t)*ε(t) 代码如下： 先编写一函数sconv用来实现连续信号的卷积： function[f,k]=sconv(f1,f2,k1,k2,p) f=conv(f1,f2); f=f*p; k0=k1(1)+k2(1); k3=length(f1)+length(f2)-2; k3=k0:p:k0+k3*p; subplot(2,2,1) plot(k1,f1) axis([0,2,-0.2,1.2]) title(f1(t)) xlabel(t) ylabel(f1(t)) subplot(2,2,2) plot(k2,f2) axis([0,2,-0.2,1.2]) title(f2(t)) xlabel(t) ylabel(f2(t)) subplot(2,2,3) plot(k3,f); axis([0,3,-0.2,2]) h=get(gca,position) h(3)=2.5*h(3); set(gca,position,h) title(f(t)=f1(t)*f2(t)) xlabel(t) ylabel(f(t)) 然后通过调用上编写的函数来实现求卷积运算。 p=0.001; k1=0:p:2; f1=heaviside(kl) k2=k1; f2=f1; [f,k]=sconv(f1,f2,k1,k2,p) （2）ε(t)*r(t) 其代码如下： 同样是通过调用sconv函数。 p=0.001; k1=0:p:2; f1=heaviside(kl) k2=k1; f2=k2; [f,k]=sconv(f1,f2,k1,k2,p) （3）e^(-3(t-2))* ε(t-2)*g1(t-0.2) 卷积图形如下 其代码如下： 同样是通过调用sconv函数，其代码如下： p=0.001; k1=0:p:8; f1=exp(-3*(k1-2)).*heaviside(k1-2) k2=-0.3:p:0.7; f2=ones(1,length(k2)); [f,k]=sconv(f1,f2,k1,k2,p) 四、实验心得体会： 通过本次实验的学习，对matlab软件有更深入的了解，把以前学过的matlab知识又巩固了一遍。 在求冲击响应，阶跃响应以及零状态响应的过程中，对其有了更深层次的了解。在求卷积的过程中，强化了对卷积的原理的理解。 同时再次感受到matlab软件是一个很强大的数学应用软件。它对于我们学好我们的专业知识将提供很大帮助，因此，我们一定要借此机会，好好学习它。