判断系统稳定性.docx

判断系统稳定性.docx

《判断系统稳定性.docx》由会员分享，可在线阅读，更多相关《判断系统稳定性.docx（21页珍藏版）》请在冰豆网上搜索。

判断系统稳定性

摘要

现今数字信号处理理论与应用已成为一门很重要的高新科学技术学科，通过功能强大的MATLAB软件与数字信号处理理论知识相互融合在一起，既使我们对数字信号处理的理论知识能够有更加深厚的解也提高了动手能力，实践并初步掌握了MATLAB的使用。

根据本次课题要求，通过使用MATLAB,方便了对系统函数的繁琐的计算，并且直观形象的用计算机进行模拟仿真，通过观察图，由图像的特征从而进一步的对系统进行形象的分析。

本课题中给出了系统函数，对其稳定性进行分析我们可以通i±MATLAB画零极图观察极点的分布，另外还可以通过MATLAB分析系统的单位阶跃响应、单位脉冲响应、幅频相频特性的图形更加具体的对系统进行分析。

关键字：

离散系统函数、MATLAB、零极点分布、系统稳定性。

一、设计原理

1•设计要求

（1）:

根据系统函数求出系统的零极点分布图并且判断系统的稳定性。

（2）:

求解系统的单位阶跃响应，并判断系统的稳定性。

（3）:

求系统的单位脉冲响应，并判斷系统的稳定性

（4）:

求出以系统频率响应，画出幅频特性和相频特性图（zp2tf,zplane,impz等）

2、系统稳定性、特性分析

进行系统分析时我主要利用MATLAB软件绘制出系统零极点的分布图、单位脉冲响应图.单位阶跃响应图等。

采用MATLAB软件进行设计时我调用了软件本身的一些函数来对课题进行绘图话口分析。

诸攻口zplane、impz、stepz、freqz等。

对系统函数的零极图而言：

极点在单位圆，则该系统稳定，极点在单位圆外，则该系统为非稳定系统。

当极点处于单位圆，系统的冲激响应曲线随着频率的增大而收敛；当极点处于单位圆上，系统的冲激响应曲线为等幅振荡；当极点处于单位圆外，系统的冲激响应曲线随着频率的增大而发散。

系统的单位阶跃响应若为有界的则系统为稳定系统。

由以上的判据配合图形对系统的稳定性进行分析，达到我们的课程要求。

系统函数H（z）的零极点分布完全决定了系统的特性，若某系统函数的零极点已知，则系统函数便可确定下来。

因此，系统函数的零极点分布对离散系统特性的分析具有非常重要意义。

通过对系统函数零极点的分析，可以分析离散系统以下几个方面的特性：

⑴系统单位样值响应h（n）的时域特性；

（2）离散系统的稳定性；

（3）离散系统的频率特性；

二、MATLAB绘图分析

MATLAB功能丰富，可扩展性强。

MATLAB软件包括基本部分和专业扩展两大部分的功能。

基本部分包括：

矩阵的运算和各种变换；代数和超越方程的求解；数据处理和傅立叶变换；数值部分等等，可以充分满足大学理工科本科的计算需要。

扩展部分称为工具箱。

它实际上是用MATLAB的基本语句辩称的各种子程序集，用于解决某一方而的专门问题，或实现某一类的新算法。

在使用MATLAB语言进行编程过程中，根据题目设计要求，需要用到得主要函数语言有clear,figure,impz,zpIane,freqz,stepz等。

clear清除变量和函数。

figure即创建图形窗口的命令。

impz绘制单位脉冲响应。

zpIane显示离散系统的零极点分布图。

freqz绘制幅频特性图。

stepz绘制单位阶跃响应图。

roots求多项式的根。

plot用于绘出函数图，plot（x,y）,其中x为横坐标，y为纵坐标，且x,y—般为一维的。

axis为人工选择坐标轴尺寸命令。

title（'加图形标題'）。

xlabel（'加X轴标记'）。

yIabeI（'加Y轴标记'）。

gridon力◎网格线。

subpIot（m,n,p）该命令将当前图形窗口分成mXn个绘图区。

（1）

%z）=

z2+5z-50

2z4-2.98z3+0.17z2+2.3418z-1.5147

1）zpIane函数求解系统函数零极点分布

对应分子多项式系数为

B=[1,5,-50];

对应分母多项式系数为

A=[2,-2.98,0.17,2.3418,-1.5147];

zpIane（B,A）;

gridon;

Iegend（'零点'，’极点'）;

titleC零极点分布图’）；

x=roots（A）;

y=roots（B）;

abs（x）;

零极点分布图为:

琴圾成0你隆I

零点rz=（-10,5）rp=（-0.9000,0.7000±0.6000iT0.9900）结合图形分析，极点都在圆，所以该系统稳定。

2）求输入为单位阶跃序列时系统的响应

B二[1,5,-50];

A=[2,一2.9&0.17,2.3418,-1.5147];

stepz（B,A,2000）;

gridon;

title（'单位阶跃系列输入的系统响应'）；单位阶跃响应图为:

单位阶跃条列输入的系统响应

o

■ssdEV

-500-1（X）0-1500-2000-Z500

020040060080010001200140016001800

nfsameles）

由图可见，该系统的单位阶跃响应曲线随着n增大最终归于有界。

因此，验证了该系统是一个稳定系统。

3）系统的单位脉冲响应

a二[1,5,-50];

b二[2,-2.9&0.17,2.341&T.5147];

impz（a,b,100）;

gridon;

title（'系统单位脉冲响应'）

单位脉冲响应图为：

承统W位脉沖H向血

n（samples）

由图可见，该系统的冲激响应曲线随着n增大而收敛。

因此，验证了该系

统是一个因果稳定系统。

4）系统相频幅频特性图及分析

cIear

a=[1,5,-50];

b=[2,-2.9&0.17,2.3418,-1.5147];

[H,w]=freqz（a,b,400,'whole'）;

Hf=abs（H）;

Hx=angIe（H）;

figured）;

pIot（w,Hf）;

titleC离散系统幅频特性曲线’）

figure

（2）plot（w,Hx）

titleC离散系统相频特性曲线’）

禹•敝癢统评：

场处KWI:

illi純

系统幅频特性图为:

3000

2500

2000

1500

500

1OOO

由图可见，该系统的幅频特性曲线是一个凹而我们可以推断该系统具有带阻滤波的特性。

系统相频特性图为:

由上我们可以看出系统在一个期间的相位变化。

z

（2）

z-0.5

1）zplane函数求解系统函数零极点分布

对应分子多项式系数为

B二⑴；

对应分母多项式系数为

A=[1,-0.5];

zpIane（B,A）;

gridon;

Iegend（*零点’J极点’）;ti11e（'零极点分布图’）;x=roots（A）;

y-roots（B）;

abs（x）;

零极点分布图为：

00.5

RealPart

2）求输入为单位阶跃序列时系统的响应

3=[1,-0.5];

b二⑴；

stepz（B,A,2000）;

gridon;

title（'单位阶跃系列输入的系统响应'）;

单位阶跃响应图为：

020040060080010001200140016001800

n（samples）

3）系统的单位脉冲响应a=[1,-0.5];

b二⑴；impz（a,b,100）;

gridon;

title（'系统单位脉冲响应'）

4）系统相频幅频特性图及分析

cIear

a=[1,-0.5];

b=[1];

[H,w]=freqz（a,b,400,'whoIe'）;

Hf=abs（H）;

Hx=angle（H）;

figured）;

pIot（w,Hf）;

titleC离散系统幅频特性曲线’）

figure

（2）plot（w,Hx）

titleC离散系统相频特性曲线’）系统幅频特性图为:

0.8

0.7

阁枚系统相频待性曲线

系统相频特性图为:

0.6

0.4

0.2

0-0.2-0.4-0.6

（3）仏⑵=亠

z-1

1）zplane函数求解系统函数零极点分布

对应分子多项式系数为

B二⑴；

对应分母多项式系数为

zpIane（B,A）;

gridon;

Iegend（'零点'，’极点'）;ti11e（'零极点分布图’）;x=roots（A）;

y二roots（B）;

abs（x）;

零极点分布图为：

•

•

O呼点

X曲点

$

■

1

•

9

•

•

■

•

・・•■

/

•

-1

-1

-0.500.5

RealPart

2）求输入为单位阶跃序列时系统的响应a二[1,T];

stepz（B,A,2000）;

gridon;

title（'单位阶跃系列输入的系统响应'）；单位阶跃响应图为:

1800

1600

1400

600

400

3）系统的单位脉冲响应

a=[1,-1];

impz（a,b,100）;

gridon;

title（'系统单位脉冲响应'）

单位脉冲响应图为：

做紆已中彳立月水”円佝应?

4）系统相频幅频特性图及分析cIear

a=[1,~1];

b=[1]

[H,w]=freqz（a,b,400,'whole'）;

Hf=abs（H）;

Hx=angle（H）;

cIffigured）;

pIot（w,Hf）;

titleC离散系统幅频特性曲线’）

figure

（2）plot（w,Hx）

titleC离散系统相频特性曲线’）

系统幅频特性图为：

离散条统幅频特性曲线

系统相频幅频特性图为:

2

1.5

1

0.5

0

-0.5

（4）

离散系统和倾特性曲线

比⑵=—

z-2

1）zplane函数求解系统函数零极点分布

对应分子多项式系数为

B二⑴；

对应分母多项式系数为

A二[1,-2];

zpIane（B,A）;

legend（*零点’J极点’）;

ti11e（'零极点分布图’）;

x=roots（A）;

y二roots（B）;

abs（x）;

零极点分布图为：

o妄点

x做总-

■

■

■

•••••••••••••••<

卜…・•・、

■

••

■

-

I

•

■

•

•

•

•

r

、-

-

•

・・・■•

•

■•

■

•

aA

-

琴极点分命閔

-1

-1-0.5

00.5

RealPart

2）求输入为单位阶跃序列时系统的响