信号与系统实验 曹玉树.docVIP

实验2 信号和系统的时域分析 一、实验目的 1、掌握利用MATLAB表示信号和对信号进行基本时域运算的方法。 2、掌握利用MATLAB对系统进行时域分析的方法。 3、掌握连续时间系统零状态响应、冲激响应和阶跃响应的求解方法。 4、掌握离散时间系统响应、单位抽样响应的方法。 5、加深对卷积积分和卷积和的理解。掌握利用计算机进行卷积积分和卷积和计算的方法。 二、实验原理 1、连续时间信号和离散时间信号的MATLAB表示与运算 1）连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用信号等时间间隔采样后的采样值近似表示的。当采样间隔足够小时，这些采样值就可以很好的近似表示出连续时间信号。 例如：t=0:0.01:10; x=sin(t); 利用plot命令即可绘制该正弦信号。 同时，还可以用符号表达式来表示信号。 例如对上面所示正弦信号：syms t; x=sin(t); 利用ezplot命令可绘制上述信号。 离散时间信号仅在一些离散时刻有定义。在MATLAB中离散时间信号用 两个向量加以表示，一个表示离散的时间点，另一个表示这些时间点上的值，例如：对于信号x[n]={-3,2,-1,2,1,-1,2,3} 可以用MATLAB表示为：n=-3:4; x=[-3 2 -1 1 -1 2 3]; 利用stem绘制波形，一般需添加“filled”选项，以绘制实心的杆状图形。 2、连续时间系统时域分析的MATLAB表示 1）连续时间系统的MATLAB表示 LTI连续系统通常有微分方程表示出，在MATLAB里可以建立模型： b=[bn,bn-1……bo]; a=[an,an-1……ao]; sys=tf(b,a); b表示x(t)的系数，a表示y（t）的系数‘ 2）连续时间系统的零状态响应 零状态响应指系统的初始状态为0，仅有输出信号所引起的响应，在MATLAB中用lsim函数表示。 例如：lsim(sys,x,t) 则绘出输入信号及响应的波形，x表示输入信号。 Y=lsim(sys,x,t) 不绘出波形，而是返回响应的数值向量‘ 3）连续时间系统的冲激响应与阶跃响应 MATLAB中利用impulse函数来求指定时间范围内，由sys描述的连续时间系统的单位冲激响应。 impulse(sys) 在默认时间范围内绘出系统冲激响应的时域波形 impulse(sys,T) 在0—T范围内绘出系统冲激响应的时域波形 impulse(systs:tp:te) 在ts—te范围内，以tp为时间间隔取样的系统冲激响应波形 函数step用于求解单位阶跃响应，调用方法同上。 3、离散时间系统时域分析的MATLAB实现 1）离散时间系统的MATLAB表示 LTI离散时间系统通常可由系统差分方程描述，在MATLAB里： b=[b0,b1……bn]; a=[a0,a-1……am]; 2）离散时间系统对任意输入的响应 MATLAB中用filter函数 y=filter(b,a,x)，x为输入序列，y为输出序列。 3）离散时间系统的单位抽样响应 MATLAB中用impz函数来求指定时间范围内，由向量b和a描述的离散时间系统的单位响应。方法与impulse函数相同。 4、卷积和与卷积积分 1）离散时间序列的卷积和 MATLAB中的conv函数可以用来计算卷积和，格式为x=conv(x1,x2). 例如：x1=ones(1,3); x2=[1 2 3 4]; x=conv(x1,x2); X= 3 6 9 7 4 2 注意卷积后时域的宽度变化。 2）连续时间信号的卷积积分 将连续时间信号分成间隔足够小的离散时间信号即可近似计算积分。具体步骤如下： 将连续信号以时间间隔进行取样，得到离散序列。 构造离散序列。 调用函数conv计算卷积积分在时的近似采样值。 构造离散序列对应的时间向量n。 三、实验内容 1、syms t; x=1+cos(pi*t/4-pi/3)+2*cos(pi*t/2-pi/4)+cos(2*pi*t); ezplot(x) 信号是周期的，周期为8。 (2) syms t; x=sin(t)+2*sin(pi*t); ezplot(x) 信号不是周期的，因为两个信号没有公周期。 n=0:20; x=2+3*sin(2*n*pi/3-pi/8); stem(n,x,filled) 信号是周期的，周期为3。 （4） n=0:20; x=cos(n*pi/6)+