本科三年级自动化专业《动态过程分析方法》深度教学设计

本科三年级自动化专业《动态过程分析方法》深度教学设计

一、教学背景与目标设定

（一）课程定位与学情分析

本课程“动态过程分析方法”是自动化专业本科三年级学生的核心专业必修课【核心基石】。在先前学习了高等数学、复变函数与积分变换、电路原理等基础课后，学生已具备必要的数学工具和一定的物理系统认知。然而，他们往往缺乏将具体物理系统抽象为数学模型、并运用统一框架分析其动态行为的系统性思维【重要转折点】。大三学生具备较强的逻辑思维和自学能力，但对工程问题的复杂性和多学科交叉特性认识尚浅，容易陷入纯数学推导而忽视物理意义和工程应用【常见误区】。因此，本教学设计旨在架设从理论到实践的桥梁，重点培养学生建立“系统观”，掌握分析动态过程的通用方法论。

（二）教学目标体系

依据布鲁姆认知目标分类及工程教育专业认证（OBE）理念，本课程设定了分层递进的教学目标：

1.知识层面（记忆与理解）：学生能准确复述动态系统、数学模型、传递函数、频率特性、稳定性等核心概念【基础】；能解释一阶、二阶系统典型环节的动态响应特性及其与极点分布的内在联系【重要】。

2.能力层面（应用与分析）：能够针对机械、电子、机电等典型工程系统【跨学科视野1】，运用机理分析法建立其微分方程模型，并推导出传递函数【核心技能】；能够熟练运用时域分析法（性能指标计算）、根轨迹法、频域分析法（Bode图、Nyquist判据）分析参数变化对系统动态性能及稳定性的影响【高频考点与难点】；能够利用MATLAB/Simulink等工具对简单动态系统进行仿真建模与特性分析【现代工程师必备】。

3.素养层面（评价与创造）：初步具备从“动态过程”角度审视工程问题的意识，能够评价不同分析方法的适用场景与优劣【高阶思维】；在小组项目（如“直流电机转速控制系统分析与设计”）中，能够综合运用所学方法，提出改进系统动态性能的方案雏形，并撰写规范的实验/分析报告【终极目标】。

（三）教学内容重构与资源整合

打破传统“数学推导先行”的枯燥模式，本设计采用“问题导向+案例贯穿+工具赋能”的策略重构内容。主线为“描述-分析-设计”三大模块。

1.描述：如何用数学模型刻画系统的动态过程？（微分方程、传递函数、状态空间）

2.分析：如何评判系统动态过程的好坏？（稳定性、准确性、快速性三大指标）

3.设计初步：如何改变系统结构或参数以获得理想的动态过程？（PID校正思想引入）

核心资源包括：自主研发的“多学科动态过程仿真案例库”（涵盖弹簧-质量-阻尼系统、RC电路网路、水温控制系统等）【热点：虚拟仿真实验】；MATLAB实时脚本及AppDesigner开发的交互式分析工具；精选国际经典教材（如Ogata《现代控制工程》）章节与前沿学术论文片段（如关于“软机器人动态建模”的综述）【跨学科视野2】。

二、教学实施过程深度设计

（一）第一阶段：思维的破冰——从“静态世界”到“动态视野”（2学时）

1.情境导入：生活与工程中的“动态之问”

课程伊始，不急于给出定义，而是展示一系列引人深思的现象视频与动画：高速列车进站时的精确停靠、无人机在阵风中的姿态调整、3D打印机喷头的精准轨迹、甚至是一杯热水自然冷却的温度变化曲线。教师连续抛出问题：“这些现象的背后，是什么量在随时间变化？它们遵循怎样的内在规律？我们如何预测并控制它们的变化过程？”【激发好奇心与求知欲】通过头脑风暴，引导学生意识到，相较于已经熟悉的“稳态”计算（如欧姆定律、静力学平衡），工程世界更本质的是充满“变化”的动态过程。

2.核心概念的建立：从具体到抽象

以一个最简单的RC电路（图1-1）和一维弹簧-质量-阻尼系统（图1-2）作为“启蒙双生子”【重要类比】。

（1）机理建模演示（教师板演）：在黑板上并行推导两个系统的微分方程。对于RC电路：基于基尔霍夫定律，得R

C

d

u

c

(

t

)

d

t

+

u

c

(

t

)

=

u

r

(

t

)

RC\frac{du_c(t)}{dt}+u_c(t)=u_r(t)

RCdtduc​(t)​+uc​(t)=ur​(t)。对于机械系统：基于牛顿第二定律，得m

d

2

y

