1、北京理工大学信号与系统实验报告 本科实验报告 实验名称： 信号与系统实验 课程名称： 数字信号处理 实验时间： 任课教师： 实验地点： 4-442 实验教师： 实验类型： 原理验证 综合设计 自主创新 学生姓名： 学号/班级： 组 号： 学 院： 信息与电子学院 同组伙伴： 专 业： 信息工程 成 绩： 实验一 信号的时域描述与运算 一、实验目的 掌握信号的MATLAB表示及其可视化方法。 掌握信号根本时域运算的MATLAB实现方法。 利用MATLAB分析常用信号，加深对信号时域特性的理解。 二、实验原理与方法 1. 连续时间信号的MATLAB表示 连续时间信号指的是在连续时间范围内有定义的信
2、号，即除了假设干个不连续点外，在任何时刻信号都有定义。在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。 从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。例如一个正弦信号可以表示如下： t=0:0.01:10; x=sin(t); 利用plot(t,x)命令可以绘制上述信号的
3、时域波形，如图1所示。 如果连续时间信号可以用表达式来描述，那么还可以采用符号表达式來表示信号。例如对于上述正弦信号，可以用符号对象表示如下： x=sin(t); ezplot(X); 利用ezplot(x)命令可以绘制上述信号的时域波形 常用的信号产生函数 函数名 功能 函数名 功能 heaviside 单位阶跃函数 rectpuls 门函数 sin 正弦函数 tripuls 三角脉冲函数 cos 余弦函数 square 周期方波 sinc sinc函数 sawtooth 周期锯齿波或三角波 exp 指数函数 2.连续时间信号的时域运算 对连续时间信号的运算包括两信号相加、相乘、微分、积分，
4、以及位移、反转、尺度变换尺度伸缩等。 1相加和相乘 信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+和“x来计算，此时要求表示两信号的向量时间范围和采样间隔相同。采用符号对象表示的两个信号，可以直接根据符号对象的运算规那么运算。 2微分和积分 对于向量表示法表示的连续时间信号，可以通过数值计算的方法计算信号的微分和积分。这里微分使用差分来近似求取的，由时间向量和采样值向量表示的连续时间信号，其微分可以通过下式求得 其中表示采样间隔。MATLAB中用diff函数来计算差分。 连续时间信号的定积分可以由MATLAB的qud函数实现，调用格式为 q
5、uad (function_name,a,b) 其中，function_name为被积函数名，a、b为积分区间。 对于符号对象表示的连续时间信号，MATLAB提供了diff函数和quad函数分别用于求微分和积分。 3.离散时间信号的MATLAB表示 离散时间信号仅在一些离散时刻有定义。在MATLAB中离散时间信号需要使用两个向量来表示，其中一个向量用于表示离散的时间点，另一个向量表示在这些时间点上的值。例如对于如下时间信号 采用MATLAB可以表示如下： n=-3:4; x=-3 2 -1 2 1 -1 2 3; stem(n,x,filled); xlabel(n); title(x(n);
6、 Stem函数用于绘制离散时间信号波形，为了与我们表示离散时间信号的习惯相同，在绘图时一般需要添加filled选项，以绘制实心的杆状图形。上述命令绘制的信号时域波形如图3所示。 4.离散时间信号的时域运算 离散时间信号的相加相乘是将两个信号对应的时间点上的值相加或相乘，可以直接使用算术运算的运算符“+和“x来计算。 离散时间信号的位移，那么可看作是将表示时间的向量平移，而表示对应时间点上的值的向量不变。 离散时间信号的反转，那么可以看作是将表示时间的向量和表示对应时间点上的值的向量以零点为基准点，一纵轴为对称轴反折，向量的反折可以利用MATLAB的fliplr函数实现。 三、实验内容 1、利用
7、MATLAB绘制以下连续时间信号的波形： 1、 t=0:0.001:10; x=(1-exp(-0.5xt).x(heaviside(t); plot(t,x) 0 1 2 3 4 5 6 7 8 9 10 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 2、 t=0:0.001:10; x=cos(pixt).x(heaviside(t)-heaviside(t-2); plot(t,x) 0 1 2 3 4 5 6 7 8 9 10 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 3、 t=0:0.001:10; x=abs
8、(t)/2.xcos(pixt).x(heaviside(t+2)-heaviside(t-2); plot(t,x) 0 1 2 3 4 5 6 7 8 9 10 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 4、 t=0:0.001:10; x=exp(-t).xsin(2xpixt).x(heaviside(t)-heaviside(t-3); plot(t,x) 0 1 2 3 4 5 6 7 8 9 10 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 2、利用MATLAB绘制以下离散时间信号的波形： 1、 x=heaviside(n-3);
9、 stem(n,x,filled) 0 1 2 3 4 5 6 7 8 9 10 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 2、 x=(-1/2).xheaviside(n); stem(n,x,filled) 3、 x=n.x(heaviside(n)-heaviside(t-5); stem(n,x,filled) 0 1 2 3 4 5 6 7 8 9 10 0 0.5 1 1.5 2 2.5 3 3.5 4 4、 x=sin(nxpi/2).xheaviside(n); stem(n,x,filled) 0 1 2 3 4 5 6 7 8 9 10 -
10、1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 3、利用MATLAB生成并绘制连续周期矩形波信号，要求周期为2，峰值为3，显示三个周期波形。 t=0:0.001:6; x=square(txpi).x3; plot(t,x) 0 1 2 3 4 5 6 -3 -2 -1 0 1 2 3 (4)、信号,用MATLAB绘出以下信号的波形 1、 t=0:0.001:10; x1=(-t+4).x(heaviside(t)-heaviside(t-4); x2=sin(2xpixt); x3=x1+x2; plot(t,x3) 0 1 2 3 4 5 6 7 8 9 1
11、0 -1 0 1 2 3 4 5 2、 t=0:0.001:10; x1=(-t+4).x(heaviside(t)-heaviside(t-4); x2=sin(2xpixt); x4=x1.xx2; plot(t,x4) 0 1 2 3 4 5 6 7 8 9 -4 -3 -2 -1 0 1 2 3 3、 syms t; x1=(4-t)x(heaviside(t)-heaviside(t-4); x1=(4+t)x(heaviside(-t)-heaviside(t+4); x5=x1+x2; ezplot(t,x5); 5、离散时间信号x(n)波形，用MATLAB绘出x(n)、x(-n
12、)、x(n+2)和x(n-2)的波形。 n=-3:4; x=0,1,2,3,3,3,3,0; n1=-fliplr(n); x1=fliplr(x); n2=n+2; n3=n-2; subplot(221); stem(n,x,filled); subplot(222); stem(n1,x1,filled); subplot(223); stem(n2,x,filled); subplot(224); stem(n3,x,filled); 6、有MATLAB编程绘制以下信号的时域波形，观察信号是否为周期信号？假设是周期信号，周期是多少？假设不是周期信号，请说明原因。 1、 syms t;
13、x=1+cos(pi/4xt-pi/3)+2xcos(pi/2xt-pi/4)+cos(2xpixt); ezplot(x,-20,20); -10 -8 -6 -4 -2 0 2 4 6 8 10 -3 -2 -1 0 1 2 3 4 5 是周期信号 T=8s 2、 x=sin(t)+2xsin(pixt); plot(t,x) -20 -15 -10 -5 0 5 10 15 20 -3 -2 -1 0 1 2 3 不是周期信号 3、 n=1:20; x=2+3xsin(2xn/3xpi-pi/8); stem(n,x,filled); 是周期信号 T=3s 4、 n=1:20; x=cos(nxpi/6)+sin(nxpi/3)+cos(nxpi/2); stem(n,x,filled); 是周期信号 T=12s 实验二 LTI系统的时域分析 一、实验目的 1、掌握利用MATLAB对系统进行时域分析的方法。 2、掌握连续时间系统零状态响应、冲激响应和阶跃响应的求解方法。 3、掌握求解离散时间系统响应、单位抽样响应的方法。 4、加深对卷积积分和卷积和的理解。掌握利用计算机进行卷积积分和卷积和计算的方法。 二、实验原理与方法 1、