《MATLAB及其在理工课程中的应用指南》课件第9章.pptx

第9章MATLAB在信号和

系统中的应用举例

9.1连续信号和系统

9.2离散信号和系统

9.3系统函数

9.4频谱及其几何意义

BACK

9.1连续信号和系统

【例9-1-1】列出单位脉冲、单位阶跃、复指数函数等连续信号的MATLAB表达式。

解：

◆建模

严格说来，MATLAB是不能表示连续信号的，因为它给出的是各个样本点的数据，只有当样本点取得很密时才可看成连续信号。所谓“密”,是相对于信号变化的快慢而言，形象地说，在相邻样本点之间的数据变化必须非常小才能看成“密”,其严格的数学定义此处不予讨论。

在编程中，先设定共同的时间坐标，然后分别列出生成三种信号的程序。

MATLAB程序

clear,t0=0;tf=5;dt=0.05;t1=1;t=[t0:dt:tf];

%(1)单位脉冲信号

%在t1(t0≤t1≤tf)处有一持续时间为dt,面积为1的脉冲信号，其余时间均为零

t=[t0:dt:tf];st=length(t);

n1=floor(t1-t0)/dt);%求t1对应的样本序号

x1=zeros(1,st);%把全部信号先初始化为零

x1(n1)=1/dt;%给出t1处幅度为1/dt的脉冲信号subplot(2,2,1),stairs(t,x1)

%绘图，注意为何用stairs而不用plot命令axis([0,5,0,1.1/dt])

%为了使脉冲顶部避开图框，改变图框坐标

%(2)单位阶跃信号

%信号从t0到tf,在t1(tO≤t1≤tf)前为0,到t1处有一跃变，以后为1

x2=[zeros(1,n1-1),ones(1,st-n1+1)];%产生阶跃信号

subplot(2,2,3),stairs(t,x2)%绘图

axis([0,5,0,1.1])

%为了使方波顶部避开图框，改变图框坐标

%(3)复数指数信号

u=-0.5;w=10;x3=exp((u+j*w)*t);

subplot(2,2,2),plot(t,real(x3))%绘图

subplot(2,2,4),plot(t,imag(x3))%绘图

◆程序运行结果

x1、x2、x3的波形如图9-1所示。注意：若要显示连续信号波形中的不连续点，用stairs命令；若要使波形光滑些，则用plot命令较好。复数指数信号x3可以分解为余弦和正弦信号，它们分别是复数信号的实部和虚部，右图中的两个衰减振荡信号就代表了这两个相位差90°的分量。

20

10

0

055

1

0.5

-1

5安电05

图9-1例9-1-1中x1、x2、x3对应的四种波形

0

0

【例9-1-2】编写求任意高阶连续常系数线性系统冲击响

应的程序。

解：这个问题在第4章4.3.5节介绍多项式函数库时已经打下基础，在第7章机械振动的例7-3-1又讨论过二阶常系数线性微分方程的解法，读者可以先看懂那些例题的解法，再看本题。任意阶次的连续线性系统可用下列线性常微分方程表述：

十…+an

=b₁

写成传递函数形式为

因此，其特性可以用系统传递函数的分子分母系数向量b和a来表示。向量b和a的长度分别为m+1和n+1,对于一切物理上可实现的系统，必有n≥m。系统的冲击响应等于传递函数的拉普拉斯反变换，问题归结为如何求出这个反变换。

其中p₁,P₂,….,pn是分母多项式的n个根，而r₁,r₂,….,r,则

是对应于这n个根的留数。一阶分式的反变换可以查表得到，所以可得到冲击响应的公式：

h(t)=rePit+r₂eP²¹+…+rePnt

如果分母系数多项式没有重根，则可以将两个多项式之比分

解成n个一阶部分分式之和，即

可见只要求出根p和留数r,线性方程的解就得到了。求根

是代数问题，当阶次高时，没有解析解。可喜的是MATLAB提供了用数值方法求根和留数的residue函数，它的调用方法如下：

[r,p]=residue(b,a)

只要给出系