离散时间信号的时域分析实验报告

Matlab离散信号时域分析实验报告1. 引言信号是信息的载体，可以通过对信号进行离散化来进行数字信号的处理和分析。

离散信号时域分析是对离散信号在时域上进行观察和分析的方法。

本实验旨在通过使用Matlab软件对离散信号进行时域分析，掌握离散信号的时域特性和基本分析方法。

2. 实验目的•了解离散信号的概念和特性；•掌握离散信号的时域分析方法；•学会通过Matlab对离散信号进行时域分析。

3. 实验原理离散信号是在时间上呈现离散的特征，可以用离散序列表示。

离散序列可以通过采样连续信号得到，也可以通过数学模型生成。

在时域分析中，通常使用的分析方法包括： - 时域图像绘制：绘制离散信号的时域图像，了解信号的振幅和波形特征； - 时域序列计算：计算离散信号的均值、方差等统计量，了解信号的基本特性；- 时域滤波：对离散信号进行滤波，去除噪声或者突发干扰。

4. 实验步骤4.1 生成离散信号首先需要生成一个离散信号序列，可以使用Matlab的随机数函数生成一个大小为N的随机序列作为离散信号。

N = 100; % 信号长度为100x = rand(1,N); % 生成随机序列4.2 时域图像绘制通过plot函数可以将离散信号在时域上绘制出来，观察信号的振幅和波形特征。

t = 1:N; % 时间序列plot(t, x);title('离散信号时域图像');xlabel('时间');ylabel('幅度');4.3 时域序列计算可以通过内置函数计算离散信号的均值、方差等统计量。

avg = mean(x); % 均值variance = var(x); % 方差4.4 时域滤波可以使用滤波器对离散信号进行滤波，去除噪声或者突发干扰。

这里以均值滤波为例，对信号进行平滑处理。

windowSize = 5; % 滑动窗口大小b = (1/windowSize)*ones(1,windowSize);a = 1;smoothed_x = filter(b, a, x);5. 实验结果与分析通过对生成的离散信号进行时域分析，得到如下结果： - 时域图像：时域图像时域图像•信号均值：0.5231•信号方差：0.0842•平滑后的信号时域图像：平滑后的时域图像平滑后的时域图像从时域图像可以观察到信号的振幅和波形特征。

实验十一 离散信号时域分析的MATLAB 实现一、实验目的1. 熟悉MATLAB 编程方法、常用语句和可视化绘图技术；2. 掌握序列时域运算的MATLAB 编程方法。

二、实验原理在用MATLAB 表示离散信号并将其可视化时，由于矩阵元素个数是有限的，因此无法表示无限长序列；另外，离散信号无法进行符号运算。

在MATLAB 中，绘制离散序列波形图的专用命令为stem( )。

其格式有：（1）stem(k,f)在图形窗口中，绘制出样值顶部为空心圆的序列f (k)波形图。

（2）stem(k,f ,’fill’)在图形窗口中，绘制出样值顶部为实心圆的序列f (k)波形图。

下面介绍离散序列的MATLAB 表示、基本运算（相加、相乘、平移、反转、尺度变换）、卷积和的实现及其图形显示方法。

