判断系统稳定性.docx
《判断系统稳定性.docx》由会员分享，可在线阅读，更多相关《判断系统稳定性.docx（22页珍藏版）》请在冰豆网上搜索。
判断系统稳定性
摘要
现今数字信号处理理论与应用已成为一门很重要的高新科学技术学科，通过功能强大的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离散系统相频特性曲线’）
系统幅频特性图为:
3000
2500
2000
1500
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）;
零极点分布图为：
2）求输入为单位阶跃序列时系统的响应
3=[1,-0.5];
b二⑴；
stepz（B,A,2000）;
gridon;
title（'单位阶跃系列输入的系统响应'）;
单位阶跃响应图为：
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];
b二⑴；
stepz（B,A,2000）;
gridon;
title（'单位阶跃系列输入的系统响应'）；
单位阶跃响应图为:
卑位卩介跃系列输入的并统响应
2000
1800
1600
1400
600
400
200
3）系统的单位脉冲响应
a=[1,-1];
b二⑴；
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）;
gridon;
legend（*零点’J极点’）;
ti11e（'零极点分布图’）;
x=roots（A）;
y二roots（B）;
abs（x）;
零极点分布图为：
o妄点
x做总-
■
■
■
•••••••••••••••<
卜…・•・、
■
••
■
-
I
•
■
•
•
•
•
r
、-
-
•
・・・■•
•
■•
■
•
aA
-
琴极点分命閔
-1
-1-0.5
00.5
RealPart
2）求输入为单位阶跃序列时系统的响应a=[1,~2];
b二⑴；
stepz（B,A,2000）;
gridon;
title（'单位阶跃系列输入的系统响应'）;
单位阶跃响应图为：
n（samples）
3）系统的单位脉冲响应a=[1,-2];
b二⑴；
impz（a,b,100）;
gridon;
title（'系统单位脉冲响应'）
单位脉冲响应图为:
系统中位脉冲II向应
4）系统相頻幅频特性图及分析
cIear
a=[1,-2];
b=[1];
[H,w]=freqz（a,b,400,'whoIe'）;
Hf=abs（H）;
Hx=angle（H）;
elf
figured）;
pIot（w,Hf）;
titleC离散系统幅频特性曲线’）
figure
（2）plot（w,Hx）
titleC离散系统相频特性曲线’）
系统幅频特性图为：
离版系銃幅:
步页僭性川I銭
■
1
1
、
•
1
28
2.6
/
/
///
、
\、
x\
-
24
ff
\
-
2.2
{
V
1
L
\
2
1.8
-
1
/
/
J
\
\
\
\
\
1
■\
-
1.6
-
f
t
/
k
\
-
14
■//
、
-
1.2
/
-//
\
\\
\
-
Z/
1
1
Au
系统相频幅频特性图为:
离散•系统朴I顾待性仙级
三、总结
我通过画零极图分析极点分布，由极点都分布在单位圆得出该系统为稳定的。
又通过作出系统的单位脉冲响应得出收敛的特性进一步验证了系统本身的稳定特性。
系统的单位阶跃响应的有界性也说明了系统的稳定性。
我还画出了系统的幅频相频特性图分析出系统具有的带阻滤波特性还有相位的变化，基本完成了课程设计的要求。
四、心得体会
通过这次数字信号处理课程设计，使我对数字信号处理这门课程有了更深刻的认识，认识到了MATLAB软件的巨大功用并能初步掌握软件的使用方法，对以后的科学研究打下了一定的总础。
课程设计过程中虽然遇到了许多问题但通过咨询老师同学，最终还是自主完成课程题目，使我更加对系统的稳定性有了深的认知这会促进我以后课程的继续学习。
通过这次课设使自己能真正的处理现实当中的相关问题，而且对信号分析与处理的基本方法有了更深一层的理解，更重要的是提高了独立分析和解决实际问题的能力，这对以后进一步学习和实验提供了宝贵的经验
参考文献
1、Matlab7.0从入门到精通人民邮电保柱等2010
2、MATLAB在数字信号处理中的应用清华大学薛年喜2003年
3、自动控制原理国防工业强2004.06
4、信号处理原理及应用林洪彬平王娜机械工业2009
年
升级会员 