第4章 习 题
A 习题(具有题解)
A 4-1 已知连续传递函数2
1()0.21
D s s s =
++,采样周期T =1s ,若分别采用
向前差分法和向后差分法将其离散化,试画出s 域和z 域对应极点的位置,并说明其稳定性。
解:
1 )S
域对应的极点为:1,21j 10
10
s =-
+
稳定
2) 向前差分法离散化:
2
12
2
2
2
1
()()|
11(0.22)0.21
(
)0.2(
)1
1
1.8 1.8
z s T
T
D z D s z z z T z T T T
T
z z -=
===--+-+-+++=
-+
z
域对应的极点为:1,2
9j 10
10
z =+
不稳定
3) 向后差分法离散化:
455
.0455.01
2.22.21
)1(
2.0)1(
1
|
)()(2
2
2
2
2
1+-=
+-=+-+-==-=
z z
z z z
z
Tz
z Tz
z s D z D Tz
z s
z
域对应的极点为:1,2
211
22
z j
=
+
稳定
(变换方法的基本练习,要求不使用MATLAB 的有关指令。) A 4-2
设连续传递函数1()0.051
D s s =
+,采样周期T =0.1s 。
(1)用突斯汀变换法求其脉冲传递函数D (z )。(2)用频率预修正突斯汀变换求其脉冲传递函数m ()D z 。(3)在转折频率20rad/s
ω
=处,分别计算()D s 、
()D z 、m ()D z 的幅值与相位,并比较之。
解:
T =0.1s 用Tustin 变换,120
1
120(1)0.5(1)
()0.051
20(1)20(1)
z s z z z D z s z z z
-=+++=
=
=
+-++
取特征频率为ω=20 rad/s 用预修正Tustin 变换得 Dm (z )=1
112.84
tan (/2)1
1
10.609(1)()0.051
(0.218)
z z s T z z z D z s z ωω--=
=+++=
=
++
在ω=20 rad/s 处D (s )的幅值与相位分别是 0.707,-45deg ;
D (z )的幅值与相位分别是
j j 1
j e
(e
)0.5(1)
0.5(1e
)0.5(1co s()sin ())T
T
T
z D z
T j T ωωωωω---==+=+=+-
幅值
j 20
(e
)
0.54T
D ωω===
相角 0
20
sin 2()
arc tan
57.28
1co s 2
j T G e ωω=∠=-=-+
D m (z )的幅值与相位分别是 0.707,-45deg ; A 4-3 设连续传递函数为2
1() 1.41
s D s s s +=
++,试用零极点匹配法使之离散化,
令1T
s =。
2
11
() 1.41
[(0.70.714)]
s s D s s s s j ++=
=
+++-±;012()(1)()()()
z z z D z k
z p z p -+=
--
1
00.368
sT
z e
e
-=== ;(0.70.714)
1,2
0.4960.7140.3750.325sT
j p e
e
j -±===∠±=±
2
(0.368)(1)()0.750.247
z z D z k
z z -+=-+;
1
()
1s D s ==;1
(10.368)(11)()1,0.397
10.750.247
z D z k k =-+===-+
所以,2
(0.3970.1444)(1)
()0.750.247
z z D z z z -+=-+
A 4-4
已知超前校正网络2()5
8
s D s s +=
+,采样周期T=0.1s ,试用突斯汀变换进
行离散化,求得其脉冲传递函数()T D z ,画出()D s 、()T D z 在0~3Hz 频段内的幅相频率特性,并比较之。
解:429
.0214
.3929.36
144555|
)()
(1
1
2
--=
--==+-=
z z z z s D z D z z T s T
1) 连续环节频率特性(见题图A4-4-1) 频率3Hz 对应于 2π36π18.85rad ?==/s
w=0:0.1:100;
[m,p]=bode(num,den,w);
subplot(211);plot(w,m),grid
subplot(212);plot(w,p),grid 2) 离散环节频率特性(见题图A4-4-2) [dm,dp]=dbode(dnum,dden,T,w); subplot(211);plot(w,dm),grid subplot(212);plot(w,dp),grid
题图A 4-4-1连续环节频率特性 题图A4-4-2离散环节频率特性 频率特性产生畸变,从离散环节频率特性中可以看见周期性。由于采样周期T =0.1较大,故使失真加大。但(0~3)Hz 低频部分类似。 A 4-5 已知连续陷波器传递函数为 1
2.011.0)
(2
2++++=s s
s s s G
1) 试用Tustin 变换方法将其离散,设采样周期T =1s ; 2) 原连续陷波器在11rad /s ω=
处频率特性幅值最小,试问Tustin 变换后,在
什么频率处幅值最小? 3) 为了使离散陷波器在s rad /11
=ω 处频率特性幅值最小,可采取什么办法。
解:1) Tustin 变换:6
.464.58.462.5)(2
1
1
2+-+-=
+-=z z
z z z G z z T s
2) 依据频率畸变公式可知,离散陷波器频率特性幅值最小的频率为 1D
2
a r c t a n 0.92732T
T ωω=
=rad/s
3) 通常可采取:
.减小采样周期T ,如取T =0.1s ,则离散陷波器频率特性幅值最小的频率为
1D 2arctan
0.9992
2
T T ωω=
=rad/s
.采用预修正Tustin 变换方法,取1
1rad /s ω=为关键频率则可保证在该频
率处离散陷波器频率特性与连续陷波器频率特性相等。 A 4-6 试用零极点匹配法求控制器()D s s a
=+的等效离散控制器,仅关注低频
段。
解:该控制器有一个零点a
s -=,没有有限极点,但其有一无限极点。依零极点
匹配法,零点a
s -=可以映射为e
aT
z -=,无限极点映射为1-=z ,因此等效离散
控制器为:
D e
()1
a T
z G z k
z --=+
根据稳态增益相等原则,可确定增益k D e
(1)(0)11
a T
z G k
G a
--===+
21e
a T
a k -=
-
故等效离散控制器为
D 2e
()1e
1
a T
a T
a z G z z ---=
-+
A 4-7 已知伺服系统被控对象的传递函数为2()(1)
G s s s =
+,串联校正装置为
0.06()0.35
0.004
s D s s +=+。采用某种合适的离散化方法,将D (s )离散为D (z ),并计算
采样周期T 分别为0.1s ,1s ,2s 时,计算机控制系统的单位阶跃响应,记录时域指标r s %,t t σ和。并说明连续域-离散化设计与采样周期T 的关系。
解:选用Tustin 变换,
T =0.1s 时
120
1
0.060.3510.3489()0.35
0.004
0.9996
z s z s z D z s z -=++-==
+-
T =1s 时 0.35980.3388
()0.996z D z z -=- T =2s 时 0.36950.3277
()0.992
z D z z -=
-
利用Simulink 进行数学仿真,可得曲线如题图A4-5所示。
t/s
c(t)
题图A 4-5单位阶跃响应
T =0.1s 时,单位阶跃响应的超调量:20.04% 峰值时间:4.7s 调节时间:8.4s
T =1s 时,单位阶跃响应的超调量:39.65% 峰值时间:4.4s 调节时间:15s T =2s 时,单位阶跃响应的超调量:67.37% 峰值时间:4.4s 调节时间:31s A 4-8 试求增量式PID 控制器(理想微分)的脉冲传递函数,设0.1c
T
T =,
I c
0.5T T =,D
c c
0.125,T T T =为临界振荡周期。
解:
p 1
2
p ()[()(1)(2)]()[]()
u k k A e k B e k C e k U z k A B z
C z
E z --?=--+-?=-+
2
1
p
2
()()()
A z
B z C
U z z U z k E z z
--+=+
I D (1//)(10.1/0.50.125/0.1) 2.45A T T T T =++=++=
D (12/) 3.5B T T =+=
D
/1.25C T T == 2
p
22.45 3.5 1.25
()z z D z k z z
-+=-
A 4-9飞行模拟转台是现代飞机飞行控制系统在地面进行仿真实验的高精度实验设备。题图A 4-9(a)是我国自行研制的三轴电动模拟转台。转台分成三个框,分别围绕各自轴转动,每轴各用一套高精度伺服系统驱动。简化后其中某一轴的伺服系统结构图如题图4-9(b)所示。所设计的控制器连续传递函数为
1()300100D s s s ??
=++ ???
2
1003001s s s ++=
试选择合适的离散化方法将其离散化,求得D (z ),并比较两个控制器的时域及频域的误差。设采样周期T =0.0005s 。
(a)
(b)
题图A 4-9 模拟转台及伺服系统结构图
解:从控制器的结构明显看出为PID 控制,所以积分项可以采用Tustin 变换:,
111()0.00025
21
1
T z z D z z z ++=
=--
微分项可以采用向后差分法:
112()100100
200000
0.0005z z D z s z
z
--===
所以得到离散后的数字控制器为
11()3000.00025
200000
1
z z D z z z
+-=++-
(1) 两个控制器的时域误差比较: 控制器结构如题图A4-9-1所示。
题图A4-9-1控制器结构图
分别加入斜坡信号、正弦信号,得到两个控制器的时域和误差曲线分别如题图A4-9-2~题图A4-9-5所示。仿真时,连续系统采用欧拉法仿真,正弦信号频率为10rad/s 。
从题图A4-9-2~题图A4-9-5可以看出,连续控制器和离散控制器针对斜坡信号、余弦信号的时域输出很接近,其误差都很小。 (2)两个控制器的频域误差比较:
1()300100D s s s ??=++ ???
2
1003001s s s ++=
2
2
11200300400300200000
()3000.00025
200000
1
z z z z D z z z
z z
+--+=++=
--
题图A 4-9-2连续与离散控制器题图A 4-9-3连续与离散控制器对斜坡信号的输出响应(已近重合)对斜坡信号的误差响应
题图A 4-9-4连续与离散控制器题图A 4-9-5连续与离散控制器对正余弦信号的输出响应(已近重合)对正余弦信号的误差响应
MATLAB仿真程序如下:
num=[100,300,1];den=[1,0];
n1=[300];d1=[1];
n2=0.00025*[1,1];d2=[1,-1];
n3=200000*[1,-1];d3=[1,0];
[n12,d12]=parallel(n1,d1,n2,d2);
[dnum,dden]=parallel(n12,d12,n3,d3);
%连续控制器频域:
bode(num,den);grid,hold on
%离散控制器频域:
T=0.0005;
dbode(dnum,dden,T)
连续控制器频域和离散控制器频域比较如题图4-7-6所示。由于采样周期较小,所以连续控制器和离散控制器频率响应特性在(0100)/rad s ω= 范围内非常
一致。
M a g n i t u d e (d B )10
10
10
10
10
10
10
P h a s e (d e g )
Bode Diagram
Frequency (rad/sec)
题图A 4-9-6 连续控制器频域和离散控制器频域
A 4-10已知连续控制器的状态方程如下
()()()()()()
x
t A x t B u t y t C x t D u t =+=+ (A4-10-1)
其拉氏变换为
1
()()
()X s sI A B U s -=- (A4-10-2)
试用Tustin 变换法进行离散化,求其离散状态方程。 解:将Tustin 变换式代入式(A4-10-2) ,可得:
1
()(1)(1)(1)()
22
A T
B T X z z I z I z I
U z -??=--++??
??
1
()()(1)()
222
A T A T
B T z I I z I
U z -?
?
=--+
+???
?
若令1
()
2
A T F I =-
;2
()2
A T F I =+
;1
2
B T G =
,上式可写为
1
121()()(1)()X z zF F z IG U z -=-+
1
1
1
1211()()()
zI F F zI I F G U z ---=-+ (A4-10-3)
依1F 、2F 定义,可以推得
1
1
1
122F F F I
---=
用该式替换式(A4-10-3)中()
zI
I +的第2项,则
1
1
1
11
1211211()()(2)()
X z zI F F zI F F F F G U z -----=-+-
进一步整理,得
1
1
21
121111()[()
2]()
X z zI F F F G F G U z ----=-+ (A4-10-4)
对式(A4-10-1) 的输出方程作z 变换,并将式(A4-10-4) 代入,则有
1
1
21
121111()[()2]()()
Y z C zI F F F G C F G U z D U z ----=-++
1
12
1
121111[()
2]()()()
C zI F F F G U z C F G
D U z ----=-++
由该式可见,这相当于一个离散状态方程的输出方程的z 变换。 令
1
1
12()()
2
2
A T A T F F F I I --==-
+
22
112()2
A T G F G I
B T
--==-
(A4-10-5)
H C
=
1
1
11()()
2
2
A T
B T E
C F G
D D C I --=+=+-
则离散系统状态方程为
(1)()()()()()
d d d d x k F x k G u k y k H x k E u k +=+=+ (A4-10-6)
所以,原连续系统状态方程式通过双线性变换离散后,其离散状态方程为式(A4-10-6),其各项矩阵如式(A4-10-5) 所示。
B 习 题
B 4-1 使用不同方法对传递函数
a
s a s D +=
)(
进行离散化近似: .前向一阶差分法 .后向一阶差分法 .双线性变换法
.预修正双线性变换法(关键频率取a
=1
ω)
.零极点匹配法。 B 4-2超前连续网络传递函数为
2
14
)(++=s s s D
试采用下述方法进行离散化,令采样周期T =0.25s 。 .前向一阶差分法 .后向一阶差分法 .双线性变换法 .预修正双线性变换法
.零阶保持器(阶跃响应不变) 法 并计算和比较各离散网络及连续网络在c 1.6rad /s
ω=处的相位及幅值。
B 4-3 已知一个简单连续陷波器为
1
1.101
)(2
2
+++=
s s
s s D
1) 用前向一阶差分近似法离散,取采样周期T =1s ,并说明该离散陷波器的稳定性。
2) 用双线性变换法求取离散陷波器,并说明此时陷波频率为多少? 3) 如若保证离散陷波器的陷波频率不变,可采取什么办法。 B 4-4已知系统结构如题图B 4-4所示,图中()()a s D s s
+=;h ()G s 为ZOH 传递
函数;01()G s s
=
。设T =0.1s 。
1) 将控制器用双线性变换法离散,试确定使系统稳定的最大a 值。
2) 试将控制器用一阶向后差分变换法离散,试确定使系统稳定的最大a 值。
题图B 4-4系统结构图
B 4-5 若离散化采用
d ()(1)(1)
()d 2u t u k u k e k t
T
+
--==
近似时称为中心差分法,试导出中心差分法替换式。
B 4-6 巴特沃斯(Butterworth) 滤波器常常用来获得锐截止阻带和平直通带频率特性的滤波器。其特性由下述幅值平方方程表示
2
2c 1
()1(/)
n
G ωωω=
+
式中n 为滤波器阶次,c ω为截止频率。若n =4,依幅值平方方程,可以得到c ω=1时的s 平面巴特沃斯(Butterworth) 滤波器的传递函数为
1
()(0.38270.9239)(0.38270.9239)
G s s j s j =
+++-
1
(0.3827
0.9239)(
0.3827
0.9239)
s j s j ?
++-+
试用零极点匹配方法求其脉冲传递函数。
B 4-7 题图B 4-7为一连续控制系统,要求将其改造为数字控制系统,设采样周期T =0.1s 。要求用零极点匹配法设计一个合适的数字控制器。离散化时要求考虑零阶保持器的影响。应比较连续控制系统与数字控制系统的单位阶跃响应。
题图B 4-7 连续控制系统结构图
B 4-8 直流电机速度伺服系统如题图B 4-8所示,采用PI 控制,对力矩干扰进行测量实现完全补偿,按连续系统方法设计,采用双线性变换方法将控制器离散,求数字控制器输出表达式。
题图B 4-8 直流电机速度伺服系统示意图