MATLAB—信号与系统中的应用

用MATLAB处理信号与系统中的问题，更加直观、方便、准确。

一. 理论知识

1.线性系统的响应

许多实际系统都可以用线性系统来模拟，如果线性系统的参数是常数，则该线性系统可以用线性常系数微分方程来描述，所以，求解线性系统的响应就转换成求线性常系数微分方程的解了！

线性系统响应对应微分方程解
对于一个线性系统，存在激励u(t)和响应y(t)，激励和响应之间的关系可以用微分方程的形式来描述：

其中，a、b为常系数；方程为n阶线性时不变（LTI）常系数微分方程。

2.微分方程的解

既然将线性系统的响应转换为微分方程的解，那么求解微分方程就是最重要的任务了！

Ⅰ.经典解

首先看一下微分方程的经典解：完全解=齐次解+特解

齐次解仅仅取决于系统本身，特解仅仅取决于激励。

Ⅱ.完全响应

系统的响应，就是完全响应，可以分为零状态响应和零输入响应，如图：

实际上，我们分析一个系统的响应，主要是从零输入响应和零状态响应入手的。

3.零输入响应

在零输入条件下，系统的响应就是零输入响应。
零输入条件下，微分方程等式右边都为0，化为齐次方程，求出特征根，就可以求出零输入响应！
当特征根都为单根时，零输入响应如下：

4.零状态响应

在零状态条件下，系统的响应就是零状态响应。
零状态条件下，初始状态都为0，即y(0-)=0…等等，此时的微分方程为非齐次方程，求出特征根，就可以求出零状态响应！
当特征根都为单根时，零状态响应如下：

利用卷积求零状态响应

5.冲激响应

冲激响应是一种特殊的零状态响应，LTI系统处于零状态下，由单位冲激信号引起的响应就是冲激响应，一般用h(t)表示，冲激响应在后面求LTI系统响应的时候十分重要！！！
看一下冲激响应的框图表示：

6.阶跃响应

阶跃响应也是一种特殊的零状态响应，是输入为阶跃函数时系统的响应：

7.卷积求零状态响应

卷积的定义：

卷积的表示：

卷积的物理含义：

卷积求零状态响应：

任意信号都可以表示为冲激序列之和，所以，任意信号作用域系统的零状态响应都可以用冲激响应表示出来，这就是卷积求零状态响应：

二.连续信号的MATLAB描述

刚才已经对求解线性系统做了一个概括，那么各种函数是怎么样在MATLAB中描述的呢？接下来就一一介绍最常用的几种函数的MATLAB描述。
要知道，利用MATLAB分析连续信号时，要对信号作离散化处理，就是取出抽样点，根据抽样点的数据画出信号的图像，所以采样的区间和采样密度就影响着信号的表示效果。

1.单位冲激信号

MATLAB的实现代码如下：
在MATLAB中，我们以持续时间为dt，面积为1的脉冲信号表示冲激函数，具体注释都已经在函数中作了描述：

冲击函数
clear
clc
t0=0;%采样开始
dt=0.01;%采样间隔
tf=5;%采样终点
t1=1;%冲击点在t1处有一个持续时间为dt，面积为1的脉冲信号
t=[t0:dt:tf];%采样点
st=length(t);
n1=floor( (t1-t0)/dt );%floor：朝负无穷方向取整数，此处为t1的序列号
x1=zeros(1,st);%zero：生成（1，st）的零矩阵，此处为初始化信号
x1(n1)=1/dt;%给出t1处的信号
subplot(2,1,1),stairs(t,x1),grid on
axis([0,5,0,30])
subplot(2,1,2),plot(t,x1),grid on
axis([0,5,0,30])

使用stairs和plot只是为了对比一下画图效果。
运行结果如示：