控制系统的阶跃响应
一、试验目的
1. 观察学习控制导流的单位阶跃相应;
2. 记录单位阶跃响应;
3. 掌握时间响应的一般方法;
二、实验
要求
对教师党员的评价套管和固井爆破片与爆破装置仓库管理基本要求三甲医院都需要复审吗
1.阻尼比与阻尼振荡频率Un对系统的影响;
2.零初值、非零初值与系统模型的关系;
3.稳态值与系统模型的关系;
4.系统零点对阶跃响应的影响;
三、试验
1. 二阶系统
10
--------------
s^2 + 2 s + 10
(1) 写程序观察记录阶跃响应曲线
程序代码为:
num=[10];
den=[1,2,10];
printsys(num,den);
step(num,den);
Damp(den);
[y,x,t]=step(num,den)
(2) Damp(den)计算系统的闭环根、阻尼比、无阻尼振荡频率并记录;
Eigenvalue Damping Freq. (rad/s)
-1.00e+000 + 3.00e+000i 3.16e-001 3.16e+000
-1.00e+000 - 3.00e+000i 3.16e-001 3.16e+000
(3) [y,x,t]=step(num,den)记录实验的峰值,峰值时间、过渡时间并与理论值比
y =
0
0.0133
0.0509
0.1096
0.1856
0.2752
0.3749
0.4809
0.5899
0.6987
0.8045
0.9049
0.9977
1.0814
1.1546
1.2165
1.2666
1.3048
1.3312
1.3463
1.3509
1.3459
1.3324
1.3116
1.2847
1.2531
1.2180
1.1807
1.1425
1.1043
1.0672
1.0321
0.9996
0.9704
0.9448
0.9233
0.9059
0.8926
0.8835
0.8783
0.8769
0.8787
0.8836
0.8910
0.9005
0.9116
0.9240
0.9371
0.9505
0.9639
0.9769
0.9892
1.0005
1.0107
1.0197
1.0272
1.0332
1.0378
1.0410
1.0427
1.0432
1.0425
1.0408
1.0381
1.0348
1.0309
1.0265
1.0219
1.0172
1.0125
1.0079
1.0036
0.9997
0.9961
0.9930
0.9904
0.9883
0.9867
0.9856
0.9850
0.9848
0.9851
0.9857
0.9867
0.9878
0.9892
0.9908
0.9924
0.9940
0.9957
0.9973
0.9988
1.0002
1.0014
1.0025
1.0034
1.0041
1.0047
1.0051
1.0053
1.0053
1.0052
1.0050
1.0047
1.0042
1.0038
1.0032
1.0027
1.0021
1.0015
1.0009
1.0004
0.9999
0.9995
0.9991
峰值
峰值时间
过渡时间并与理论值比
2.修改参数分别实现
程序代码为:
n0=10;d0=[1;2;10];step(n0,d0); ξ=0.316
hold on;
n1=n0;d1=[1;6.32;10];step(n1,d1); ξ=1
n2=n0;d2=[1;612.64;10];step(n2,d2);
n0/n0 =
10
--------------
s^2 + 2 s + 10
n1/n1 =
10
-----------------
s^2 + 6.32 s + 10
n2/n12=
10
------------------
s^2 + 12.64 s + 10
3. 修改参数分别实现
Wn1= 0.5Wn0
Wn2= 0.5Wn1的响应;
4. 改变参数与原系统比较
C1(s) =
10
--------------
s^2 + 2 s + 10
C11(s) =
2 s + 10
--------------
s^2 + 2 s + 10
C12(s) =
s^2 + 0.5 s + 10
----------------
s^2 + 2 s + 10
C13(s)=
s^2 + 0.5 s
--------------
s^2 + 2 s + 10
C14(s) =
s
--------------
s^2 + 2 s + 10
程序代码为:
n0=10;d0=[1,2,10];
printsys(n0,d0);
step(n0,d0);
hold on;
n1=[2,10];d1=[1, 2,10];
printsys(n1,d1);
step(n1,d1);
hold on;
n2=[1,0.5,10];d2=[1,2,10];
printsys(n2,d2);
step(n2,d2);
hold on;
n3=[1,0.5,0];d3=[1,2,10];
printsys(n3,d3);
step(n3,d3);
hold on;
n4=[1,0];d4=[1,2,10];
printsys(n4,d4);
step(n4,d4);