计算机控制系统高金源版课后答案.docx

计算机控制系统高金源版课后答案

135
2020年4月19日
文档仅供参考，不当之处，请联系改正。

第1章 习 题
B 习 题
B1-1 举例说明2-3个你熟悉的计算机控制系统，并说明与常规连续模拟控制系统相比的优点。
B1-2 利用计算机及接口技术的知识，提出一个用同一台计算机控制多个被控参量的分时巡回控制方案。
B1-3 题图B1-3是一典型模拟式火炮位置控制系统的原理结构图。由雷达测出目标的高低角、方位角和斜距，信号经滤波后，由模拟式计算机计算出伺服系统高低角和方位角的控制指令，分别加到炮身的高低角和方位角伺服系统，使炮身跟踪指令信号。为了改进系统的动态和稳态特性，高低角和方位角伺服系统各自采用了有源串联校正网络和测速反馈校正，同时利用逻辑电路实现系统工作状态的控制(如偏差过大时可断开主反馈，实现最大速度控制，当偏差小于一定值后实现精确位置控制)。试将其改造为计算机控制系统，画出系统原理结构图。

题图B1-3典型模拟式火炮位置控制系统的原理结构图
B1-4水位高度控制系统如题图B.1-4所示。水箱水位高度指令由W1 电位计指令电压ur确定，水位实际高度h由浮子测量，并转换为电位计W2 的输出电压uh。用水量Q1 为系统干扰。当指令高度给定后，系统保持给定水位，如打开放水管路后，水位下降，系统将控制电机，打开进水阀门，向水箱供水，最终保持水箱水位为指令水位。试把该系统改造为计算机控制系统。画出原理示意图及系统结构图。

题图B1-4 水箱水位控制系统原理示意图
B1-5 题图B1-5为一机械手控制系统示意图。将其控制器改造为计算机实现，试画出系统示意图及控制系统结构图。

题图B1-5机械手控制系统示意图
B1-6题图B1-6为仓库大门自动控制系统示意图。试将其改造为计算机控制系统，画出系统示意图。

题图B1-6 仓库大门自动控制系统示意图
B1-7车床进给伺服系统示意图如题图B1-7所示。电动机经过齿轮减速机构带动丝杠转动，进而使工作台面实现直线运动。该系统为了改进系统性能，利用测速电机实现测速反馈。试将该系统改造为计算机控制系统，画出系统示意图。

题图B1-7车床进给伺服系统示意图
B1-8 现代飞机普遍采用数字式自动驾驶仪稳定飞机的俯仰角、滚转角和航向角。连续模拟式控制系统结构示意图如题图B1-8所示。图中所有传感器、舵机及指令信号均为连续模拟信号。试把该系统改造为计算机控制系统，画出系统结构图。

题图B1-8 飞机连续模拟式姿态角控制系统结构示意图

第2章 习 题
A 习题（具有题解）
A 2-1 下述信号被理想采样开关采样，采样周期为T，试写出采样信号的表示式。

解：
1) ；2) ；
3)
A 2-2 已知f(t) 的拉氏变换式F(s) ，试求采样信号的拉氏变换式F* (s)(写成闭合形式) 。

解：
1) 首先进行拉氏反变换，得；

因为 , ，(依等比级数公式)
类似，，，因此有

A 2-3 试分别画出及其采样信号的幅频曲线(设采样周期T＝0.1s)。
解：连续函数的频率特性函数为：。
连续幅频曲线能够用如下MATLAB程序绘图：
step=0.1;
Wmax=100;
w2=-Wmax;
y2=5*abs(1/(10+w2*i));
W=[w2]; Y=[y2];
for w=-Wmax:step:Wmax
y=5*abs(1/(10+w*i));
W=[W,w]; Y=[Y,y];
end
plot(W,Y); axis([-Wmax Wmax 0 0.6])
grid
结果如题图A 2-3-1所示。

题图A 2-3-1
该函数的采样信号幅频谱数学表示式为

