2022年信号与系统实验报告实验一信号与系统的时域分析.pdf
《2022年信号与系统实验报告实验一信号与系统的时域分析.pdf》由会员分享,可在线阅读,更多相关《2022年信号与系统实验报告实验一信号与系统的时域分析.pdf(20页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、信号与系统实验报告实验一信号与系统的时域分析实验一信号与系统的时域分析一、实验目的1、熟悉与掌握常用的用于信号与系统时域仿真分析的MA TLAB 函数 ; 2、掌握连续时间与离散时间信号的MATLAB产生 ,掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MATLAB 编程 ; 3、 牢固掌握系统的单位冲激响应的概念,掌握 LTI 系统的卷积表达式及其物理意义,掌握卷积的计算方法、卷积的基本性质; 4、 掌握利用 MATLAB 计算卷积的编程方法,并利用所编写的MATLAB 程序验证卷积的常用基本性质 ; 掌握 MATLAB 描述 LTI 系统的常用方法及有关函数,并学会
2、利用MATLAB求解 LTI 系统响应 ,绘制相应曲线。基本要求:掌握用MATLAB描述连续时间信号与离散时间信号的方法,能够编写MATLAB程序 ,实现各种信号的时域变换与运算,并且以图形的方式再现各种信号的波形。掌握线性时不变连续系统的时域数学模型用MA TLAB 描述的方法 ,掌握卷积运算、 线性常系数微分方程的求解编程。二、实验原理信号 (Signal) 一般都就是随某一个或某几个独立变量的变化而变化的,例如 ,温度、压力、声音 ,还有股票市场的日收盘指数等,这些信号都就是随时间的变化而变化的,还有一些信号,例如在研究地球结构时,地下某处的密度就就是随着海拔高度的变化而变化的。一幅图片
3、中的每一个象素点的位置取决于两个坐标轴,即横轴与纵轴,因此 ,图像信号具有两个或两个以上的独立变量。在信号与系统课程中,我们只关注这种只有一个独立变量 (Independent variable) 的信号 ,并且把这个独立变量统称为时间变量 (Time variable) ,不管这个独立变量就是否就是时间变量。在自然界中 ,大多数信号的时间变量都就是连续变化的,因此这种信号被称为连续时间信号(Continuous-Time Signals)或模拟信号 (Analog Signals) ,例如前面提到的温度、压力与声音信号就就是连续时间信号的例子。但就是,还有一些信号的独立时间变量就是离散变化的
4、,这种信号称为离散时间信号。前面提到的股票市场的日收盘指数,由于相邻两个交易日的日收盘指数相隔24 小时 ,这意味着日收盘指数的时间变量就是不连续的,因此日收盘指数就是离散时间信号。而系统则用于对信号进行运算或处理,或者从信号中提取有用的信息,或者滤出信号中某些无用的成分 ,如滤波 ,从而产生人们所希望的新的信号。系统通常就是由若干部件或单元组成的一个整体(Entity) 。系统可分为很多不同的类型,例如 ,根据系统所处理的信号的不同,系统可分为 连续时间系统(Continuous-time system)与离散时间系统(Discrete-time system),根精品资料 - - - 欢迎
5、下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 1 页,共 20 页 - - - - - - - - - - 信号与系统实验报告实验一信号与系统的时域分析据系统所具有的不同性质,系统又可分为因果系统 (Causal system)与非因果系统(Noncausal system)、稳定系统 (Stable system) 与不稳定系统 (Unstable system) 、线性系统 (Linear system)与 非 线 性 系 统 (Nonlinear system) 、 时 变 系 统 (Time-variant system
6、) 与 时 不 变 系 统(Time-invariant system)等等。然而 ,在信号与系统与数字信号处理中,我们所分析的系统只就是所谓的线性时不变系统,这种系统同时满足两个重要的基本性质,那就就是线性性与时不变性,通常称为线性时不变(LTI) 系统。1、 信号的时域表示方法1、1 将信号表示成独立时间变量的函数例如x(t)=sin( t) 与 xn=n(0 、5)nun 分别表示一个连续时间信号与一个离散时间信号。在 MATLAB 中有许多内部函数,可以直接完成信号的这种表达,例如 : sin():正弦信号cos():余弦信号exp():指数信号1、2 用信号的波形图来描述信号用函数曲
7、线表示一个信号,图 1、1 就就是一个连续时间信号与一个离散时间信号的波形图。图 1、1 连续时间信号与离散时间信号的波形图1、3 将信号用一个数据序列来表示对于离散时间信号,还可以表示成一个数的序列,例如 : xn= 、 、 、 、, 0、1, 1、1, -1、2, 0, 1、3, 、 n=0 在信号与系统与数字信号处理课程中,上述三种信号的描述方法就是经常要使用的。2 用 MATLAB仿真连续时间信号与离散时间信号在 MATLAB中,无论就是连续时间信号还就是离散时间信号,MATLAB都就是用一个数精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归
8、纳 - - - - - - - - - -第 2 页,共 20 页 - - - - - - - - - - 信号与系统实验报告实验一信号与系统的时域分析字序列来表示信号,这个数字序列在MATLAB中叫做向量 (vector)。通常的情况下,需要与时间变量相对应。如前所述 ,MA TLAB有很多内部数学函数可以用来产生这样的数字序列,例如sin()、cos()、exp()等函数可以直接产生一个按照正弦、余弦或指数规律变化的数字序列。2、1 连续时间信号的仿真程序 Program1_1 就是用MATLAB对一个正弦信号进行仿真的程序,请仔细阅读该程序 ,并在计算机上运行,观察所得图形。% Prog
9、ram1_1 % This program is used to generate a sinusoidal signal and draw its plot clear, % Clear all variables close all, % Close all figure windows dt = 0、01; % Specify the step of time variable t = -2:dt:2; % Specify the interval of time x = sin(2*pi*t); % Generate the signal plot(t,x) % Open a figu
10、re window and draw the plot of x(t) title(Sinusoidal signal x(t) xlabel(Time t (sec) 常用的图形控制函数axis(xmin,xmax,ymin,ymax):图型显示区域控制函数,其中 xmin 为横轴的显示起点,xmax为横轴的显示终点,ymin 为纵轴的显示起点,ymax 为纵轴的显示终点。有时 ,为了使图形具有可读性,需要在所绘制的图形中,加上一些网格线来反映信号的幅度大小。 MATLAB中的 grid on/grid off可以实现在您的图形中加网格线。grid on:在图形中加网格线。grid off:
11、 取消图形中的网格线。x = input( Type in signal x(t) in closed form: ) 在信号与系统课程中,单位阶跃信号u(t) 与单位冲激信号 (t) 就是二个非常有用的信号。它们的定义如下0, 0)(1)(ttdttt1、1(a) 0, 00, 1)(tttu1、1(b) 这里分别给出相应的简单的产生单位冲激信号与单位阶跃信号的扩展函数。产生单位冲激信号的扩展函数为: function y = delta(t) dt = 0、01; y = (u(t)-u(t-dt)/dt; 产生单位阶跃信号的扩展函数为: % Unit step function func
12、tion y = u(t) y = (t=0); % y = 1 for t 0, else y = 0 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 3 页,共 20 页 - - - - - - - - - - 信号与系统实验报告实验一信号与系统的时域分析请将这二个MATLAB函数分别以delta 与 u 为文件名保存在work 文件夹中 ,以后 ,就可以像教材中的方法使用单位冲激信号 (t) 与单位阶跃信号u(t)。2、2 离散时间信号的仿真程序 Program1_2 用来产生离散时间信号xn=si
13、n(0 、2n) 。% Program1_2 % This program is used to generate a discrete-time sinusoidal signal and draw its plot clear, % Clear all variables close all, % Close all figure windows n = -10:10; % Specify the interval of time x = sin(0、2*pi*n); % Generate the signal stem (n,x) % Open a figure window and d
14、raw the plot of xn title (Sinusoidal signal xn) xlabel (Time index n) 请仔细阅读该程序,比较程序Program1_1 与 Program1_2 中的不同之处 ,以便自己编程时能够正确使用这种方法方针连续时间信号与离散时间信号。程序 Program1_3 用来仿真下面形式的离散时间信号: xn= 、 、 、 、, 0、1, 1、1, -1、2, 0, 1、3, 、 n=0 % Program1_3 % This program is used to generate a discrete-time sequence % and
15、 draw its plot clear, % Clear all variables close all, % Close all figure windows n = -5:5; % Specify the interval of time, the number of points of n is 11、x = 0, 0, 0, 0, 0 、1, 1、1, -1、2, 0, 1、3, 0, 0; % Generate the signal stem(n,x,、 ) % Open a figure window and draw the plot of xn grid on, title
16、(A discrete-time sequence xn) xlabel (Time index n) 由于在程序的stem(n,x,、) 语句中加有 、选项 ,因此绘制的图形中每根棒条线的顶端就是一个实心点。如果需要在序列的前后补较多的零的话,可以利用函数zeros(),其语法为 : zeros(1, N): 圆括号中的1 与 N 表示该函数将产生一个一行N 列的矩阵 ,矩阵中的所有元素均为零。利用这个矩阵与序列xn 进行组合 ,从而得到一个长度与n 相等的向量。例如 ,当xn= 0 、1, 1、1, -1、 2, 0, 1 、3 时,为了得到程序Program1_3 中的序列, n=0 可
17、以用这个MATLAB语句 x = zeros(1,4) x zeros(1, 2) 来实现。用这种方法编写的程序如下: % Program1_4 % This program is used to generate a discrete-time sinusoidal signal % and draw its plot clear, % Clear all variables 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 4 页,共 20 页 - - - - - - - - - - 信号与系统实验报告
18、实验一信号与系统的时域分析close all, % Close all figure windows n = -5:5; % Specify the interval of time x = zeros(1,4), 0 、1, 1、1, -1、2, 0, 1、3, zeros(1,2); % Generate the sequence stem (n,x,、) % Open a figure window and draw the plot of xn grid on, title (A discrete-time sequence xn) xlabel (Time index n) 离散时间
19、单位阶跃信号un定义为0, 00, 1nnnu1、2 离散时间单位阶跃信号un除了也可以直接用前面给出的扩展函数来产生,还可以利用MATLAB内部函数ones(1,N) 来实现。 这个函数类似于zeros(1,N),所不同的就是它产生的矩阵的所有元素都为1。值得注意的就是,利用ones(1,N) 来实现的单位阶跃序列并不就是真正的单位阶跃序列,而就是一个长度为N 单位门 (Gate)序列 ,也就就是un-un-N 。 但就是在一个有限的图形窗口中,我们瞧到的还就是一个单位阶跃序列。在绘制信号的波形图时,有时我们需要将若干个图形绘制在图一个图形窗口中,这就需要使用 MA TLAB 的图形分割函数
20、subplot(), 其用法就是在绘图函数stem或 plot 之前 ,使用图形分割函数subplot(n1,n2,n3), 其中的参数n1,n2 与 n3 的含义就是 ,该函数将把一个图形窗口分割成 n1xn2 个子图 ,即将绘制的图形将绘制在第n3 个子图中。2、3信号的时域变换2、3、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
21、秒之后的结果。在 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 % This program is used to implement the time-shift operation % on a continuous-time signal and to obtain its time-shifted versions % a
22、nd to draw their plots 、clear,close all, t = -5:0、01:5; x = exp(-0、5*t) 、*u(t); % Generate the original signal x(t) x1 = exp(-0、5*(t+2) 、*u(t+2); % Shift x(t) to the left by 2 second to get x1(t) x2 = exp(-0、5*(t-2) 、*u(t-2); % Shift x(t) to the right by 2 second to get x2(t) subplot(311) plot(t,x)
23、% Plot x(t) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 5 页,共 20 页 - - - - - - - - - - 信号与系统实验报告实验一信号与系统的时域分析grid on, title (Original signal x(t) subplot (312) plot (t,x1) % Plot x1(t) grid on, title (Left shifted version of x(t) subplot (313) plot (t,x2) % Plot x2(t) grid o
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2022 信号 系统 实验 报告 时域 分析
限制150内