简介这是一份面向机械工程及相关专业学生和初学者的《机械工程控制基础》课程PPT源自三峡大学机械与材料学院方子帆教授的课堂讲义聚焦控制理论的基本概念、系统工作原理与组成并通过恒温箱温度控制、钢铁轧制等案例讲解自动控制与反馈控制的核心逻辑。资源包仅含1个PPT文件整体大小约9.07MB内容紧凑适合用于课程预习、期末复习或教学备课。目前已有199人浏览学习。PPT以第一章绪论及基本概念为主线包含控制的定义、被控对象与控制器等术语解读、人工/自动恒温控制动态过程、系统原理方块图、反馈控制思想等页面直观重点突出能帮助读者快速建立控制系统知识框架理解传递函数、反馈等核心概念在实际工程中的作用也为后续学习现代控制理论打下基础。1. 机械工程控制基础PPT的实质不是排版是把控制过程讲顺做机械工程控制基础PPT的人通常是从课件正文开始动手的先找一版现成的讲稿把公式和波形图贴上去再调整配色。但真正的难点不在排版而在内容编排——控制理论的每个结论都依赖前一步的定义。比如讲二阶系统时如果没有先把传递函数和特征参数 ξ、ωn 讲透后面画 Bode 图、讲相位裕度时就只能靠“背结论”一旦被问“为什么相位裕度大就更稳定”就会卡住。这篇文章从建模、时域/频域分析、稳定性判据到仿真出图把机械工程控制基础课件里最常见的五十多页内容整理成一条可以直接复现的路径。适合准备新课件的高校教师、研究生助教以及需要把控制基础讲给同事听的企业内训师。文中代码按可运行优先所有图都能直接导出成 PPT 能用的 PNG。2. 机械工程控制基础建模从质量-弹簧-阻尼到传递函数2.1 为什么控制基础PPT的第一页永远是系统的数学模型机械工程控制基础这门课研究的是线性定常系统教学主线是“建模—分析—校正”。如果建模这一步跳过后面所有仿真都只是装饰。机械系统里最典型的对象是质量-弹簧-阻尼系统它用三个元件浓缩了惯性、弹性、能量耗散三种物理特征几乎所有高阶机械系统都能在一定频率范围内折合成这个形式。按牛顿第二定律写出运动方程m·x¨(t) c·ẋ(t) k·x(t) F(t)。对等式两边做拉普拉斯变换并假设零初始条件得到传递函数G(s) X(s) / F(s) 1 / (m·s² c·s k)这里做拉氏变换的理由要在 PPT 里讲清楚它把微积分运算变成代数运算把“求解微分方程”变成“整理分式多项式”。学生只有认同这一步才会接受后面所有方块图变换。常见做法是在这一页放一张物理示意图旁边放对应的方块图让学生看到“物理模型→数学模型→方块图”三者是一一对应的。2.2 用control.tf建立二阶系统的最小代码给机械工程控制基础 PPT 配仿真图不需要从零写求解器。Python 的 control 库已经把传递函数建模、时域响应、Bode 图、根轨迹这些课程主线工具都封装好了。先装依赖pip install control numpy matplotlib下面的代码建立上面提到的质量-弹簧-阻尼系统并打印传递函数import control as ct m 1.0 # 质量, kg c 2.0 # 阻尼系数, N·s/m k 20.0 # 弹簧刚度, N/m sys ct.tf([1], [m, c, k]) print(sys)代码里ct.tf(分子多项式, 分母多项式)的两个参数都按 s 的降幂排列。分子是常数 1对应输入力到输出位移的增益分母依次是 m、c、k正是上面微分方程左侧系数。打印结果会是1 ------------- s^2 2 s 20这个最小代码可以作为 PPT 里“学生动手环节”的第一步。很多教材习惯把二阶系统写成 K/(T²s² 2ξTs 1) 的形式两套参数之间的换算是考试常考点ωn √(k/m)ξ c / (2√(mk))。我一般会在这一页用旁注标出换算关系而不是单独再讲一遍因为换算本身不难难的是让学生理解 ξ 为什么叫阻尼比。2.3 建模章节PPT里的参数表该列什么这一章学生扣分大多发生在单位换算和符号混淆上所以参数表比大段文字有用。表格要包含参数、符号、单位、物理含义和在传递函数中的位置参数符号国际单位对应传递函数位置质量mkgs² 项系数阻尼系数cN·s/ms 项系数弹簧刚度kN/m常数项系数无阻尼自然频率ωnrad/s√(k/m)阻尼比ξ无单位c / (2√(mk))表格放在这一章最后一页后面讲时域指标时还要回来引用。注意不要把“阻尼系数”和“阻尼比”写在同一行这两个量一个带单位一个不带学生最容易在计算调节时间时混用。PPT 里的建议是这一章控制在 8 页以内因为后面时域、频域、稳定性三章才是需要大量时间展开的部分。3. 机械工程控制基础的时域与频域指标用Python算出能写进PPT的数3.1 时域指标超调量、调节时间、稳态误差是什么关系时域分析是机械工程控制基础里信息密度最高的一节。对欠阻尼二阶系统0 ξ 1单位阶跃响应有四个核心指标超调量、峰值时间、调节时间、稳态误差。课程里给的解析式是超调量 Mp exp(-πξ / √(1-ξ²)) 峰值时间 tp π / (ωn·√(1-ξ²)) 调节时间 ts ≈ 4 / (ξ·ωn) ±2% 准则三个公式之间有明确关联ξ 决定超调量和调节时间ωn 只决定响应快慢。讲课时最怕学生误以为“ ωn 越大超调越大”实际上 ωn 只改变时间轴缩放不改变超调百分比。我在 PPT 上会并排画两幅阶跃响应图一组固定 ξ、改变 ωn另一组固定 ωn、改变 ξ。两幅图一对比结论一目了然比任何文字描述都省时间。3.2 用step_info一键输出指标并核对公式手算公式容易抄错课堂演示时直接用 control 库输出指标。继续用前面二阶系统加一行step_infoimport control as ct sys ct.tf([20], [1, 4, 20]) # ωn²20, 2ξωn4对应 ξ≈0.447 info ct.step_info(sys) print(info)输出字典里包含RiseTime、SettlingTime、Overshoot、PeakTime、SteadyStateValue等键。用这组数据可以现场核对ξ 4/(2√20) ≈ 0.447代入超调公式约 20.5%而info[Overshoot]会给出 20.4% 左右误差来自数值仿真步长。这个对照过程本身就是很好的教学演示。需要注意step_info里的Overshoot带百分号含义但没有单位PPT 上标注时写“约 20%”而不是“0.2”RiseTime默认是响应从 10% 升到 90% 的时间调节时间默认用 ±2% 准则和很多教材一致。如果教材用的是 ±5% 准则对应系数是 3/(ξ·ωn)讲课前先看自己用的教材版本避免两套数对不上。3.3 频域Bode图和相位裕度一张图兼顾稳定性和快速性频域分析是机械工程控制基础从时域思维转向频域思维的关卡。Bode 图把系统在不同频率下的幅值增益和相位滞后画在同一幅图里相位裕度 PM 是剪切频率处相位与 -180° 的差值。工程经验是 PM 在 30° 到 60° 之间系统既有足够的稳定裕度又有可接受的响应速度。PM 过小阶跃响应会明显振荡PM 过大系统响应会变得迟钝。用 control 库算相位裕度只需要一行import control as ct sys ct.tf([20], [1, 4, 20]) gm, pm, wgm, wpm ct.margin(sys) print(f增益裕度 gm {gm:.3f}, 相位裕度 pm {pm:.2f}°) print(f相位穿越频率 wgm {wgm:.3f} rad/s, 剪切频率 wpm {wpm:.3f} rad/s)ct.margin返回四个值分别对应增益裕度、相位裕度、相位穿越频率、剪切频率。增益裕度 gm 通常用 dB 表示时写作 20·log10(gm)很多教材直接用分贝值PPT 上要统一标注单位。这一节讲课顺序我建议是先放 Bode 图标出 PM 的位置再放一页时域阶跃响应让学生直观看到 PM 与超调量的对应关系。4. 机械工程控制基础稳定性判据劳斯表与根轨迹的边界怎么算4.1 劳斯判据表格里最容易出错的全零行稳定性是机械工程控制基础课程承上启下的环节。特征方程次数超过三阶时直接求根已经不方便劳斯判据通过系数表格判断右半平面根的个数是 PPT 上必须详细演示的代数方法。学生上手最容易错的是两种特殊情况某一列首元素为零、整行为零。整行为零说明系统存在对称于原点的根比如共轭虚根需要回到上一行构造辅助方程继续。下面这段代码实现劳斯表第一列的计算适合课堂现场演示四阶以下系统import numpy as np def routh_first_column(coeffs): # coeffs: 特征多项式系数从最高次到常数项 n len(coeffs) - 1 # 系统阶次 m (n 1) // 2 # 每行长度 arr np.zeros((n 1, m)) arr[0, :len(coeffs[0::2])] coeffs[0::2] arr[1, :len(coeffs[1::2])] coeffs[1::2] for i in range(2, n 1): a arr[i-2, 0] b arr[i-1, 0] if abs(b) 1e-10: b 1e-10 # 首列为零时用小量代替 for j in range(m - 1): arr[i, j] (a * arr[i-1, j1] - arr[i-2, j1] * b) / b return arr[:, 0] # 例s^3 2s^2 3s 4 print(routh_first_column([1, 2, 3, 4]))提示这段代码用 1e-10 代替零元素能处理大多数课堂例题。真正出现全零行时需要先用辅助方程求导再填表代码不会自动处理讲课时要把这个边界条件讲清楚。第一列符号变化次数就是右半平面根的个数。运行上面的特征多项式第一列为 [1, 2, -1, -4]符号变化一次说明系统有一个不稳定的闭环极点。这正是劳斯判据的教学价值不求解特征方程也能判断稳定性边界。4.2 根轨迹K值范围的三种确认方式根轨迹法讲的是“参数从 0 变化到无穷时闭环极点怎么走”。课程选用的参数一般是开环增益 K根轨迹与虚轴的交点就是稳定边界。用 control 库画根轨迹import control as ct import matplotlib.pyplot as plt # 开环传函G(s) K / (s(s2)(s6)) sys ct.tf([1], [1, 8, 12, 0]) roots, gains ct.root_locus(sys, plotTrue) plt.show()root_locus返回每一支轨迹的根位置和对应的增益数组plotTrue时直接出图。从图上可以看到K 从 0 增大时三个极点分别从 0、-2、-6 出发其中两支在某个 K 值处穿过虚轴进入右半平面临界 K 值就是系统稳定的上限。工程上常用三种方式互相验证从根轨迹图目测交点增益用劳斯判据求临界 K 的解析值再用第 3 章的margin函数检查闭环系统的相位裕度。三种结果应一致不一致时优先怀疑坐标轴单位或特征多项式写错。PPT 上把三种结果并列展示比单讲一种方法更有说服力。4.3 三种稳定性判据怎么分配讲义位置很多老师把劳斯判据、奈奎斯特判据、根轨迹法当成并列的三节讲学生到期末容易搞混。我建议在 PPT 里加一张选用对照表把每个判据的信息来源和结论类型写清楚判据需要的信息输出结论适合出现在课程的哪个阶段劳斯判据特征方程系数是否稳定、右半平面根个数时域分析之后频域分析之前奈奎斯特判据开环频率响应曲线闭环稳定的条件与稳定裕度频域分析章节收尾衔接校正根轨迹法开环零极点、参数变化范围参数边界与动态性能变化趋势参数设计与校正章节这张表放在三章内容全部讲完之后作为回顾页。学生如果能自己填出第三列说明已经理解了三个判据分别回答什么问题、需要什么输入。如果填不出来问题多半出在“不知道为什么要用这个判据”上而不是计算过程。5. 机械工程控制基础配图实战一张能直接放进PPT的仿真图5.1 一条命令得到课程标准的阶跃响应图很多讲义的仿真图是从教材截图里“借用”的清晰度不够且坐标风格不统一。自己生成一套风格统一的图工作量其实很低。先用 rcParams 统一全局字体和网格再画阶跃响应import control as ct import matplotlib.pyplot as plt plt.rcParams.update({ font.size: 12, font.family: Noto Sans CJK SC, # Windows 可换 Microsoft YaHei axes.grid: True, grid.alpha: 0.25, }) sys ct.tf([5], [1, 2, 5]) t, y ct.step_response(sys, T8) fig, ax plt.subplots(figsize(7, 4.5)) ax.plot(t, y, linewidth1.8, color#0072B2) ax.set_xlabel(时间 t (s)) ax.set_ylabel(输出 x(t)) ax.set_title(单位阶跃响应) fig.savefig(step_response.png, dpi300, bbox_inchestight)这段代码里值得解释的参数有三个figsize(7,4.5)控制图幅比例接近 PPT 内容区的横向半幅插入后不需要再裁剪dpi300适合打印如果只投屏 150 就够文件体积更小bbox_inchestight会自动裁掉空白边避免插入 PPT 后周围留白太多。grid.alpha调低网格透明度打印时不会抢主曲线的视觉权重。5.2 字体、线宽和DPI打印和投屏两套配置不同使用场景对图片参数的要求差异很大。我整理了一张常用配置表场景DPI字体线宽说明打印纸质讲义300衬线体宋体/SimSun1.8黑白打印时靠线型区分曲线不靠颜色课堂投屏150无衬线体微软雅黑/Noto2.0远处能看清坐标文字导出后粘贴到 Word200跟随正文样式1.5不要用透明背景Word 打印易出问题中文字体在不同系统上名称不一致Windows 推荐Microsoft YaHeimacOS 推荐PingFang SCLinux 推荐Noto Sans CJK SC。如果在 rcParams 里写死字体名换电脑跑脚本时会报找不到字体我一般把字体配置单独放在脚本顶部方便每台机器改一处。5.3 两种配色方案能避免把图打成糊仿真图的配色第一个原则是不要只看颜色还要看线型和符号。第二个原则是考虑打印场景很多老师会把 PPT 打成一页六张的讲义彩色图变成灰度后靠红色和绿色区分的两条曲线会完全分不开。推荐两套配色浅底模板用色盲友好顺序#0072B2蓝、#D55E00橙、#009E73绿深色模板则把主曲线改成白色或浅灰其他曲线用较亮的辅助色。多曲线时同时设置linestyle用实线、虚线、点画线区分曲线这样即使打印成黑白也分得清。具体做法plt.rcParams[axes.prop_cycle] plt.cycler( color[#0072B2, #D55E00, #009E73] )设置之后同一幅图里连续三根 plot 会自动套用三色循环。超过三条曲线时建议拆分成子图而不是在一张图里堆六种颜色。6. 机械工程控制基础PPT的收尾技巧公式、例题与动画节奏6.1 用PPT动画让根轨迹“动”起来根轨迹是参数变化的过程静态图很难体现“K 增大闭环极点沿分支移动”的动态感。不推荐插入外部动画视频因为视频字体和格式容易在教室电脑上出问题。常见做法是把根轨迹按 K 值分割成三到五帧静态图用 PowerPoint 自带的“擦除”动画依次播放。具体操作先导出一张完整根轨迹图复制三份分别只显示 K 在 0-10、10-30、30-100 区间内的分支片段然后给这三份都设置“上一动画之后”。这样播放时曲线会像被一笔画出来一样逐段展开而且每一帧都可以停下来讲解。课堂控制的关键是每页只做一个动画动作不要把“擦除”和“浮入”叠加在同一张图上。6.2 公式编辑器与符号约定保持一致机械工程控制基础课件里公式密度高统一使用 Office 自带的公式编辑器快捷键 Alt不要用截图贴公式。截图公式在缩放时模糊而且无法像文本公式一样被 CtrlF 搜索。符号约定方面建议全篇统一时间变量 t 用斜体拉氏变量 s 用小写斜体ωn 和 ξ 作为带下标的希腊字母用正体。每个公式给出后单独一行写“变量说明”只解释本页新出现的变量。这样学生复习时不需要往前翻十页找符号定义。如果团队共用课件可以在第一页放一张符号简表后续所有页面按这个表执行能省掉大量校对时间。6.3 用Excel随机参数组快速生成课后练习最后一个技巧非常实用但很少人用用 Excel 的 RANDBETWEEN 函数快速生成一组课后练习参数。做法是在表格里生成两列随机数——阻尼比 ξ 取 0.3 到 0.8ωn 取 2 到 10然后用前面公式自动计算超调量、峰值时间、调节时间最后把公式结果粘贴成数值避免学生打开文件时随机数刷新导致答案变化。每次上课前改一次随机数种子Excel 里可以通过重新输入实现几分钟就能得到一套新参数。把这个 Excel 模板和绘图脚本放在同一个目录每学期开学改一次参数范围就能在十分钟内得到全套练习、答案和配图。本文还有配套的精品资源点击获取
简介这是一份面向机械工程及相关专业学生和初学者的《机械工程控制基础》课程PPT源自三峡大学机械与材料学院方子帆教授的课堂讲义聚焦控制理论的基本概念、系统工作原理与组成并通过恒温箱温度控制、钢铁轧制等案例讲解自动控制与反馈控制的核心逻辑。资源包仅含1个PPT文件整体大小约9.07MB内容紧凑适合用于课程预习、期末复习或教学备课。目前已有199人浏览学习。PPT以第一章绪论及基本概念为主线包含控制的定义、被控对象与控制器等术语解读、人工/自动恒温控制动态过程、系统原理方块图、反馈控制思想等页面直观重点突出能帮助读者快速建立控制系统知识框架理解传递函数、反馈等核心概念在实际工程中的作用也为后续学习现代控制理论打下基础。1. 机械工程控制基础PPT的实质不是排版是把控制过程讲顺做机械工程控制基础PPT的人通常是从课件正文开始动手的先找一版现成的讲稿把公式和波形图贴上去再调整配色。但真正的难点不在排版而在内容编排——控制理论的每个结论都依赖前一步的定义。比如讲二阶系统时如果没有先把传递函数和特征参数 ξ、ωn 讲透后面画 Bode 图、讲相位裕度时就只能靠“背结论”一旦被问“为什么相位裕度大就更稳定”就会卡住。这篇文章从建模、时域/频域分析、稳定性判据到仿真出图把机械工程控制基础课件里最常见的五十多页内容整理成一条可以直接复现的路径。适合准备新课件的高校教师、研究生助教以及需要把控制基础讲给同事听的企业内训师。文中代码按可运行优先所有图都能直接导出成 PPT 能用的 PNG。2. 机械工程控制基础建模从质量-弹簧-阻尼到传递函数2.1 为什么控制基础PPT的第一页永远是系统的数学模型机械工程控制基础这门课研究的是线性定常系统教学主线是“建模—分析—校正”。如果建模这一步跳过后面所有仿真都只是装饰。机械系统里最典型的对象是质量-弹簧-阻尼系统它用三个元件浓缩了惯性、弹性、能量耗散三种物理特征几乎所有高阶机械系统都能在一定频率范围内折合成这个形式。按牛顿第二定律写出运动方程m·x¨(t) c·ẋ(t) k·x(t) F(t)。对等式两边做拉普拉斯变换并假设零初始条件得到传递函数G(s) X(s) / F(s) 1 / (m·s² c·s k)这里做拉氏变换的理由要在 PPT 里讲清楚它把微积分运算变成代数运算把“求解微分方程”变成“整理分式多项式”。学生只有认同这一步才会接受后面所有方块图变换。常见做法是在这一页放一张物理示意图旁边放对应的方块图让学生看到“物理模型→数学模型→方块图”三者是一一对应的。2.2 用control.tf建立二阶系统的最小代码给机械工程控制基础 PPT 配仿真图不需要从零写求解器。Python 的 control 库已经把传递函数建模、时域响应、Bode 图、根轨迹这些课程主线工具都封装好了。先装依赖pip install control numpy matplotlib下面的代码建立上面提到的质量-弹簧-阻尼系统并打印传递函数import control as ct m 1.0 # 质量, kg c 2.0 # 阻尼系数, N·s/m k 20.0 # 弹簧刚度, N/m sys ct.tf([1], [m, c, k]) print(sys)代码里ct.tf(分子多项式, 分母多项式)的两个参数都按 s 的降幂排列。分子是常数 1对应输入力到输出位移的增益分母依次是 m、c、k正是上面微分方程左侧系数。打印结果会是1 ------------- s^2 2 s 20这个最小代码可以作为 PPT 里“学生动手环节”的第一步。很多教材习惯把二阶系统写成 K/(T²s² 2ξTs 1) 的形式两套参数之间的换算是考试常考点ωn √(k/m)ξ c / (2√(mk))。我一般会在这一页用旁注标出换算关系而不是单独再讲一遍因为换算本身不难难的是让学生理解 ξ 为什么叫阻尼比。2.3 建模章节PPT里的参数表该列什么这一章学生扣分大多发生在单位换算和符号混淆上所以参数表比大段文字有用。表格要包含参数、符号、单位、物理含义和在传递函数中的位置参数符号国际单位对应传递函数位置质量mkgs² 项系数阻尼系数cN·s/ms 项系数弹簧刚度kN/m常数项系数无阻尼自然频率ωnrad/s√(k/m)阻尼比ξ无单位c / (2√(mk))表格放在这一章最后一页后面讲时域指标时还要回来引用。注意不要把“阻尼系数”和“阻尼比”写在同一行这两个量一个带单位一个不带学生最容易在计算调节时间时混用。PPT 里的建议是这一章控制在 8 页以内因为后面时域、频域、稳定性三章才是需要大量时间展开的部分。3. 机械工程控制基础的时域与频域指标用Python算出能写进PPT的数3.1 时域指标超调量、调节时间、稳态误差是什么关系时域分析是机械工程控制基础里信息密度最高的一节。对欠阻尼二阶系统0 ξ 1单位阶跃响应有四个核心指标超调量、峰值时间、调节时间、稳态误差。课程里给的解析式是超调量 Mp exp(-πξ / √(1-ξ²)) 峰值时间 tp π / (ωn·√(1-ξ²)) 调节时间 ts ≈ 4 / (ξ·ωn) ±2% 准则三个公式之间有明确关联ξ 决定超调量和调节时间ωn 只决定响应快慢。讲课时最怕学生误以为“ ωn 越大超调越大”实际上 ωn 只改变时间轴缩放不改变超调百分比。我在 PPT 上会并排画两幅阶跃响应图一组固定 ξ、改变 ωn另一组固定 ωn、改变 ξ。两幅图一对比结论一目了然比任何文字描述都省时间。3.2 用step_info一键输出指标并核对公式手算公式容易抄错课堂演示时直接用 control 库输出指标。继续用前面二阶系统加一行step_infoimport control as ct sys ct.tf([20], [1, 4, 20]) # ωn²20, 2ξωn4对应 ξ≈0.447 info ct.step_info(sys) print(info)输出字典里包含RiseTime、SettlingTime、Overshoot、PeakTime、SteadyStateValue等键。用这组数据可以现场核对ξ 4/(2√20) ≈ 0.447代入超调公式约 20.5%而info[Overshoot]会给出 20.4% 左右误差来自数值仿真步长。这个对照过程本身就是很好的教学演示。需要注意step_info里的Overshoot带百分号含义但没有单位PPT 上标注时写“约 20%”而不是“0.2”RiseTime默认是响应从 10% 升到 90% 的时间调节时间默认用 ±2% 准则和很多教材一致。如果教材用的是 ±5% 准则对应系数是 3/(ξ·ωn)讲课前先看自己用的教材版本避免两套数对不上。3.3 频域Bode图和相位裕度一张图兼顾稳定性和快速性频域分析是机械工程控制基础从时域思维转向频域思维的关卡。Bode 图把系统在不同频率下的幅值增益和相位滞后画在同一幅图里相位裕度 PM 是剪切频率处相位与 -180° 的差值。工程经验是 PM 在 30° 到 60° 之间系统既有足够的稳定裕度又有可接受的响应速度。PM 过小阶跃响应会明显振荡PM 过大系统响应会变得迟钝。用 control 库算相位裕度只需要一行import control as ct sys ct.tf([20], [1, 4, 20]) gm, pm, wgm, wpm ct.margin(sys) print(f增益裕度 gm {gm:.3f}, 相位裕度 pm {pm:.2f}°) print(f相位穿越频率 wgm {wgm:.3f} rad/s, 剪切频率 wpm {wpm:.3f} rad/s)ct.margin返回四个值分别对应增益裕度、相位裕度、相位穿越频率、剪切频率。增益裕度 gm 通常用 dB 表示时写作 20·log10(gm)很多教材直接用分贝值PPT 上要统一标注单位。这一节讲课顺序我建议是先放 Bode 图标出 PM 的位置再放一页时域阶跃响应让学生直观看到 PM 与超调量的对应关系。4. 机械工程控制基础稳定性判据劳斯表与根轨迹的边界怎么算4.1 劳斯判据表格里最容易出错的全零行稳定性是机械工程控制基础课程承上启下的环节。特征方程次数超过三阶时直接求根已经不方便劳斯判据通过系数表格判断右半平面根的个数是 PPT 上必须详细演示的代数方法。学生上手最容易错的是两种特殊情况某一列首元素为零、整行为零。整行为零说明系统存在对称于原点的根比如共轭虚根需要回到上一行构造辅助方程继续。下面这段代码实现劳斯表第一列的计算适合课堂现场演示四阶以下系统import numpy as np def routh_first_column(coeffs): # coeffs: 特征多项式系数从最高次到常数项 n len(coeffs) - 1 # 系统阶次 m (n 1) // 2 # 每行长度 arr np.zeros((n 1, m)) arr[0, :len(coeffs[0::2])] coeffs[0::2] arr[1, :len(coeffs[1::2])] coeffs[1::2] for i in range(2, n 1): a arr[i-2, 0] b arr[i-1, 0] if abs(b) 1e-10: b 1e-10 # 首列为零时用小量代替 for j in range(m - 1): arr[i, j] (a * arr[i-1, j1] - arr[i-2, j1] * b) / b return arr[:, 0] # 例s^3 2s^2 3s 4 print(routh_first_column([1, 2, 3, 4]))提示这段代码用 1e-10 代替零元素能处理大多数课堂例题。真正出现全零行时需要先用辅助方程求导再填表代码不会自动处理讲课时要把这个边界条件讲清楚。第一列符号变化次数就是右半平面根的个数。运行上面的特征多项式第一列为 [1, 2, -1, -4]符号变化一次说明系统有一个不稳定的闭环极点。这正是劳斯判据的教学价值不求解特征方程也能判断稳定性边界。4.2 根轨迹K值范围的三种确认方式根轨迹法讲的是“参数从 0 变化到无穷时闭环极点怎么走”。课程选用的参数一般是开环增益 K根轨迹与虚轴的交点就是稳定边界。用 control 库画根轨迹import control as ct import matplotlib.pyplot as plt # 开环传函G(s) K / (s(s2)(s6)) sys ct.tf([1], [1, 8, 12, 0]) roots, gains ct.root_locus(sys, plotTrue) plt.show()root_locus返回每一支轨迹的根位置和对应的增益数组plotTrue时直接出图。从图上可以看到K 从 0 增大时三个极点分别从 0、-2、-6 出发其中两支在某个 K 值处穿过虚轴进入右半平面临界 K 值就是系统稳定的上限。工程上常用三种方式互相验证从根轨迹图目测交点增益用劳斯判据求临界 K 的解析值再用第 3 章的margin函数检查闭环系统的相位裕度。三种结果应一致不一致时优先怀疑坐标轴单位或特征多项式写错。PPT 上把三种结果并列展示比单讲一种方法更有说服力。4.3 三种稳定性判据怎么分配讲义位置很多老师把劳斯判据、奈奎斯特判据、根轨迹法当成并列的三节讲学生到期末容易搞混。我建议在 PPT 里加一张选用对照表把每个判据的信息来源和结论类型写清楚判据需要的信息输出结论适合出现在课程的哪个阶段劳斯判据特征方程系数是否稳定、右半平面根个数时域分析之后频域分析之前奈奎斯特判据开环频率响应曲线闭环稳定的条件与稳定裕度频域分析章节收尾衔接校正根轨迹法开环零极点、参数变化范围参数边界与动态性能变化趋势参数设计与校正章节这张表放在三章内容全部讲完之后作为回顾页。学生如果能自己填出第三列说明已经理解了三个判据分别回答什么问题、需要什么输入。如果填不出来问题多半出在“不知道为什么要用这个判据”上而不是计算过程。5. 机械工程控制基础配图实战一张能直接放进PPT的仿真图5.1 一条命令得到课程标准的阶跃响应图很多讲义的仿真图是从教材截图里“借用”的清晰度不够且坐标风格不统一。自己生成一套风格统一的图工作量其实很低。先用 rcParams 统一全局字体和网格再画阶跃响应import control as ct import matplotlib.pyplot as plt plt.rcParams.update({ font.size: 12, font.family: Noto Sans CJK SC, # Windows 可换 Microsoft YaHei axes.grid: True, grid.alpha: 0.25, }) sys ct.tf([5], [1, 2, 5]) t, y ct.step_response(sys, T8) fig, ax plt.subplots(figsize(7, 4.5)) ax.plot(t, y, linewidth1.8, color#0072B2) ax.set_xlabel(时间 t (s)) ax.set_ylabel(输出 x(t)) ax.set_title(单位阶跃响应) fig.savefig(step_response.png, dpi300, bbox_inchestight)这段代码里值得解释的参数有三个figsize(7,4.5)控制图幅比例接近 PPT 内容区的横向半幅插入后不需要再裁剪dpi300适合打印如果只投屏 150 就够文件体积更小bbox_inchestight会自动裁掉空白边避免插入 PPT 后周围留白太多。grid.alpha调低网格透明度打印时不会抢主曲线的视觉权重。5.2 字体、线宽和DPI打印和投屏两套配置不同使用场景对图片参数的要求差异很大。我整理了一张常用配置表场景DPI字体线宽说明打印纸质讲义300衬线体宋体/SimSun1.8黑白打印时靠线型区分曲线不靠颜色课堂投屏150无衬线体微软雅黑/Noto2.0远处能看清坐标文字导出后粘贴到 Word200跟随正文样式1.5不要用透明背景Word 打印易出问题中文字体在不同系统上名称不一致Windows 推荐Microsoft YaHeimacOS 推荐PingFang SCLinux 推荐Noto Sans CJK SC。如果在 rcParams 里写死字体名换电脑跑脚本时会报找不到字体我一般把字体配置单独放在脚本顶部方便每台机器改一处。5.3 两种配色方案能避免把图打成糊仿真图的配色第一个原则是不要只看颜色还要看线型和符号。第二个原则是考虑打印场景很多老师会把 PPT 打成一页六张的讲义彩色图变成灰度后靠红色和绿色区分的两条曲线会完全分不开。推荐两套配色浅底模板用色盲友好顺序#0072B2蓝、#D55E00橙、#009E73绿深色模板则把主曲线改成白色或浅灰其他曲线用较亮的辅助色。多曲线时同时设置linestyle用实线、虚线、点画线区分曲线这样即使打印成黑白也分得清。具体做法plt.rcParams[axes.prop_cycle] plt.cycler( color[#0072B2, #D55E00, #009E73] )设置之后同一幅图里连续三根 plot 会自动套用三色循环。超过三条曲线时建议拆分成子图而不是在一张图里堆六种颜色。6. 机械工程控制基础PPT的收尾技巧公式、例题与动画节奏6.1 用PPT动画让根轨迹“动”起来根轨迹是参数变化的过程静态图很难体现“K 增大闭环极点沿分支移动”的动态感。不推荐插入外部动画视频因为视频字体和格式容易在教室电脑上出问题。常见做法是把根轨迹按 K 值分割成三到五帧静态图用 PowerPoint 自带的“擦除”动画依次播放。具体操作先导出一张完整根轨迹图复制三份分别只显示 K 在 0-10、10-30、30-100 区间内的分支片段然后给这三份都设置“上一动画之后”。这样播放时曲线会像被一笔画出来一样逐段展开而且每一帧都可以停下来讲解。课堂控制的关键是每页只做一个动画动作不要把“擦除”和“浮入”叠加在同一张图上。6.2 公式编辑器与符号约定保持一致机械工程控制基础课件里公式密度高统一使用 Office 自带的公式编辑器快捷键 Alt不要用截图贴公式。截图公式在缩放时模糊而且无法像文本公式一样被 CtrlF 搜索。符号约定方面建议全篇统一时间变量 t 用斜体拉氏变量 s 用小写斜体ωn 和 ξ 作为带下标的希腊字母用正体。每个公式给出后单独一行写“变量说明”只解释本页新出现的变量。这样学生复习时不需要往前翻十页找符号定义。如果团队共用课件可以在第一页放一张符号简表后续所有页面按这个表执行能省掉大量校对时间。6.3 用Excel随机参数组快速生成课后练习最后一个技巧非常实用但很少人用用 Excel 的 RANDBETWEEN 函数快速生成一组课后练习参数。做法是在表格里生成两列随机数——阻尼比 ξ 取 0.3 到 0.8ωn 取 2 到 10然后用前面公式自动计算超调量、峰值时间、调节时间最后把公式结果粘贴成数值避免学生打开文件时随机数刷新导致答案变化。每次上课前改一次随机数种子Excel 里可以通过重新输入实现几分钟就能得到一套新参数。把这个 Excel 模板和绘图脚本放在同一个目录每学期开学改一次参数范围就能在十分钟内得到全套练习、答案和配图。本文还有配套的精品资源点击获取