显然，采用的项数越大，则计算得到的值越逼近于实际值。这里采用来进行计算。
采样幅频曲线能够用如下MATLAB程序绘图：
T=0.1; %采样周期
ws=2*pi/T; %采样频率
num=50; %每个采样周期的计算点数
step=ws/num; %计算步长
Wmax=150; %画图显示的频率范围
GW=4*Wmax; %计算的频率范围
g0=(1/T)*5*abs(1/(1+10*GW*i)); G00=[g0];
g0=(1/T)*5*abs(1/(10+(GW+ws)*i)); G11=[g0];
g0=(1/T)*5*abs(1/(10+(GW-ws)*i)); G12=[g0];
g0=(1/T)*5*abs(1/(10+(GW+2*ws)*i)); G21=[g0];
g0=(1/T)*5*abs(1/(10+(GW-2*ws)*i)); G22=[g0];
g0=(1/T)*5*abs(1/(10+(GW+3*ws)*i)); G31=[g0];
g0=(1/T)*5*abs(1/(10+(GW-3*ws)*i)); G32=[g0];
g0=(1/T)*5*abs(1/(10+(GW+4*ws)*i)); G41=[g0];
g0=(1/T)*5*abs(1/(10+(GW-4*ws)*i)); G42=[g0];
g0=(1/T)*5*abs(1/(10+(GW+5*ws)*i)); G51=[g0];
g0=(1/T)*5*abs(1/(10+(GW-5*ws)*i)); G52=[g0];
g0=(1/T)*5*abs(1/(10+(GW+6*ws)*i)); G61=[g0];
g0=(1/T)*5*abs(1/(10+(GW-6*ws)*i)); G62=[g0];
g0=(1/T)*5*abs(1/(10+(GW+7*ws)*i)); G71=[g0];
g0=(1/T)*5*abs(1/(10+(GW-7*ws)*i)); G72=[g0];
g0=(1/T)*5*abs(1/(10+(GW+8*ws)*i)); G81=[g0];
g0=(1/T)*5*abs(1/(10+(GW-8*ws)*i)); G82=[g0];
g0=(1/T)*5*abs(1/(10+(GW+9*ws)*i)); G91=[g0];
g0=(1/T)*5*abs(1/(10+(GW-9*ws)*i)); G92=[g0];
其余类似，最后可得，结果如题图A 2-3-2所示。

题图A 2-3-2
A 2-4 若数字计算机的输入信号为，试根据采样定理选择合理的采样周期T。设信号中的最高频率为定义为。
解： ；；
因此有

由此可得
依采样定理得： rad/s
A 2-5 已知信号x＝，试画出该信号的频谱曲线以及它经过采样器和理想滤波器以后的信号频谱。设采样器的采样频率分别为4w1，1.5w1，和w1 3种情况。解释本题结果。
解：的频谱为脉冲，如题图A 2-5-1所示。
当采样频率时，采样频谱如题图A 2-5-1所示。由于满足采样定理，经过理想滤波器后，能够不失真恢复原连续信号。(见题图A 2-5-2)
当采样频率时，采样频谱如题图A 2-5-1所示。由于不满足采样定理，采样频率发生折叠，当经过理想滤波器后，只保留了折叠后的低频信号，其频率为。(见题图A 2-5-2)
当采样频率时，采样频谱如题图A 2-5-1所示。由于不满足采样定理，采样频率发生折叠，折叠后的低频信号位于处，当经过理想滤波器后，只保留了折叠后的低频信号，其频率为，即直流信号。(见题图A 2-5-2)

题图A 2-5-1

题图A 2-5-2
A 2-6 已知信号x＝，经过采样频率的采样器以后。又由零阶保持器恢复成连续信号，试画出恢复以后信号的频域和时域曲线；当时，情况又如何?比较结果。
解：见题图A 2-6。

题图A 2-6
结果表明，当采样频率较低时，零阶保持器输出阶梯较大，高频分量较大。
A 2-7 已知信号试求各采样信号的x(kT)及y(kT)，并说明由此结果所得结论。
解： ；
；
；
。
；

结果表明，不满足采样定理，高频信号将变为低频信号。
A 2-8 试证明ZOH传递函数中的s＝0不是Gh (s)的极点，而中，只有一个单极点s＝0。
解：
可见只有一个s=0极点。
表明分母上实际不存在积分环节。
A 2-9 对一信号进行采样，信号频谱如题图A 2-9所示，其中感兴趣的频率范围为(0~ω1 ) ，已知干扰信号频率ωf =5ω1，试讨论采样周期及前置滤波器的选择。