(

t

)

d

t

2

+

f

d

y

(

t

)

d

t

+

k

y

(

t

)

=

F

(

t

)

m\frac{d^2y(t)}{dt^2}+f\frac{dy(t)}{dt}+ky(t)=F(t)

mdt2d2y(t)​+fdtdy(t)​+ky(t)=F(t)。

（2）类比与归纳（师生互动）：引导学生对比两个方程的结构相似性——均包含输入、输出、以及输出的各阶导数。教师由此引出【核心基石1】——动态系统的本质：其输出响应不仅取决于当前输入，还与系统的“记忆”（即过去的状态，由微分方程的阶次体现）有关。进而给出线性时不变（LTI）系统的定义及其重要特性：齐次性与叠加性。

（3）数学工具的引入：为什么需要拉普拉斯变换？通过解上述微分方程的繁琐过程（回顾经典解法：齐次解+特解，需由初始条件确定积分常数），展示其在处理复杂输入或高阶系统时的局限性。此时，引入拉普拉斯变换这一“降维打击”利器【核心工具】。简要回顾其定义、典型函数的变换、微分定理（重中之重），并立刻应用到RC电路中：在零初始条件下，将微分方程转化为代数方程：R

C

s

U

c

(

s

)

+

U

c

(

s

)

=

U

r

(

s

)

RCsU_c(s)+U_c(s)=U_r(s)

RCsUc​(s)+Uc​(s)=Ur​(s)。

3.【核心基石2】传递函数：系统的“基因指纹”

由上述代数方程，直接定义传递函数G

(

s

)

=

U

c

(

s

)

U

r

(

s

)

=

1

R

C

s

+

1

G(s)=\frac{U_c(s)}{U_r(s)}=\frac{1}{RCs+1}

G(s)=Ur​(s)Uc​(s)​=RCs+11​。强调其核心含义：系统固有属性的频域表征，与输入无关。它浓缩了系统内部结构的所有信息，是连接输入与输出的“黑箱”内核。学生第一次体会到，一个复杂的物理过程，其动态本质可以被一个简洁的有理分式所概括【思维震撼】。布置课后思考题：传递函数的分母多项式（特征方程）为什么是系统“记忆”和“生命力”的源泉？【预习伏笔】

（二）第二阶段：动态行为的时域揭示——系统在时间轴上的“心电图”（6学时）

本阶段聚焦于在时间域内直接观察系统对典型输入信号的响应，这是评价动态过程最直观的方式。

1.典型输入信号的工程意义

阐述为什么要研究单位阶跃、单位脉冲、单位斜坡等信号【重要】。它们分别对应着工程中的突然变化、瞬间冲击和匀速跟踪场景，是对实际复杂输入的一种理想化抽象和“压力测试”。

2.一阶系统的动态过程分析【基础与重点】

（1）数学模型与特征：以G

(

s

)

=

K

T

s

+

1

G(s)=\frac{K}{Ts+1}

G(s)=Ts+1K​为统一模型。明确K为稳态增益，T为时间常数。时间常数T是表征一阶系统惯性或热容量的核心参数，其大小决定了系统响应速度的唯一因素。

（2）阶跃响应分析（教师精讲+MATLAB演示）：推导输出c

(

t

)

=

K

(

1

−

e

−

t

/

T

)

c(t)=K(1-e^{-t/T})

c(t)=K(1−e−t/T)。绘制响应曲线，逐一定义并解释动态性能指标：【高频考点】

1.3.时间常数T：响应从0上升到稳态值的63.2%所需的时间。几何意义是响应曲线在t=0处切线与稳态值交点对应的时刻。

2.4.调节时间t

s

t_s

ts​：响应进入并保持在稳态值±2%或±5%误差带内所需的最短时间。估算公式t

s

≈

3

T

t_s\approx3T

ts​≈3T（5%误差带）或4

T

4T

4T（2%误差带）。

3.5.稳态值：K

K

K。没有超调量和峰值时间。

（3）脉冲响应分析：推导c

(

t

)

=

K

T

e

−

t

/

T

c(t)=\frac{K}{T}e^{-t/T}

c(t)=TK​e−t/T。指出脉冲响应实质上是系统传递函数的拉普拉斯逆变换，包含了系统的全部动态信息【难点理解】。脉冲响应与阶跃响应之间存在积分关系。

（4）重要性质总结：【核心基石3】一阶系统响应单调上升，无振荡，其动态过程完全由T决定。T越小，响应越快。

（5）交互式实验：学生分组在MATLAB中修改T值，实时观察响应曲线的变化，亲身体验参数对动态过程的影响。

