基础 1

数学与电路基础

学完这一章你应该能

  • 熟练地在代数形式和极坐标形式之间转换复数,知道「乘法 = 模相乘、辐角相加」;
  • 看到 \(e^{-t/T}\)、\(e^{-\sigma t}\sin\omega t\) 就能想象出曲线的样子;
  • 求二次方程的复数根,知道「方程的根决定运动模式」;
  • 解简单的一阶线性微分方程;
  • 会用分贝和对数坐标;
  • 会用一阶泰勒展开做线性近似;
  • 用「虚短、虚断」推导运放电路。

前置知识:高中数学,以及一点微积分(会求导和简单积分)。这一章是后面所有章节的工具,用到哪里会回头链接到这里。

先测一测

下面 5 道题都能很快做对的话,可以直接跳到拉普拉斯变换。做错了哪题,就去看对应的小节。

1. 把 \(z=-1-j\) 写成极坐标形式。(→ 1 节)

\(|z|=\sqrt2\approx1.414\),辐角在第三象限:\(\angle z=-135^\circ\)。所以 \(z=1.414\angle-135^\circ\)。如果你写的是 45°,说明掉进了「只算 arctan(b/a)」的坑。

2. \(e^{-t/2}\) 衰减到初值的 5% 大约需要多久?(→ 2 节)

时间常数 \(T=2\) s,衰减到 5% 大约要 \(3T=6\) s。

3. 求 \(s^2+4s+13=0\) 的根。(→ 3 节)

\(s=\dfrac{-4\pm\sqrt{16-52}}{2}=-2\pm j3\)。

4. 幅值放大 100 倍是多少 dB?缩小到 1/2 呢?(→ 5 节)

\(20\lg100=40\) dB;\(20\lg0.5\approx-6\) dB。

5. 理想运放反相放大器,输入电阻 10 kΩ,反馈电阻 100 kΩ,放大倍数是多少?(→ 7 节)

\(-R_f/R_i=-10\),即反相放大 10 倍。

1. 复数与复平面

控制理论里到处是复数:拉氏变换的变量 \(s=\sigma+j\omega\) 是复数,系统的极点、零点是复数,频率特性 \(G(j\omega)\) 也是复数。这一节是全书用得最多的数学工具。

代数形式与复平面

虚数单位 \(j\) 满足 \(j^2=-1\)(数学课上写作 \(i\),工程里用 \(j\),避免和电流 \(i\) 混淆)。复数写成

$$z=a+jb,\qquad a=\operatorname{Re}z\ (\text{实部}),\quad b=\operatorname{Im}z\ (\text{虚部})$$

把实部当横坐标、虚部当纵坐标,每个复数就对应平面上的一个点,也可以看成从原点指向该点的向量。这个平面叫复平面,控制里叫 s 平面。

实轴 虚轴 z = a + jb ab r = |z| θ σ jω p s s − p ∠(s−p)

模、辐角与极坐标形式

易错点:辐角要看象限。\(\arctan(b/a)\) 只能给出 −90°~90° 之间的角。\(z=-1-j\) 和 \(z=1+j\) 的 \(b/a\) 都是 1,但前者在第三象限,辐角是 −135°(或 225°),不是 45°。编程时用 atan2(b, a),它会自动处理象限。

欧拉公式与指数形式

$$e^{j\theta}=\cos\theta+j\sin\theta\qquad\Longrightarrow\qquad z=re^{j\theta}$$

由它还可以得到 \(\cos\theta=\dfrac{e^{j\theta}+e^{-j\theta}}{2}\),\(\sin\theta=\dfrac{e^{j\theta}-e^{-j\theta}}{2j}\),拉氏变换那一章会用到。

四则运算

「乘法 = 模相乘、辐角相加」是整本笔记最重要的数学事实之一。频率特性 \(G(j\omega)=\dfrac{K}{(j\omega-p_1)(j\omega-p_2)}\) 的模等于 \(K\) 除以各因子的模,相角等于各因子相角之和取负。波特图能「叠加着画」、根轨迹的相角条件,都来自这一条。

例 1 计算 \(\dfrac{-1+j2}{4j}\)(拉氏变换那一章会遇到)。

分子:\(|-1+j2|=\sqrt5=2.236\),在第二象限,辐角 \(180^\circ-\arctan2=116.57^\circ\)。分母:\(4j=4\angle90^\circ\)。所以

$$\frac{-1+j2}{4j}=\frac{2.236}{4}\angle(116.57^\circ-90^\circ)=0.559\angle26.57^\circ$$

