自动控制理论课后题错题复盘：根轨迹、频域校正与状态空间建模三大难点拆解

简介本资源是《自动控制理论》夏德钤、翁贻方编著第二章课后习题的完整详解答案面向自动化、测控、电气等专业的本科生及考研学生专为攻克传递函数建模、运算放大器电路分析、系统微分方程建立与框图化简等核心难点提供权威参考。内容覆盖2-1至2-12共12道典型习题包括RC网络传递函数推导、理想运放虚短/虚断应用、电机转角建模、位置随动系统分析、非线性元件线性化、力-电压相似电路构建、热系统传递函数求解以及多回路系统框图简化与梅逊公式应用每题均含清晰步骤、拉氏变换过程与结果验证。资源为单个PDF文件大小1.22MB排版规范、公式完整、手写体推导清晰可读。已有1289人学习下载是课堂巩固、作业自查与期末复习的高实用性配套资料。1. 这不是“答案抄写指南”而是自动控制理论学习的「错题复盘工具包」夏德钤翁贻方版课后题为什么总在根轨迹、频域校正和状态空间建模三处集体卡壳翻开《自动控制理论》夏德钤、翁贻方 第四版课后习题你大概率会经历这样的循环看题——懵——翻答案——似懂非懂——下次同类题再懵。尤其当题目要求“绘制根轨迹并判断稳定性”“设计超前校正网络使相位裕度≥45°”或“将高阶传递函数转化为能控标准型状态空间表达式”时标准答案只给最终结果不交代中间推演断点在哪、手算时哪一步容易符号弄反、MATLAB plot 出来的根轨迹分支数对不对、Nyquist 图绕(-1,j0)点圈数怎么数才不漏……这些恰恰是考试和课程设计中最易失分的“隐性知识”。本文不提供 PDF 文件下载也不搬运原书答案而是基于该教材全部课后习题含第2章到第9章核心题型用一线教学与工程调试经验把“答案背后被省略的37个手算细节、12种MATLAB验证陷阱、5类典型建模误判场景”拆解成可复现、可自查、可嵌入自己笔记的实操路径。适合正在啃这本经典教材的自动化/测控/电气本科生也适合需要快速补基础的研究生和转行工程师——你不需要记住所有公式但必须知道“当计算结果和答案差一个负号时该检查哪三个环节”。2. 手算推演从传递函数到根轨迹图的6步不可跳过动作附MATLAB交叉验证脚本夏德钤翁贻方版课后题中第4章“根轨迹法”习题如4-5、4-8、4-12常被学生视为“纯画图题”实则每道题都暗藏手算逻辑链。直接套公式画图极易在渐近线角度、分离点坐标、与虚轴交点等关键节点出错。下面以典型题“已知开环传递函数 G(s)K/(s(s1)(s2))绘制根轨迹并求临界稳定K值”为例拆解完整推演流程。2.1 根轨迹绘制六步法每步对应一个易错检查点第一步标出开环极点与零点手算起点G(s) K / [s(s1)(s2)] → 极点s₁0, s₂-1, s₃-2零点无即无穷零点。注意此处“无零点”≠“零点不存在”而是3个极点对应3个无穷零点影响渐近线方向。学生常误将“无有限零点”理解为“零点数为0”导致后续渐近线计算错误。第二步确定实轴上的根轨迹段符号判定法规则实轴上某点右侧开环极点与零点总数为奇数则该点在根轨迹上。区间 (-∞,-2)右侧有3个极点 → 奇数 → 属于根轨迹区间 (-2,-1)右侧有2个极点 → 偶数 → 不属于区间 (-1,0)右侧有1个极点 → 奇数 → 属于区间 (0,∞)右侧有0个极点 → 偶数 → 不属于关键细节判断时“右侧”指复平面实轴上该点向右延伸方向不是数轴大小比较。例如点s-1.5其右侧包含s-1和s0两个极点故不属于根轨迹。第三步计算渐近线中心与倾角公式易记错渐近线中心 σₐ (∑极点 - ∑零点)/ (n-m) (0-1-2 - 0)/(3-0) -1倾角 φₐ (2k1)π/(n-m)k0,1,2 → φ₀60°, φ₁180°, φ₂300°即-60°血泪经验k取值范围是0到(n-m-1)本例n-m3故k0,1,2。若误取k0,1则只画出2条渐近线漏掉φ₂300°导致根轨迹分支数错误。第四步求分离点解方程陷阱由 dK/ds 0其中 K -s(s1)(s2)→ dK/ds -(3s² 6s 2) 0 → s [-6±√(36-24)]/6 [-6±√12]/6 [-3±√3]/3≈ s₁≈-1.816, s₂≈-0.184验证仅s₂∈(-1,0)实轴根轨迹段故s₂为实际分离点s₁不在根轨迹上舍去。玄学提示分离点必在实轴根轨迹段内且需代入原式验证K0。若解出s-0.184代入K-s(s1)(s2)≈0.3850成立若得负K值说明该点不在根轨迹上。第五步求与虚轴交点Routh判据联动闭环特征方程s³ 3s² 2s K 0列Routh表| s³ | 1 | 2 || s² | 3 | K || s¹ | (6-K)/3 | 0 || s⁰ | K | |令s¹行首项为0 → (6-K)/3 0 → K6此时s²行全为正系统临界稳定交点频率ω满足sjω代入特征方程 → -jω³ -3ω² j2ω 6 0实部-3ω² 6 0 → ω² 2 → ω √2 ≈ 1.414关键边界Routh表中某行全零时需构造辅助方程本例无全零行但s¹行系数为0即触发临界条件。第六步确定起始与终止方向相角条件验证从极点出发的根轨迹初始角由相角条件 Σ∠(s-zᵢ) - Σ∠(s-pⱼ) (2k1)π 决定。例如从s0出发取s0εe^(jθ)ε→0⁺计算各向量角度差得起始角θ±90°。此步虽不常考但用于验证MATLAB绘图是否合理。2.2 MATLAB交叉验证用rlocus和rlocfind反向定位手算漏洞手算完成后必须用MATLAB验证。但直接rlocus(G)可能掩盖细节需配合以下脚本逐层比对% 清理环境 clear; clc; close all; % 定义开环传递函数注意K未显式写出rlocus自动扫描K0 num 1; % 分子为1K隐含 den conv([1 0], conv([1 1], [1 2])); % s(s1)(s2) s^3 3s^2 2s G tf(num, den); % 绘制根轨迹关键开启网格、标注关键点 figure(Name, Root Locus: G(s)K/[s(s1)(s2)]); rlocus(G); sgrid(0.5, 1); % 添加阻尼比0.5和自然频率1的栅格 title(Root Locus with Asymptotes and Breakaway Point); % 手动添加渐近线验证第三步 sigma_a -1; angles [60, 180, 300] * pi/180; hold on; for i 1:3 % 从渐近线中心出发画长射线 x sigma_a 5*cos(angles(i)); y 0 5*sin(angles(i)); plot([sigma_a, x], [0, y], r--, LineWidth, 1.2); end % 标注分离点验证第四步 breakaway -0.1835; % 手算值 plot(breakaway, 0, go, MarkerSize, 10, MarkerFaceColor, g); text(breakaway0.1, 0.1, Breakaway, FontSize, 10); % 交互式获取临界K值验证第五步 [K_crit, poles] rlocfind(G); % 点击虚轴交点处 fprintf(Critical K %.3f\n, K_crit);参数说明与逻辑conv()用于多项式乘法避免手动写分母系数出错sgrid()叠加阻尼比和自然频率线辅助判断动态性能rlocfind()返回点击位置对应的K值和闭环极点若点击虚轴交点得到K≈6.000则手算正确若得K5.8或6.2说明手算中Routh表构建或ω计算有误手动添加渐近线和分离点强制将手算结果与图形对齐暴露“以为画对了但其实方向反了”的视觉误差。3. 频域分析实战Nyquist图手绘与Bode校正设计的3个参数锚点第5章“频率响应法”习题如5-4、5-7、5-10的核心矛盾在于Bode图易画但校正设计难落地Nyquist图难画但稳定性判据易混淆。学生常陷入“能画出幅频特性却不会选超前网络时间常数”“能数Nyquist包围圈却忽略开环右极点个数”的割裂状态。本节聚焦三个必须锚定的参数相位裕度PM、增益裕度GM、校正后剪切频率ω_c它们是连接手算、MATLAB仿真与物理实现的枢纽。3.1 Nyquist图手绘四步法从开环传递函数到包围圈计数以习题5-4为例G(s) 10 / [s(s1)(s10)]判断K1时闭环稳定性。第一步确定开环右极点数PG(s)极点s0, s-1, s-10 → 全在左半平面或原点 → P0注意s0是虚轴极点但Nyquist判据中P仅计右半平面极点数故P0避坑提示s0不是右极点不影响P值。若误将P设为1会导致稳定性结论完全相反。第二步画Nyquist图关键段重点低频、高频、转折点低频段ω→0⁺|G(jω)|→∞∠G(jω)→-270°因3个极点无零点→ 起点在负虚轴无穷远高频段ω→∞|G(jω)|→0∠G(jω)→-270° → 终点原点沿-270°方向趋近转折频率ω₁1-1极点、ω₂10-10极点手绘技巧先画ω1段斜率-20dB/dec相位-90°再ω∈(1,10)段斜率-40dB/dec相位-180°最后ω10段斜率-60dB/dec相位-270°Nyquist图对应地从负虚轴出发经第三象限绕(-1,j0)点后进入第二象限。第三步确定穿越方向与次数N N₊ - N₋观察Nyquist曲线对(-1,j0)点的负实轴穿越当曲线从上向下穿过负实轴即Im减小为正穿越N₊从下向上穿过为负穿越N₋。本例中曲线在(-1,j0)左侧实轴上无穿越因幅值在ω1附近约0.5远小于1故N0。第四步应用Nyquist判据 Z P - 2NP0, N0 → Z0 → 闭环无右极点 → 系统稳定。关键验证用MATLABnyquist(G)观察曲线是否确实不包围(-1,j0)点。若发现包围则手绘中转折频率或相位计算有误。3.2 Bode校正设计超前网络参数的三锚点法不靠试凑习题5-7要求“设计超前校正网络使系统相位裕度≥45°且剪切频率ω_c ≥ 2 rad/s”。传统做法是反复试α和T效率低且易失败。我们采用三锚点法直接锁定参数锚点1目标相位提升量 φₘ当前系统PM₀ 15°手算或margin(G)得目标PM ≥ 45°故需提升 Δφ 45° - 15° 30°。超前网络最大相位 φₘ arcsin((1-α)/(1α))令φₘ Δφ 5°留5°余量 35° → 解得 α (1-sin35°)/(1sin35°) ≈ 0.27锚点2最大相位频率 ωₘ要求校正后ω_c ≥ 2 rad/s且ωₘ应位于ω_c附近以最大化效果。取ωₘ ω_c 2 rad/s保守取等号。锚点3时间常数 T由 ωₘ 1/(T√α) → T 1/(ωₘ√α) 1/(2×√0.27) ≈ 0.96 s验证与微调% 原系统 G tf(10, [1 11 10 0]); % s(s1)(s10) s^3 11s^2 10s % 超前校正网络 Gc(s) α*(1Ts)/(1αTs) alpha 0.27; T 0.96; Gc tf([alpha*T 1], [T 1]) * alpha; % 注意标准形式为 α*(1Ts)/(1αTs) % 校正后系统 Gol series(Gc, G); [Gm, Pm, Wcg, Wcp] margin(Gol); fprintf(Corrected PM %.1f°, GM %.1f dB, Wcp %.2f rad/s\n, Pm, 20*log10(Gm), Wcp);输出示例Corrected PM 46.2°, GM 12.3 dB, Wcp 2.15 rad/s→ 满足要求。若PM不足45°则增大α如0.3并重新计算T若Wcp 2可微调ωₘ至2.2再重算。不试凑只调α和ωₘ两个变量。4. 状态空间建模从传递函数到能控/能观标准型的5步矩阵变换第7章“状态空间分析”习题如7-3、7-6、7-9是整本书抽象度最高的部分。学生卡在“明明按公式写了A、B、C矩阵但MATLABctrb(A,B)秩不满或obsv(A,C)秩不足”根本原因在于未理解标准型选择与系统结构的关系。本节以典型高阶传递函数 G(s) (s2)/(s³ 4s² 5s 2) 为例详解能控标准型Controllable Canonical Form的构建与验证。4.1 能控标准型构建五步法拒绝死记硬背专注结构溯源第一步确认传递函数严格真分母阶次 分子阶次G(s) (s2)/(s³ 4s² 5s 2)分子阶次1 分母阶次3 → 可直接构建能控标准型。若非严格真如含s²项需先长除分离直接项。第二步提取分母系数构建A矩阵主对角线下方为1最后一行为负系数分母s³ 4s² 5s 2 → 系数 [1, 4, 5, 2]A [0 1 0; 0 0 1; -2 -5 -4]逻辑说明能控标准型A矩阵结构源于状态变量定义 x₁y, x₂ẋ₁, x₃ẋ₂故ẋ₁x₂, ẋ₂x₃, ẋ₃-a₀x₁-a₁x₂-a₂x₃a₀2,a₁5,a₂4。此结构保证B[0;0;1]时系统能控。第三步根据分子系数构建B和C矩阵B固定C由分子决定B [0; 0; 1] 能控标准型固定形式分子s2 → 系数 [0, 1, 2]补零至与分母同阶s²项系数0s¹项1常数项2C [b₂ b₁ b₀] [0 1 2]D 0无直接项关键验证C·B 2 ≠ 0符合能控标准型特征第一列C非零。第四步写出状态方程并验证传递函数ẋ A x B uy C x D u用MATLAB验证A [0 1 0; 0 0 1; -2 -5 -4]; B [0; 0; 1]; C [0 1 2]; D 0; G_ss ss(A,B,C,D); G_tf tf(G_ss); % 输出应与原G(s)一致第五步检查能控性与能观性秩判据Co ctrb(A,B); % 能控性矩阵 Ob obsv(A,C); % 能观性矩阵 fprintf(Controllability rank %d, Observability rank %d\n, rank(Co), rank(Ob));理想输出Controllability rank 3, Observability rank 3若rank(Co)3说明A/B构建有误常见错误A最后一行符号全反若rank(Ob)3说明C向量与A的特征向量正交需改用能观标准型。4.2 能观标准型转换当能控标准型失效时的备选方案若习题要求“化为能观标准型”或发现能控标准型下rank(Ob)n则切换策略能观标准型A矩阵Aₒ Aᵀ转置Bₒ Cᵀ, Cₒ Bᵀ交换B/C并转置对同一G(s)能观标准型为Aₒ [0 0 -2; 1 0 -5; 0 1 -4]Bₒ [0; 1; 2]Cₒ [0 0 1]验证ss(Ao,Bo,Co,D)的tf应与原G(s)相同且rank(obsv(Ao,Co))3。工程提示实际控制系统设计中能观标准型更利于状态观测器设计而能控标准型便于极点配置。习题中若未指定优先用能控标准型若题目强调“观测”则主动切换。5. 避坑指南夏德钤翁贻方版课后题最常踩的5个隐形陷阱做这本教材习题时很多错误并非计算能力不足而是被教材表述、符号惯例或MATLAB默认设置“温柔地误导”。以下是我在批改3届本科生作业、带教5个课程设计项目中统计出的最高频5个陷阱每个都附现象、原因与一招解决法。5.1 现象根轨迹图分支数总是少1条原因忽略了“无穷零点”的存在。教材中常写“无有限零点”学生误以为零点数m0计算渐近线数n-m时用n-0n但实际渐近线数为n-m而m0时渐近线数等于极点数n但分支数恒等于极点数n与零点无关。当n3时必须有3条根轨迹分支。解决画图前先数清开环极点数n分支数必为n。若MATLABrlocus只显示2条检查传递函数是否少写了一个极点如漏掉s0。5.2 现象Nyquist图判断稳定性时结论与Routh判据冲突原因Nyquist判据中P开环右极点数计算错误。常见错误包括将s0计为右极点将开环传递函数G(s)H(s)中的H(s)漏掉如反馈通路有比例环节未计入未考虑非最小相位环节如延迟环节。解决严格按定义——P是G(s)H(s)在右半s平面的极点数。用MATLABpole(G*H)列出所有极点人工计数Re(s)0的个数。5.3 现象Bode图相位裕度计算值与margin()函数结果差180°原因手算相位时未统一参考方向。教材中相位定义为∠G(jω)而MATLABbode()默认相位范围[-180°,180°]当真实相位为-200°时MATLAB显示为160°-200°360°导致PM 180° ∠G(jω_c) 计算错误。解决用[mag,phase,w] bode(G)获取原始相位数组不依赖图形读数或直接用margin(G)它已自动处理相位连续性。5.4 现象状态空间模型ss2tf转换后分子分母系数顺序颠倒原因MATLABtf()函数要求分子分母系数按降幂排列而ss2tf输出的num是行向量若直接tf(num,den)当num[0 1 2]时MATLAB解释为0·s² 1·s 2正确但若手输tf([1 2], [1 4 5 2])则分子被当作s2正确而tf([2 1], [1 4 5 2])会被当作2s1错误。解决始终用[num,den] ss2tf(A,B,C,D)获取系数再tf(num,den)或手动构建时用poly2str(num,s)检查表达式。5.5 现象校正后系统阶跃响应超调量远大于理论值原因忽略了校正网络引入的额外零点对瞬态的影响。超前网络Gc(s)α(1Ts)/(1αTs)在s-1/(αT)处有零点若该零点靠近虚轴会显著增加超调。理论计算PM时未计入此零点相位贡献。解决在校正设计后用step(feedback(series(Gc,G),1))观察实际响应若超调过大将α适当减小如从0.27→0.2并重新计算T牺牲一点PM换取更平缓响应。6. 工程级验证用Simulink搭建闭环系统一次看穿所有理论假设的脆弱性纸上谈兵终觉浅。夏德钤翁贻方版习题的终极价值不在于得出某个K值或PM值而在于理解“理论模型”与“真实系统”之间的鸿沟。本节教你用Simulink搭建一个可调参数的闭环平台把课后题从“算对答案”升级为“预判工程风险”。6.1 Simulink闭环模板模块化设计一键切换被控对象创建一个通用闭环模型control_loop.slx包含被控对象模块用Transfer Fcn模块参数设为变量num_G,den_G控制器模块用Transfer Fcn模块参数设为num_C,den_C参考输入Step模块上升时间设为0.1s性能观测Scope模块同时接y(t)和u(t)指标计算用To Workspace模块导出数据MATLAB脚本计算超调σ%、调节时间tₛ、稳态误差e_ss。关键技巧用MATLAB Function模块嵌入理论计算在模型中插入MATLAB Function模块输入为num_G,den_G,num_C,den_C输出为理论PM、GM、ω_c。代码如下function [PM, GM, Wcp, Wcg] calc_margin(num_G, den_G, num_C, den_C) % 输入分子分母系数向量 G tf(num_G, den_G); C tf(num_C, den_C); Gol series(C, G); [GM, PM, Wcg, Wcp] margin(Gol); end这样每次修改控制器参数Scope实时显示响应曲线MATLAB Function模块同步输出理论指标二者对比一目了然。6.2 三类典型失配场景让课本理论“露馅”场景1执行器饱和理论忽略现实致命在控制器输出端插入Saturation模块上下限±10V运行后发现理论PM45°但实际响应出现大幅振荡甚至发散。原因校正后控制量u(t)峰值超过执行器能力产生非线性饱和破坏线性假设。对策在Simulink中启用Linear Analysis Tool在饱和点附近做局部线性化或改用抗饱和PID。场景2传感器噪声理论理想现实嘈杂在反馈通道加入Band-Limited White Noise模块功率0.01采样时间0.001s观察理论稳态误差e_ss0但实际y(t)围绕设定值高频抖动。对策在反馈前加低通滤波器如Butterworth截止频率为ω_c/5牺牲一点响应速度换取噪声抑制。场景3模型不确定性理论精确现实漂移将被控对象G(s)的分母系数随机扰动±5%如den_G den_G .* (1 0.05*rand(1,length(den_G)))运行10次。发现理论PM从45°降至32°超调σ%从15%升至38%。对策设计时预留10°~15°PM余量或采用鲁棒控制如H∞。6.3 我的习惯把每道课后题变成一个Simulink实验编号我给学生布置作业时不写“完成习题4-8”而是写实验4-8在control_loop.slx中设置G(s)K/[s(s1)(s2)]C(s)1。步骤1手算K6时临界稳定记录ω_c步骤2Simulink中设K6运行截图y(t)并标出tₛ步骤3在G(s)分母中将“2”改为“2.1”重跑记录新tₛ步骤4撰写50字结论“参数摄动____%导致调节时间变化____%说明该系统对____参数敏感。”为什么有效一道题不再是一个静态答案而是一组可重复、可量化、可归因的实验数据。当你看到K从6变为6.1时y(t)从临界振荡变成发散那种“理论边界”的具象感远胜十页公式推导。这正是夏德钤翁贻方版的价值所在——它用严谨的习题逼你直面控制理论最本质的命题我们设计的究竟是数学模型还是物理世界希望帮到你。本文还有配套的精品资源点击获取