题型④:频域分析(奈奎斯特判稳 + Bode 图求裕度)
写在最前面
你已经做出了三类题:求传递函数、时域分析、根轨迹。这一个模块,我们进第 5 章——频域分析。
先给你打个预防针,也给你吃颗定心丸:
- 预防针:第 5 章是整门课里"换了个世界观"的一章。前四章你一直在看"输入一个阶跃,输出怎么随时间跑";这一章突然改成"输入正弦波,看输出的大小和快慢怎么变"。这个转弯是第 5 章最大的坎,大部分人卡的不是计算,是"这到底在干嘛"。我会用一整节从零把这个坎讲平。
- 定心丸:一旦转过这个弯,这类大题的动手部分是全书最套路化的之一——五步走,每步都是死流程。考点分析里学长的原话是"这种题出的比较基础,也是按套路做题"(来源见文件顶部)。
目标一句话:这个模块走完,下面那道 20 分的大题,你自己拿笔能写出来。
篇幅长,是因为频域的名词特别多(幅频、相频、分贝、对数坐标、奈氏图、包围、穿越、裕度……),我一个都不打算让它裸奔。真正的解题动作只有五步。
1. 先看目标:今天要拿下的这道题
北方工业大学《控制理论基础》阶段考试卷(2) 二、(20分): 已知线性系统开环传递函数为 $\(G(s)=\frac{4}{s(s+1)^2}\)$ (1) 绘制系统的开环幅相曲线,计算并标注与坐标轴相交点的值; (2) 根据绘制的 Nyquist 曲线,判断该反馈控制系统的稳定性。注意增补曲线的绘制。
关于这道题的来源,我要说实话(免得你把它当成考研真题去记):
- 它来自学校的阶段考试卷(期末性质),不是北方工业大学 851 考研真题的原题。文件路径见本文顶部。
- 我为什么挑它:因为它题干是纯文字——给你一个传递函数,图要你自己画。频域这类题很多是"给你一张图,从图上读参数",而我们手上的真题扫描件里,那些图在转成文字的过程中大量丢失了。挑纯文字题,能保证你看到的题是完整的、能真做的。
- 它和考研真题是不是一路货?是。 考点分析里白纸黑字写着 851 的常考大题包含"频域分析:已知开环传递函数绘制 bode/Nyquist 图,判断稳定性"(来源见顶部)。2025 年考研第六大题就是同一个套路,只不过把数字换成了抽象字母(那道题原文在 ../01_结构化/真题/北方工业大学真题__北方工业大学2025年真题-BUFF控制考研_扫描版_p004.md,是"单位负反馈开环传函 \(G(s)=\dfrac{k(\tau s+1)}{s^2(Ts+1)}\),画出 Nyquist 曲线图,用奈氏判据判稳")。先用数字版把套路练熟,抽象版只是把数字换成字母,我在最后一节会点破它多出来的那一层难点。
现在你可能连"幅相曲线"四个字都不认识。正常。看完第 3、4 节回头再看这道题,它会变小。
2. 这类题到底在考什么(一句话)
给你一个开环传递函数,让你画出它的奈奎斯特曲线(幅相曲线),然后判断闭环系统稳不稳定。
这里有四个词你八成还不认识,先各给一句大白话(后面第 3 节再展开):
- 开环传递函数:信号从输入走一圈到反馈回来、但还没接上那个减号时,整条通路的传递函数,常记作 \(G(s)H(s)\) 或 \(G(s)\)。你可以想成"把反馈线剪断,量一量这条链子的加工配方"。
- 闭环系统:把反馈线接回去、真正在工作的那个系统。我们真正关心稳不稳的是闭环。
- 奈奎斯特曲线 / 幅相曲线 / 极坐标图 / Nyquist 图:这四个词在考研里指的是同一个东西(不同教材叫法不同,别被吓到)。它是一条画在平面上的曲线,第 3 节会讲清它每个点是什么意思。
- 判稳:回答"闭环系统稳定还是不稳定"。稳定的意思在第 3 章你已经见过——受到扰动后能自己回到平衡,而不是越晃越大。
【你可能会以为】 第 3 章已经有劳斯判据能判稳、第 4 章根轨迹也能判稳,为什么还要再学一个奈氏判据? 但其实,它们的适用场合和"附赠品"不一样,因为: ① 奈氏判据只需要开环的东西(开环传递函数,甚至只需要一张实验测出来的曲线),不需要去解闭环特征方程的根。工程上系统的数学模型常常不精确,但频率特性可以拿仪器实测出来——给设备输入不同频率的正弦波,量输出,曲线就有了。这是频域法最大的实用价值(PPT《第5章频域分析h-1》开篇讲的第一条优点就是"频率特性具有明确的物理意义,可以用实验的方法……",此处为大白话转述,非逐字)。 ② 劳斯判据只回答"稳/不稳"这个是非题;奈氏判据顺带能告诉你"离不稳还有多远"(就是本模块副线要讲的相角裕度、幅值裕度)。这个"还有多远"是第 6 章校正设计的地基。 ③ 含延迟环节 \(e^{-\tau s}\) 这类系统,劳斯判据根本用不了(它不是多项式),奈氏判据照用不误。
3. 最小基础:只补做这道题真正用得上的东西
第 5 章还有一大堆内容(典型环节的频率特性、Bode 图详细画法、闭环频率特性、谐振峰值、带宽……)。做这道题用不到那些,这一节我一个都不补,用到哪补哪。(这不是"略掉了不重要的",而是这道题确实不需要;Bode 图那一块我在第 7 节副线里会补它需要的最小量。)
3.1 最大的坎:频域分析到底在干嘛
这一小节不讲公式,只讲"世界观"。这是全章最该慢慢看的地方。
前四章你研究系统的方式是:给它一个信号(阶跃、脉冲),看输出随时间怎么跑。 这叫时域分析——"域"就是"看问题的场地","时域"就是"在时间这条轴上看"。
频域分析换了个玩法:我给系统输入一个正弦波,而且是各种不同频率的正弦波,一个一个试过去,看输出变成什么样。
【补课】正弦波是什么、频率是什么 正弦波就是 \(r(t)=A\sin(\omega t)\) 这种上下来回均匀摆动的信号(把它想成秋千、心电图那种规则波浪)。 - \(A\):幅值,波浪的高度(摆多大)。 - \(\omega\)(读作"欧米伽"):角频率,单位是 rad/s(弧度每秒),表示摆得多快。\(\omega\) 大 = 抖得快 = 高频;\(\omega\) 小 = 晃得慢 = 低频。 - \(t\):时间。 你完全不用会推正弦的性质,只要接受"正弦波有两个特征:多高(幅值)、多快(频率)"就够了。
对线性定常系统(就是我们这门课研究的这种系统:线性 = 输入放大 2 倍输出也放大 2 倍;定常 = 系统本身的参数不随时间变),有一条极其漂亮的事实:
【补课】"稳态":刚给信号的头一段时间,输出会乱晃一阵(叫暂态),晃过去之后进入的那个稳定重复的状态,叫稳态。这里说的输出,指的是晃过去以后的那段。
也就是说,输入 \(A\sin(\omega t)\),稳态输出长这样:
- 频率没变:还是 \(\omega\)。系统改不了频率,这是线性系统的特权。
- 高度变了:乘上了一个倍数 \(A(\omega)\)。这个倍数叫 幅频特性——"输出幅值 ÷ 输入幅值"。它是 \(\omega\) 的函数,意思是频率不同,放大倍数不同。
- 时间上错开了:多了一个 \(\varphi(\omega)\)(读作"斐"),叫 相频特性——输出比输入超前或滞后了多少角度。也是 \(\omega\) 的函数。\(\varphi<0\) 表示输出落后于输入(最常见),\(\varphi>0\) 表示超前。
用生活例子把这两个量坐实:想象你推秋千。
- 你慢慢地来回推(低频),秋千跟得上你,幅度大,而且几乎跟你同步 → \(A(\omega)\) 大、\(\varphi(\omega)\) 接近 0。
- 你飞快地来回推(高频),秋千根本来不及反应,只能小幅抖动,而且明显"慢半拍" → \(A(\omega)\) 小、\(\varphi(\omega)\) 负得厉害。
这就是频域分析的全部直觉:同一个系统,对不同快慢的信号,反应的"力度"和"拖延"是不一样的。把每个频率下的"力度"和"拖延"都记下来,就完整刻画了这个系统。
【难点·这个坎难在哪】 难在你会问:"我关心的是阶跃响应,你给我讲正弦干嘛?" 正确理解:任何一个实际信号(包括阶跃、包括各种干扰),在数学上都能拆成一大堆不同频率的正弦波叠加(这叫傅里叶的思想,本模块不展开推导,考研也不考这个推导)。所以只要你知道系统对每个频率的正弦怎么反应,你就等于知道了它对任何信号怎么反应。频域不是另起炉灶,它是把同一个系统换个角度看,信息一点没少。 这一段第一遍看不懂很正常,卡点通常在"为什么正弦这么特殊"。你现在可以先按"这是这门课的一个约定玩法"接受下来,不影响你把题做对;等第 6 章校正时你会自己感觉到它的好用。
3.2 频率特性怎么算:一个替换动作 \(s \to j\omega\)
上面说的 \(A(\omega)\) 和 \(\varphi(\omega)\),怎么从传递函数算出来?规则简单到过分:
这是 PPT《第5章频域分析h-1》给的定义(幻灯片对应文字为"……为系统的频率特性",此处按教材做了大白话展开)。
【补课】\(j\) 是什么 \(j\) 是虚数单位,满足 \(j^2=-1\)。工科用 \(j\),数学书里写 \(i\),是同一个东西(电工怕跟电流 \(i\) 混,才改用 \(j\))。 一个复数长这样:\(z=a+jb\)。\(a\) 叫实部,\(b\) 叫虚部。你可以把它当成平面上的一个点:横坐标 \(a\)、纵坐标 \(b\)。 你现在只需要三条运算规矩: ① \(j\cdot j=-1\);② 复数加减 = 实部跟实部加、虚部跟虚部加;③ 复数乘法照普通乘法展开,遇到 \(j^2\) 就换成 \(-1\)。 想从零补复数,看 预备知识。
一个复数有两个"长相"描述:
- 模(读作 mó,写作 \(|z|\)):这个点离原点多远,\(|a+jb|=\sqrt{a^2+b^2}\)(就是勾股定理)。
- 幅角(写作 \(\angle z\)):从正实轴(横轴正方向)转到这个点的方向,转过多少度。第一象限里 \(\angle z=\arctan\frac{b}{a}\)。
【补课】\(\arctan\) 是什么:\(\tan\theta\) 是"对边比邻边";\(\arctan x\) 就是反过来问"正切值等于 \(x\) 的那个角是多少度"。例如 \(\arctan 1=45^\circ\)、\(\arctan 0=0^\circ\)、\(\arctan\infty=90^\circ\)。考试常用的就这几个特殊值,以及"\(x\) 越大角越大、上限 \(90^\circ\)"这个趋势。25 年起考场禁用计算器,所以真题里的数字一般都凑成特殊角,或者只要求你写表达式。
于是:
这两行正是 PPT《第5章频域分析h-1》幻灯片上写的 \(A(\omega)=|G(j\omega)|\)、\(\varphi(\omega)=\angle G(j\omega)\)(公式一致,文字说明为我补的)。
用大白话再说一遍这两行公式:把 \(s\) 换成 \(j\omega\) 之后,你手里是一个复数(它随 \(\omega\) 变);这个复数的长度就是放大倍数,这个复数的方向角就是相位滞后。 一个复数,同时装下了"多大"和"多拖"两条信息——这就是为什么频率特性非要用复数不可。
一个能省一半计算量的实用公式(考研里天天用):
- 大白话:求模时,分子分母各自的模照乘照除;求角时,分子的角相加、分母的角相减。
- \(N\) 指分子上的因子,\(D\) 指分母上的因子。
- 用处:\(G(j\omega)=\dfrac{4}{j\omega(j\omega+1)^2}\) 这种式子,不用先通分成 \(a+jb\),直接一块一块算模和角,快得多。
3.3 奈奎斯特曲线(幅相曲线)是什么:一条"画出来的档案"
现在你手上有 \(G(j\omega)\):每给一个 \(\omega\),它就是一个复数,也就是平面上的一个点。
那么让 \(\omega\) 从 \(0\) 一直增加到 \(+\infty\),这些点连起来,就画出一条曲线。这条曲线就是奈奎斯特曲线(= 幅相曲线 = 极坐标图 = Nyquist 图)。
这张平面图的坐标:横轴叫 Re(实轴,Real),纵轴叫 Im(虚轴,Imaginary)。曲线上每个点:
- 它离原点的距离 = 该频率下的放大倍数 \(A(\omega)\);
- 它相对于正实轴的方向角 = 该频率下的相位 \(\varphi(\omega)\);
- 曲线上一般要标一个箭头,表示 \(\omega\) 增大的方向。这个箭头非常重要,判稳全靠它。
【你可能会以为】 奈氏图的横轴是频率 \(\omega\)。 但其实不是。 横轴是实部、纵轴是虚部,\(\omega\) 藏在曲线上,它是"沿着曲线走"的那个参数,不占任何一根坐标轴。这是零基础最常见的误解,一旦搞混,整张图就画不对了。 你可以这样理解:\(\omega\) 像时间,曲线像一只蚂蚁爬过的轨迹——纸上看不见时间,但箭头告诉你它往哪边爬。
3.4 为什么盯着 \((-1,j0)\) 这个点
奈氏判据的所有戏,都围着复平面上一个点转:\((-1,\,j0)\),即实部 \(=-1\)、虚部 \(=0\) 的那个点(横轴上 \(-1\) 的位置)。
为什么偏偏是它? 讲透这一层,你就不用死记了:
闭环系统的稳定性,由闭环特征方程 \(1+G(s)H(s)=0\) 的根决定(这条你在第 3、4 章见过:特征方程的根就是闭环极点,根只要有一个跑到复平面右半边,系统就不稳定)。
把 \(1+G(s)H(s)=0\) 挪一下:
而 \(-1\) 这个数,写成复数就是 \(-1+j0\),也就是点 \((-1,j0)\)。
大白话:开环特性一旦等于 \(-1\),意味着"信号绕一圈回来,大小一模一样(模=1)、方向正好反了 \(180^\circ\)",再经过反馈处那个减号一减,就变成了正反馈自我循环——系统卡在临界状态,再多一点就发散。所以 \((-1,j0)\) 是临界点,曲线跟它的相对位置(绕不绕它、从它哪边过)就决定了稳定性。
【难点】 这里最容易糊的是:曲线画的是开环的 \(G(j\omega)H(j\omega)\),判的却是闭环稳不稳。 正确理解:两者的桥梁就是上面那个式子——闭环特征方程 \(1+GH=0\) 里含着开环的 \(GH\)。奈氏判据的全部本事,就是从开环曲线绕 \((-1,j0)\) 的圈数,反推出闭环有几个根在右半平面。所以千万别说成"曲线包围 \((-1,j0)\) 说明开环不稳定"——那是错的,开环稳不稳是另一回事(\(P\) 就是描述它的)。
3.5 奈氏判据:三个符号必须讲清
判据有两种等价写法,考研两种都可能用,我都给你,并且每个符号都交代。
写法一(圈数版,PPT《第5章频域分析h-2》幻灯片18的原意):
- \(Z\):闭环传递函数在 s 右半平面的极点个数。\(Z=0\) 就是闭环稳定;\(Z\) 是几,就有几个不稳定的根。这是我们最终要回答的东西。
- "s 右半平面"= 复平面上实部大于 0 的那一半。根落在那儿意味着响应里有 \(e^{+at}\) 这种越来越大的成分 → 发散 → 不稳定(这一条第 3 章讲过)。
- \(P\):开环传递函数在 s 右半平面的极点个数。这个你直接从开环表达式上数:把分母因式分解,看有几个根实部为正。
- 注意:\(s=0\)(原点上的极点,也就是积分环节)不算右半平面极点,它在虚轴上,单独处理(见 3.6 的"增补")。
- \(R\):奈氏曲线逆时针包围 \((-1,j0)\) 点的圈数。逆时针记正,顺时针记负(PPT 幻灯片18原话即"逆时针包围为正,顺时针包围为负")。
写法二(穿越版,考研做题最常用):
- \(N^+\)(正穿越次数):曲线从下往上穿过实轴上 \((-1,j0)\) 左边那一段(即从 \(-\infty\) 到 \(-1\) 这一段)的次数。
- \(N^-\)(负穿越次数):曲线从上往下穿过那一段的次数。
- \(N\):两者之差。
- 判稳还是看 \(Z=0\) 则闭环稳定。
等一下,"正穿越"到底是从下往上还是从上往下?这里必须小心。 PPT《第5章频域分析h-2》幻灯片27 的文字是:"正穿越:逆时针包含(-1,j0)点时,必然是从上到下穿越,伴随着相角的增大;负穿越:顺时针包含(-1,j0)点时,必然是从下到上穿越,伴随着相角的减小。" 这段 PPT 文字与国内主流教材(胡寿松《自动控制原理》)的表述方向相反——教材是"负穿越:相角减小,曲线自上而下穿过 \(-1\) 左侧负实轴"。我判断这是 PPT 文字在制作时把"上/下"写反了(它自己后面幻灯片36的算例 \(N^-=1,\ N=-1,\ Z=P-2N=2\) 不稳定,与教材口径一致)。 给你的安全做法(记这一条就够,不会错): $\(\boxed{\text{相角"增大"(往逆时针方向转)}=\text{正穿越};\ \text{相角"减小"(往顺时针方向转)}=\text{负穿越}}\)$ 相角增减这个判断永远不会被上下之争坑到,而且和"逆时针包围为正"完全自洽。 这一条我标
【存疑·待核】——存疑的是 PPT 那句话的上下方向,不是判据本身;判据本身(正穿越对应逆时针/相角增大)是确定的。做题时你按上面框里那句走。
两种写法为什么等价? 因为"逆时针绕一整圈"必然要从 \((-1,j0)\) 左侧的负实轴上穿过去两次(上去一次、下来一次)。所以 1 圈 = 2 次净穿越,\(R=2N\)。半闭合曲线(只画 \(\omega:0\to+\infty\) 那一半)时数圈数不方便,数穿越次数才方便——这就是考研为什么爱用 \(Z=P-2N\)。
【补课】"最小相位系统":开环传递函数的所有零点和极点都在 s 左半平面(且不含延迟环节)的系统。对最小相位系统,\(P=0\),判据立刻简化成: $\(Z=-2N,\quad\text{所以}\ \boxed{N=0\ (\text{正负穿越次数相等,含都为0})\Longleftrightarrow \text{闭环稳定}}\)$ 这正是 PPT 幻灯片28 那句"最小相位系统,正负穿越之差=0时,系统稳定"。851 的题绝大多数是最小相位系统,即 \(P=0\),所以实战中你多半只需要判断"曲线有没有绕过/穿过 \((-1,j0)\) 左边那段"。 名字里"最小相位"的由来:同样的幅频特性下,这种系统的相位滞后是所有可能系统里最小的。这个由来考研不考,知道有这么回事即可。
3.6 "增补"是什么(题目里点名要求的那一步)
题目第(2)问原文写着"注意增补曲线的绘制"。这五个字不解释清楚,这道题拿不到满分。
问题出在哪:我们这道题的 \(G(s)=\dfrac{4}{s(s+1)^2}\) 分母有一个 \(s\),也就是在原点 \(s=0\) 处有一个极点(叫积分环节,物理上表示系统里有"累加/积分"作用,比如电机的角度是转速的积分)。
后果:当 \(\omega\to 0^+\) 时,\(|G(j\omega)|=\left|\frac{4}{j\omega(\cdots)}\right|\to\infty\)——曲线的起点跑到无穷远去了,这条曲线是"断"的、不闭合的。而奈氏判据要数"包围圈数",曲线不闭合就没法数圈(PPT 幻灯片19 说的"s 平面的闭合曲线通过原点时,其映射不再封闭,这时判据不成立"就是这个意思)。
解决办法(这就是"增补"):人为在无穷远处补一段圆弧,把曲线接成闭合的。规则是:
- \(\nu\)(读作"纽"):分母里 \(s\) 的幂次,也叫系统型别(\(\nu=0\) 叫 0 型,\(\nu=1\) 叫 I 型,\(\nu=2\) 叫 II 型)。本题 \(G=\frac{4}{s(s+1)^2}\),分母的 \(s\) 是一次方,所以 \(\nu=1\),I 型系统。
- 顺时针补 \(\nu\cdot 90^\circ\):本题就是补 \(90^\circ\)。这段弧在无穷远处,从正虚轴方向(\(+90^\circ\))顺时针转到 \(\omega=0^+\) 处曲线的出发方向(\(-90^\circ\))。
- 这段弧整个待在右半平面(实部为正的一侧),所以它不会穿过 \((-1,j0)\) 左边的负实轴 → 它对 \(N\) 没有贡献。
- PPT《第5章频域分析h-2》幻灯片24/25 写的是"从 \(\omega=0^-\) 到 \(0^+\) 以无穷大半径顺时针补;由于对称性可以补一半……一般只需绘制半闭合开环幅相曲线"。我这里给的是只画 \(\omega:0^+\to+\infty\) 半条曲线时的实用版本(考研标准做法),不是逐字照抄。
【难点】为什么补的是顺时针、不是逆时针? 一句话原因:这段弧是在 s 平面上用一个小半圆绕过原点极点映射过来的,绕行方向定死了它的转向。完整的复变函数论证(幅角原理)PPT 幻灯片10-17 有,考研不会让你证它,只会让你用它。所以这里我明写:这一步的"为什么"我不展开,因为它需要复变函数的幅角原理,属于超纲的推导;你需要掌握的是"\(\nu\) 个原点极点 → 顺时针补 \(\nu\cdot 90^\circ\)"这条操作规则,以及"这段弧在右半平面、不影响穿越计数"这个结论。
4. 套路:五步拿下(记死这五步)
- 第一步 · 把 \(s\) 换成 \(j\omega\),写出 \(G(j\omega)\)。 顺手数两件事:\(P\)(开环右半平面极点数)和 \(\nu\)(原点极点数/型别)。
- 第二步 · 定起点和终点。 算 \(\omega\to 0^+\) 时的模和相角(起点)、\(\omega\to+\infty\) 时的模和相角(终点)。这两头定了,曲线的大致走向就定了。
- 第三步 · 求与实轴、虚轴的交点。 令虚部 \(=0\) 解出 \(\omega\),回代求实部 → 与实轴的交点(这是本类题的核心得分点,题目明确要求"计算并标注")。同理令实部 \(=0\) 求与虚轴的交点(有些题不要求)。
- 第四步 · 画图 + 增补。 把起点、终点、交点标上,连成曲线,标 \(\omega\) 增大的箭头;若 \(\nu\ge 1\),在无穷远处顺时针补 \(\nu\cdot 90^\circ\) 圆弧。
- 第五步 · 用判据算 \(Z\),下结论。 数穿越次数 \(N\),代入 \(Z=P-2N\);\(Z=0\) 稳定,否则不稳定并说明有几个右半平面闭环极点。
每一步为什么这么做,下面例题里边做边说。
5. 例题手把手:就用开场那道 20 分大题
已知开环传递函数 \(G(s)=\dfrac{4}{s(s+1)^2}\)。(1) 画开环幅相曲线,标出与坐标轴交点;(2) 判断闭环稳定性,注意增补。
第一步:换 \(s\to j\omega\),数 \(P\) 和 \(\nu\)
- 数 \(P\):开环极点就是分母 \(s(s+1)^2\) 的根,即 \(s=0\)(一个)、\(s=-1\)(两重)。\(s=-1\) 在左半平面,\(s=0\) 在虚轴上——右半平面一个都没有,所以 \(\boxed{P=0}\)。(顺带:零点没有,极点都不在右半平面,所以这是个最小相位系统。)
- 数 \(\nu\):分母含 \(s^1\),所以 \(\boxed{\nu=1}\),I 型系统,后面必须增补 \(90^\circ\)。
为什么先做这步:\(P\) 是判据公式里的输入量,\(\nu\) 决定要不要增补。这两个数不先定死,后面全是白算。
第二步:定起点(\(\omega\to0^+\))和终点(\(\omega\to+\infty\))
用 3.2 那条"模照乘除、角照加减"的规矩,先把幅频、相频写出来:
- 拆解给你看:\(|j\omega|=\sqrt{0^2+\omega^2}=\omega\)(取 \(\omega>0\));\(|j\omega+1|=|1+j\omega|=\sqrt{1^2+\omega^2}\),因为是平方所以整体再平方。
- 拆解:分子是正实数 4,方向就是正实轴,角度 \(0^\circ\);\(j\omega\) 是纯虚数、在正虚轴上,角度 \(90^\circ\),在分母所以减掉;\(1+j\omega\) 的角是 \(\arctan\frac{\omega}{1}=\arctan\omega\),平方就是角度乘 2,也在分母所以减掉。
- (这一行 \(\varphi(\omega)=-90^\circ-2\arctan\omega\) 和原卷手写解答里的 "\(\varphi(\omega)=-90^\circ-2\tan^{-1}\omega\)" 一致;\(\tan^{-1}\) 是 \(\arctan\) 的另一种写法,同一个东西。)
现在代两头:
- 起点 \(\omega\to0^+\):\(A\to\dfrac{4}{0\cdot 1}\to+\infty\);\(\varphi\to-90^\circ-2\arctan 0=-90^\circ\)。 → 起点在无穷远处,方向是 \(-90^\circ\)(正下方,负虚轴方向)。
- 终点 \(\omega\to+\infty\):\(A\to\dfrac{4}{\infty}\to 0\);\(\varphi\to-90^\circ-2\times 90^\circ=-270^\circ\)。 → 终点回到原点,进入方向是 \(-270^\circ\)(等价于 \(+90^\circ\),即从正上方切进原点)。
(原卷手写解答写的"起于 \(-90^\circ\),终于 \(-270^\circ\)"就是这个意思,我们自己算出来了。)
为什么先定两头:奈氏图不要求你逐点描 100 个点,考研只要定性画对形状。而形状是由"从哪来、到哪去、中间穿没穿轴"决定的。两头一定,你就知道曲线要从第三/四象限的无穷远,一路扫到原点。
第三步:求与实轴的交点(本题的核心得分点)
"与实轴相交"的数学含义:虚部 \(=0\)。所以要先把 \(G(j\omega)\) 化成 \(a+jb\) 的形式。
先算分母: $\(j\omega(j\omega+1)^2\)$ 先展开 \((j\omega+1)^2=(j\omega)^2+2j\omega+1=-\omega^2+2j\omega+1=(1-\omega^2)+2j\omega\) (用了 \(j^2=-1\)。)
再乘 \(j\omega\): $\(j\omega\big[(1-\omega^2)+2j\omega\big]=j\omega(1-\omega^2)+2j^2\omega^2=-2\omega^2+j\,\omega(1-\omega^2)\)$
所以 $\(G(j\omega)=\frac{4}{-2\omega^2+j\,\omega(1-\omega^2)}\)$
分母有 \(j\),不能直接看虚部,要"有理化"——分子分母同乘分母的共轭(共轭就是把虚部的符号取反:\(a+jb\) 的共轭是 \(a-jb\);这么做的好处是 \((a+jb)(a-jb)=a^2+b^2\) 变成实数):
分母现在是个正实数,记作 \(M=4\omega^4+\omega^2(1-\omega^2)^2\)。于是:
令虚部 \(=0\):\(-4\omega(1-\omega^2)=0\) → \(\omega=0\) 或 \(\omega^2=1\) 即 \(\omega=1\)(\(\omega\) 取正)。\(\omega=0\) 是起点(在无穷远,不算交点),所以真正的交点在
回代求实部(直接用原式更快):\(\omega=1\) 时分母 \(=-2(1)^2+j\cdot 1\cdot(1-1)=-2\),所以
(这与原卷手写的"与负实轴交点 \((-2,0j)\)"一致。)
与虚轴的交点:令实部 \(\dfrac{-8\omega^2}{M}=0\) → 只有 \(\omega=0\),而 \(\omega=0\) 时曲线在无穷远。所以除了起点跑向无穷远之外,曲线不与虚轴相交(答题时这么写一句即可,别空着)。
为什么这步最重要:因为判稳的全部信息就是"曲线跟 \((-1,j0)\) 的相对位置"。而曲线跟负实轴的交点在 \(-2\) —— \(-2\) 比 \(-1\) 更靠左,这就意味着曲线是从 \((-1,j0)\) 的外侧扫过去的,它把这个点圈进去了。答案基本已经出来了。
第四步:画图 + 增补
考场上你要画一张草图。我用文字把它精确描述给你(这里不引用任何图片文件,因为我们手上真题扫描件里的这张图不保证清晰完整;文字描述足以让你自己画出来):
Im (虚轴)
|
| ↖ 终点:ω→∞ 时曲线从上方切入原点
| (方向 -270°,即正上方)
----(-2,j0)-------+------------------- Re (实轴)
● (-1,j0) 原点
| ×
| ↑ 临界点
曲线在 ω=1 处穿过负实轴于 -2
|
| ω 增大方向:从下往上(相角从 -90° 单调减到 -270°,
| 即顺时针转;曲线由第三象限下方无穷远处上来)
↓
起点:ω→0+ 时在无穷远,方向 -90°(正下方)
用语言把形状说清楚(考场按这个画):
- 从下方无穷远处出发(\(\omega\to0^+\),相角 \(-90^\circ\),即负虚轴方向的无穷远);
- 随着 \(\omega\) 增大,相角从 \(-90^\circ\) 单调地减小到 \(-270^\circ\)(因为 \(\varphi=-90^\circ-2\arctan\omega\),\(\arctan\omega\) 随 \(\omega\) 单调增),模从 \(\infty\) 单调减到 \(0\);
- 途中在 \(\omega=1\) 处穿过负实轴于 \(-2\) 点;
- 最后从正上方(\(-270^\circ\) 方向)切入原点;
- 增补:\(\nu=1\),在无穷远处从 \(\omega=0^+\) 那一点(方向 \(-90^\circ\))逆着 \(\omega\) 方向往回、以无穷大半径顺时针补 \(90^\circ\) 圆弧,即这段弧扫过 \(-90^\circ\to 0^\circ\to +90^\circ\) 这个范围(整段都在右半平面),把曲线接闭合。
【难点·画图时最容易错的地方】 很多人把这条曲线画成"从右边过来",或者忘了标箭头。 记住两条保命规则:① 起点方向 \(=-\nu\cdot 90^\circ\)(本题 \(\nu=1\),起点 \(-90^\circ\));② 终点方向 \(=-(n-m)\times 90^\circ\),其中 \(n\) 是分母最高次、\(m\) 是分子最高次(本题分母 \(s(s+1)^2\) 展开最高次 \(n=3\),分子是常数 \(m=0\),所以终点 \(=-270^\circ\),和我们算的一致)。这两条规则能让你不算极限也能秒定两头。
第五步:判据算 \(Z\),下结论
数穿越次数:我们关心的是曲线穿过实轴上 \((-\infty,-1)\) 这一段(即 \((-1,j0)\) 左边)的情况。
- 唯一的穿越发生在 \(\omega=1\)、位置 \(-2\) —— \(-2\) 确实在 \(-1\) 的左边,算数。
- 穿越方向:相角 \(\varphi\) 从 \(-90^\circ\) 减小到 \(-270^\circ\),是减小的 → 按 3.5 那个框:相角减小 = 负穿越。
- (对照几何图像:曲线此处是从下半平面走到上半平面,也就是"自下而上"穿过 \(-1\) 左侧的负实轴,这与教材"负穿越"的图像一致。)
- 增补的那段无穷大圆弧全在右半平面,不穿越,不计。
所以: $\(N^+=0,\quad N^-=1,\quad N=N^+-N^-=-1\)$
代入判据: $\(Z=P-2N=0-2\times(-1)=2\)$
结论:\(Z=2\ne 0\),闭环系统不稳定,并且有 2 个闭环极点位于 s 右半平面。
(用另一种写法验算:曲线顺时针包围 \((-1,j0)\) 一圈,所以逆时针圈数 \(R=-1\),\(Z=P-R=0-(-1)\cdot 2\)……注意这里要小心:\(Z=P-R\) 用的是完整曲线(\(\omega:-\infty\to+\infty\))的圈数,完整曲线顺时针绕 1 整圈对应 \(R=-2\)?不要在这里绕。考研统一用 \(Z=P-2N\) 这一套,把半闭合曲线的穿越数数清楚即可,不要两套混着用——混用是这类题最大的翻车点。原卷手写解答那行 "\(Z=P-R=0-2(0-1)=2\)" 其实写的就是 \(Z=P-2N\) 的算法,结果 \(Z=2\) 与我们一致。)
这道 20 分的大题做完了。 答案: - (1) 幅相曲线起于无穷远(\(-90^\circ\)方向)、终于原点(\(-270^\circ\)方向),在 \(\omega=1\) 处与负实轴交于 \((-2,j0)\),与虚轴无交点;\(\nu=1\) 需在无穷远处顺时针增补 \(90^\circ\) 圆弧。 - (2) \(P=0\),\(N=-1\),\(Z=P-2N=2\),闭环不稳定,有 2 个右半平面闭环极点。
停一下,回头看第 1 节那道题。 你刚刚把它做出来了。
6. 你来做:变式题(带详细答案)
练习题(自编,仿真题难度——特此说明:下面这道题是我按同等难度自己编的,不是任何一年的真题原题,别当真题记):
已知单位负反馈系统的开环传递函数为 $\(G(s)=\frac{1}{s(s+1)^2}\)$ (1) 求幅频特性 \(A(\omega)\) 和相频特性 \(\varphi(\omega)\),定出幅相曲线的起点和终点; (2) 求曲线与负实轴的交点; (3) 用奈氏判据判断闭环系统的稳定性。 (4) 进阶(想挑战再做):若把分子的 1 换成 \(K>0\),即 \(G(s)=\dfrac{K}{s(s+1)^2}\),求使闭环系统稳定的 \(K\) 取值范围。
注意:这道题和例题只差一个系数(4 变成 1)。这不是偷懒,是我故意的——你会看到结论翻了个个儿,从而亲眼看见"增益 \(K\) 怎么把一个系统从稳定推向不稳定"。这正是第 6 章校正的全部动机。
自己先做。卡住了对照回看:
- 不会把 \(s\) 换成 \(j\omega\)、不会算模和角? → 回看 3.2 节
- 不知道起点终点怎么定? → 回看 第 5 节第二步,或用那两条"保命规则"(起点 \(-\nu\cdot90^\circ\)、终点 \(-(n-m)\cdot90^\circ\))
- 不会求交点(通分、共轭那一步)? → 回看 第 5 节第三步
- 不会用判据、不知道 \(P\) 怎么数? → 回看 3.5 节