例 2 \(s\) 和 \(p\) 都是复数时,\(s-p\) 就是从 \(p\) 指向 \(s\) 的向量(上图右):它的模是两点之间的距离,辐角是这个向量相对正实轴方向的角度。第4章的根轨迹就是在 s 平面上量这些向量的长度和角度。

2. 指数函数、正弦与时间常数

指数函数:衰减还是发散

\(e^{at}\):\(a<0\) 时衰减到 0,\(a>0\) 时越来越大,\(a=0\) 时保持不变。控制里常写成 \(e^{-t/T}\),\(T>0\) 叫时间常数。

\(t\)\(T\)\(2T\)\(3T\)\(4T\)\(5T\)
\(e^{-t/T}\)(剩下多少)36.8%13.5%5.0%1.8%0.7%
\(1-e^{-t/T}\)(完成多少)63.2%86.5%95.0%98.2%99.3%

经验:过了 3~4 个时间常数,指数过程就基本结束了。这就是一阶系统调节时间 \(t_s=3T\sim4T\) 的来源。

正弦量

\(A\sin(\omega t+\varphi)\):\(A\) 为幅值,\(\omega\) 为角频率(rad/s),\(\varphi\) 为初相位。角频率和频率 \(f\)(Hz)、周期 \(T_0\) 的关系:\(\omega=2\pi f=2\pi/T_0\)。\(\varphi>0\) 叫超前,\(\varphi<0\) 叫滞后。

合并公式(拉氏反变换时常用):

$$a\cos\omega t+b\sin\omega t=R\cos(\omega t-\phi),\qquad R=\sqrt{a^2+b^2},\quad \phi=\operatorname{atan2}(b,a)$$

例:\(\cos2t-0.5\sin2t\):\(R=\sqrt{1+0.25}=1.118\),\(\phi=\operatorname{atan2}(-0.5,1)=-26.57^\circ\),所以等于 \(1.118\cos(2t+26.57^\circ)\)。

衰减振荡

\(e^{-\sigma t}\sin\omega t\)(\(\sigma>0\))是一个振幅按 \(e^{-\sigma t}\) 缩小的正弦波:\(\pm e^{-\sigma t}\) 叫它的包络线。\(\sigma\) 决定衰减快慢,\(\omega\) 决定振荡快慢。后面会看到,它对应 s 平面上的一对复数 \(-\sigma\pm j\omega\):实部管衰减,虚部管振荡。

3. 多项式与方程的根

二次方程与复数根

$$as^2+bs+c=0\quad\Longrightarrow\quad s=\frac{-b\pm\sqrt{b^2-4ac}}{2a}$$

判别式 \(\Delta=b^2-4ac<0\) 时,根是一对共轭复数:\(s=-\dfrac{b}{2a}\pm j\dfrac{\sqrt{4ac-b^2}}{2a}\)。例:\(s^2+2s+5=0\Rightarrow s=-1\pm j2\)。

实系数多项式的复数根一定成对共轭出现,所以 s 平面上的极点图总是关于实轴对称。

根与因式

如果 \(s_1,\ldots,s_n\) 是 \(a_ns^n+\cdots+a_0=0\) 的根,那么多项式可以写成 \(a_n(s-s_1)(s-s_2)\cdots(s-s_n)\)。反过来,展开因式就得到系数,比如

$$(s+1)(s+3)=s^2+4s+3,\qquad s(s+1)(s+3)=s^3+4s^2+3s$$

根与系数的关系(韦达定理):对首项系数为 1 的多项式 \(s^n+a_{n-1}s^{n-1}+\cdots+a_0\),所有根之和 \(=-a_{n-1}\),所有根之积 \(=(-1)^na_0\)。第4章根轨迹的「根之和不变」就是这条。

三次以上的方程手算求根很难。控制理论的一大半内容(劳斯判据、根轨迹、频域法)都在回答:不解方程,怎样判断根在哪里。实际工作中求根直接用计算机,比如 Python 的 numpy.roots([1, 4, 3, 12])。

4. 微分方程初步

含有未知函数的导数的方程叫微分方程。控制里自变量是时间 \(t\),常把 \(\dfrac{dy}{dt}\) 写成 \(\dot y\),\(\dfrac{d^2y}{dt^2}\) 写成 \(\ddot y\)。本课程处理的几乎都是线性常系数微分方程:

$$a_n y^{(n)}+\cdots+a_1\dot y+a_0y=(\text{输入及其导数})$$

解的结构:全解 = 齐次方程的通解(自由运动,由系统本身和初始条件决定)+ 一个特解(强迫运动,由输入决定)。

