实验一连续时间信号的时域基本运算一、实验目的（1）掌握连续时间信号时域运算的基本方法；（2）掌握相关函数的调用格式及作用；（3）掌握连续信号的基本运算。

二、实验原理信号的基本运算包括信号的相加（减）和相乘(除).信号的时域变换包括信号的平移、翻转、倒相以及尺度变换。

（1）加减： f(t)=f1(t)±f2(t)（2）乘: f(t)=f1(t)×f2(t)（3）延时或平移：f(t)→(t-t0) t0>0时右移；t0<0时左移（4）翻转： f(t)→f(-t)（5）尺度变换：f(t)→ f(at) |a|>1时尺度缩小；|a|<1时尺度放大；a<0时，尺度翻转。

（6）标量乘法：f(t)→af(t)（7）倒相： f(t)→-f(t)（8）微分： f(t)→df(t)/dt（9）积分： f(t)→∫t -∞f(t)d(t)三、涉及的MATLAB函数及其运算1、stepfun函数功能：产生一个阶跃信号。

调用格式：Stepfun(t,t0)其中，t是时间区间，在该区间内阶跃信号一定会产生；t是信号发生从0到1跳跃的时刻。

2、diff函数调用格式：diff (f) : 求函数f对预设独立变数的一次微分值。

diff (f, ’t’) : 求函数f对独立变数t的一次微分值。

3、int函数调用格式:Int（f）: 函数f对预设独立变数的积分值。

Int(f,’t’): 函数f对独立变数t的积分值。

4、heaviside函数Heaviside（t）:产生没有移位的阶跃信号。

Heaviside（t-k）:产生向右平移K单位的阶跃信号。

四、实验内容与方法1、验证性实验（1）移位实现连续信号的移位，即f(t-t0),或者f(t+t0);常数t0>0。

MATLAB程序clear allt=0:0.0001:2y=sin(2*pi*(t));y1=sin(2*pi*(t-0.2));plot(t,y,'-',t,y1,'--');ylabel('f(t)');xlabel('t');title('信号的移位')运行结果：（2）翻转（信号的翻转就是将信号的波形以纵轴为对称轴翻转180度，将信号f(t)的自变量t替换为-t即可得到其翻转信号。

）MATLAB程序：clear all;t=0:0.02:1;t1=-1:0.02:0;g1=3*t;g2=3*(-t1);plot(t,g1,'--',t1,g2);grid on;xlabel('t');ylabel('g(t)');title('信号的反折');运行结果;（3）尺度变换将信号f(t)中的自变量t替换为at。

