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