一阶例子:RC 充电

\(RC\dot u+u=U\),\(u(0)=0\)(\(U\) 为常数)。

  1. 齐次方程 \(RC\dot u+u=0\) 的通解:\(u_h=Ce^{-t/RC}\)(代回去验证一下);
  2. 特解:输入是常数,猜 \(u_p=\) 常数,代入得 \(u_p=U\);
  3. 全解 \(u=U+Ce^{-t/RC}\),由 \(u(0)=0\) 得 \(C=-U\):
    $$u(t)=U\left(1-e^{-t/RC}\right)$$

电容电压按时间常数 \(T=RC\) 指数上升,这就是第3章「一阶系统」的原型。

二阶齐次方程:特征方程

对 \(\ddot y+a_1\dot y+a_0y=0\),设 \(y=e^{\lambda t}\) 代入,得到特征方程 \(\lambda^2+a_1\lambda+a_0=0\)。按根的情况:

特征根通解运动形式
两个不同实根 \(\lambda_1,\lambda_2\)\(C_1e^{\lambda_1t}+C_2e^{\lambda_2t}\)两个指数的叠加
重根 \(\lambda\)\((C_1+C_2t)e^{\lambda t}\)指数,带一个 \(t\) 因子
共轭复根 \(\sigma\pm j\omega\)\(e^{\sigma t}(C_1\cos\omega t+C_2\sin\omega t)\)振荡,实部决定衰减或发散

例:\(\ddot y+3\dot y+2y=0\),\(\lambda=-1,-2\),\(y=C_1e^{-t}+C_2e^{-2t}\);\(\ddot y+2\dot y+5y=0\),\(\lambda=-1\pm j2\),\(y=e^{-t}(C_1\cos2t+C_2\sin2t)\)。

特征根决定了系统的固有运动模式:实部为负就衰减,实部为正就发散,有虚部就振荡。后面的「极点」就是这些特征根,「稳定性」就是问「所有特征根的实部是否都为负」。拉氏变换会给出一套更系统、能带上初始条件和任意输入的解法。

5. 对数与分贝

\(\lg x=\log_{10}x\)。基本性质:\(\lg(ab)=\lg a+\lg b\),\(\lg(a/b)=\lg a-\lg b\),\(\lg a^n=n\lg a\)。几个要记的值:\(\lg1=0\),\(\lg10=1\),\(\lg2\approx0.301\)。

分贝:幅值比 \(A\) 用分贝表示为 \(L=20\lg A\) dB。

幅值比 \(A\)0.10.50.707123.1610100
分贝−20−6−306102040

6. 泰勒展开与线性近似

在点 \(x_0\) 附近,光滑函数可以用它的切线近似:

$$f(x)\approx f(x_0)+f'(x_0)(x-x_0)\qquad\text{即}\qquad \Delta y\approx f'(x_0)\,\Delta x$$

例:\(\sin\theta\approx\theta\)(\(|\theta|<14^\circ\) 时误差小于 1%);\(e^x\approx1+x\)。

控制里的用途:真实系统多少都有非线性,只要它在某个工作点附近小范围变化,就可以这样线性化,然后用线性理论分析(第2章)。例:水箱出水流量 \(Q=k\sqrt H\),在液位 \(H_0\) 附近

$$\Delta Q\approx\frac{k}{2\sqrt{H_0}}\,\Delta H$$

系数 \(\dfrac{k}{2\sqrt{H_0}}\) 随工作点变化:液位越高,同样的液位变化引起的流量变化越小。

两个变量时:\(f(x,y)\approx f(x_0,y_0)+\dfrac{\partial f}{\partial x}\Delta x+\dfrac{\partial f}{\partial y}\Delta y\)。

7. 电路基础:R、L、C 与运放

三种元件

元件电压—电流关系特点复阻抗(第2章用)
电阻 \(R\)\(u=Ri\)没有记忆\(R\)
电容 \(C\)\(i=C\dfrac{du}{dt}\)电压不能突变,储存电场能\(\dfrac{1}{Cs}\)
电感 \(L\)\(u=L\dfrac{di}{dt}\)电流不能突变,储存磁场能\(Ls\)

两条定律:KCL,流入一个节点的电流之和为零;KVL,沿任一回路的电压之和为零。两个电阻串联分压:\(u_{\text{out}}=\dfrac{R_2}{R_1+R_2}u_{\text{in}}\)。

理想运放的两条规则

运放接成负反馈时,可以用两条规则分析:

uᵢ Rᵢ A(虚地) Rf −+ uₒ i →