MATLAB程序：clear all;t=0:0.001:1;a=2;y=sin(2*pi*t);y1=sin(2*a*pi);subplot(211)plot(t,y);ylabel('y(t)');xlabel('t');title('尺度变换');subplot(212)plot(t,y1);ylabel('y1(t)');xlabel('t');运行结果：（4）综合MATLAB程序syms tf=sym('(t/2+1)*(heaviside(t+2)-heaviside(t-2))')subplot(2,3,1);ezplot(f,[-3,3]);title('f');y1=subs(f,t,t+2);subplot(2,3,2);ezplot(y1,[-5,1]);title('y1');y2=subs(f,t,t-2);subplot(2,3,3);ezplot(y2,[-1,5]);title('y2');y3=subs(f,t,-t);subplot(2,3,4);ezplot(y1,[-3,3]);title('y3');y4=subs(f,t,2*t);subplot(2,3,5);ezplot(y1,[-2;2]);title('y4');y5=-f;subplot(2,3,6);ezplot(y1,[-3,3]);title('y5');表达式f =(t/2+1)*(heaviside(t+2)-heaviside(t-2))运行结果：2、程序设计实验syms tf=sym('(-2*t/3+2/3)*(heaviside(t+2)-heaviside(t-1))')subplot(2,4,1);ezplot(f,[-3,3]);title('f');y1=subs(f,t,-t);subplot(2,4,2);ezplot(y1,[-3,3]);title('y1');y2=subs(f,t,-t+2);subplot(2,4,3);ezplot(y2,[-1,5]);title('y2');y3=subs(f,t,-t-2);subplot(2,4,4);ezplot(y1,[-2,4]);title('y3');y4=subs(f,t,2*t);subplot(2,4,5);ezplot(y1,[-2,3]);title('y4');y5=subs(f,t,1/2*t);subplot(2,4,6);ezplot(y1,[-2,4]);title('y5');y6=subs(f,t,t-2);subplot(2,4,7);ezplot(y1,[-2,3]);title('y6');y7=subs(f,t,-1/2*t+1);subplot(2,4,8);ezplot(y1,[-4,8]);title('y7运行结果：五、实验总结-通过本次实验学会了连续时间信号时域运算的基本方法，掌握了stepfun函数，diff函数、int函数、heaviside函数调用格式以及作用，并且掌握了连续信号的基本运算。

实验二LTI系统的响应一、实验目的1. 熟悉连续时间系统的单位冲激响应、阶跃响应的意义及求解方法2. 熟悉连续（离散）时间系统在任意信号激励下响应的求解方法3. 熟悉应用MATLAB实现求解系统响应的方法二、实验原理1.连续时间系统对于连续的LTI系统，当系统输入为f(t)，输出为y(t)，则输入与输出之间满足如下的线性常系数微分方程：，当系统输入为单位冲激信号δ(t)时产生的零状态响应称为系统的单位冲激响应，用h(t)表示。

若输入为单位阶跃信号ε(t)时，系统产生的零状态响应则称为系统的单位阶跃响应，记为g(t)，如下图所示。

系统的单位冲激响应h(t)包含了系统的固有特性，它是由系统本身的结构及参数所决定的，与系统的输入无关。

我们只要知道了系统的冲激响应，即可求得系统在不同激励下产生的响应。

因此，求解系统的冲激响应h(t)对我们进行连续系统的分析具有非常重要的意义。

在MATLAB中有专门用于求解连续系统冲激响应和阶跃响应，并绘制其时域波形的函数impulse( ) 和step( )。

如果系统输入为f(t)，冲激响应为h(t)，系统的零状态响应为y(t)，则有：y(t)h(t)f(t)。

若已知系统的输入信号及初始状态，我们便可以用微分方程的经典时域求解方法，求出系统的响应。

但是对于高阶系统，手工计算这一问题的过程非常困难和繁琐。

在MATLAB中，应用lsim( )函数很容易就能对上述微分方程所描述的系统的响应进行仿真，求出系统在任意激励信号作用下的响应。

lsim( )函数不仅能够求出连续系统在指定的任意时间范围内系统响应的数值解，而且还能同时绘制出系统响应的时域波形图。

以上各函数的调用格式如下：⑴impulse( ) 函数函数impulse( )将绘制出由向量a和b所表示的连续系统在指定时间范围内的单位冲激响应h(t)的时域波形图，并能求出指定时间范围内冲激响应的数值解。

impulse(b,a) 以默认方式绘出由向量a和b所定义的连续系统的冲激响应的时域波形。

impulse(b,a ,t0) 绘出由向量a和b所定义的连续系统在0 ~ t0时间范围内冲激响应的时域波形。

impulse(b,a,t1:p:t2) 绘出由向量a和b所定义的连续系统在t1 ~ t2时间范围内，并且以时间间隔p均匀取样的冲激响应的时域波形。

y=impulse(b,a,t1:p:t2) 只求出由向量a和b所定义的连续系统在t1 ~ t2时间范围内，并且以时间间隔p均匀取样的冲激响应的数值解，但不绘出其相应波形。

⑵step( ) 函数函数step( )将绘制出由向量a和b所表示的连续系统的阶跃响应，在指定的时间范围内的波形图，并且求出数值解。

和impulse( )函数一样，step( )也有如下四种调用格式：step( b,a)step(b,a,t0)step(b,a,t1:p:t2)y=step(b,a,t1:p:t2)上述调用格式的功能和impulse( )函数完全相同，所不同只是所绘制（求解）的是系统的阶跃响应g(t)，而不是冲激响应h(t)。

⑶lsim( )函数根据系统有无初始状态，lsim( )函数有如下两种调用格式：①系统无初态时，调用lsim( )函数可求出系统的零状态响应，其格式如下：lsim(b,a,x,t)绘出由向量a和b所定义的连续系统在输入为x和t所定义的信号时，系统零状态响应的时域仿真波形，且时间范围与输入信号相同。

其中x和t是表示输入信号的行向量，t为表示输入信号时间范围的向量，x则是输入信号对应于向量t所定义的时间点上的取样值。

y=lsim(b,a,x,t) 与前面的impulse 和step函数类似，该调用格式并不绘制出系统的零状态响应曲线，而只是求出与向量t定义的时间范围相一致的系统零状态响应的数值解。

②系统有初始状态时，调用lsim( )函数可求出系统的全响应，格式如下：lsim(A,B,C,D,e,t,X0)绘出由系数矩阵A,B,C,D所定义的连续时间系统在输入为e和t所定义的信号时，系统输出函数的全响应的时域仿真波形。

t为表示输入信号时间范围的向量，e则是输入信号e(t)对应于向量t所定义的时间点上的取样值，X0表示系统状态变量X=[x1,x2,…..xn]'在t=0时刻的初值。