实验一：用MATLAB产生时域离散信号

实验1 利用matlab进行系统的时域分析

一． 实验目的:

1． 了解离散时间序列卷积与的matlab实现;

2． 利用卷积与求解系统的零状态响应;

二． 实验原理:

1． 连续时间系统零状态响应的求解

连续时间LTI系统以常系数微分方程描述,系统的零状态响应可通过求解初始状态为零的微分方程得到。在MATLAB中,控制系统工具箱提供了一个用于求解零初始状态微分方程数值解的函数lsim。其调用方式为

y= lsim( sys,x,t)

式中t表示计算系统响应的抽样点向量,x就是系统输入信号向量,sys就是连续时间LTI系统模型,用来表示微分方程、差分方程、状态方程。在求解微分方程时,微分方程的连续时间LTI系统模型sys要借助tf函数获得,其调用方式为

sys= tf(b,a)

式中b与a分别为微分方程右端与左端各项的系数向量。例如对3阶微分方程

+++=+++

可用

a=[ a3, a2, a1, a0];b=[b3 ,b2, b1,b0]; sys=tf( b,a)

获得连续时间LTI模型。注意微分方程中为零的系数一定要写入向量a与b中。

【例2-1】描述某力学系统中物体位移y(t)与外力f(t)的关系为

++y(t)=x(t)

物体质量m=l kg,弹簧的弹性系数ks= 100 N/m,物体与地面的摩擦系数fd=2 N·s/m,系统的初始储能为零,若外力x(t)就是振幅为10、周期为1的正弦信号,求物体的位移y(t)。

解:由已知条件,系统的输入信号为x(t)=10sin(2πt),系统的微分方程为

++100y(t)=x(t)

计算物体位移y(t)的MATLAB程序如下:

%program2_1微分方程求解

ts=0;te=5;dt=0、01;

sys=tf([1],[1 2 100]);

实验2离散系统的时域分析

实验目的：加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。

实验原理：离散系统

][nx][nyDiscrete-timesystme

其输入、输出关系可用以下差分方程描述：

MkkNkkknxpknyd00][][

输入信号分解为冲激信号，

mmnmxnx][][][。

记系统单位冲激响应

][][nhn，

则系统响应为如下的卷积计算式：

mmnhmxnhnxny][][][][][

当Nkdk,...2,1,0时，h[n]是有限长度的（n：[0，M]），称系统为FIR系统；反之，称系统为IIR系统。

在MATLAB中，可以用函数y=filter(p,d,x)实现差分方程的仿真，也可以用函数 y=conv(x,h)计算卷积，用y=impz(p,d,N)求系统的冲激响应。

实验内容和要求：

1、以下程序中分别使用conv和filter函数计算h和x的卷积y和y1，运行程序，并分析y和y1是否有差别，为什么要使用x[n]补零后的x1来产生y1；具体分析当h[n]有i个值，x[n]有j个值，使用filter完成卷积功能，需要如何补零？

% Program P2_7

clf;

h = [3 2 1 -2 1 0 -4 0 3]; % impulse response

x = [1 -2 3 -4 3 2 1]; % input sequence y = 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);

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

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

一、实验目的

本实验旨在通过MATLAB软件，对离散时间信号进行时域分析，包括信号的显示、基本运算（如加法、减法、乘法、反转等）、以及频域变换（如傅里叶变换）等，以加深对离散时间信号处理的基本概念和原理的理解。

二、实验原理

离散时间信号是在时间轴上离散分布的信号，其数学表示为离散时间函数。与连续时间信号不同，离散时间信号只能在特定的时间点取值。离散时间信号的时域分析是研究信号的基本属性，包括幅度、时间、频率等。通过时域分析，我们可以对信号进行各种基本运算和变换，以提取有用的信息。

三、实验步骤

1. 信号生成：首先，我们使用MATLAB生成两组简单的离散时间信号，一组为正弦波，另一组为方波。我们将这些信号存储在数组中，以便后续分析和显示。

2. 信号显示：利用MATLAB的绘图功能，将生成的信号在时域中显示出来。这样，我们可以直观地观察信号的基本属性，包括幅度和时间关系。

3. 基本运算：对生成的信号进行基本运算，包括加法、减法、乘法、反转等。将这些运算的结果存储在新的数组中，并绘制出运算后的信号波形。

4. 傅里叶变换：使用MATLAB的FFT（快速傅里叶变换）函数，将信号从时域变换到频域。我们可以得到信号的频谱，进而分析信号的频率属性。

5. 结果分析：对上述步骤得到的结果进行分析，包括比较基本运算前后的信号波形变化，以及傅里叶变换前后的频谱差异等。

四、实验结果 1. 信号显示：通过绘制图形，我们观察到正弦波和方波在时域中的波形特点。正弦波呈现周期性的波形，方波则呈现明显的阶跃特性。

2. 基本运算：通过对比基本运算前后的信号波形图，我们可以观察到信号经过加法、减法、乘法、反转等运算后，其波形发生相应的变化。例如，两个信号相加后，其幅度和时间与原信号不同。反转信号则使得波形在时间轴上反向。

3. 傅里叶变换：通过FFT变换，我们将时域中的正弦波和方波转换到频域。正弦波的频谱显示其频率为单一的直流分量，方波的频谱则显示其主要频率分量是直流分量和若干奇数倍的谐波分量。这表明离散时间信号的傅里叶变换可以揭示信号的频率属性。

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

离散时间信号是一种离散的信号形式，其具有离散的时间间隔。这种信号在数字信号处理中得到了广泛的应用。时域分析是分析信号的一种方法，它通常包括分析信号的幅度、相位、频率等参数，并从中获得信号的特征。在本实验中，我们将探讨离散时间信号的时域分析方法。

1.实验目的

• 了解离散时间信号的基本概念和性质。

• 熟悉MATLAB软件的使用，理解信号处理工具箱的使用方法。

2.实验原理

离散时间信号是一种在离散时间点上定义的数列。它通常用序列来表示，序列的元素是按照一定的时间间隔离散采样得到的。离散时间信号的采样频率通常表示为Fs，单位是赫兹。

离散时间信号可以写成如下的形式：

x(n) = [x(0),x(1),x(2),...,x(N-1)]

其中，n表示离散时间点的下标，N表示离散时间信号的长度。

• 幅度分析：指分析离散时间信号的振幅大小。

离散时间信号的幅度、相位、频率的分析通常使用傅里叶变换、离散傅里叶变换等变换方法来实现。

3.实验步骤

3.1 生成离散时间信号

使用MATLAB编写程序，生成一个离散时间信号。例如，我们可以生成一个正弦信号：

t = 0:0.01:1;

x = sin(2*pi*100*t);

其中，t表示时间向量，x表示正弦信号。

将信号进行离散化，得到离散时间信号：

其中，fs表示采样频率，n表示采样时间点，xn表示采样后的信号。 使用MATLAB的plot函数，绘制离散时间信号的时域图像。

figure(1);

plot(n, xn);

xlabel('Time');

ylabel('Amplitude');

其中，figure(1)表示创建一个新的窗口，用于显示图像。xlabel和ylabel用于设置图像的横轴和纵轴标签。