6.二阶系统的动态过程分析【核心与高频高点】

二阶系统是控制工程中最常见、最具代表性的系统，其动态过程远比一阶丰富，能呈现振荡、衰减等多种形态。

（1）标准型与特征参数：引出二阶系统闭环传递函数的标准型G

(

s

)

=

ω

n

2

s

2

+

2

ζ

ω

n

s

+

ω

n

2

G(s)=\frac{\omega_n^2}{s^2+2\zeta\omega_ns+\omega_n^2}

G(s)=s2+2ζωn​s+ωn2​ωn2​​。定义【重要参数】：无阻尼自然振荡频率ω

n

\omega_n

ωn​和阻尼比ζ

\zeta

ζ。强调这两个参数是决定二阶系统动态行为的“基因”。

（2）极点分布与阻尼比的关系：建立特征方程s

2

+

2

ζ

ω

n

s

+

ω

n

2

=

0

s^2+2\zeta\omega_ns+\omega_n^2=0

s2+2ζωn​s+ωn2​=0，求解极点s

1

,

2

=

−

ζ

ω

n

±

ω

n

ζ

2

−

1

s_{1,2}=-\zeta\omega_n\pm\omega_n\sqrt{\zeta^2-1}

s1,2​=−ζωn​±ωn​ζ2−1<pathd="M263,681c0.7,0,18,39.7,52,119

c34,79.3,68.167,158.7,102.5,238c34.3,79.3,51.8,119.3,52.5,120

c340,-704.7,510.7,-1060.3,512,-1067

l0-0

c4.7,-7.3,11,-11,19,-11

H40000v40H1012.3

s-271.3,567,-271.3,567c-38.7,80.7,-84,175,-136,283c-52,108,-89.167,185.3,-111.5,232

c-22.3,46.7,-33.8,70.3,-34.5,71c-4.7,4.7,-12.3,7,-23,7s-12,-1,-12,-1

s-109,-253,-109,-253c-72.7,-168,-109.3,-252,-110,-252c-10.7,8,-22,16.7,-34,26

c-22,17.3,-33.3,26,-34,26s-26,-26,-26,-26s76,-59,76,-59s76,-60,76,-60z

M100180h400000v40h-400000z">

​。根据ζ

\zeta

ζ的取值范围，将系统分为五种情况：【必考分类】

1.7.过阻尼(ζ

>

1

\zeta>1

ζ>1)：两个不等的负实根。

2.8.临界阻尼(ζ

=

1

\zeta=1

ζ=1)：两个相等的负实根。

3.9.欠阻尼(0

<

ζ

<

1

0<\zeta<1

0<ζ<1)：一对具有负实部的共轭复根。

4.10.无阻尼(ζ

=

0

\zeta=0

ζ=0)：一对纯虚根。

5.11.负阻尼(ζ

<

0

\zeta<0

ζ<0)：极点位于右半平面。

（3）欠阻尼二阶系统的阶跃响应分析（深度剖析）：

推导时域响应表达式（通过拉普拉斯逆变换）：c

(

t

)

=

1

−

e

−

ζ

ω

n

t

1

−

ζ

2

sin

⁡

(

ω

d

t

+

β

)

c(t)=1-\frac{e^{-\zeta\omega_nt}}{\sqrt{1-\zeta^2}}\sin(\omega_dt+\beta)

c(t)=1−1−ζ2<pathd="M95,702

c-2.7,0,-7.17,-2.7,-13.5,-8c-5.8,-5.3,-9.5,-10,-9.5,-14

c0,-2,0.3,-3.3,1,-4c1.3,-2.7,23.83,-20.7,67.5,-54

c44.2,-33.3,65.8,-50.3,66.5,-51c1.3,-1.3,3,-2,5,-2c4.7,0,8.7,3.3,12,10

s173,378,173,378c0.7,0,35.3,-71,104,-213c68.7,-142,137.5,-285,206.5,-429

c69,-144,104.5,-217.7,106.5,-221

l0-0

c5.3,-9.3,12,-14,20,-14

H400000v40H845.2724

s-225.272,467,-225.272,467s-235,486,-235,486c-2.7,4.7,-9,7,-19,7

c-6,0,-10,-1,-12,-3s-194,-422,-194,-422s-65,47,-65,47z

M83480h400000v40h-400000z">

​e−ζωn​t​sin(ωd​t+β)，其中ω

d

=

ω

n

1

−

ζ

2

\omega_d=\omega_n\sqrt{1-\zeta^2}

ωd​=ωn​1−ζ2<pathd="M263,681c0.7,0,18,39.7,52,119