1.单位序列δ(k ) 单位序列的定义：⎩⎨⎧≠==0,00,1)(k k k δ下面为绘制δ(k-k 0)波形图的子程序：function impseq(k1,k2,k0) %单位序列δ(k-k0),k0为时移量k=k1:k2; %k1,k2为序列的起止序列号 n=length(k); x=zeros(1,n);x(1,k0-k1+1)=1; %在k0时刻信号赋值为1 stem(k,x,'fill') axis([k1,k2,0,1.1])title('单位序列d(k-k0)')输入如下命令，则可获得单位序列δ(k-3)的波形图，如图11-1所示。

impseq(-1,5,3)图11-12.单位阶跃序列ε(k ) 单位序列的定义：⎩⎨⎧<≥=0,00,1)(k k k ε下面为绘制ε(k-k 0)波形图的MATLAB 子程序。

function stepseq(k1,k2,k0) %单位阶跃序列,k0为时移量 k=k1:k0-1; %k1,k2为序列的起止序列号 kk=length(k);x=zeros(1,kk); %k0前信号赋值为0 stem(k,x,'fill') %绘出k1~k0-1的波形(0值) hold onn=k0:k2; nn=length(n);x=ones(1,nn); %k0后信号赋值为1 stem(n,x,'fill') %绘出k0~k2的波形(1值) hold offaxis([k1,k2,0,1.1]) title('单位阶跃序列')运行如下命令，则可获得单位序列ε(k-3)的波形图，如图11-2所示。

第1篇一、实验目的1. 理解时域离散信号的基本概念和特性。

2. 掌握时域离散信号的表示方法。

3. 熟悉常用时域离散信号的产生方法。

4. 掌握时域离散信号的基本运算方法。

5. 通过MATLAB软件进行时域离散信号的仿真分析。

二、实验原理时域离散信号是指在时间轴上取离散值的一类信号。

这类信号在时间上不连续，但在数值上可以取到任意值。

时域离散信号在数字信号处理领域有着广泛的应用，如通信、图像处理、语音处理等。

时域离散信号的基本表示方法有：1. 序列表示法：用数学符号表示离散信号，如 \( x[n] \) 表示离散时间信号。

2. 图形表示法：用图形表示离散信号，如用折线图表示序列。

3. 时域波形图表示法：用波形图表示离散信号，如用MATLAB软件生成的波形图。

常用时域离散信号的产生方法包括：1. 单位阶跃信号：表示信号在某个时刻发生突变。

2. 单位冲激信号：表示信号在某个时刻发生瞬时脉冲。

3. 正弦信号：表示信号在时间上呈现正弦波形。

4. 矩形脉冲信号：表示信号在时间上呈现矩形波形。

时域离散信号的基本运算方法包括：1. 加法：将两个离散信号相加。

2. 乘法：将两个离散信号相乘。

3. 卷积：将一个离散信号与另一个离散信号的移位序列进行乘法运算。

4. 反褶：将离散信号沿时间轴翻转。

三、实验内容1. 实验一：时域离散信号的表示方法（1）使用序列表示法表示以下信号：- 单位阶跃信号：\( u[n] \)- 单位冲激信号：\( \delta[n] \)- 正弦信号：\( \sin(2\pi f_0 n) \)- 矩形脉冲信号：\( \text{rect}(n) \)（2）使用图形表示法绘制以上信号。

2. 实验二：时域离散信号的产生方法（1）使用MATLAB软件生成以下信号：- 单位阶跃信号- 单位冲激信号- 正弦信号（频率为1Hz）- 矩形脉冲信号（宽度为2）（2）观察并分析信号的波形。

3. 实验三：时域离散信号的基本运算（1）使用MATLAB软件对以下信号进行加法运算：- \( u[n] \)- \( \sin(2\pi f_0 n) \)（2）使用MATLAB软件对以下信号进行乘法运算：- \( u[n] \)- \( \sin(2\pi f_0 n) \)（3）使用MATLAB软件对以下信号进行卷积运算：- \( u[n] \)- \( \sin(2\pi f_0 n) \)（4）使用MATLAB软件对以下信号进行反褶运算：- \( u[n] \)4. 实验四：时域离散信号的仿真分析（1）使用MATLAB软件对以下系统进行时域分析：- 系统函数：\( H(z) = \frac{1}{1 - 0.5z^{-1}} \)（2）观察并分析系统的单位冲激响应。

一、实验目的1. 理解数字信号处理的基本概念和原理。

2. 掌握离散时间信号的基本运算和变换方法。

3. 熟悉数字滤波器的设计和实现。

4. 培养实验操作能力和数据分析能力。

二、实验原理数字信号处理（Digital Signal Processing，DSP）是利用计算机对信号进行采样、量化、处理和分析的一种技术。

本实验主要涉及以下内容：1. 离散时间信号：离散时间信号是指时间上离散的信号，通常用序列表示。

2. 离散时间系统的时域分析：分析离散时间系统的时域特性，如稳定性、因果性、线性等。

3. 离散时间信号的变换：包括离散时间傅里叶变换（DTFT）、离散傅里叶变换（DFT）和快速傅里叶变换（FFT）等。

4. 数字滤波器：设计、实现和分析数字滤波器，如低通、高通、带通、带阻滤波器等。

三、实验内容1. 离散时间信号的时域运算（1）实验目的：掌握离散时间信号的时域运算方法。

（2）实验步骤：a. 使用MATLAB生成两个离散时间信号；b. 进行时域运算，如加、减、乘、除等；c. 绘制运算结果的时域波形图。

2. 离散时间信号的变换（1）实验目的：掌握离散时间信号的变换方法。

（2）实验步骤：a. 使用MATLAB生成一个离散时间信号；b. 进行DTFT、DFT和FFT变换；c. 绘制变换结果的频域波形图。

3. 数字滤波器的设计和实现（1）实验目的：掌握数字滤波器的设计和实现方法。

（2）实验步骤：a. 设计一个低通滤波器，如巴特沃斯滤波器、切比雪夫滤波器等；b. 使用MATLAB实现滤波器；c. 使用MATLAB对滤波器进行时域和频域分析。

4. 数字滤波器的应用（1）实验目的：掌握数字滤波器的应用。

（2）实验步骤：a. 采集一段语音信号；b. 使用数字滤波器对语音信号进行降噪处理；c. 比较降噪前后的语音信号，分析滤波器的效果。

四、实验结果与分析1. 离散时间信号的时域运算实验结果显示，通过MATLAB可以方便地进行离散时间信号的时域运算，并绘制出运算结果的时域波形图。

离散时间信号与系统的时域分析实验报告报告⼆：⼀、设计题⽬1．绘制信号)()(1k k f δ=和)2()(2-=k k f δ的波形2.绘制直流信号)()(1k k f ε=和)2(2-=k f ε的波形3绘制信号)()(6k G k f =的波形⼆实验⽬的1．掌握⽤MATLAB 绘制离散时间信号（序列）波形图的基本原理。

2．掌握⽤MATLAB 绘制典型的离散时间信号（序列）。

3．通过对离散信号波形的绘制与观察，加深理解离散信号的基本特性。

三、设计原理离散时间信号(也称为离放序列)是指在时间上的取值是离散的，只在⼀些离放的瞬间才有定义的,⽽在其他时间没有定义，简称离放信号(也称为离散序列) 序列的离散时间间隔是等间隔(均匀)的，取时间间隔为T.以f(kT)表⽰该离散序列,k 为整数(k=0,±1.±2，...)。

为了简便，取T=1.则f(kT)简记为f(k)， k 表⽰各函数值在序列中出现的序号。

序列f(k)的数学表达式可以写成闭合形式,也可逐⼀列出f(k)的值。

通常,把对应某序号K0的序列值称为序列的第K0个样点的“样点值”。

四、设计的过程及仿真1clear all; close all; clc;k1=-4;k2=4;k=k1:k2;n1=0;n2=2;f1=[(k-n1)==0];f2=[(k-n2)==0];subplot(1,2,1)stem(k,f1,'fill','-k','linewidth',2);xlabel('k');ylabel('f_1(k)');title('δ(k)')axis([k1,k2,-0.1,1.1]);subplot(1,2,2)stem(k,f2,'filled','-k','linewidth',2);ylabel('f_2(k)');title('δ(k-2)')axis([k1,k2,-0.1,1.1]);程序运⾏后，仿真绘制的结果如图所⽰:2c lear all; close all; clc;k1=-2;k2=8;k=k1:k2;n1=0;n2=2; %阶跃序列开始出现的位置f1=[(k-n1)>=0]; f2=[(k-n2)>=0];subplot(1,2,1)stem(k,f1,'fill','-k','linewidth',2);xlabel('k');ylabel('f_1(k)');title('ε(k)')axis([k1,k2+0.2,-0.1,1.1])subplot(1,2,2)stem(k,f2,'filled','-k','linewidth',2);xlabel('k');ylabel('f_2(k)');title('ε(k-2)')axis([k1,k2+0.2,-0.1,1.1]);程序运⾏后，仿真绘制的结果如图所⽰:3clear all; close all; clc;k1=-2;k2=7;k=k1:k2; %建⽴时间序列n1=0;n2=6; f1=[(k-n1)>=0];f2=[(k-n2)>=0];f=f1-f2;stem(k,f,'fill','-k','linewidth',2);xlabel('k');ylabel('f(k)');title('G_6(k)')axis([k1,k2,-0.1,1.1]);程序运⾏后，仿真绘制的结果如图所⽰:五、设计的结论及收获实现了⽤matlab绘制离散时间信号, 通过对离散信号波形的绘制与观察，加深理解离散信号的基本特性。

离散时间信号的时域分析实验报告实验名称：离散时间信号的时域分析⼀、实验⽬的1.学会⽤MATLAB在时域中产⽣⼀些基本的离散时间信号，并对这些信号进⾏⼀些基本的运算。

2.熟悉MATLAB中产⽣信号和绘制信号的基本命令。

⼆、实验内容1.在载波信号xH[n]和调制信号xL[n]采⽤不同频率、不同调制指数m的情况下，运⾏程序P1.6，以产⽣振幅调制信号y[n]。

2.编写matlab程序，以产⽣图1.1和图1.2所⽰的⽅波和锯齿波序列，并将序列绘制出来。

三、主要算法与程序1.n=0:100;m=0.6;fH=0.2;fL=0.02;xH=sin(2*pi*fH*n);xL=sin(2*pi*fL*n);y=(1+m*xL).*xH;stem(n,y);grid;xlabel('时间序列');ylabel('振幅');通过改变m,fH和fL来产⽣不同情况下的振幅调制信号。

2.画出图⼆：n=0:1:30;y=3*square(n*pi/5,60);stem(n,y),grid onaxis([0,30,-4,4]);xlabel('时间序号n');ylabel('振幅');为画出图三，将占空⽐由图⼆的60改为30。

画出图四：n=0:1:50;y=2*sawtooth(n*pi/10,1);stem(n,y),grid onaxis([0,50,-2,2]);xlabel('时间序号n');ylabel('振幅');为画出图五，将图四中从-1到1的范围由1改为0.5。

四、实验结果与分析图⼀确定了数值：m=0.6，fH=0.2，fL=0.02，绘出图像。

图⼆图三图四图五五、实验⼩结通过这次实验，我熟悉MATLAB中产⽣信号和绘制信号的基本命令，学会⽤MATLAB在时域中产⽣⼀些基本的离散时间信号，并对这些信号进⾏⼀些基本的运算。

信号与系统实验报告一、信号的时域基本运算1．连续时间信号的时域基本运算两实验之一实验分析：输出信号值就等于两输入信号相加（乘）。

由于b=2，故平移量为2时，实际是右移1，符合平移性质。

两实验之二心得体会：时域中的基本运算具有连续性，当输入信号为连续时，输出信号也为连续。

平移，伸缩变化都会导致输出结果相对应的平移伸缩。

2.离散时间信号的时域基本运算两实验之一实验分析：输出信号的值是对应输入信号在每个n值所对应的运算值，当进行拉伸变化后，n值数量不会变，但范围会拉伸所输入的拉伸系数。

两实验之二心得体会：离散时间信号可以看做对连续时间信号的采样，而得到的输出信号值，也可以看成是连续信号所得之后的采样值。

二、连续信号卷积与系统的时域分析1.连续信号卷积积分两实验之一实验分析：当两相互卷积函数为冲激函数时，所卷积得到的也是一个冲激函数，且该函数的冲激t值为函数x，函数y冲激t值之和。

两实验之二心得体会：连续卷积函数每个t值所对应的卷积和可以看成其中一个在k值取得的函数与另外一个函数相乘得到的一个分量函数，并一直移动k值直至最后，最后累和出来的最终函数便是所得到的卷积函数。

3.RC电路时域积分两实验之一实验分析：全响应结果正好等于零状态响应与零输入响应之和。

两实验之二心得体会：具体学习了零状态，零输入，全响应过程的状态及变化，与之前所学的电路知识联系在一起了。

三、离散信号卷积与系统的时域分析1.离散信号卷积求和两实验之一实验分析：输出结果的n值是输入结果的k号与另一个n-k的累和两实验之二心得体会：直观地观察到卷积和的产生，可以看成连续卷积的采样形式，从这个方面去想，更能深入地理解卷积以及采样的知识。

2.离散差分方程求解两实验之一实验分析：其零状态响应序列为0 0 4 5 7.5，零输入响应序列为2 4 5 5.5 5.75，全状态响应序列为2 4 9 10.5 13.25，即全状态=零输入+零状态。

两实验之二心得体会：求差分方程时，可以根据全状态响应是由零输入输入以及零状态相加所得，分开来求，同时也加深了自己对差分方程的求解问题的理解。

武汉工程大学数字信号处理实验报告二专业班级：14级通信03班学生姓名：秦重双学号：1404201114实验时间：2017年5月3日实验地点：4B315指导老师: 杨述斌实验一离散时间信号的分析实验一、实验目的①认识常用的各种信号,理解其数学表达式和波形表示。

②掌握在计算机中生成及绘制数值信号波形的方法。

③掌握序列的简单运算及计算机实现与作用。

④理解离散时间傅里叶变换、Z变换及它们的性质和信号的频域特性。

二、实验设备计算机,MATLAB语言环境。

三、实验基础理论1、序列的相关概念离散时间信号用一个称为样本的数字序列来表示。

一般用｛x［n］｝表示，其中自变量n的取值范围是﹣∞到﹢∞之间的整数。

为了表示方便，序列通常直接用x[n]表示。

离散时间信号可以是一个有限长序列，也可以是一个无限长序列。

有限长（也称为有限时宽)序列仅定义在有限的时间间隔中：﹣∞≤N1 ≤N2 ≤+∝。

