控制工程中波特图法综合调节器的方程组求解问题
嗨，我来帮你拆解这个方程组的求解过程，其实核心是用三角函数的差角公式简化计算，一步步来就很清晰了：
首先，咱们利用第一个方程的比值关系减少变量数量。已知 $\frac{\omega_{lb1}}{\omega_{ln1}} = 12.6$，咱们令 $k = 12.6$，那么可以把 $\omega_{lb1}$ 表示成 $\omega_{lb1} = k \cdot \omega_{ln1}$，这样第二个方程里就只剩一个未知量了。
接下来看第二个方程：$\arctan \left(\frac{2}{\omega_{lb1}} \right) - \arctan \left( \frac{2}{\omega_{ln1}} \right) = -6^\circ$，把 $\omega_{lb1} = k \cdot \omega_{ln1}$ 代入进去，就变成：
$$\arctan \left(\frac{2}{k \cdot \omega_{ln1}} \right) - \arctan \left( \frac{2}{\omega_{ln1}} \right) = -6^\circ$$
为了更简洁，咱们再设 $x = \frac{2}{\omega_{ln1}}$，那么上式可以改写为：
$$\arctan\left(\frac{x}{k}\right) - \arctan(x) = -6^\circ$$
这里要用到反正切的差角公式：$\arctan(A) - \arctan(B) = \arctan\left( \frac{A - B}{1 + A \cdot B} \right)$，不过要注意结果是负角度，所以两边取正切的话，$\tan(-6^\circ) = -\tan(6^\circ)$。
把 $A = \frac{x}{k}$、$B = x$ 代入公式，计算分子分母：
- 分子：$\frac{x}{k} - x = x \cdot \frac{1 - k}{k}$
- 分母：$1 + \frac{x}{k} \cdot x = \frac{k + x^2}{k}$
分子除以分母后，$k$ 会约掉，得到：
$$\frac{x(1 - k)}{k + x^2} = -\tan(6^\circ)$$
代入 $k=12.6$，$1 - k = -11.6$，两边的负号可以抵消，再把 $\tan(6^\circ) \approx 0.1051$ 代入，就得到：
$$\frac{11.6x}{12.6 + x^2} = 0.1051$$
整理成标准的二次方程形式：
$$0.1051x^2 - 11.6x + 0.1051 \times 12.6 = 0$$
计算常数项：$0.1051 \times 12.6 \approx 1.324$，所以方程变为：
$$0.1051x^2 - 11.6x + 1.324 = 0$$
接下来用二次方程求根公式 $x = \frac{-b \pm \sqrt{b^2 - 4ac}}{2a}$（这里 $a=0.1051$，$b=-11.6$，$c=1.324$），计算判别式：
$$\sqrt{(-11.6)^2 - 4 \times 0.1051 \times 1.324} \approx \sqrt{134.56 - 0.557} \approx 11.576$$
得到两个根：
- $x_1 = \frac{11.6 + 11.576}{2 \times 0.1051} \approx 110.26$
- $x_2 = \frac{11.6 - 11.576}{2 \times 0.1051} \approx 0.114$
现在要筛选合理的根：因为咱们是在做波特图的调节器综合（滞后/超前补偿），这些频率都是低频段的数值，所以 $x = \frac{2}{\omega_{ln1}}$ 不能太小（否则 $\omega_{ln1}$ 会太大，不符合低频设计）。所以取 $x_1$：
- $\omega_{ln1} = \frac{2}{x_1} \approx \frac{2}{110.26} \approx 0.0181 \text{ rad/sec}$
- $\omega_{lb1} = 12.6 \times \omega_{ln1} \approx 12.6 \times 0.0181 \approx 0.228 \text{ rad/sec}$
题目里给的解（0.016和0.2）是近似值，应该是计算时对 $\tan(6^\circ)$ 或中间步骤做了更粗略的四舍五入，本质上是一致的。
备注：内容来源于stack exchange，提问作者FlatAssembler