题图 A 2-9
解：依采样定理要求，为使采样信号不失真，要求采样频率应满足；另外，对干扰频率来说，为使其不进入感兴趣的频率范围内，要求，因此，要求。因此有2种情况：
1) 如果，那么干扰信号并不会与数据信号相混叠，干扰可经过数字滤波器滤掉；
2) 采用抗混叠前置滤波器进行滤除，则采样频率取。如要求干扰信号在信号频率处衰减20 dB，那么一个n阶滤波器的最大衰减率为，因此为了到达在十倍频程内衰老20 dB，应取。
A 2-10 用z变换法求解下列差分方程。
（1），已知输入信号，初始条件。
（2），已知初始条件。
（3），已知初始条件。
求。
解：
(1) 对差分方程进行z变换，得
，因此，，
z反变换，得
(2) 对差分方程进行z变换，得
，，

；；
；。

z反变换，
A 2-11 已知以下离散系统的差分方程，求系统的脉冲传递函数。
（1）；
（2）且初始条件为零。
解：
(1) 对差分方程进行z变换，得

A 2-12 试列出题图A 2-12所示计算机控制系统的状态方程和输出方程。图中
。

题图A 2-12 题A 2-12系统框图

解：1）被控对象离散化：
=
依串行法写状态方程：

2) 控制器离散化

状态方程为

3) 闭环系统方程

1 0]
A 2-13试用表示题图A 2-13所列系统的输出，指出哪些系统能够写出输出对输入的脉冲传递函数，哪些不能写出。

题图A 2-13题A 2-13所示系统
解：
(1) 不能 ；(2) 能(输出加虚拟开关) ；
(3) 能(输出加虚拟开关) ；(4) 不能 ；
(5) 能 ；(6) 不能
A 2-14 试分别求如题图A 2-14所示的两个系统的阶跃响应采样序列，并比较其结果可得什么结论 (设T=1s) 。

题图A 2-14 系统方块图
解：(a)
；
经过长除法，得

(b) ；

经过长除法，得
比较可见，加入零阶保持器后，系统响应升起较慢，振荡性加强，稳定性差。
A 2-15热蒸汽加热系统如题图A 2-15（a）所示。进气阀门开度由线圈控制的铁心带动。水箱内水温由热电偶检测。系统方块图如题图A 2-15（b）所示。若，T=0.2s，试求闭环传递函数、单位阶跃响应和稳态值。

(b)
题图A 2-15 题A 2-15加热系统结构图
解：

A 2-16 题图A 2-16（a）是以太阳能作动力的“逗留者号”火星漫游车，由地球上发出的路径控制信号能对该装置实施摇控，控制系统结构如图（b）所示，其中为干扰（如岩石）信号。控制系统的主要任务就是保证漫游车对斜坡输入信号具有较好的动态跟踪性能，并对干扰信号具有较好的抑制能力。若令数字控制器和增益，试求输出对输入信号及干扰信号的输出表示式 (设T=0.1s) 。
解：

(a)
（b）
题图A 2-16 火星漫游车控制系统
A 2-17 气体成分控制系统如题图A 2-17（a）所示。其中阀门开度由线圈控制的铁心位移控制。培育室内二氧化碳含量由气体分析仪测定，气体分析仪是一个时滞环节。系统动态结构图如题图A 2-17（b）所示。若采样周期，试求闭环传递函数。令K=1，D(z)=1。

题图A 2-17 题A 2-17气体成分控制系统
解：
；
其中
因此，

若采样周期，则有
A 2-18 已知某经过零阶保持器采样的连续系统由下述差分方程描述，如若可能，试确定它所对应的连续时间系统。
1）
2)
解：
1) 该差分方程能够转换为下述离散状态方程

相对应的连续时间系统的状态方程为

上述两方程应有下述关系
；
由此可求得 ；
，因此，
2) 该差分方程能够转换为下述离散状态方程

类似，其连续系统应满足，但该式无解，故没有对应的连续系统存在。可见，并不是所有离散系统都能找到相对应的连续时间系统。能够证明，只有离散系统矩阵G没有在负实轴上特征根时才存在对应的连续时间系统。
A 2-19 已知数字控制器的脉冲传递函数为

试求其频率特性，并画出其幅相频率特性曲线。
解：

依该式即可画出当的幅相特性曲线。分析可见，当时，虚部趋于，而实部趋于常数。当时，虚部等于0，而实部等于。幅相特性曲线的大致形状如题图A 2-19所示。

题图A 2-19 题A 2-19幅相特性曲线
A 2-20 采样系统如题图A 2-20所示，输入信号为，试求采样系统输出c(t) ，式中是信号与采样时刻相角差。
解：输入信号采样后，得