c34,79.3,68.167,158.7,102.5,238c34.3,79.3,51.8,119.3,52.5,120

c340,-704.7,510.7,-1060.3,512,-1067

l0-0

c4.7,-7.3,11,-11,19,-11

H40000v40H1012.3

s-271.3,567,-271.3,567c-38.7,80.7,-84,175,-136,283c-52,108,-89.167,185.3,-111.5,232

c-22.3,46.7,-33.8,70.3,-34.5,71c-4.7,4.7,-12.3,7,-23,7s-12,-1,-12,-1

s-109,-253,-109,-253c-72.7,-168,-109.3,-252,-110,-252c-10.7,8,-22,16.7,-34,26

c-22,17.3,-33.3,26,-34,26s-26,-26,-26,-26s76,-59,76,-59s76,-60,76,-60z

M100180h400000v40h-400000z">

​为阻尼振荡频率，β

=

arccos

⁡

ζ

\beta=\arccos\zeta

β=arccosζ。

结合推导过程和响应曲线，逐一定义并推导各动态性能指标与ζ

\zeta

ζ、ω

n

\omega_n

ωn​的定量关系【高频考点与难点】：

6.12.上升时间t

r

t_r

tr​：响应从终值的10%上升到90%（或从0第一次上升到终值）所需的时间。近似有t

r

≈

π

−

β

ω

d

t_r\approx\frac{\pi-\beta}{\omega_d}

tr​≈ωd​π−β​。

7.13.峰值时间t

p

t_p

tp​：响应达到第一个峰值所需的时间。严格推导得t

p

=

π

ω

d

t_p=\frac{\pi}{\omega_d}

tp​=ωd​π​。可见ω

d

\omega_d

ωd​越大，t

p

t_p

tp​越小。

8.14.最大超调量M

p

M_p

Mp​：响应最大峰值超出稳态值的百分比。M

p

=

e

−

π

ζ

/

1

−

ζ

2

×

100

%

M_p=e^{-\pi\zeta/\sqrt{1-\zeta^2}}\times100\%

Mp​=e−πζ/1−ζ2<pathd="M95,702

c-2.7,0,-7.17,-2.7,-13.5,-8c-5.8,-5.3,-9.5,-10,-9.5,-14

c0,-2,0.3,-3.3,1,-4c1.3,-2.7,23.83,-20.7,67.5,-54

c44.2,-33.3,65.8,-50.3,66.5,-51c1.3,-1.3,3,-2,5,-2c4.7,0,8.7,3.3,12,10

s173,378,173,378c0.7,0,35.3,-71,104,-213c68.7,-142,137.5,-285,206.5,-429

c69,-144,104.5,-217.7,106.5,-221

l0-0

c5.3,-9.3,12,-14,20,-14

H400000v40H845.2724

s-225.272,467,-225.272,467s-235,486,-235,486c-2.7,4.7,-9,7,-19,7

c-6,0,-10,-1,-12,-3s-194,-422,-194,-422s-65,47,-65,47z

M83480h400000v40h-400000z">

​×100%。这是一个【极其重要的公式】，它唯一地由ζ

\zeta

ζ决定，ζ

\zeta

ζ越大，M

p

M_p

Mp​越小。典型值：ζ

=

0.707

\zeta=0.707

ζ=0.707时，M

p

≈

4.3

%

M_p\approx4.3\%

Mp​≈4.3%。

9.15.调节时间t

s

t_s

ts​：估算公式t

s

≈

3.5

ζ

ω

n

t_s\approx\frac{3.5}{\zeta\omega_n}

ts​≈ζωn​3.5​（2%误差带）。可见乘积ζ

ω

n

\zeta\omega_n

ζωn​（即极点实部的绝对值）决定了响应衰减的速度。

（4）临界阻尼与过阻尼分析：简要介绍这两种情况下响应无超调、单调上升的特点，并与欠阻尼情况进行对比。

（5）重要结论总结：【核心枢纽】二阶系统的动态过程是“快”（由ω

n

\omega_n

ωn​影响t

r

,

t

p

t_r,t_p

tr​,tp​）与“稳”（由ζ

\zeta

ζ影响M

p

M_p

Mp​）的矛盾统一体。设计者需要根据工程要求在两者间进行折衷。

（6）MATLAB/Simulink深度仿真实验（2学时）：

10.16.任务1：构建标准二阶系统仿真模型。固定ω

n

=

10

\omega_n=10

ωn​=10，分别设置ζ

=

0.2

,

0.5

,

0.707

,

1.0