有限长序列的长度或时宽为N=N1 -N2+1。

满足x[n+kN］=x[n]（对于所有n）的序列称为周期为N的周期序列,其中N取任意正整数；k取任意整数；2、常见序列常见序列有单位取样值信号、单位阶跃序列、矩形序列、斜变序列、单边指数序列、正弦序列、复指数序列等。

3、序列的基本运算序列的基本运算有加法、乘法、倒置（反转）、移位、尺度变换、卷积等。

4、离散傅里叶变换的相关概念5、Z变换的相关概念四．实验内容与步骤1、知识准备认真复习以上基础理论,理解本实验所用到的实验原理。

2、离散时间信号（序列）的产生利用MATLAB语言编程和绘制单位样值信号、单位阶跃序列、指数序列、正弦序列及随机离散信号的波形，以加深对离散信号时域表示的理解。

①单位取样值信号Matlab程序x=0;y=1;stem(x，y）；title（'单位样值’）；axis([—2,2，0,1]）；②单位阶跃序列Matlab程序n0=0;n1=—5;n2=5；n=[n1:n2］;x=［（n—n0）>=0];stem(n，x）；xlabel（'n'）；ylabel（’x(n)’）；title(’单位阶跃序列’）;③指数序列、正弦序列Matlab程序n=［0：10］;x=(1/3)。

matlab离散信号时域分析实验报告Matlab离散信号时域分析实验报告引言：离散信号时域分析是数字信号处理中的重要内容，通过对信号在时域上的分析，可以了解信号的特征和性质。

本实验使用Matlab软件进行离散信号的时域分析，通过实验数据的采集和处理，探索信号的频率、幅度、相位等重要参数。

实验目的：1. 了解离散信号的基本概念和性质；2. 掌握Matlab软件在离散信号时域分析中的应用；3. 分析离散信号的频谱特性和时域波形。

实验步骤：1. 信号采集与导入首先，我们需要采集一段离散信号的数据，并将其导入Matlab中进行分析。

在实验中，我们选择了一个简单的正弦信号作为实验对象。

通过Matlab中的数据采集工具，我们可以方便地获取该信号的采样数据，并导入到Matlab中。

2. 时域波形绘制在Matlab中，我们可以使用plot函数绘制离散信号的时域波形。

通过对信号的采样数据进行绘制，我们可以直观地观察到信号的变化规律。

同时，我们还可以通过调整绘图参数，比如线型、颜色等，使得波形图更加美观。

3. 信号频谱分析离散信号的频谱分析是了解信号频率特性的重要手段。

在Matlab中，我们可以使用fft函数对信号进行频谱分析。

通过对信号的采样数据进行傅里叶变换，我们可以得到信号的频谱图。

频谱图可以清晰地展示信号的频率分布情况，对于分析信号的频率成分非常有帮助。

4. 信号幅度与相位分析除了频率特性，离散信号的幅度和相位也是需要关注的重要参数。

在Matlab中，我们可以使用abs函数计算信号的幅度，使用angle函数计算信号的相位。

通过对信号的采样数据进行计算，我们可以得到信号的幅度和相位信息。

这些信息对于了解信号的时域特性非常有帮助。

实验结果与分析：通过以上实验步骤，我们得到了离散信号的时域波形、频谱特性、幅度和相位信息。

通过观察实验结果，我们可以发现信号的频率成分、幅度变化以及相位差异等重要特征。

结论：通过本次实验，我们深入了解了离散信号的时域分析方法，并通过Matlab软件进行了实际操作。

《数字信号处理》 实验报告学院 专业 电子信息工程 班级 姓名 学号 时间实验一 时域离散信号与系统分析一、实验目的1、熟悉连续信号经理想采样后的频谱变化关系，加深对时域采样定理的理解。

2、熟悉时域离散系统的时域特性，利用卷积方法观察分析系统的时域特性。

3、学会离散信号及系统响应的频域分析。

4、学会时域离散信号的MATLAB 编程和绘图。

5、学会利用MATLAB 进行时域离散系统的频率特性分析。

二、实验内容1、序列的产生（用Matlab 编程实现下列序列（数组），并用stem 语句绘出杆图。

（要求标注横轴、纵轴和标题）（1）. 单位脉冲序列x(n)=δ（n ） （2）. 矩形序列x(n)=R N (n) ,N=10nδ(n )nR N (n )图1.1 单位脉冲序列 图1.2 矩形序列(3) . x(n)=e （0.8＋3j ）n ； n 取0－15。

4n|x (n )|201321111053 陈闽焜n<x (n )/R a d图1.3 复指数序列的 模 图1.4 复指数序列的 相角（4）. x(n)=3cos （0. 25πn ＋0.3π）＋2sin （0.125πn ＋0.2π） n 取0－15。

ny (n )图1.4 复合正弦实数序列（5）. 把第（3）小题的复指数x(n)周期化，周期20点，延拓3个周期。

4m|y (m )|201321111053 陈闽焜图1.5 第（3）的20点周期延拓杆图（6）. 假设x(n)= [1，－3，2，3，-2 ]， 编程产生以下序列并绘出杆图：y(n) y(n)= x(n)－2x(n+1)+x(n-1)+x(n-3)；5201321111053 陈闽焜图1.6 y(n)序列杆图（7）、编一个用户自定义matlab 函数，名为stepshf （n0，n1，n2）实现单位阶跃序列u[n －n1]。

其中位移点数n1在起点n0和终点n2之间任意可选。

自选3个入口参数产生杆图。

离散信号与系统的时域分析实验报告1. 引言离散信号与系统是数字信号处理中的重要基础知识，它涉及信号的采样、量化和表示，以及离散系统的描述和分析。

本实验通过对离散信号在时域下的分析，旨在加深对离散信号与系统的理解。

在实验中，我们将学习如何采样和显示离散信号，并通过时域分析方法分析信号的特性。

2. 实验步骤2.1 信号的采样与显示首先，我们需要准备一个模拟信号源，例如函数发生器，来产生一个连续时间域的模拟信号。

通过设置函数发生器的频率和振幅，我们可以产生不同的信号。

接下来，我们需要使用一个采样器来对模拟信号进行采样，将其转化为离散时间域的信号。

使用合适的采样率，我们可以准确地获取模拟信号的离散样本。

最后，我们将采样后的信号通过合适的显示设备进行显示，以便观察和分析。

2.2 信号的观察与分析在实验中，我们可以选择不同类型的模拟信号，例如正弦波、方波或脉冲信号。

通过观察采样后的离散信号，我们可以观察到信号的周期性、频率、振幅等特性。

通过对不同频率和振幅的信号进行采样，我们可以进一步研究信号与采样率之间的关系，例如采样定理等。

2.3 信号的变换与滤波在实验中，我们可以尝试对采样后的离散信号进行变换和滤波。

例如，在频域下对信号进行离散傅里叶变换（DFT），我们可以将时域信号转换为频域信号，以便观察信号的频谱特性。

通过对频谱进行分析，我们可以观察到信号的频率成分和能量分布情况。

此外，我们还可以尝试使用不同的数字滤波器对离散信号进行滤波，以提取感兴趣的频率成分或去除噪声等。

3. 实验结果与分析通过实验，我们可以得到许多有关离散信号与系统的有趣结果。

例如，在观察信号的采样过程中，我们可以发现信号频率大于采样率的一半时，会发生混叠现象，即信号的频谱会发生重叠，导致采样后的信号失真。

而当信号频率小于采样率的一半时，可以还原原始信号。

此外，我们还可以观察到在频域下，正弦波信号为离散频谱，而方波信号则有更多的频率成分。

4. 结论通过本实验，我们对离散信号与系统的时域分析有了更深入的理解。

实验二离散时间信号的时域分析1.实验目的（1）学习MATLAB软件及其在信号处理中的应用，加深对常用离散时间信号的理解。

（2）利用MATLAB产生常见离散时间信号及其图形的显示，进行简单运算。

（3）熟悉MATLAB对离散信号的处理及其应用。

2.实验原理离散时间信号是时间为离散变量的信号。

其函数值在时间上是不连续的“序列”。

（1）单位抽样序列如果序列在时间轴上面有K个单位的延迟，则可以得到，即：该序列可以用MATLAB中的zeros函数来实现。

（2）正弦序列可以利用sin函数来产生。

（3）指数序列在MATLAB中通过：和来实现。

3.实验内容及其步骤（1）复习有关离散时间信号的有关内容。

（2）通过程序实现上述几种信号的产生，并进行简单的运算操作。

单位抽样序列参考：% Generation of a Unit Sample Sequenceclf;% Generate a vector from -10 to 20n = -10:20;% Generate the unit sample sequenceu = [zeros(1,10) 1 zeros(1,20)];% Plot the unit sample sequencestem(n,u);xlabel('Time index n');ylabel('Amplitude');title('Unit Sample Sequence');axis([-10 20 0 1.2]);如果序列在时间轴上面有K个单位的延迟，则可以得到，即：，通过程序来实现如下所示结果。

正弦序列参考：% Generation of a sinusoidal sequencen = 0:40; f = 0.1;phase = 0; A = 1.5;arg = 2*pi*f*n - phase; x = A*cos(arg);clf; % Clear old graphstem(n,x); % Plot the generated sequenceaxis([0 40 -2 2]); grid;title('Sinusoidal Sequence'); xlabel('Time index n');ylabel('Amplitude'); axis;指数序列参考：% Generation of a real exponential sequenceclf; n = 0:35; a = 1.2; K = 0.2;x = K*a.^n; stem(n,x);xlabel('Time index n'); ylabel('Amplitude');（3）加深对离散时间信号及其特性的理解，对于离散信号能进行基本的运算（例如信号加、乘、延迟等等），并且绘出其图形。

实验2 离散系统的时域分析一、实验目的1、熟悉并掌握离散系统的差分方程表示法；2、加深对冲激响应和卷积分析方法的理解。

二、实验原理在时域中，离散时间系统对输入信号或者延迟信号进行运算处理，生成具有所需特性的输出信号，具体框图如下:其输入、输出关系可用以下差分方程描述：输入信号分解为冲激信号,记系统单位冲激响应，则系统响应为如下的卷积计算式：当时，h[n]是有限长度的（），称系统为FIR系统；反之，称系统为IIR系统。

三、实验内容1、用MATLAB求系统响应1）卷积的实现线性移不变系统可由它的单位脉冲响应来表征。

若已知了单位脉冲响应和系统激励就可通过卷积运算来求取系统响应，即程序：x=input(‘Type in the input sequence=’); %输入xh=input(‘Type in the impulse response sequence=’); %输入hy=conv(x,h); % 对x，h进行卷积N=length(y)-1; %求出N的值n=0:1:N; %n从0开始，间隔为1的取值取到N为止disp(‘output sequence=’); disp(y); %输出ystem(n,y); %画出n为横轴，y为纵轴的离散图xlabel(‘Time index n’); ylable(‘Amplitude’); % 规定x轴y 轴的标签输入为：x=[-2 0 1 -1 3]h=[1 2 0 -1]图形:2）单位脉冲响应的求取线性时不变因果系统可用MATLAB的函数filter来仿真y=filter(b,a,x);其中，x和y是长度相等的两个矢量。

矢量x表示激励，矢量a，b 表示系统函数形式滤波器的分子和分母系数，得到的响应为矢量y。

例如计算以下系统的单位脉冲响应y(n)+0.7y(n-1)-0.45y(y-2)-0.6y(y-3)=0.8x(n)-0.44x(n-1)+0.36x(n-2)+0.02x(n-3)程序:N=input(‘Desired impuse response length=’);b=input(‘Type in the vector b=’);a=input(‘Type in the vector a=’);x=[1 zeros(1,N-1)];y=filter(b,a,x);k=0:1:N-1;stem(k,y);xlabel(’Time index n’); ylable(‘Amplitude’);输入：N=41b=[0.8 -0.44 0.36 0.02]a=[1 0.7 -0.45 -0.6]图形:2、以下程序中分别使用conv和filter函数计算h和x的卷积y和y1，运行程序，并分析y和y1是否有差别，为什么要使用x[n]补零后的x1来产生y1；具体分析当h[n]有i个值，x[n]有j个值，使用filter完成卷积功能，需要如何补零？程序:clf;h = [3 2 1 -2 1 0 -4 0 3]; %impulse responsex = [1 -2 3 -4 3 2 1]; %input sequencey = conv(h,x);n = 0:14;subplot(2,1,1);stem(n,y);xlabel('Time index n'); ylabel('Amplitude');title('Output Obtained by Convolution'); grid;x1 = [x zeros(1,8)];y1 = filter(h,1,x1);subplot(2,1,2);stem(n,y1);xlabel('Time index n'); ylabel('Amplitude');title('Output Generated by Filtering'); grid;图形:因为在y=filter（b,a,x）中，利用给定矢量a和b对x中的数据进行滤波，结果放入y矢量中，y与x长度要相等，所以要使用x[n]补零后的x1来产生y1。

matlab离散信号时域分析实验报告实验目的：本实验旨在通过使用Matlab对离散信号进行时域分析，探究离散信号的特性和变化规律，加深对信号处理理论的理解，提高实际应用能力。

实验仪器和材料：Matlab软件实验步骤：1. 生成离散信号：首先，我们使用Matlab生成一个离散信号，可以是正弦信号、方波信号等。

通过调整频率、幅度等参数，可以得到不同特性的信号。

2. 时域分析：接下来，我们对生成的离散信号进行时域分析，包括信号的幅度、相位、周期等特性进行分析，通过Matlab提供的函数和工具进行计算和可视化展示。

3. 变换分析：除了时域分析外，我们还可以对离散信号进行变换分析，如傅里叶变换、离散傅里叶变换等，通过观察频谱图和功率谱图等来分析信号的频率成分和能量分布情况。

4. 实验结果分析：最后，根据实验结果进行分析，总结离散信号的特性和变化规律，对信号处理理论进行深入理解。

实验结果：通过实验，我们得到了生成的离散信号的时域特性、频域特性等数据和图表，并对其进行了分析和总结。

我们发现不同频率、幅度的离散信号具有不同的时域特性和频域特性，这为我们理解信号处理理论提供了直观的实验数据和实例。

实验结论：通过本次实验，我们深入了解了Matlab对离散信号进行时域分析的方法和步骤，加深了对信号处理理论的理解，提高了实际应用能力。

同时，我们也发现了离散信号的特性和变化规律，为进一步的研究和应用提供了基础。

实验心得：本次实验让我对离散信号的时域分析有了更深入的理解，也提高了我在Matlab软件上的操作能力。

通过实验，我对信号处理理论有了更直观的认识，为今后的学习和研究打下了坚实的基础。

希望能够通过更多的实验和学习，不断提升自己在信号处理领域的能力和水平。

实验报告课程名称：数字信号处理授课班级：学号：姓名：指导老师：实验一离散时间信号及系统的时域分析实验类别：基础性实验1实验目的：(1）了解MA TLAB 程序设计语言的基本特点，熟悉MA TLAB软件运行环境。

(2）熟悉MA TLAB中产生信号和绘制信号的基本命令，学会用MA TLAB在时域中产生一些基本的离散时间信号，并对这些信号进行一些基本的运算。

(3）通过MA TLAB仿真一些简单的离散时间系统，并研究它们的时域特性。

(4）通过MA TLAB进行卷积运算，利用卷积方法观察分析系统的时域特性。

2. 实验报告要求●简述实验原理及目的。

●结合实验中所得给定典型序列幅频特性曲线，与理论结果比较，并分析说明误差产生的原因以及用FFT作谱分析时有关参数的选择方法。

●记录调试运行情况及所遇问题的解决方法。

3.实验内容：思考题：9.2.1 运行程序P9.2.1，哪个参数控制该序列的增长或衰减：哪个参数控制该序列的振幅？若需产生实指数序列，应对程序作何修改？9.2.2运行程序P9.2.1，该序列的频率是多少？怎样改变它？哪个参数控制该序列的相位？哪个参数可以控制该序列的振幅？该序列的周期是多少？9.2.3 运行程序P9.2.3，对加权输入得到的y(n)与在相同权系数下输出y1(n)和y2(n)相加得到的yt(n)进行比较，这两个序列是否相等？该系统是线性系统吗？9.2.4 假定另一个系统为y(n)=x(n)x(n-1)修改程序，计算这个系统的输出序列y1(n),y2(n)和y(n)。

比较有y(n)和yt(n)。

这两个序列是否相等？该系统是线性系统吗？(提高部分)9.2.5运行程序P9.2.4，并比较输出序列y(n)和yd(n-10)。

这两个序列之间有什么关系？该系统是时不变系统吗？9.2.6 考虑另一个系统：修改程序,以仿真上面的系统并确定该系统是否为时不变系统。

（选做）n = 0:40; D = 10;a = 3.0;b = -2;x = a*cos(2*pi*0.1*n) + b*cos(2*pi*0.4*n);[x1,n1]=sigmult(n,n,x,n)[x2,n2]=sigshift(x,n,1)[y,ny1]= sigadd(x1,n1,x2,n2)[y1,ny11]= sigshift(y,ny1,D)[sx,sn]= sigshift(x,n,D)[sx1,sn1]=sigmult(n,n, sx,sn)[sx2,sn2]=sigshift(sx,sn,1)[y2,ny2]= sigadd(sx1,sn1,sx2,sn2)D= sigadd(y1,ny11,y2,ny22)六、实验心得体会：实验时间批阅老师实验成绩实验二 FFT 实现数字滤波实验类别：提高性实验 1.实验目的(1) 通过这一实验，加深理解FFT 在实现数字滤波（或快速卷积）中的重要作用，更好的利用FFT 进行数字信号处理。