题图A 2-20 采样系统频率特性的测试
脉冲序列函数的傅里叶级数展开，能够写成

其中

因此

因此，

进入计算机的信号包括基频信号和各次旁频信号，计算机输出也同样包括上述信号。但后续环节F(s) 一般是低通网络，由于频带限制，高频被滤除。
1) 如果测试频率较低，，能够认为输出信号即为基频信号：

2) 如果测试频率，依采样频谱分析可知， k=1旁频与基频相重叠，因此

由于频率较高，常被系统滤除，因此
输入为
考虑到此时，因此输出为

=
因为 ， 因此

频率特性为
可见此时c(t) 、与起始相角有关。题图A 2-20-1说明这种情况(图中设T=1s) 。
．当，即测试信号与采样开关同步时，，=0，如题图A 2-20-1 (a) 所示。
．当，即测试信号与采样开关不同步时，如题图A 2-20-1 (b) 所示()。
．时，即测试信号与采样开关相移时时，如题图A 2-20-1 (c) 所示。
上述结果表明，采样系统是一种特殊的时变系统，它的输出与采样时刻有关。

（a） (b)

(c)
题图A 2-20-1 不同起始相角时时域响应曲线
3）当，但非常接近，会产生另一种频率干涉现象
假定，此时
+高频部分(此部分被滤除)

这即为一种差拍现象，高频信号被一低频信号进行幅值调制。
若，测试频率，则

经过后续环节后，幅值有衰减，但形状不变。

B 习 题
B 2-1 若已知的采样信号拉氏变换，试问时，=？，并就所得结果进行说明。
B 2-2 若，试由此证明，均为的极点(m为正整数)，并说明的零点与零点的关系。
B 2-3若连续信号的频谱如题图B 2-3所示，若采样频率分别为
时，试画出采样信号的频谱。

题图B 2-3连续信号的频谱
B 2-4 若信号被理想采样开关采样，并经过零阶保持器，试画出零阶保持器输出信号的频谱。假定分别大于和小于奈奎斯特频率。
B 2-5 若加到采样-零阶保持器上，采样周期，试求
1) 该保持器在w=3rad/s处有一输出分量，试求它的幅值与相位；
2) 对w=15rad/s、w=27rad/s，重复上述计算。
B 2-6 已知采样周期T=0.5s，试问在系统截止频率wc=2rad/s处，零阶保持器所产生的相移为多少？若使零阶保持器所产生的相移为-5o ，试问应取多大的采样周期。
B 2-7 已知连续信号x(t)=sin(w1t ) ，ws =4w1 ，试画出题图B 2-7上A、B、C点的波形图。

题图B 2-7 采样-保持示意图
B 2-8 已知连续信号，采样频率，试说明该信号采样又经过零阶保持器后，恢复为一直流信号。
B 2-9用相机拍一个转轮的图片，转轮上标有标记(如题图B 2-9所示) ，转轮转动频率为，照相机快门开关频率为，试讨论
时，从相机所拍图像上看到的情况。

题图B 2-9 转轮示意图
B 2-10 已知一连续信号为被采样，其采样频率为，试求表示采样信号的频率（令）。
B 2-11 一阶保持器在数学仿真中常有应用，试推导一阶保持器的传递函数。
B 2-12已知采样系统的脉冲传递函数为

试证明
并用该式求取
的值。
B 2-13车床进给伺服系统如题图B 2-13（a）所示。电动机经过齿轮减速机构带动丝杠转动，进而使工作台面实现直线运动。该系统为了改进系统性能，利用测速电机实现测速反馈。试将该系统改造为计算机控制系统。连续系统的结构框图如题图B 2-13（b）所示。若，试求数字闭环系统传递函数。令T=0.1s，K1 =Kx =1，K2 =0.1，Km =40，a=2。

题图B 2-13 习题B 2-13车床进给伺服系统
B 2-14 已知连续传递函数，试求取G(z)=Z[]，并讨论其零点随采样周期的变化情况。
B 2-15已知连续传递函数，如采用零阶保持器时，试求取其脉冲传递函数，并确定当采样周期为多大时，其零点均在单位圆内。
B 2-16若开环传递函数为，试绘制连续系统奈奎斯特图及带零阶保持器和不带零阶保持器离散系统的奈奎斯特图，设采样周期T=0.2s。

