北京理工大学实验1信号的时域描述与及运算

《北京理工大学实验1信号的时域描述与及运算》由会员分享，可在线阅读，更多相关《北京理工大学实验1信号的时域描述与及运算（133页珍藏版）》请在金锄头文库上搜索。

1、实验1信号的时域描述与及运算一、实验目的1.掌握信号的MATLAB表示及其可视化方法。2.掌握信号基本时域运算的MATLAB实现方法。3.利 用MATLAB分析常用信号，加深对信号时域特性的理解。二、实验原理与方法1.连续时间信号的MATLAB表示连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出

2、连续时间信号，这种表示方法称为向量表示法。表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。例如一个正弦信号可以表示如下：t=0:0.01:10;x=sin(t);利 用plot(t,x)命令可以绘制上述信号的时域波形，如 图1所示。如果连续时间信号可以用表达式来描述，则还可以采用符号表达式来表示信号。例如对于上述正弦信号，可以用符号对象表示如下：x=sin(t);ezplot(X);利用ezplot(x)命令可以绘制上述信号的时域波形10.80.60.40.20Time(seconds)图1利用向量表示连续时间信号-1

3、0.50-0.5sin(t)24图 2 利 用 符 号 对 象 表 示 连 续 时 间 信 号0t6常用的信号产生函数函数名功能函 数名功能heaviside单 位 阶 跃 函数rectpuls门函数sin正弦函数tripuls三 角 脉 冲 函数cos余弦函数square周期方波sinesine函数sawtooth周 期 锯 齿 波或三角波exp指数函数2.连续时间信号的时域运算对连续时间信号的运算包括两信号相加、相乘、微分、积 分，以及位移、反转、尺度变换（尺度伸缩）等。1 ）相加和相乘信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和

4、“*”来 计 算，此时要求表示两信号的向量时间范围和采样间隔相同。采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。2）微分和积分对于向量表示法表示的连续时间信号，可以通过数值计算的方法计算信号的微分和积分。这里微分使用差分来近似求 取 的，由时间向量 乙,，2,，八 和采样值向量 尤 1，%2,表示的连续时间信号，其微分可以通过下式求得其中，表示采样间隔。MATLAB中 用 d iff函数来计算差分4+1 f。连续时间信号的定积分可以由MATLAB的 qud函数实现，调用格式为quad（function_name,a,b）其 中，function_name为被积函数名，a、b

5、为积分区间。对于符号对象表示的连续时间信号，MATLAB提供了diff函数和quad函数分别用于求微分和积分。3.离散时间信号的MATLAB表示离散时间信号仅在一些离散时刻有定义。在MATLAB中离散时间信号需要使用两个向量来表示，其中一个向量用于表示离散的时间点，另一个向量表示在这些时间点上的值。例如对于如下时间信号x()=-3,2,-1,2,1,-1,2,3)采 用MATLAB可以表示如下：n=-3:4;x=-3 2-1 2 1 -1 2 3;stem(n,x,filled);xlabel(n);title(x(n);Stem函数用于绘制离散时间信号波形，为了与我们表示离散时间信号的习惯相

6、同，在绘图时一般需要添加filled选 项，以绘制实心的杆状图形。上述命令绘制的信号时域波形如图3所示。x(n)图3离散时间信号示例4.离散时间信号的时域运算离散时间信号的相加相乘是将两个信号对应的时间点上的值相加或相乘，可以直接使用算术运算的运算符“+”和“*”来计算。离散时间信号的位移，则可看作是将表示时间的向量平移，而表示对应时间点上的值的向量不变。离散时间信号的反转，则可以看作是将表示时间的向量和表示对应时间点上的值的向量以零点为基准点，一纵轴为对 称 轴 反 折，向量的反折可以利用MATLAB的fliplr函数实现。三、实验内容(1)利 用 MATLAB绘制下列连续时间信号波形。加)

