北理工-信号与系统实验报告.docxVIP

实验1信号的时域描述与运算（基础实验）

实验目的

掌握信号的MATLAB表示及可视化方法。

掌握信号基本时域运算的MATLAB实现方法。

利用MATLAB分析常用信号，加深对信号时域特性的理解。

实验原理与方法

连续时间信号的MATLAB表示

连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。

从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用信号等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法成为向量表示法。表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。

连续时间信号的时域运算

对连续时间信号的运算包括两个信号的相加、相乘、微分、积分、以及移位、反转、尺度变换（尺度伸缩）等。

相加和相乘

信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算，此时要求表示两信号的向量时间范围和采样间隔相同。采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。

微分和积分

对于向量表示法表示的连续时间信号，可以通过数值计算的方法计算信号的微分和积分。这里微分是用差分来近似求取的，由时间向量[t1,t2,…,tn]和采样值向量[x1,x2,…,xn]表示的连续时间信号，其微分可以通过下式来求得

x

其中?t表示采样间隔。MATLAB中用diff函数来计算差分x

连续时间信号的定积分可以由MATLAB的quad函数实现，调用格式为

quad

其中，function_name为被积函数名，a,b为被积区间。

对于符号对象表示的连续时间信号，MATLAB提供了diff函数和quad函数分别用于求微分和积分。

移位、反转和尺度变换

信号的移位：信号x（t）的自变量为（t-t0），表示x（t）波形在t轴上整体移动，当t0

信号的反转：信号x（t）的自变量t更换为-t，x（t）的波形相当于的波形以t=0为轴反转过来。

信号的尺度变换：信号x（t）的自变量t更换为at，x（at）表示信号压缩a1或拉伸a1。

离散时间信号的MATLAB表示

离散时间信号仅在一些离散时刻有定义。在MATLAB中离散时间信号需要使用两个向量来表示，其中一个向量用于表示离散的时间点，另一个向量表示在这些时间点上的值。

离散时间信号的时域运算

离散时间信号的相加相乘是将两个信号对应的实践点上的值相加或相乘，可以直接使用算术运算的运算符“+”或“*”来计算。

离散时间信号的位移，则可看作是将表时间的向量平移，而表示对应时间点上的值的向量不变。

离散时间信号的反转，则可看做是将表示时间的向量和表示对应时间点上的值的向量以零点为基准点，以纵轴为对称轴反折，向量的反折可以利用MATLAB的fliplr函数实现。

实验内容

利用MATLAB绘制下列连续时间信号波形：

x

实验过程与实验结果：

在编辑器中编写程序代码如下：

symst;

x=(1-exp(-0.5*t))*heaviside(t);

ezplot(x,[-0.5,10]);

得到的时域图形如下：

x

实验过程与实验结果：

在编辑器中编写程序代码如下：

symst;

x=cos(pi*t)*(heaviside(t)-heaviside(t-2));

ezplot(t,x);

得到的时域图形如下：

x

实验过程与实验结果：

在编辑器中编写程序代码如下：

symst;

x=abs(t)/2*cos(pi*t)*(heaviside(t+2)-heaviside(t-2));

ezplot(x,[-3,3]);

得到的时域图形如下：

x

实验过程与实验结果：

在编辑器中编写程序代码如下：

symst;

x=exp(-t)*sin(2*pi*t)*(heaviside(t)-heaviside(t-3));

ezplot(t,x);

得到时域图形如下：

利用MATLAB绘制下列离散时间信号波形：

x

实验过程与实验结果：

在编辑器中编写程序代码如下：

n=0:10;

x=heaviside(n-3);

stem(n,x,filled);

得到的时域图形如下：

x

实验过程与实验结果：

在编辑器中编写程序代码如下：

n=0:10;

x=(-1/2).^n.*heaviside(n);

stem(n,x,filled);

得到时域图形如下：

x

实验过程与实验结果：

在编辑器中编写程序代码如下：

n=0:10;

x=n.*(heaviside(n)-hea