B 2-17一般，直流电动机可用下述连续传递函数或状态空间模型描述

式中为电机转角，U为电机控制电压。若令，试确定
1) 经过零阶保持器采样时，系统的离散状态空间模型；
2) 脉冲传递函数；
3) 输入与输出的差分方程；
4) 脉冲传递函数极点与零点随采样周期变化的关系。
B 2-18试用级数展开法求题图B 2-18系统离散状态方程，并画出结构图。

题图B 2-18 系统结构图
B 2-19试推导下述连续系统相对应的具有零阶保持器的离散状态方程。(T=1s)
1)
2)
B 2-20 试证明题图B 2-20(a) 表示近似微分。试证明题图B 2-20(b) 表示为一种积分器(一般称之为无延迟数字积分器)，即
若假设。
试证明题图B 2-20(c) 表示为另一种积分器(一般称之为有延迟数字积分器)，即
若假设。

题图B 2-20 近似微分及两种数字积分器结构图

第3章 习 题
A 习 题（具有题解）
A3-1 平面上有3对极点，分别为
，试求在z平面上相应极点的位置。
解:
1）对
；
2)
3)
映射结果见题图 A3-1。

题图 A3-1 题A3-1映射结果
A3-2 已知平面上实轴平行线上点的位置（）如题图A3-2（a）和（b）所示，试分别画出映射到z平面上点的位置。

题图A3-2 题A3-2图
解：依据进行判断。
(1) 题图A3-2 -1(a) ：
Ai 各点均映射在z平面单位圆内正实轴上同一点。
Bi 各点均映射在z平面单位圆内正实轴上同一点，但更靠近z=1点。
Ci各点均映射在z平面单位圆外正实轴上同一点。
(2) 题图A3-2 -1(b) ：
Ai 各点均映射在z平面单位圆内负实轴上同一点。
Bi 各点均映射在z平面单位圆内负实轴上同一点，但更靠近z=-1点。
Ci各点均映射在z平面单位圆外负实轴上同一点。
(a)（b）
题图A3-2-1
A3-3 已知z平面上的点，试求其映射至平面上的位置，设采样周期。
解：因为，因此有，

，因此有

A3-4 已知平面上封闭曲线如题图A3-4所示（①→②→③→④→⑤→①），试画出映射至z 平面的封闭曲线。

题图A3-4 习题A3-4图
解：如题图A3-4-1所示。图形对横轴是对称的。

题图A3-4-1 习题A3-4图解答

A3-5 已知离散系统闭环特征方程分别为
（1）
（2）
（3），试判断其稳定性。
解：
(2)依2阶系统稳定条件，
系统稳定。
(3)3阶系统，依所提供的朱利稳定判据程序：
%
clear
a=[1 2 1.31 0.28]; % 生成特征多项式系数数组
n=length(a); % 求数组维数
b=a; c(1)=a(1);
for i=1:n-1
p=b(1:n-i+1); % 取n-1维不为0数组
if abs(p(1))<10^(-10) break
else an=p(n-i+1)/p(1); % 计算朱利判据第一行系数
end
pp=fliplr(p); % 翻转数组
b=p-pp*an; % 计算第二行及与第一行之差
c(i+1)=b(1); % 取第一个数
end
c % 给出朱利判据系数
运行该程序，结果为：
c = 1.0000 0.9216 0.3112 0.0141，所有参数大于0，系统稳定。
利用matlab求特征根可得:
P=[1 2 1.31 0.28] ;
rs=roots(P)'
rs = -0.8000 -0.7000 -0.5000
即系统的极点为：
其模值均小于1，系统稳定。
A3-6 已知系统的结构图如题图A3-6所示，其中，输入，试用稳态误差系数法求稳态误差，并分析误差系数与的关系。

题图A3-6
解：

可见加入该信号，稳态误差为1，且与采样周期无关。
A3-7 汽车行驶速度控制系统的结构图如题图A3-7所示。设，试判断干扰力矩为单位阶跃时所产生的稳态误差（依图直接判断）。若，求使系统稳定的值范围。若该系统为连续系统时，结果又如何。比较说明之。

题图A3-7 习题A3-7汽车行驶速度控制系统的结构图
解：1) 从图中可见，稳态时为对消Mf的干扰，综合点处误差

如折算到速度v，则
2)
=

因此，
3) 若为连续系统，由于闭环系统为2阶系统，故有。