推导反相放大器:同相端接地,由虚短,反相端 A 点电压也为 0(叫虚地)。流过 \(R_i\) 的电流 \(i=u_i/R_i\);由虚断,这个电流全部流过 \(R_f\),所以 \(u_o=0-iR_f\):

$$u_o=-\frac{R_f}{R_i}u_i\qquad\text{推广到阻抗:}\quad \frac{U_o(s)}{U_i(s)}=-\frac{Z_f(s)}{Z_i(s)}$$

把反馈电阻换成电容(\(Z_f=\frac{1}{Cs}\)),就得到积分器 \(-\dfrac{1}{RCs}\),即 \(u_o=-\dfrac{1}{RC}\int u_i\,dt\)。换不同的元件可以搭出各种典型环节,实验课的模拟实验箱就是这个原理(第2章 2.3 节)。

本章小结

  • 复数:\(z=a+jb=r\angle\theta=re^{j\theta}\);乘法模相乘、辐角相加;辐角要看象限;\(s-p\) 是从 \(p\) 指向 \(s\) 的向量。
  • 指数:过 3~4 个时间常数基本结束;\(e^{-\sigma t}\sin\omega t\) 的实部管衰减,虚部管振荡。
  • 方程的根决定运动模式;实系数方程的复根成对共轭;根之和等于 \(-a_{n-1}\)。
  • 线性微分方程的解 = 自由运动 + 强迫运动;特征根就是后面的极点。
  • 分贝 \(20\lg A\):×10 是 +20 dB,×2 约 +6 dB,×0.707 约 −3 dB;−20 dB/dec 表示 \(A\propto1/\omega\)。
  • 线性近似 \(\Delta y\approx f'(x_0)\Delta x\),系数和工作点有关。
  • 运放负反馈:虚短、虚断;反相电路 \(-Z_f/Z_i\)。

练习

1. 把 \(z=-2+j2\) 写成极坐标形式。

\(|z|=2\sqrt2\approx2.83\),第二象限,\(\angle z=135^\circ\),所以 \(z=2.83\angle135^\circ\)。

2. 计算 \(\dfrac{1+j}{1-j}\)。

极坐标:\(\dfrac{\sqrt2\angle45^\circ}{\sqrt2\angle-45^\circ}=1\angle90^\circ=j\)。或者上下同乘 \(1+j\):\(\dfrac{(1+j)^2}{2}=\dfrac{2j}{2}=j\)。

3. 设 \(G(j\omega)=\dfrac{1}{j\omega(j\omega+1)}\),求 \(\omega=1\) 时的模和相角。

\(|j|=1\),\(|j+1|=\sqrt2\),所以 \(|G|=\dfrac{1}{\sqrt2}\approx0.707\)。相角 \(=-(90^\circ+45^\circ)=-135^\circ\)。这就是第5章算频率特性的方法。

4. 解 \(\dot y+2y=4\),\(y(0)=0\)。

齐次通解 \(Ce^{-2t}\),特解 \(y_p=2\),由初始条件 \(C=-2\):\(y(t)=2(1-e^{-2t})\)。时间常数 0.5 s,约 1.5~2 s 后基本稳定在 2。

5. 方程 \(\ddot y+4\dot y+4y=0\) 的通解是什么?它会振荡吗?

特征方程 \(\lambda^2+4\lambda+4=0\),重根 \(\lambda=-2\),\(y=(C_1+C_2t)e^{-2t}\)。没有虚部,不振荡;这是第3章说的「临界阻尼」。

6. 一个环节在 ω=1 rad/s 时幅值为 20 dB,幅频斜率为 −20 dB/dec,ω=100 rad/s 时幅值是多少?

从 1 到 100 是两个十倍频程,下降 40 dB,所以是 −20 dB,即幅值比 0.1。

7. 在 \(x_0=1\) 附近线性化 \(y=x^3\)。

\(f'(1)=3\),\(\Delta y\approx3\Delta x\),即 \(y\approx1+3(x-1)\)。在 \(x=1.1\) 时近似值 1.3,精确值 1.331。

8. 运放反相电路,输入支路是电阻 \(R\),反馈支路是电阻 \(R_f\) 与电容 \(C\) 并联,传递函数是什么?

\(Z_f=R_f\parallel\dfrac{1}{Cs}=\dfrac{R_f}{R_fCs+1}\),所以 \(\dfrac{U_o}{U_i}=-\dfrac{R_f/R}{R_fCs+1}\):增益 \(R_f/R\)、时间常数 \(R_fC\) 的惯性环节(带负号)。实验一里的一阶惯性环节就是这样搭的。