,

2.0

\zeta=0.2,0.5,0.707,1.0,2.0

ζ=0.2,0.5,0.707,1.0,2.0，运行仿真并记录阶跃响应曲线。在同一张图上对比，直观感受阻尼比对响应形态（超调、振荡次数）的决定性影响。

11.17.任务2：固定ζ

=

0.5

\zeta=0.5

ζ=0.5，分别设置ω

n

=

5

,

10

,

20

\omega_n=5,10,20

ωn​=5,10,20，观察响应曲线。对比发现ω

n

\omega_n

ωn​增大时，上升时间和峰值时间显著缩短，但超调量保持不变（验证了理论）。

12.18.任务3（探究性）：引入附加零点或极点（高阶系统简化铺垫），观察其对二阶系统动态性能的“污染”效应，例如闭环零点会增加超调，闭环非主导极点会延缓响应【高阶思维训练】。

13.19.要求：每组提交一份简明的仿真实验报告，包含图形和基于数据的分析结论。

（三）第三阶段：动态过程的根轨迹透视——参数变化的全景地图（4学时）

在掌握了时域分析方法后，学生面临新问题：当系统某个参数（如增益K）变化时，闭环极点如何运动？动态性能又将如何随之变化？根轨迹法提供了一幅清晰的全景图【核心工具2】。

1.根轨迹的基本概念

定义：当系统开环传递函数的某一参数（通常为开环增益K）从0变化到∞时，闭环极点在s平面上描绘出的轨迹。强调其【核心价值】：

1.2.避免了每改变一次K就求解一次特征方程的繁琐计算。

2.3.直观展示了K变化对系统稳定性、动态性能（通过极点的位置判断）的全貌影响。

4.绘制根轨迹的幅值条件和相角条件【基础】

从闭环特征方程1

+

G

(

s

)

H

(

s

)

=

0

1+G(s)H(s)=0

1+G(s)H(s)=0推导出G

(

s

)

H

(

s

)

=

−

1

G(s)H(s)=-1

G(s)H(s)=−1。由此引出两个等价条件：

1.5.幅值条件：∣

G

(

s

)

H

(

s

)

∣

=

1

|G(s)H(s)|=1

∣G(s)H(s)∣=1（用于确定给定极点对应的K值）

2.6.相角条件：∠

G

(

s

)

H

(

s

)

=

(

2

k

+

1

)

π

\angleG(s)H(s)=(2k+1)\pi

∠G(s)H(s)=(2k+1)π（用于寻找所有在根轨迹上的点）

重点强调，相角条件是根轨迹的“灵魂”，是决定s平面上一个点是否在根轨迹上的充要条件；而幅值条件主要用于“定速”。

7.常规根轨迹的绘制法则（8条基本法则精讲）【高频考点】

结合一个具体示例系统G

(

s

)

H

(

s

)

=

K

s

(

s

+

1

)

(

s

+

2

)

G(s)H(s)=\frac{K}{s(s+1)(s+2)}

G(s)H(s)=s(s+1)(s+2)K​进行逐步图解演示，讲解并应用主要法则：

（1）起点与终点：起于开环极点，终于开环零点或无穷远。

（2）分支数与对称性：分支数等于闭环极点数，且对称于实轴。

（3）实轴上的根轨迹：实轴上某线段右侧的开环零、极点个数之和为奇数时，该线段是根轨迹的一部分【核心法则，需举例反复练习】。

（4）渐近线：当n

>

m

n>m

n>m时，指向无穷远的根轨迹沿着渐近线行进。计算渐近线与实轴的交点σ

a

\sigma_a

σa​和夹角φ

a

\varphi_a

φa​。

（5）分离点与会合点：根轨迹在实轴上相遇后分离的点。通过求解方程d

K

d

s

=

0

\frac{dK}{ds}=0

dsdK​=0或∑

1

s

−

p

i

=

∑

1

s

−

z

i

\sum\frac{1}{s-p_i}=\sum\frac{1}{s-z_i}

∑s−pi​1​=∑s−zi​1​来获得候选点。

（6）与虚轴的交点：应用劳斯稳定判据或令s

=

j

ω

s=j\omega

s=jω代入特征方程求解。该交点对应的K值即为临界稳定增益K

m

K_m

Km​，是【系统稳定性的重要边界】。

（7）出射角与入射角：对于复数极点或零点，计算根轨迹离开或进入时的角度。

（8）闭环极点的确定：给定K值，利用幅值条件和试探法在根轨迹上找到对应的闭环极点。

8.根轨迹分析系统动态性能【核心应用】