A3-8 已知单位反馈离散系统开环传递函数为

试求使系统稳定k与T的关系式。
解：闭环系统特征方程为

依2阶系统稳定条件，有
，此条件成立。

由此解得 ，因此可得
进一步，由第3个条件，可得 ，即
A3-9 试确定题图A3-9所示系统使系统稳定的值范围，令采样周期趋于0，值又如何？若将该系统作为连续系统，结果又如何？对上述结果进行讨论。

题图A3-9 习题离散系统结构图

解：1)

2) 当T趋于0时，上式的极限值为
3) 若为连续系统，则特征方程为
为使系统稳定要求 (k+2) >0；故
结论：1) 离散系统稳定性比连续系统差，稳定增益范围小；
2) T趋于0时，系统并不等于连续系统，按采样系统计算k范围较小。
A3-10 给定系统如题图A3-10所示，设指令输入，扰动输入，，当，系统的稳态误差如何？

题图A3-10
解：；
1) 首先判稳定性
；系统稳定。
2) 系统为0型系统；因此输入信号稳态误差
3) 求干扰引起的输出
；
取D(z)=k=2，
稳态值即为稳态误差，
4) 总误差
A3-11 写出开环脉冲传递函数的脉冲响应表示式。
解：

；
;

-----------------------------------------------------------
--------------------------------------------------------------
A3-12 如题图A2-16所示的火星漫游车控制系统，若分别为及，试确定使系统稳定的值范围。
解：1) T=0.1s

依可得
可得
该式成立。
2) T=1s

依，可得
得k>0 ；，得k<11.42
最后可知，系统稳定要求。

A3-13 双关节机械臂如题图A3-13（a）所示。简化后系统结构图如题图A3-13（b）所示。若，试画出连续系统及采样周期及开环对数频率特性曲线。

题图A3-13 习题A3-13双关节机械臂结构图
解：

T=0.1 开环对数频率特性曲线如题图A3-13-1所示。
T=1 开环对数频率特性曲线如题图A3-13-2所示。
MATLAB的计算程序如下
num1=[0.004 0.0045];den1=[1 -1.723 0.74];
w=logspace(-2,2);
[m1,p1]=dbode(num1,den1,0.1,w);
semilogx(w,20*log10(m1))
grid
xlabel('Fraquency(rad/s)'),ylabel('Magnitude(dB)'
semilogx(w,p1)
>> grid
>> xlabel('Fraquency(rad/s)'),ylabel('phase(deg)')

题图A3-13-1开环对数频率特性曲线(T=0.1秒)

题图A3-13-2开环对数频率特性曲线(T=1秒)

B 习 题
B3-1 已知z平面复极点，试求相应s平面极点的阻尼比及无阻尼自然频率。
B3-2 试确定使开环传递函数为

单位负反馈闭环系统稳定的K值。
B3-3 试确定题图B3-3所示系统的稳定性和单位阶跃输入时的稳态输出值。其中
1 (比例控制器) ；
2 (积分控制器) 。

题图B3-3系统结构图
B3-4 试求B3-3所示系统在斜坡输入时的稳态误差。
B3-5 题图B3-5为水位高度控制系统略图。电机经过减速器控制N个阀门的开度，水箱高度为h(t) ，水箱底面积为A，进水量为(--电机转角)，出水量，因此，水箱中水位高度由下述方程描述

因此有
对该系统，根据已给参数，可知

直流电机的传递函数为

驱动电机的功率放大器系数ka =50；电位计的传递系数ks =1；减速比i=100。
1) 若=1，T=0.05，试求使系统稳定的最大阀门数N；
2) 如考虑A/D的转换误差为5%，试求系统保持水位高度的稳态误差。

题图B3-5 水箱控制系统原理示意图
B3-6 微机控制的直流电机速度控制系统如题图B3-6所示。其中V，，，p=100脉冲/周。设采样周期T=0.1s。试求使系统稳定的值以及时，系统单位阶跃响应特性及稳态值。

题图B3-6 直流电机速度控制系统示意图
B3-7 数字飞船控制系统如题图B3-7 所示。若采样周期T=0.264s，，，试推导系统开环及闭环传递函数，并求使系统稳定的临界值。

