《北理工信号及系统实验报告》由会员分享，可在线阅读，更多相关《北理工信号及系统实验报告（90页珍藏版）》请在金锄头文库上搜索。
1、实验1 信号的时域描述与运算2实验2 LTI系统的时域分析18实验3 信号的频域分析32实验4 LTI系统的频域分析52实验5 连续时间系统的复频域分析64实验6 离散时间系统的z域分析79实验1 信号的时域描述与运算一、实验目的1.掌握信号的MATLAB表示及其可视化方法；2.掌握信号基本时域运算的MATLAB实现方法；3.利用MATLAB分析常用信号，加深对信号时域特性的理解；二、实验原理与方法1.连续时间信号的MATLAB表示连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号
2、对象表示法。从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。例如一个正弦信号可以表示如下： t=0:0.01:10; x=sin(t);利用plot(t,x)命令可以绘制上述信号的时域波形。如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。例如对于上述正弦信号，可以用符号对象表示如下：
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为积分区间。对于符号对象表示的连续时间信号，MATLAB提供了diff函数和quad函
5、数分别用于求微分和积分。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选项，以绘制实心的杆状图形。4.离散时间信号的时域运算离散时间信号的相加相乘是将两个信
6、号对应的时间点上的值相加或相乘，可以直接使用算术运算的运算符“+”和“*”来计算。离散时间信号的位移，则可看作是将表示时间的向量平移，而表示对应时间点上的值的向量不变。离散时间信号的反转，则可以看作是将表示时间的向量和表示对应时间点上的值的向量以零点为基准点，一纵轴为对称轴反折，向量的反折可以利用MATLAB的fliplr函数实现。三、实验内容1.利用MATLAB绘制下列连续时间信号波形：(1) 程序代码：%one_1_1.mt=0:0.01:10 ;x=1-exp(-0.5*t) ;plot(t,x);xlabel(t(s);title(x(t);grid on;波形图：(2) 程序代码：%
7、one_1_2.mt=0:0.01:2;x=cos(pi*t);plot(t,x);xlabel(t(s);title(x(t);grid on;波形图：(3) 程序代码：%one_1_3.mt=-2:0.01:2;x=0.5*abs(t).*cos(pi*t);plot(t,x);xlabel(t(s);title(x(t);grid on;波形图：(4) 程序代码：%one_1_4.mt=0:0.01:3;x=exp(-t).*sin(2*pi*t);plot(t,x);xlabel(t(s);title(x(t);grid on;波形图：2.利用MATLAB绘制下列离散时间信号波形：(1
8、) 程序代码：%one_2_1.mn=0:7;x=0 0 0 1 1 1 1 1;stem(n,x,filled);xlabel(t(s);title(x(t);波形图：(2) 程序代码：%one_2_2.mn=-2:5;x=0 0 1 -0.5 0.25 -0.125 0.0625 -0.03125;stem(n,x,filled);xlabel(t(s);title(x(t);波形图：(3) 程序代码：%one_2_3.mn=-2:7;x=0 0 0 1 2 3 4 0 0 0;stem(n,x,filled);xlabel(t(s);title(x(t);波形图：(4) 程序代码：%on
9、e_2_4.mn=0:7;x=sin(0.5*pi*n);stem(n,x,filled);xlabel(t(s);title(x(t);波形图：3.利用MATLAB生成并绘制连续周期矩形波信号，要求周期为2，峰值为3，显示三个周期的波形。程序代码：%one_3.m%t=0:0.01:10;t=linspace(0,6*pi);x=3*square(t);plot(t/pi,x);xlabel(t(s);title(x(t);grid on;波形图：4.已知信号，及信号，用MATLAB绘出下列信号的波形。(1) (2) (3) (4) 程序代码：%one_4.mt=-4:0.01:4;x1=(
10、4-t).*(heaviside(t)-heaviside(t-4);x2=sin(2*pi*t); subplot(221); % x3(t)=x1(t)+x2(t)x3=x1+x2;plot(t,x3);xlabel(t(s);title(x_3(t)=x_1(t)+x_2(t);grid on; subplot(222); % x4(t)=x1(t)*x2(t)x4=x1.*x2;plot(t,x4);xlabel(t(s);title(x_4(t)=x_1(t)*x_2(t);grid on; subplot(223); % x5(t)=x1(-t)+x1(t)t1=-fliplr(t
11、);x11=fliplr(x1); % x11(t)=x1(-t)x5=x11+x1;plot(t1,x5);xlabel(t(s);title(x_5(t)=x_1(-t)*x_1(t);grid on; subplot(224); % x6(t)=x2(t)+x3(t-1)t2=1-t;x31=fliplr(x3); % x31(t)=x3(t-1)x6=x31.*x2;plot(t2,x6);xlabel(t(s);title(x_6(t)=x_2(t)*x_3(t-1);grid on;波形图：5.已知离散时间信号，用MATLAB绘出的波形。程序代码：%one_5.mn=-3:3;x=
12、0 1 2 3 3 3 3;subplot(221); % x(t)stem(n,x,filled);xlabel(n);title(x(n);grid on;subplot(222); % x2(t)=x(-n)n2=-fliplr(n);x2=fliplr(x);stem(n2,x2,filled);xlabel(n);title(x(-n);grid on;subplot(223); % x3(t)=x(n+2)stem(n-2,x,filled);xlabel(n);title(x(n+2);grid on;subplot(224); % x4(t)=x(n-2)stem(n+2,x,
13、filled);xlabel(n);title(x(n-2);grid on;波形图：6.用MATLAB编程绘制下列信号的时域波形，观察信号是否为周期信号？若是周期信号，周期是多少？若不是周期信号，请说明原因。(1) 程序代码：%one_6_1.msyms t;x=1+cos(pi*t/4-pi/3)+2*cos(pi*t/2-pi/4)+cos(2*pi*t);ezplot(x,-10,10);xlabel(t(s);title(x(t);grid on;波形图：所以原信号是周期信号，周期为8。(2) 程序代码：%one_6_2.msyms t;x=sin(t)+2*sin(pi*t);ezplot(x,-20,20);xlabel(t(s);title(x(t);grid on;波形图：所以原信号不是周期信号。原函数由一个周期为2和一个周期为2的函数组成，2和2没有最小公倍数，所以没有周期。(3) 程序代码：%one_6_3.mn=-10:10;x=2+3*sin(2*n*pi/3-pi/8);stem(n,x);xlabel(n);title(x(n);grid on;波形图：所以原信号是周期函数，周期为3。(4) 程序代码：%on