（1）稳定性判断：只要根轨迹全部位于s左半平面，系统对所有K稳定；一旦根轨迹越过虚轴进入右半平面，系统便不稳定。从根轨迹图上一目了然。

（2）动态性能分区：在根轨迹图上划分区域。

1.9.稳定区域（左半平面）：其中靠近虚轴的部分为“主导极点区”，主要决定动态响应。

2.10.不稳定区域（右半平面）。

3.11.等阻尼比线：过原点与负实轴夹角为β

\beta

β（ζ

=

cos

⁡

β

\zeta=\cos\beta

ζ=cosβ）的射线。根轨迹与某条等阻尼比线的交点，即为具有该阻尼比的闭环极点，由此可确定对应的K值和系统超调量。

4.12.等自然频率线：以原点为圆心，ω

n

\omega_n

ωn​为半径的圆弧。用于确定系统的响应速度。

（3）案例分析：对于示例系统，分析当K从小变大时，系统动态过程的变化规律：

5.13.0

<

K

<

K

1

0<K<K_1

0<K<K1​（分离点对应增益）：三个负实极点，过阻尼或临界阻尼响应，无超调。

6.14.K

1

<

K

<

K

m

K_1<K<K_m

K1​<K<Km​：一对共轭复极点+一个负实极点，欠阻尼响应，且随着K增加，共轭复极点向虚轴移动，阻尼比ζ

\zeta

ζ减小，超调量M

p

M_p

Mp​增大，调节时间先减后增（需具体计算）。

7.15.K

=

K

m

K=K_m

K=Km​：系统处于等幅振荡边界，临界稳定。

8.16.K

>

K

m

K>K_m

K>Km​：闭环极点进入右半平面，系统不稳定。

（4）MATLAB辅助根轨迹分析实验（1学时）：

9.17.任务：使用MATLAB的rlocus命令绘制上述示例系统的根轨迹图。

10.18.交互操作：使用rlocfind命令在根轨迹图上选取特定点，获取该点的增益K和对应的所有闭环极点。学生可以动态拖动滑块改变K值，观察极点位置和对应阶跃响应曲线的实时联动变化【热点：虚拟仿真交互】。这种“所见即所得”的方式，能深刻建立参数K、极点位置、动态性能三者之间的内在联系。

（四）第四阶段：动态过程的频域阐释——从另一个维度看世界（6学时）

时域分析直观，但求解高阶系统响应复杂，且难以直接从物理上指导系统校正。频域分析法利用系统对正弦信号的稳态响应来分析动态过程，具有实验便捷、物理概念清晰的特点【核心工具3】。

1.频率特性的基本概念

（1）定义与获取：对于稳定LTI系统，输入正弦信号r

(

t

)

=

A

sin

⁡

ω

t

r(t)=A\sin\omegat

r(t)=Asinωt，其稳态输出也是同频率的正弦信号c

s

(

t

)

=

A

∣

G

(

j

ω

)

∣

sin

⁡

(

ω

t

+

∠

G

(

j

ω

)

)

c_s(t)=A|G(j\omega)|\sin(\omegat+\angleG(j\omega))

cs​(t)=A∣G(jω)∣sin(ωt+∠G(jω))。定义频率特性G

(

j

ω

)

=

∣

G

(

j

ω

)

∣

e

j

∠

G

(

j

ω

)

G(j\omega)=|G(j\omega)|e^{j\angleG(j\omega)}

G(jω)=∣G(jω)∣ej∠G(jω)，它是传递函数s

=

j

ω

s=j\omega

s=jω时的特例。强调频率特性描述的是系统对不同频率正弦信号的“幅值放大/衰减能力”和“相位滞后/超前能力”【核心思想】。

（2）物理意义：演示一个低通滤波器的实例，输入低频正弦波可以顺利通过（幅值增益≈1，相位滞后小），输入高频正弦波则被严重衰减（幅值增益≈0，相位滞后大）。让学生建立“系统对频率有选择性”的直观印象。

2.频率特性的图示方法

（1）极坐标图（Nyquist图）：介绍绘制方法（以ω

\omega

ω为参变量，将G

(

j

ω

)

G(j\omega)

G(jω)的实部和虚部画在复平面上）。重点分析典型环节（比例、积分、惯性、一阶微分、振荡、二阶微分）的Nyquist图特点。

（2）对数坐标图（Bode图）【重中之重】：讲解其构成——对数幅频特性图（L

(

ω

)

=

20

lg

⁡

∣

G

(

j

ω

)

∣

L(\omega)=20\lg|G(j\omega)|