7、=(1-/)MATLAB程序如下：clear,close all,dt=0.01;t=-2:dt:2;x=(1-exp(-0.5*t).*heaviside(t);plot(t,x)title(x(t)xlabel(t)波形图如下：0.7x(t)0.6/0.5/0.4-0.3 r-0.2 0.1-Q 匕 一 r r-r 一:-2-1.5-1-0.5 0 0.5 1 1.5 2tQ x(f)=cos(m)w(r)-(r-2)MATLAB程序如下：clear,close all,dt=0.01;t=-4:dt:4;x=(cos(pi*t).*(heaviside(t)-heaviside(t-2)

8、;plot(t,x)title(x(t)xlabel(t)波形图如下：.M f)省 cos(砌必+2)-5 2)MATLAB程序如下：clear,close all,dt=0.01;t=-2:dt:2;x=(abs(t)/2).*(cos(pi*t).*(heaviside(t+2)-heaviside(t-2);plot(t,x)title(x(t)xlabel(t)波 形 图 如 下：x(r)=e sin(2m)(f)-u(t-3)MATLAB程 序 如 下：clear,close all,dt=0.01;t=-6:dt:6;x=(exp(-1*t).*(sin(2*pi*t).*(hea

9、viside(t)-heaviside(t-3);plot(t,x)title(x(t)xlabelft)波形图如下:(2)利 用MATLAB绘制下列离散时间信号波形 x(n)=w(n-3)MATLAB程序如下:clear,close all,n=-10.9999:10.9999;x=heaviside(n-3);stem(n,x)title(xn)xlabel(n)波形图如下：。尤()=(一 1 /2)(九)MATLAB程序如下:clear,close all,n=-10.9999:10.9999;x=(-0.5).An).*heaviside(n);stem(n,x)title fxn)x

10、label(n)波形图如下：x()=nu(n)-u(n-5)MATLAB程序如下:clear,close all,n=-10.9999:10.9999;x=n.*(heaviside(n)-heaviside(n-5);stem(n,x)title(xn)xlabel(n)波 形 图 如 下：4.543.532.521.510.50 xn x(n)=sin(/2)w(n)MATLAB程序如下:clear,close all,n=-10.9999:10.9999;x=(sin(n*pi/2).*heaviside(n);stem(n,x)title(xln1)xlabel(n1)波形图如下：10

11、.80.60.40.20-0.2-0.4-0.6-0.8-1-15-10-5 0 5 10 15n(3 )利 用MATLAB生成并绘制连续周期矩形波信号，要求周期为2,峰值为3,显示三个周期的波形。利 用MATLAB提供的square函数来生成方波信号xnMATLAB程序如下:clear,close all,t=0:0.01:6;x=square(pi*t).*3;plot(t,x);title(x(t)(4)已知信号否，及信号工2(，)=sin(2 R)，用MATLAB绘出下列信号的波形。由教材上石的波形可知不为一个三角波的右半部分。x3(t)=xl(t)+x2(t)MATLAB程序如下：c

12、lear,close all,t=-10:0.01:10;x1=4.*tripuls(t,8).*heaviside(t);x2=sin(2*pi*t);x3=x1+x2;plot(t,x3);xlabel(t);title(x(t);波形图如下：x(t)52x4(/)=x,(r)xx2(r)MATLAB程序如下：clear,close all,t=-10:0.01:10;x1=4.*tripuls(t,8).*heaviside(t);x2=sin(2*pi*t);x3=x1.*(x2);plot(t,x3);xlabel(t);title(x(t)；波 形 图 如 下：3 天=M(T)+X

13、|)MATLAB程 序 如 下：clear,close all,t=-10:0.01:10;x1=4.*tripuls(t,8).*heaviside(t);x2=4.*tripuls(t,8).*heaviside(-t);x3=x1+x2;plot(t,x3);xlabel(t);title(x(t);波形图如下:Xb(t)=X2(t)XX3(t-1)MATLAB程序如下：clear,close all,t=-10:0.01:10;x1=4.*tripuls(t-1,8).*heaviside(t-1);x2=sin(2*pi*t);x3=x1+sin(2*pi.*(t-1);x6=x2.

14、*(x3);plot(t,x6);xlabel(t);title(x(t);波形图如下：(5)已 知 离 散 时 间 信 号 ()，用M ATLAB绘出尤(),x(+2)和一 2)的 波 开 乡。由教材上x()的波形可知、()=12?，3,3,3 x()MATLAB程序如下：clear,close all,n=-10.9999:10.9999;x=heaviside(n+2)+heaviside(n+1)+heaviside(n);stem(n,x);xlabel(n);title(xn);波形图如下：。x(-n)MATLAB程序如下：clear,close all,n=-10.9999:10

15、.9999;x=heaviside(-n+2)+heaviside(-n+1)+heaviside(-n);stem(n,x);xlabel(n);title(xn);波形图如下：3x(+2)MATLAB程序如下：clear,close all,n=-10.9999:10.9999;x=heaviside(n+4)+heaviside(n+3)+heaviside(n+2);stem(n,x);xlabel(n);title(xn);波形图如下：x(n-2)MATLAB程序如下:clear,close all,n=-10.9999:10,9999;x=heaviside(n)+heavisid

16、e(n-1)+heaviside(n-2);stem(n,x);xlabel(n);title(xn);波形图如下：(6)用MATLAB编程绘制下列信号的时域波形，观察信号是否为周期信号？若是周期信号，周期是多少？若不是周期信 号，请说明原因。MATLAB程序如下：clear,close all,t=-10:0.01:10;x=1+cos(0.25*pi*t-pi/3)+2.*cos(0.5*pi*t-0.25*pi)+cos(2*pi*t)；plot(t,x);xlabel(t);title(x(t);波形图如下：x(t)5-10-8-6-4-2 0 2 4 6 8 10t该信号是周期信号，周期为8。x(f)=sin(f)+2sin(m)MATLAB程序如下:clear,close all,t=-10:0.01:10;x=sin(t)+2.*sin(pi*t);plot(t,x);xlabel(t);title(x(t);波形图如下:不是周期信号，此函数由一个周期为2TT和一个周期为2的函数组成，2rr和2没有最小公倍数，所以没有周期。0)x()=2 +3 sin(等 -令MATLA