《北京理工大学 信号及系统实验报告》由会员分享，可在线阅读，更多相关《北京理工大学 信号及系统实验报告（83页珍藏版）》请在金锄头文库上搜索。
1、 本科实验报告实验名称： 信号与系统实验 课程名称：数字信号处理实验时间：任课教师：范哲意实验地点：4-442实验教师：何冰松, 范哲意实验类型： 原理验证 综合设计 自主创新学生姓名：学号/班级：组 号：学 院：信息与电子学院同组搭档：专 业：信息工程成 绩：实验一 信号的时域描述与运算一、实验目的 掌握信号的MATLAB表示及其可视化方法。掌握信号基本时域运算的MATLAB实现方法。利用MATLAB分析常用信号，加深对信号时域特性的理解。二、实验原理与方法1. 连续时间信号的MATLAB表示连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。在M
2、ATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。例如一个正弦信号可以表示如下： t=0:0.01:10; x=sin(t);利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。如果连续时间信号可以用表达式来描述，则还
3、可以采用符号表达式來表示信号。例如对于上述正弦信号，可以用符号对象表示如下： x=sin(t); ezplot(X);利用ezplot(x)命令可以绘制上述信号的时域波形常用的信号产生函数 函数名 功能 函数名 功能heaviside单位阶跃函数rectpuls门函数sin正弦函数tripuls三角脉冲函数cos余弦函数square周期方波sincsinc函数sawtooth周期锯齿波或三角波exp指数函数2.连续时间信号的时域运算 对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换（尺度伸缩）等。1）相加和相乘信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个
4、采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算，此时要求表示两信号的向量时间范围和采样间隔相同。采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。2）微分和积分对于向量表示法表示的连续时间信号，可以通过数值计算的方法计算信号的微分和积分。这里微分使用差分来近似求取的，由时间向量和采样值向量表示的连续时间信号，其微分可以通过下式求得其中表示采样间隔。MATLAB中用diff函数来计算差分。连续时间信号的定积分可以由MATLAB的qud函数实现，调用格式为quad (function_name,a,b)其中，function_name为被积函数名，a、b为积分区间。对
5、于符号对象表示的连续时间信号，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);Stem函数用于绘制离散时间信号波形，为了与我们表示离散时间信号的习惯相同，在绘图时一般需要添加filled选项，以绘制
6、实心的杆状图形。上述命令绘制的信号时域波形如图3所示。4.离散时间信号的时域运算 离散时间信号的相加相乘是将两个信号对应的时间点上的值相加或相乘，可以直接使用算术运算的运算符“+”和“*”来计算。离散时间信号的位移，则可看作是将表示时间的向量平移，而表示对应时间点上的值的向量不变。离散时间信号的反转，则可以看作是将表示时间的向量和表示对应时间点上的值的向量以零点为基准点，一纵轴为对称轴反折，向量的反折可以利用MATLAB的fliplr函数实现。三、实验内容（1）、利用MATLAB绘制下列连续时间信号的波形：1、 t=0:0.001:10; x=(1-exp(-0.5*t).*(heavisid
7、e(t); plot(t,x)01234567891000.10.20.30.40.50.60.70.80.912、 t=0:0.001:10; x=cos(pi*t).*(heaviside(t)-heaviside(t-2); plot(t,x)012345678910-1-0.8-0.6-0.4-0.200.20.40.60.83、 t=0:0.001:10; x=abs(t)/2.*cos(pi*t).*(heaviside(t+2)-heaviside(t-2); plot(t,x)012345678910-0.6-0.4-0.200.20.40.60.814、 t=0:0.001:
8、10; x=exp(-t).*sin(2*pi*t).*(heaviside(t)-heaviside(t-3); plot(t,x)012345678910-0.6-0.4-0.200.20.40.60.8（2）、利用MATLAB绘制下列离散时间信号的波形：1、 x=heaviside(n-3);stem(n,x,filled) 01234567891000.10.20.30.40.50.60.70.80.92、 x=(-1/2).n.*heaviside(n); stem(n,x,filled)3、 x=n.*(heaviside(n)-heaviside(t-5); stem(n,x,f
9、illed)01234567891000.511.522.533.544、 x=sin(n*pi/2).*heaviside(n); stem(n,x,filled)012345678910-1-0.8-0.6-0.4-0.200.20.40.60.8（3）、利用MATLAB生成并绘制连续周期矩形波信号，要求周期为2，峰值为3，显示三个周期波形。 t=0:0.001:6; x=square(t*pi).*3; plot(t,x)0123456-3-2-10123(4)、已知信号,用MATLAB绘出下列信号的波形1、 t=0:0.001:10; x1=(-t+4).*(heaviside(t)-
10、heaviside(t-4); x2=sin(2*pi*t);x3=x1+x2;plot(t,x3)012345678910-10123452、 t=0:0.001:10; x1=(-t+4).*(heaviside(t)-heaviside(t-4); x2=sin(2*pi*t); x4=x1.*x2;plot(t,x4)0123456789-4-3-2-101233、syms t;x1=(4-t)*(heaviside(t)-heaviside(t-4);x1=(4+t)*(heaviside(-t)-heaviside(t+4);x5=x1+x2;ezplot(t,x5);5）、已知离
11、散时间信号x(n)波形（），用MATLAB绘出x(n)、x(-n)、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编程绘制下列信号的时域波形，观察信号是否为周期信号？若是周期信号，周期是多少？若不是周期信号，
12、请说明原因。1、syms t;x=1+cos(pi/4*t-pi/3)+2*cos(pi/2*t-pi/4)+cos(2*pi*t);ezplot(x,-20,20);-10-8-6-4-20246810-3-2-1012345是周期信号 T=8s2、 x=sin(t)+2*sin(pi*t); plot(t,x)-20-15-10-505101520-3-2-10123不是周期信号3、n=1:20;x=2+3*sin(2*n/3*pi-pi/8);stem(n,x,filled);是周期信号 T=3s4、n=1:20;x=cos(n*pi/6)+sin(n*pi/3)+cos(n*pi/2);stem(n,x,filled); 是周期信号 T=12s实验二 LTI系统的时域分析一、实验目的1、掌握利用MATLAB对系统进行时域分析的方法。2、掌握连续时间系统零状态响应、冲激响应和阶跃响应的求解方法。3、掌握求解离散时间系统响应、单位抽样响应的方法。4、加深对卷积积分和卷积和的理解。掌握利用计算机进行卷积积分和卷积和计算的方法。 二、实验原理与方法1、连续时间系统时域分析的MATLAB实现（1)、连续
