信号与系统matlab实验指导书.docx
- 文档编号:16359824
- 上传时间:2023-07-12
- 格式:DOCX
- 页数:57
- 大小:306.31KB
信号与系统matlab实验指导书.docx
《信号与系统matlab实验指导书.docx》由会员分享,可在线阅读,更多相关《信号与系统matlab实验指导书.docx(57页珍藏版)》请在冰点文库上搜索。
信号与系统matlab实验指导书
实验一基本信号的波形
一、实验目的:
1.掌握matlab软件的基本操作。
2.熟悉matlab的基本命令的使用。
3.掌握用matlab绘出基本信号。
二、实验原理:
1.信号的时域表示方法
1.1将信号表示成独立时间变量的函数
例如x(t)=sin(ωt)和x[n]=n(0.5)nu[n]
分别表示一个连续时间信号和一个离散时间信号。
在MATLAB中有许多内部函数,可以直接完成信号的这种表达,例如:
sin():
正弦信号
cos():
余弦信号
exp():
指数信号
1.2用信号的波形图来描述信号
用函数曲线表示一个信号,图1.1就是一个连续时间信号和一个离散时间信号的波形图。
图1.1连续时间信号与离散时间信号的波形图
1.3将信号用一个数据序列来表示
对于离散时间信号,还可以表示成一个数的序列,例如:
x[n]={....,0.1,1.1,-1.2,0,1.3,….}
↑n=0
在《信号与系统》和《数字信号处理》课程中,上述三种信号的描述方法是经常要使用的。
2用MATLAB仿真连续时间信号和离散时间信号
在MATLAB中,无论是连续时间信号还是离散时间信号,MATLAB都是用一个数字序列来表示信号,这个数字序列在MATLAB中叫做向量(vector)。
通常的情况下,需要与时间变量相对应。
如前所述,MATLAB有很多内部数学函数可以用来产生这样的数字序列,例如sin()、cos()、exp()等函数可以直接产生一个按照正弦、余弦或指数规律变化的数字序列。
2.1连续时间信号的仿真
程序Program1_1是用MATLAB对一个正弦信号进行仿真的程序,请仔细阅读该程序,并在计算机上运行,观察所得图形。
%Program1_1
%Thisprogramisusedtogenerateasinusoidalsignalanddrawitsplot
clear,%Clearallvariables
closeall,%Closeallfigurewindows
dt=0.01;%Specifythestepoftimevariable
t=-2:
dt:
2;%Specifytheintervaloftime
x=sin(2*pi*t);%Generatethesignal
plot(t,x)%Openafigurewindowanddrawtheplotofx(t)
title('Sinusoidalsignalx(t)')
xlabel('Timet(sec)')
常用的图形控制函数
axis([xmin,xmax,ymin,ymax]):
图型显示区域控制函数,其中xmin为横轴的显示起点,xmax为横轴的显示终点,ymin为纵轴的显示起点,ymax为纵轴的显示终点。
有时,为了使图形具有可读性,需要在所绘制的图形中,加上一些网格线来反映信号的幅度大小。
MATLAB中的gridon/gridoff可以实现在你的图形中加网格线。
gridon%在图形中加网格线。
gridoff%取消图形中的网格线。
holdon%图形显示窗口原来的图像保持。
holdoff%关闭图形保持功能。
figure%打开新的显示窗口。
x=input(‘Typeinsignalx(t)inclosedform:
’)%接收键盘输入。
在绘制信号的波形图时,有时我们需要将若干个图形绘制在同一个图形窗口中,这就需要使用MATLAB的图形分割函数subplot(),其用法是在绘图函数stem或plot之前,使用图形分割函数subplot(n1,n2,n3),其中的参数n1,n2和n3的含义是,该函数将把一个图形窗口分割成n1×n2个子图,即将绘制的图形将绘制在第n3个子图中。
在《信号与系统》课程中,单位阶跃信号u(t)和单位冲激信号δ(t)是二个非常有用的信号。
它们的定义如下
1.1(a)
1.1(b)
这里分别给出相应的简单的产生单位冲激信号和单位阶跃信号的扩展函数。
产生单位冲激信号的扩展函数为:
functiony=delta(t)
dt=0.01;
y=(u(t)-u(t-dt))/dt;
产生单位阶跃信号的扩展函数为:
functiony=u(t)
y=(t>=0);%y=1fort>0,elsey=0
请将这二个MATLAB函数分别以delta和u为文件名保存在work文件夹中,以后,就可以像教材中的方法使用单位冲激信号δ(t)和单位阶跃信号u(t)。
单位冲激信号:
x=-10:
0.1:
10;
y=(x==0);%就是个条件判断,只有x=0的时候,y才为“1”
plot(x,y);
如果是单位脉冲信号(离散的),改成stem(x,y);就行了
t=-5:
0.01:
5;
y=(t==0);
subplot(121);
plot(t,y,'r');
n=-5:
5;
x=(n==0);
subplot(122);
>>stem(n,x);
>>x=-2:
0.01:
2;
>>y=100*(x==0);
>>stairs(x,y)
练习产生正弦信号(sin、cos)、指数信号(exp)、符号函数(sign)、抽样函数(sinc)、矩形脉冲(rectpuls(t,width))、三角脉冲(tripuls(t,width,skew))
2.2离散时间信号的仿真
程序Program1_2用来产生离散时间信号x[n]=sin(0.2πn)。
%Program1_2
%Thisprogramisusedtogenerateadiscrete-timesinusoidalsignalanddrawitsplot
clear,%Clearallvariables
closeall,%Closeallfigurewindows
n=-10:
10;%Specifytheintervaloftime
x=sin(0.2*pi*n);%Generatethesignal
stem(n,x)%Openafigurewindowanddrawtheplotofx[n]
title('Sinusoidalsignalx[n]')
xlabel('Timeindexn')
请仔细阅读该程序,比较程序Program1_1和Program1_2中的不同之处,以便自己编程时能够正确使用这种方法方针连续时间信号和离散时间信号。
程序Program1_3用来仿真下面形式的离散时间信号:
x[n]={....,0.1,1.1,-1.2,0,1.3,….}
↑n=0
%Program1_3
%Thisprogramisusedtogenerateadiscrete-timesequence
%anddrawitsplot
clear,%Clearallvariables
closeall,%Closeallfigurewindows
n=-5:
5;%Specifytheintervaloftime,thenumberofpointsofnis11.
x=[0,0,0,0,0.1,1.1,-1.2,0,1.3,0,0];%Generatethesignal
stem(n,x,'.')%Openafigurewindowanddrawtheplotofx[n]
gridon,
title('Adiscrete-timesequencex[n]')
xlabel('Timeindexn')
由于在程序的stem(n,x,'.')语句中加有'.'选项,因此绘制的图形中每根棒条线的顶端是一个实心点。
如果需要在序列的前后补较多的零的话,可以利用函数zeros(),其语法为:
zeros(1,N):
圆括号中的1和N表示该函数将产生一个一行N列的矩阵,矩阵中的所有元素均为零。
利用这个矩阵与序列x[n]进行组合,从而得到一个长度与n相等的向量。
例如,当x[n]={0.1,1.1,-1.2,0,1.3}时,为了得到程序Program1_3中的序列,
↑n=0
可以用这个MATLAB语句x=[zeros(1,4)xzeros(1,2)]来实现。
用这种方法编写的程序如下:
%Program1_4
%Thisprogramisusedtogenerateadiscrete-timesinusoidalsignal
%anddrawitsplot
clear,%Clearallvariables
closeall,%Closeallfigurewindows
n=-5:
5;%Specifytheintervaloftime
x=[zeros(1,4),0.1,1.1,-1.2,0,1.3,zeros(1,2)];%Generatethesequence
stem(n,x,'.')%Openafigurewindowanddrawtheplotofx[n]
gridon,
title('Adiscrete-timesequencex[n]')
xlabel('Timeindexn')
离散时间单位阶跃信号u[n]定义为
1.2
离散时间单位阶跃信号u[n]除了也可以直接用前面给出的扩展函数来产生,还可以利用MATLAB内部函数ones(1,N)来实现。
这个函数类似于zeros(1,N),所不同的是它产生的矩阵的所有元素都为1。
单位冲激序列:
n=-10:
10;
y=[zeros(1,10),1,zeros(1,10)];
stem(n,y)
单位阶跃序列:
n=-10:
10;
y=[zeros(1,10),ones(1,11)];
stem(n,y)
值得注意的是,利用ones(1,N)来实现的单位阶跃序列并不是真正的单位阶跃序列,而是一个长度为N单位门(Gate)序列,也就是u[n]-u[n-N]。
但是在一个有限的图形窗口中,我们看到的还是一个单位阶跃序列。
三、实验内容及步骤
Q1-1:
修改程序Program1_1,将dt改为0.2,再执行该程序,保存图形,看看所得图形的效果如何?
dt=0.01时的信号波形dt=0.2时的信号波形
此处粘贴图形此处粘贴图形
这两幅图形有什么区别,哪一幅图形看起来与实际信号波形更像?
答:
Q1-2:
修改程序Program1_1,并以Q1_2为文件名存盘,产生实指数信号x(t)=e-2t。
要求在图形中加上网格线,并使用函数axis()控制图形的时间范围在0~2秒之间。
然后执行该程序,保存所的图形。
修改Program1_1后得到的程序Q1_2如下:
信号x(t)=e-2t的波形图
Q1-3:
将实验原理中所给的单位冲激信号和单位阶跃信号的函数文件在MATLAB文件编辑器中编写好,并分别以以文件名delta和u存入work文件夹中以便于使用。
抄写函数文件delta如下:
抄写函数文件u如下:
Q1-4:
修改程序Program1_4,并以Q1_4为文件名存盘,利用axis()函数,将图形窗口的横坐标范围改为-2≤n≤5,纵坐标范围改为-1.5≤x≤1.5。
修改Program1_4后得到的程序Q1_4如下:
信号的波形图
此处粘贴图形
实验二连续系统的时域分析
一、实验目的
1.掌握信号的时域变换。
2.信号的时域运算。
3.掌握卷积的计算。
4.掌握微分方程的求解。
二、实验原理
1信号的时域变换
1.1信号的时移
信号的时移可用下面的数学表达式来描述:
设一个连续时间信号为x(t),它的时移y(t)表示为:
y(t)=x(t-t0)1.3
其中,t0为位移量。
若t0为正数,则y(t)等于将x(t)右移t0秒之后的结果。
反之,若t0为负数,则y(t)等于将x(t)左移t0秒之后的结果。
在MATLAB中,时移运算与数学上习惯表达方法完全相同。
程序Program1_5对给定一个连续时间信号x(t)=e-0.5tu(t),对它分别左移2秒钟和右移2秒钟得到信号x1(t)=e-0.5(t+2)u(t+2)和x2(t)=e-0.5(t-2)u(t-2)。
%Program1_5
%Thisprogramisusedtoimplementthetime-shiftoperation
%onacontinuous-timesignalandtoobtainitstime-shiftedversions
%andtodrawtheirplots.
clear,closeall,
t=-5:
0.01:
5;
x=exp(-0.5*t).*u(t);%Generatetheoriginalsignalx(t)
x1=exp(-0.5*(t+2)).*u(t+2);%Shiftx(t)totheleftby2secondtogetx1(t)
x2=exp(-0.5*(t-2)).*u(t-2);%Shiftx(t)totherightby2secondtogetx2(t)
subplot(311)
plot(t,x)%Plotx(t)
gridon,
title('Originalsignalx(t)')
subplot(312)
plot(t,x1)%Plotx1(t)
gridon,
title('Leftshiftedversionofx(t)')
subplot(313)
plot(t,x2)%Plotx2(t)
gridon,
title('Rightshiftedversionofx(t)')
xlabel('Timet(sec)')
1.2信号的时域反褶
对一个信号x[n]的反褶运算在数学上表示为
y[n]=x[-n]1.4
这种反褶运算,用MATLAB实现起来也是非常简单的。
有多种方法可以实现信号的反褶运算。
方法一,修改绘图函数plot(t,x)和stem(n,x)中的时间变量t和n,即用-t和-n替代原来的t和n,这样绘制出来的图形,看起来就是原信号经时域反褶后的版本。
方法二,直接利用原信号与其反褶信号的数学关系式来实现。
这种方法最符合信号反褶运算的实际意义。
1.3信号的时域尺度变换
信号x(t)的时域尺度变换在数学描述为
y(t)=x(at),1.5
其中a为任意常数。
根据a的不同取值,这种时域尺度变换对信号x(t)具有非常不同的影响。
当a=1时,y(t)=x(t);
当a=-1时,y(t)=x(-t),即y(t)可以通过将x(t)反褶运算而得到;
当a>1时,y(t)=x(at),y(t)是将x(t)在时间轴上的压缩而得到;
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 信号 系统 matlab 实验 指导书