系统的卷积响应.docxVIP

山东轻工业学院实验报告 课程名称：信号与系统教程 指导教师：贺爱玲 实验F1期：201L11.3 院(系)： 理学院 专业班级：信息09?2 实验地点：机电楼A座 学生姓名： 学号： 200903021076 同组人：± 实验项目名称： 实验二：系统的卷积响应 一、实验目的： 1、 已知有系统: y⑴ + 7y⑴ +12 曲)=4xf(t) +16x(0 求该系统的冲激相应、阶跃响应。 并求当 兀⑴=2eT乜⑴ 时的零状态响应。 2、 求以下卷积: ￡(r)*g(r) w(r)*厂⑴ €-3(~2)￡(/_2)乜卫一 0 2) 若需取时间间隔，全取为()?()() 1 二、实验器材： 安装有matlab软件的计算机一台。 三、实验内容： 1、求冲激相应，阶跃响应及零状态响应。 为求该系统的冲击响应与阶跃响应，现在m文件编辑窗口屮，写入一下程序: 代码如下： b二[4J6]; 圧[1,7,12]; subplot( 1,2,1) impulse(b,a) titleC冲激响应*) subplot(l,2,2) step(b,a) title邙介跃响应) 通过编译连接，可得到冲激相应打阶跃响应图形如下: 4 0 0.5 1 1.5 Time (sec) 冲激响应 5 3 3 5 2 5 ? ? 2 1 ① PB一-duj 5 o o 0.5 1 1.5 2 Time (sec) 阶跃响应 同样通过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、求以下各卷积: ￡ (t)* E (t) 1 0.5 10.5 1 0.5 0 ■ ■ ■ 0 0.5 1 1.5 2 0 0.5 1 1.5 2 f(t)=f1(t)*f2(t) 代码如下： 先编写一函数sconv用来实现连续信号的卷积: function ff,k]=sconv(fl ,f2,kl ,k2,p) f=conv(fl,f2); f=f*p; k0=kl(l)+k2(l); k3=length(f 1 )+length(f2)-2; k3二k():p:k()+k3*p; subplot(2,2,l) plot(kl,fl) axis(l0,2,-0.2,1.2J) titleCflCt)1) xlabcl(t) ylabelCfUt)*) subplot(2,2,2) plot(k2,f2) axis([0,2,-0.2,1.2]) titlc(T2(t)J xlabel(tr) ylabeRW1) subplot(2,2,3) plot(k3,f); axis([0,3,-0.2,2]) h=get(gca,position) h(3)=2.5*h(3); sct(gca,position,h) title(f(t)二fl(t)*f2(t)) xlabel(t) ylabelCfCt)1) 然后通过调用上编写的函数来实现求卷积运算。 p=0.001; kl=0:p:2; fl=heaviside(kl) k2=kl; f2=fl; [f,k]=sconv(fl ,f2,k 1 ,k2,p) ￡ (t)*r(t) f1(t) 1 0.5 0 2 f(t)=f1 (t)*f2(t) 0 0.5 1 t 其代码如下： 同样是通过调用sconv函数。 p=().()()l; kl=0:p:2; fl=heaviside(kl) k2=kl; f2=k2; [f,k]=sconv(fl ,f2,k 1 ,k2,p) eA(-3(t-2))* e (t-2)*gl(t-0.2) 卷积图形如下 f1 (t) f2(t) 0-0.2 0 0.2 0.4 0.6tf(t)=f1(t)*f2(t)t 0 -0.2 0 0.2 0.4 0.6 t f(t)=f1(t)*f2(t) t 其代码如下： 同样是通过调用sconv函数，其代码如下： p=0.001; kl=0:p:8; fl =exp(-3*(k 1 -2)).*heaviside(k 1 -2) k2=().3:p:0.7; f2=oncs( 1 ,length(k2)); [f,k]=sconv(fl ,f2,k 1 ,k2,p) 四、实验心得体会： 通过本次实验的学习，对matlab软件冇更深入的了解，把以前学过的matlab 识乂巩 固了一遍。 在求冲击响应，阶跃响应以及零状态响应的过程中，对其有了更深层次的了解。在求卷 积的过程中，强化了对卷积的原理的理解。 同时再次感受到matlab软件是一个很