西华大学《现代控制理论》实验报告样本(02) (1).doc

西华大学实验报告

西华大学实验报告
第组
.
西华大学实验报告（计算机类）

开课学院及实验室：机械工程专业实验中心实验时间：2014年4月16日
学生姓名学号成绩学生所在学院机械工程与自动化学院年级/专业/班课程名称现代控制理论课程代码1202389实验项目名称控制系统的稳定性分析及仿真项目代码指导教师肖继学项目学分一、实验目的
了解利用Mathlab实现控制系统稳定性分析及仿真的方法。

二、内容与设计思想
（1）实验内容
分析系统的稳定性，其中，
，，，

（2）实验思想
可以通过稳定性概念或者线性定常系统稳定性的判据来分析。

三、使用环境
1、WindowsXP
2、Mathlab5.0或以上版本的Mathlab

四、核心代码及调试过程

A=[1,0,1,1;2,-3,0,0;1,0,-2,0;4,-1,-2,-4];
B=[0,;0;1;2];
C=[3,0,1,0];
D=0;
flag1=0;
flag2=0;
[z,p,k]=ss2zp(A,B,C,D,1);
sys=ss(A,B,C,D)
pzmap(sys)
z,p,k
n=length(A);
%李亚谱诺夫第一方法?
fori=1:n
ifreal(p(i))>0flag1=1;
end
end
flag1
ifflag1==1
disp('由李亚谱诺夫第一法判定，系统不稳定');
elsedisp('由李亚谱诺夫第二法判定，系统稳定')
end

%李亚谱诺夫第二方法?
Q=eye(4,4);%Q=I?
P=lyap(A,Q);%求解矩阵P?
fori=1:n
det(P(1:i,1:i));
if(det(P(1:i,1:i))<=0)
flag2=1;
end
end
flag2
ifflag2==1
disp('由李亚谱诺夫第二法判定，系统不稳定');
elsedisp('由李亚谱诺夫第二法判定，系统稳定')
end
五、总结

六、附录
a=
x1x2x3x4
x11011
x22-300
x310-20
x44-1-2-4

b=
u1
x10
x20
x31
x42

c=
x1x2x3x4
y13010

d=
u1
y10

Continuous-timemodel.

z=

-10.9247
-1.2169
-2.8585

p=

1.7913
-5.0000
-2.0000
-2.7913

k=

1.0000

flag1=

1

由李亚谱诺夫第一法判定，系统不稳定

flag2=

1

由李亚谱诺夫第二法判定，系统不稳定