信号实验三.doc

等级： 实验报告课程名称： MATLAB应用基础 姓 名： 游 金 涛 学 号： 200808002304 院 （系）： 计算机与通信工程系 专业、班级 ： 　 通信0803班 任课教师： 李 玲 香 精品.实验三、 连续时间系统的时域分析一、实验目的 1.深刻理解卷积运算，掌握离散线性卷积、连续线性卷积的计算方法； 2.加深对线性时不变系统中零状态响应概念的理解，掌握其求解方法；3.掌握给定连续系统的冲激响应和阶跃响应 二、实验原理 1. 线性时不变 (LTI) 连续时间系统用常系数线性微分方程进行描述，系统的零状态响应就是在系统初始状态为零条件下微分方程的解MATLAB控制系统工具箱提供了一个lsim函数来求解连续时间系统的零状态响应设系统方程为： 该方程左边、右边的系数向量分别为，，所对应的系统模型sys可借助MATLAB中的tf函数得到：sys = tf(b, a) . 这样，系统的零状态响应为：y = lsim(sys, f, t) ，其中f是输入信号向量，t是与f对应的时间变量。

2. 连续系统的冲激响应、阶跃响应分别是输入信号为和所对应的零状态响应MATLAB控制系统工具箱专门提供了两个函数求解连续系统的冲激响应和阶跃响应 冲激响应：y = impulse(sys, t) ;阶跃响应：y = step(sys, t) .其中sys, t的含义同上精品. 3. 卷积是信号与系统中一个最基本、也是最重要的概念之一在时域中，对于LTI 连续时间系统，其零状态响应等于输入信号与系统冲激响应的卷积；而利用卷积定理，这种关系又对应频域中的乘积我们可以用离散卷积来代替连续卷积，只要取样时间间隔足够小时，就可得到满意的效果 MATLAB信号处理工具箱提供了一个计算两个离散序列卷积和的函数conv设向量a、b代表待卷积的两个序列，则c = conv(a, b)就是a与b卷积后得到的新序列 我们知道两个序列卷积以后，一般而言所得新序列的时间范围、序列长度都会发生变化例如设f1(n)长度为5，－3≤n≤1；f2(n)长度为7，2≤n≤8；则卷积后得到的新序列长度为11，－1≤n≤9但是用conv函数求出卷积后没有给出新序列所对应的时间变量为此，我们在下面的程序示例中给出了一个函数文件dconv，它在完成conv函数功能的同时，还产生了一个对应新序列的时间变量。

4. 对于连续卷积 令（为整数），则 (*)由(*)式，连续卷积积分可由离散卷积和近似代替，只要取样时间间隔足够小，就可以得到高精度卷积积分的数值计算在示例3中给出了一个函数文件cconv来完成该功能精品.三、实验内容与步骤 1. 已知系统的微分方程为：计算信号通过系统的零状态响应、冲激响应和阶跃响应，并画出相应的图形 程序代码及截图如下：a = [1 2 2];b = [1 ];sys = tf(b, a);td = 0.0001;t = 0 : td : 10;f = exp(-t);y = lsim(sys, f, t);subplot(131);plot(t,y);subplot(132);impulse(sys,t);subplot(133);step(sys,t);2. 用MATLAB计算如下两序列的卷积和，绘出它们的时域波形 精品. 程序代码及截图如下：f1=[ 1 2 1];n1=[0 1 2];f2=ones(1,4);n2=-2:1;f=conv(f1,f2)n0=n1(1)+n2(1);n3=length(n1)+length(n2)-1;n=n0:n0+n3-1;subplot(131); stem(n1,f1);subplot(132);stem(n2,f2);subplot(133);stem(n,f); 运行结果：f =1 3 4 4 3 13.编程实现如下图所示的两个波形；并利用cconv函数计算这两个信号的卷积、画出卷积后的波形。

精品.程序代码及截图如下：1. %函数文件dconv.mfunction n=dconv(n1,n2,p)n0=n1(1)+n2(1);n3=length(n1)+length(n2)-2;n=n0:p:n0+n3*p;2. %函数文件cconv.mfunction [f, k]=cconv(f1,f2,k1,k2,p)f1=2;f2=1;k1=-1:0.01:1;k2=-2:0.01:2;p=0.01;f=conv(f1,f2); f=f*p;k0=k1(1)+k2(1); k3=length(f1)+length(f2)-2; k=k0:p:k0+k3*p; k=dconv(k1,k2,p)subplot(3,3,1)plot(k1,f1) title('f1(t)')xlabel('t')ylabel('f1(t)')subplot(3,3,4)plot(k2,f2) 精品.title('f2(t)')xlabel('t')ylabel('f2(t)')subplot(3,3,7)plot(k,f); 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)') 四、思考题 1. 在“实验内容与步骤”（1），零状态响应和阶跃响应是否相同？为什么？ 答：不同；阶跃响应是阶跃信号作用系统的零状态响应。

2. 两序列进行卷积后得到新的序列，说明新序列在时域长度、时域区间上与原来两序列的关系答：新序列的时间范围、序列长度都会发生变化例如f1(n)长度为5，－3≤n≤1；f2(n)长度为7，2≤n≤8；则卷积后得到的新序列长度为11，区间－1≤n≤9精品.精品.如有侵权请联系告知删除，感谢你们的配合！精品。