L(ω)=20lg∣G(jω)∣，单位分贝dB）和对数相频特性图。阐述使用对数坐标的巨大优点：

1.3.将幅值的乘除运算转化为加减运算。

2.4.可以方便地用渐近线绘制近似的幅频特性曲线。

3.5.能够在一张图上涵盖从低频到高频很宽的频率范围。

分步骤详细讲解典型环节的Bode图绘制方法，特别是惯性环节的转折频率ω

=

1

/

T

\omega=1/T

ω=1/T、振荡环节的谐振峰（与ζ

\zeta

ζ相关）等【高频考点】。

6.乃奎斯特稳定判据（NyquistCriterion）【难点与核心考点】

（1）引入辅助函数与映射定理：从开环频率特性G

(

s

)

H

(

s

)

G(s)H(s)

G(s)H(s)判断闭环系统的稳定性，是频域分析的精髓。简要介绍幅角原理的基本思想（绕行s平面封闭曲线，映射到F

(

s

)

=

1

+

G

(

s

)

H

(

s

)

F(s)=1+G(s)H(s)

F(s)=1+G(s)H(s)平面上的曲线包围原点的次数，进而关联到G

(

s

)

H

(

s

)

G(s)H(s)

G(s)H(s)平面上的曲线包围(-1,j0)点的次数）。

（2）Nyquist判据的标准表述：若开环系统稳定（P=0），则闭环系统稳定的充要条件是开环频率特性G

(

j

ω

)

H

(

j

ω

)

G(j\omega)H(j\omega)

G(jω)H(jω)曲线（ω

\omega

ω从0

→

∞

0\to\infty

0→∞）不包围(-1,j0)点；若包围，则闭环不稳定。若开环不稳定（有P个右半平面极点），则闭环系统稳定的充要条件是开环频率特性曲线逆时针包围(-1,j0)点P/2次（按常规正负方向规定）。

（3）穿越的概念与半Nyquist曲线：引入正负穿越的概念，将判据简化，便于在Bode图上应用。

（4）稳定裕度（StabilityMargins）【工程应用关键】：

1.7.定义：衡量一个稳定系统距离不稳定边界的程度，是系统“稳健性”的重要指标。

2.8.相角裕度（PM,PhaseMargin）：在增益穿越频率ω

c

\omega_c

ωc​（∣

G

(

j

ω

)

H

(

j

ω

)

∣

=

1

|G(j\omega)H(j\omega)|=1

∣G(jω)H(jω)∣=1）处，使系统达到临界稳定所需增加的相位滞后量。公式P

M

=

180

∘

+

∠

G

(

j

ω

c

)

H

(

j

ω

c

)

PM=180^\circ+\angleG(j\omega_c)H(j\omega_c)

PM=180∘+∠G(jωc​)H(jωc​)。通常要求PM在30

∘

→

60

∘

30^\circ\to60^\circ

30∘→60∘之间【工程经验值】。

3.9.幅值裕度（GM,GainMargin）：在相位穿越频率ω

g

\omega_g

ωg​（∠

G

(

j

ω

g

)

H

(

j

ω

g

)

=

−

180

∘

\angleG(j\omega_g)H(j\omega_g)=-180^\circ

∠G(jωg​)H(jωg​)=−180∘）处，使系统达到临界稳定所需增加的开环增益倍数。G

M

=

1

/

∣

G

(

j

ω

g

)

H

(

j

ω

g

)

∣

GM=1/|G(j\omega_g)H(j\omega_g)|

GM=1/∣G(jωg​)H(jωg​)∣，或用分贝表示G

M

d

B

=

20

lg

⁡

G

M

GM_{dB}=20\lgGM

GMdB​=20lgGM。通常要求GM>6dB。

（5）Bode图与稳定裕度的关系：引导学生掌握在Bode图上直接读取相角裕度和幅值裕度的方法，这对于工程设计和系统调试具有极高的实用价值【核心技能】。

10.频域性能指标与时域性能的关联【高阶综合】

这是打通频域和时域的桥梁，也是本阶段的【思维高潮】。

（1）闭环频域指标：谐振峰值M

r

M_r

Mr​和谐振频率ω

r

\omega_r

ωr​、带宽ω

b

\omega_b

ωb​。

（2）二阶系统中频域与时域指标的定量关系（推导并总结）：

1.11.M

r

M_r

Mr​与ζ

\zeta

ζ的关系：对于欠阻尼二阶系统，M

r

=

1

2

ζ

1

−

ζ

2

M_r=\frac{1}{2\zeta\sqrt{1-\zeta^2}}

