连续时间的时域分析.docVIP

实验二 连续时间系统的时域分析 一、实验目的： 1、掌握用Matlab进行卷积运算的数值方法和解析方法,加深对卷积积分的理解。 2、学习利用Matlab实现LTI系统的冲激响应、阶跃响应和零状态响应 二、原理说明： 1、卷积的计算 借助MATLAB的内部函数conv()可以很容易地完成两个信号的卷积积分运算。其语法为： y = conv(x,h)。其中x和h分别是两个作卷积运算的信号，y为卷积结果。 为了正确地运用这个函数计算卷积，这里有必要对conv(x,h)做一个详细说明。conv(x,h)函数实际上是完成两个多项式的乘法运算。例如，两个多项式p1和p2分别为： 这两个多项式在MATLAB中是用它们的系数构成一个行向量来表示的，如果用x来表示多项式p1，h表示多项式p2，则x和h分别为： x=[1 2 3 4]，h=[4 3 2 1] 在MATLAB命令窗口依次输入 x = [1 2 3 4]; h = [4 3 2 1]; y=conv(x,h) 在屏幕上得到的显示结果： y = 4 11 20 30 20 11 4 这表明，多项式p1和p2的卷积为： 正如前所述，用MATLAB处理连续时间信号时，独立时间变量t的变化步长应该是很小的，假定用符号dt表示时间变化步长，那么，用函数conv（）作两个信号的卷积积分时，应该在这个函数之前乘以时间步长方能得到正确的结果。也就是说，正确的语句形式应为：y=dt*conv(x,h)。 对于定义在不同时间段的两个时限信号x（t），，和。如果用y(t)来表示它们的卷积结果，则y(t)的持续时间范围要比x(t)或h(t)要长，其时间范围为t0+t2 ≤ t ≤ t1+t3。这个特点很重要，利用这个特点，在处理信号在时间上的位置时，可以很容易地将信号的函数值与时间轴的位置和长度关系保持一致性。 ①由于MATLAB运算的特点，两个连续信号f1(t)、f2(t)的卷积f(t)=f1(t)*f2(t)，用MATLAB 实现的过程应为： （1）将连续信号f1(t)、f2(t)以时间间隔Δ进行取样，得离散序列f1(k )、f2(k )； （2）构造f1(k )、f2(k )与相对应的时间向量k1和k2； （3）调用conv()函数计算卷积积分f(t)的近似向量f(k )； （4）构造f(k )对应的时间向量k。 下面是利用MATLAB实现连续信号卷积运算的通用函数sconv()，它在计算出卷积积分近 似值的同时，还绘出f(t)的波形图。 function [f,k]=sconv(f1,f2,k1,k2,p) %计算连续信号卷积积分f(t)=f1(t)*f2(t) %f:卷积积分f(t)对应的非零值向量 %k:f(t)的对应时间向量 %f1,f2:f1(t),f2(t)的非零样值向量 %k1,k2:f1(t),f2(t)的对应时间向量 %p:取样时间间隔 f=conv(f1,f2); %计算序列f1,f2的卷积和f 例1：已知编制一个m文件，绘出f1(t)、f2(t)和f(t)=f1(t)* f2(t)的波形。 解：实现所要求运算的m文件如下，运算结果如图所示。 clear;clc; p=0.01; k1=0:p:2 f1=exp(-k1); k2=0:p:3; f2=ones(1,length(k2)); [f,k]=sconv(f1,f2,k1,k2,p) 例题2：根据给定的两个连续时间信号x(t) = t[u(t)-u(t-1)]和h(t) = u(t)-u(t-1)，编写 程序，完成这两个信号的卷积运算，并绘制它们的波形图。 解：% Program1_6 % This program computes the convolution of two continuou-time signals clear;close all;%close all是关闭所有图形窗口 t0 = -2; t1 = 4; dt = 0.01; t = t0:dt:t1;%-2:0.01:4 x =stepfun(t,0)-stepfun(t,1); h = t.*(stepfun(t,0)-stepfun(t,1)); y = dt*conv(x,h); % Compute the convolution of x(t) and h(t) subplot(221) plot(t,x), grid on, title(Signal x(t)), axis([t0,t1,-0.2,1.2]) subplot(222) plot(t,h), grid on, title(Signal h(t)), axis([t0,t1,-0.2,1.2]) subplot(212) t = 2*t0:dt:2