题图B3-7 数字飞船控制系统
B3-8 已知单位负反馈闭环系统传递函数为
T=1s
试求开环传递函数，并绘制Bode图，求相位、增益稳定裕度。
B3-9 试求题A.2-15中热蒸汽加热系统的相位、幅值稳定裕度及单位阶跃响应特性和稳态误差。令分别为1、10，采样周期T分别为0.2s及1s。
B3-10 若开环传递函数为，试绘制连续系统奈奎斯特图及带零阶保持器和不带零阶保持器离散系统的奈奎斯特图，设采样周期T=0.2s。
B3-11 试求题图B3-11(a) 及(b) 所示系统干扰所引起的稳态误差表示式，并说明为减少干扰所引起的稳态误差应如何选取系统参数。

题图B3.11 题B3-11系统结构图
B3-12 系统结构如题图B3-12(a)所示。其中控制算法为

其中为延迟时间。
1) 试求当分别为0及T时，使闭环系统稳定的控制器增益值的大小。
2) 如连续系统如题图B3-12(b)所示，试求当分别为0及T时，使闭环系统稳定的控制器增益值的大小，并与上述结果进行比较。

题图B3.12 题B3-12系统结构图

第4章 习 题
A 习题（具有题解）
A 4-1 已知连续传递函数，采样周期T=1s，若分别采用向前差分法和向后差分法将其离散化，试画出s域和z域对应极点的位置，并说明其稳定性。
解：
1 )S域对应的极点为： 稳定
2) 向前差分法离散化：

z域对应的极点为： 不稳定
3) 向后差分法离散化：

z域对应的极点为： 稳定
(变换方法的基本练习，要求不使用MATLAB的有关指令。)
A 4-2 设连续传递函数，采样周期T=0.1s。
（1）用突斯汀变换法求其脉冲传递函数D(z)。（2）用频率预修正突斯汀变换求其脉冲传递函数。（3）在转折频率处，分别计算、、的幅值与相位，并比较之。
解：
T=0.1s用Tustin变换，
取特征频率为=20 rad/s 用预修正Tustin变换得
Dm（z）=
在=20 rad/s处D（s）的幅值与相位分别是 0.707，-45deg；
D（z）的幅值与相位分别是

幅值
相角
Dm（z）的幅值与相位分别是 0.707，-45deg；
A 4-3 设连续传递函数为，试用零极点匹配法使之离散化，令。
;
；
；
；

因此，
A 4-4 已知超前校正网络，采样周期T=0.1s，试用突斯汀变换进行离散化，求得其脉冲传递函数，画出、在0~3Hz频段内的幅相频率特性，并比较之。
解：
1) 连续环节频率特性(见题图A4-4-1)
频率3Hz对应于 /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低频部分类似。
A4-5 已知连续陷波器传递函数为

1) 试用Tustin变换方法将其离散，设采样周期T=1s；
2) 原连续陷波器在 处频率特性幅值最小，试问Tustin变换后，在什么频率处幅值最小？
3) 为了使离散陷波器在 处频率特性幅值最小，可采取什么办法。
解：1) Tustin变换：
2) 依据频率畸变公式可知，离散陷波器频率特性幅值最小的频率为
rad/s
3) 一般可采取：
．减小采样周期T，如取T=0.1s，则离散陷波器频率特性幅值最小的频率为
rad/s
．采用预修正Tustin变换方法，取为关键频率则可保证在该频率处离散陷波器频率特性与连续陷波器频率特性相等。
A 4-6 试用零极点匹配法求控制器的等效离散控制器，仅关注低频段。
解：该控制器有一个零点，没有有限极点，但其有一无限极点。依零极点匹配法，零点能够映射为，无限极点映射为，因此等效离散控制器为：

根据稳态增益相等原则，可确定增益k

故等效离散控制器为

A 4-7 已知伺服系统被控对象的传递函数为，串联校正装置为。采用某种合适的离散化方法，将D(s)离散为D(z)，并计算采样周期T分别为0.1s，1s，2s时，计算机控制系统的单位阶跃响应，记录时域指标。并说明连续域-离散化设计与采样周期T的关系。
解：选用Tustin变换，
T=0.1s时
T=1s时
T=2s时
利用Simulink进行数学仿真，可得曲线如题图A4-5所示。

题图A 4-5单位阶跃响应
T=0.1s时，单位阶跃响应的超调量：20.04% 峰值时间：4.7s 调节时间：8.4s
T=1s时，单位阶跃响应的超调量：39.65% 峰值时间：4.4s 调节时间：15s
T=2s时，单位阶跃响应的超调量：67.37% 峰值时