Mr​=2ζ1−ζ2<pathd="M95,702

c-2.7,0,-7.17,-2.7,-13.5,-8c-5.8,-5.3,-9.5,-10,-9.5,-14

c0,-2,0.3,-3.3,1,-4c1.3,-2.7,23.83,-20.7,67.5,-54

c44.2,-33.3,65.8,-50.3,66.5,-51c1.3,-1.3,3,-2,5,-2c4.7,0,8.7,3.3,12,10

s173,378,173,378c0.7,0,35.3,-71,104,-213c68.7,-142,137.5,-285,206.5,-429

c69,-144,104.5,-217.7,106.5,-221

l0-0

c5.3,-9.3,12,-14,20,-14

H400000v40H845.2724

s-225.272,467,-225.272,467s-235,486,-235,486c-2.7,4.7,-9,7,-19,7

c-6,0,-10,-1,-12,-3s-194,-422,-194,-422s-65,47,-65,47z

M83480h400000v40h-400000z">

​1​(ζ

≤

0.707

\zeta\le0.707

ζ≤0.707)。由此可知，M

r

M_r

Mr​越大（谐振峰越高），超调量M

p

M_p

Mp​也越大，两者正相关。

2.12.ω

r

\omega_r

ωr​与ω

n

\omega_n

ωn​、ζ

\zeta

ζ的关系：ω

r

=

ω

n

1

−

2

ζ

2

\omega_r=\omega_n\sqrt{1-2\zeta^2}

ωr​=ωn​1−2ζ2<pathd="M263,681c0.7,0,18,39.7,52,119

c34,79.3,68.167,158.7,102.5,238c34.3,79.3,51.8,119.3,52.5,120

c340,-704.7,510.7,-1060.3,512,-1067

l0-0

c4.7,-7.3,11,-11,19,-11

H40000v40H1012.3

s-271.3,567,-271.3,567c-38.7,80.7,-84,175,-136,283c-52,108,-89.167,185.3,-111.5,232

c-22.3,46.7,-33.8,70.3,-34.5,71c-4.7,4.7,-12.3,7,-23,7s-12,-1,-12,-1

s-109,-253,-109,-253c-72.7,-168,-109.3,-252,-110,-252c-10.7,8,-22,16.7,-34,26

c-22,17.3,-33.3,26,-34,26s-26,-26,-26,-26s76,-59,76,-59s76,-60,76,-60z

M100180h400000v40h-400000z">

​。

3.13.带宽ω

b

\omega_b

ωb​与响应速度：带宽越宽，系统复现快速变化信号的能力越强，即上升时间t

r

t_r

tr​越短，响应越快。

4.14.相角裕度PM与阻尼比ζ

\zeta

ζ的近似关系：对于高阶系统，常用经验公式P

M

≈

100

ζ

PM\approx100\zeta

PM≈100ζ（ζ

≤

0.7

\zeta\le0.7

ζ≤0.7）进行快速估算。

（3）案例分析：给定一个系统的Bode图，要求估算其时域动态性能（超调量、调节时间等），并与直接时域仿真的结果进行对比验证，使学生深刻体会频域分析方法的工程实用价值。

（五）第五阶段：动态过程分析方法的综合应用——以“直流电机转速控制系统”为载体的项目式学习（4学时+课外）

本阶段旨在将前序学习的碎片化知识整合起来，应用于一个真实的工程对象，实现能力的升华。

1.项目发布：直流电机转速控制系统建模、分析与校正初步

提供一个实际的直流电机（如Maxon品牌某型号）的参数（电枢电阻、电感、力矩常数、反电势常数、转子惯量、粘性摩擦系数等）。

2.项目实施流程（小组协作，教师引导）

（1）机理建模【跨学科整合】：

1.3.电网络方程：u

a

(

t

)

=

R

a

i

a

(

t

)

+

L

a

d

i

a

(

t

)

d

t

+

e

b

(

t

)

u_a(t)=R_ai_a(t)+L_a\frac{di_a(t)}{dt}+e_b(t)

ua​(t)=Ra​ia​(t)+La​dtdia​(t)​+eb​(t)。

2.4.机械运动方程：J

d

ω

(

t

)

d

t

+

b

ω

(

t

)

=

T

e

(

t

)

=

K

t

i

a

(

t

)

J\frac{d\omega(t)}{dt}+b\omega(t)=T_e(t)=K_ti_a(t)

Jdtdω(t)​+bω(t)=Te​(t)=Kt​ia​(t)。

3.5.反电势方程：e

b

(

t

)

=

K

e

ω

(

t

)

e_