串联控制环路仿真解析
在环路设计时，我们经常使用串联控制环路，如下图所示：
需要调试内环、外环，甚至2个以上的环路。机器在工作过程中同时采集和控制多个不同的物理量，比如电压、电流、温度、频率等等，而这些物理量又常常相互影响和制约。而串联控制环路是一种简单有效的处理方法，它的工作原理是将外环的输出作为内环的输入，内环的输出作为整个系统的输出，从而使所有被控物理量工作在稳定状态。
内外环路的相互影响
一个双环的控制系统，等效于将内环的闭环传递函数串联在外环的开环传递函数上，那么必然在增益和相位上对外环造成影响。假设内环和外环的开环传递函数为：
\[\begin{align}
C_{inner}(s) &= \frac{s+100}{s} \\
C_{outer}(s) &= \frac{s+500}{s}
\end{align}
\]
在 Python 中，我们对这两个环路以及总的传递函数进行 Bode 图绘制：
import numpy as np
import matplotlib.pyplot as plt
from scipy import signal
def bode_plot(tf_list, f_min, f_max, ax1=None, ax2=None, colors=None, linestyles=None, labels=None, title=""):
"""
生成伯德图 - 支持多个传递函数
参数:
tf_list: 传输函数列表或单个传输函数
f_min: 最小频率 (Hz)
f_max: 最大频率 (Hz)
ax1: 幅度图轴对象 (可选)
ax2: 相位图轴对象 (可选)
colors: 颜色列表 (可选)
linestyles: 线条样式列表 (可选)
labels: 标签列表 (可选)
title: 图形标题 (可选)
"""
# 如果输入是单个传递函数，转换为列表
if not isinstance(tf_list, (list, tuple)):
tf_list = [tf_list]
# 设置默认值
if colors is None:
colors = plt.cm.tab10(np.linspace(0, 1, len(tf_list)))
if linestyles is None:
linestyles = ['-'] * len(tf_list)
if labels is None:
labels = [f'TF {i+1}' for i in range(len(tf_list))]
# 确保列表长度一致
if len(colors) < len(tf_list):
colors = colors * (len(tf_list) // len(colors) + 1)
if len(linestyles) < len(tf_list):
linestyles = linestyles * (len(tf_list) // len(linestyles) + 1)
if len(labels) < len(tf_list):
labels = labels + [f'TF {i+1}' for i in range(len(labels), len(tf_list))]
# 如果没有提供轴对象，则创建新的图形
if ax1 is None or ax2 is None:
fig, (ax1, ax2) = plt.subplots(2, 1, figsize=(10, 8))
show_plot = True
else:
show_plot = False
# 为每个传递函数生成Bode图
for i, tf in enumerate(tf_list):
frequencies = np.logspace(np.log10(f_min), np.log10(f_max), 1000)
omega = 2 * np.pi * frequencies
w, mag, phase = signal.bode(tf, omega)
# 绘图参数
color = colors[i] if i < len(colors) else None
linestyle = linestyles[i] if i < len(linestyles) else '-'
label = labels[i] if i < len(labels) else f'TF {i+1}'
# 幅度图
ax1.semilogx(w, mag, color=color, linestyle=linestyle, label=label)
ax1.set_ylabel('幅度 (dB)')
ax1.grid(True, which="both", ls="-", alpha=0.3)
# 相位图
ax2.semilogx(w, phase, color=color, linestyle=linestyle, label=label)
ax2.set_xlabel('角频率 (Rad/s)')
ax2.set_ylabel('相位 (度)')
ax2.grid(True, which="both", ls="-", alpha=0.3)
# 设置标题和图例
if show_plot:
if title:
plt.suptitle(title)
ax1.legend()
ax2.legend()
plt.tight_layout()
plt.show()
return ax1, ax2
def series_transfer_functions(tf1, tf2):
"""
将两个传递函数串联
参数:
tf1: 第一个传递函数 (TransferFunction对象)
tf2: 第二个传递函数 (TransferFunction对象)
返回:
串联后的传递函数
"""
# 串联传递函数：H(s) = H1(s) * H2(s)
# 分子：num = conv(num1, num2)
# 分母：den = conv(den1, den2)
num_result = np.convolve(tf1.num, tf2.num)
den_result = np.convolve(tf1.den, tf2.den)
return signal.TransferFunction(num_result, den_result)
def safe_vector_add_right(vec1, vec2):
"""Right-aligned vector addition using basic loops"""
len1 = len(vec1)
len2 = len(vec2)
if len1 >= len2:
# vec1 is longer
result = []
# Copy the beginning of vec1 that doesn't overlap
for i in range(len1 - len2):
result.append(vec1[i])
# Add the overlapping parts
for i in range(len2):
result.append(vec1[len1 - len2 + i] + vec2[i])
else:
# vec2 is longer
result = []
# Copy the beginning of vec2 that doesn't overlap
for i in range(len2 - len1):
result.append(vec2[i])
# Add the overlapping parts
for i in range(len1):
result.append(vec1[i] + vec2[len2 - len1 + i])
return result
def close_loop_transfer_function(tf):
"""
计算闭环传输函数
参数:
tf: 传递函数 (TransferFunction对象)
返回:
闭环传输函数
"""
return signal.TransferFunction(tf.num, safe_vector_add_right(tf.den, tf.num))
if __name__ == "__main__":
# 输入参数
f_min = 1
f_max = 1e6
# 外环传递函数
gain = 1
zeros = [-100]
poles = [0]
num1 = np.array(np.poly(zeros))*gain
den1 = np.poly(poles)
tf_outside = signal.TransferFunction(num1, den1)
# 内环传递函数
gain = 1
zeros = [-500]
poles = [0]
num2 = np.array(np.poly(zeros))*gain
den2 = np.poly(poles)
tf_inner = signal.TransferFunction(num2, den2)
# 内环闭环传输函数
close_tf = close_loop_transfer_function(tf_inner)
# 串联后的开环传递函数
series_tf = series_transfer_functions(tf_outside, close_tf)
# 传递函数列表
tf_list = [tf_outside, tf_inner, series_tf, close_tf]
labels = ['C_outside', 'C_inner', 'C_series_open', 'C_inner_close']
colors = ['blue', 'red', 'green', 'black']
linestyles = ['-', '--', '-.', ':']
# 绘制所有传递函数
bode_plot(tf_list, f_min, f_max, colors=colors, linestyles=linestyles,
labels=labels, title="串联控制环路仿真")
运行 Python，输出以下结果：
从结果中可以看出，内环控制器的闭环传递函数在一些频率段，增益低于 0dB，相位产生延迟。在与外环控制器串联后，这些负增益、相位延迟叠加到了串联后的开环传递函数上。
原因在于，假如将内环传递函数写为
\[\begin{align}
& C_{inner}(s) = \frac{X(s)}{Y(s)} \\
& C_{close}(s) = \frac{C_{inner}(s)}{1+C_{inner}(s)} = \frac{X(s)}{X(s)+Y(s)} \\
& \Rightarrow |C_{close}(s)| = |X(s)| - |X(s)+Y(s)| \\
& \Rightarrow \begin{cases}
|C_{close}(s)| \geq |X(s)| - |X(s)| - |Y(s)| = -|Y(s)| \\
|C_{close}(s)| \leq |X(s) - X(s) - Y(s)| = |Y(s)|
\end{cases} \\
& \angle C_{close}(s) = \tan^{-1}\left(\frac{\mathrm{Im}(X(s))}{\mathrm{Re}(X(s))}\right)-\tan^{-1}\left(\frac{\mathrm{Im}(X(s)+Y(s))}{\mathrm{Re}(X(s)+Y(s))}\right)
\end{align}
\]
从上面的推导可以看出，在一个传递函数闭环后，增益在\([-|Y(s)|,|Y(s)|]\)范围内产生增强或者衰减，相位会产生超前或者滞后，外环控制器与这些偏移叠加后可能存在不确定、不可控的情况。
解决方法也很简单，在内环传递函数上再乘以一个较大的增益，使它的模远大于1，以上问题在很大程度上得到改善：
\[\begin{align}
& C_{inner}(s) = \frac{AX(s)}{Y(s)} , |AX(s)| \gg |Y(s)| \\
& C_{close}(s) = \frac{AX(s)}{AX(s)+Y(s)} \\
& \Rightarrow |C_{close}(s)| = |AX(s)| - |AX(s)+Y(s)| \approx |AX(s)| - |AX(s)| = 0 \\
& \Rightarrow \angle C_{close}(s) = \tan^{-1}\left(\frac{\mathrm{Im}(AX(s))}{\mathrm{Re}(AX(s))}\right)-\tan^{-1}\left(\frac{\mathrm{Im}(AX(s)+Y(s))}{\mathrm{Re}(AX(s)+Y(s))}\right) \\
& \approx \tan^{-1}\left(\frac{\mathrm{Im}(AX(s))}{\mathrm{Re}(AX(s))}\right)-\tan^{-1}\left(\frac{\mathrm{Im}(AX(s))}{\mathrm{Re}(AX(s))}\right) = 0
\end{align}
\]
修改代码重新测试：
# 内环传递函数
gain = 100
zeros = [-500]
poles = [0]
num2 = np.array(np.poly(zeros))*gain
den2 = np.poly(poles)
tf_inner = signal.TransferFunction(num2, den2)
运行结果如下：
结果显示内环闭环传递函数幅度值为 0dB，相位值为 \(0^\circ\)，开环传递函数在串联前后波形几乎重合了。
变增益系统的分析
考虑一种情况，假设内环是一个 Buck 电路的输出电压控制环路，其中包含 Buck 的传递函数：
\[\begin{align}
G_{vd}(s) = \frac{\hat{v}_o}{\hat{d}} = \frac{RV_o}{LRCs^2+Ls+R}
\end{align}
\]
和一个 PI 控制器：
\[\begin{align}
G_c(s) = \frac{s+a}{s}
\end{align}
\]
其中，\(R=10\Omega,C=300\mu F,L=1mH, a=900\) 。
外环可能输出变化的内环设定值，即变化的输出电压设定值。这意味着在稳态下，内环传递函数将工作在动态的增益 \(V_o\) 下。这时，我们可以将这个 \(V_o\) 分离出来，用一个变量 \(K\) 来表示，这个环路就成为一个变增益系统：
它的特征方程为：
\[\begin{align}
1+KG'_{vd}(s)G_c(s) = 0
\end{align}
\]
其中，\(G'_{vd}\) 是 \(G_{vd}\) 分离 \(V_o\) 后的传递函数：
\[\begin{align}
G'_{vd}(s) = \frac{R}{LRCs^2+Ls+R}
\end{align}
\]
将特征方程变形为：
\[\begin{align}
G'_{vd}(s)G_c(s) = -\frac{1}{K}
\end{align}
\]
在 Python 中分别画出 \(G'_{vd}(s)G_c(s)\) 和 \(-\frac{1}{K}\) 在复平面的轨迹，根据 Nyquist 判据，如果 \(G'_{vd}(s)G_c(s)\) 的轨迹没有包围 \(-\frac{1}{K}\) 的轨迹，那么系统是稳定的。
import numpy as np
import matplotlib.pyplot as plt
from control import TransferFunction, nyquist_plot, series
from scipy import signal
def variable_gain_function(ax, A_range=None):
"""
绘制非线性环节的 -1/N(A) 轨迹
参数:
ax : Matplotlib坐标轴对象
A_range : tuple, optional
幅值范围，格式为 (min_A, max_A)。如果为None，则自动缩放
"""
# 如果没有指定A_range，则设置默认范围
if A_range is None:
A_range = (0.1, 10)
# 生成增益K的范围（这里用A表示增益）
K_values = np.linspace(A_range[0], A_range[1], 100)
# 计算-1/K的实部和虚部
neg_inv_K_real = [-1/K for K in K_values]
neg_inv_K_imag = [0 for _ in K_values] # 对于线性增益，虚部为0
# 在给定的坐标轴上绘制-1/K轨迹
ax.plot(neg_inv_K_real, neg_inv_K_imag, 'r-', linewidth=2, label='-1/K 轨迹')
# 标记一些特殊的K值点
K_special = np.linspace(A_range[0], A_range[1], 5)
for K in K_special:
ax.plot(-1/K, 0, 'ro') # 用红点标记特殊点
ax.text(-1/K, 0.1, f'K={K:.1f}', ha='center', va='bottom')
# 添加标题和标签
ax.set_xlabel('实轴')
ax.set_ylabel('虚轴')
ax.grid(True, which='both', linestyle=':', alpha=0.6)
ax.legend()
# 设置坐标轴比例一致
ax.axis('equal')
return ax
def generate_nyquist_plot(tf, omega_range=None, plot_title="Nyquist Plot",
variable_gain=False, A_range=None):
"""
基于给定的传递函数生成奈奎斯特图
参数:
tf : TransferFunction
控制系统传递函数
omega_range : tuple, optional
频率范围，格式为 (min_freq, max_freq)。如果为None，则自动缩放
plot_title : str
图形标题
variable_gain : bool
是否包含变增益 -1/K 的轨迹
返回:
fig: Matplotlib图形对象
ax: 坐标轴对象
"""
# 创建图形
fig, ax = plt.subplots(figsize=(10, 8))
# 生成奈奎斯特图
if omega_range:
omega = np.logspace(np.log10(omega_range[0]), np.log10(omega_range[1]), 1000)
nyquist_plot(tf, omega=omega, ax=ax, label="传递函数的 Nyquist 图")
else:
nyquist_plot(tf, ax=ax, label="传递函数的 Nyquist 图")
# 添加标题和网格
ax.set_title(plot_title)
ax.grid(True)
# 如果需要，添加非线性环节的 -1/K 轨迹
if variable_gain:
variable_gain_function(ax, A_range)
ax.legend()
return fig, ax
# 示例用法:
if __name__ == "__main__":
# 传递函数: G(s) = 10/(3e-6*s^2 + 1e-3*s + 10)
num1 = [10]
den1 = [3e-6, 1e-3, 10]
tf1 = TransferFunction(num1, den1)
# 传递函数: H(s) = (s+900)/s
num2 = [1, 900]
den2 = [1, 0]
tf2 = TransferFunction(num2, den2)
tf = series(tf1, tf2)
A_range = (0.1, 1)
# 绘制包含非线性环节的奈奎斯特图
fig2, ax2 = generate_nyquist_plot(tf,
plot_title="变增益系统 Nyquist 图",
variable_gain=True, A_range=A_range)
plt.show()
从输出结果可以看出，要想使这个系统稳定，\(K\) 值必须小于 0.6，才能保证 \(-\frac{1}{K}\) 的轨迹处于 \(G'_{vd}(s)G_c(s)\) 的外侧。
总结
本文通过仿真实例深入分析了串联控制环路的设计原理和关键问题。主要内容包括：
-
内外环路相互影响机制：内环的闭环特性会直接影响外环的开环传递函数，表现为增益衰减和相位滞后，这对整个系统的稳定性构成挑战。
-
解决方案验证：通过提高内环增益的方法，有效解决了内外环相互干扰的问题，使得内环闭环传递函数近似为单位增益系统，从而减少了对其他环路的影响。
-
变增益系统分析方法：针对实际应用中常见的变增益系统（如Buck电路输出电压控制），提出了基于Nyquist图的稳定性分析方法，通过绘制\(G'_{vd}(s)G_c(s)\)轨迹和\(-1/K\)轨迹，直观判断系统在不同增益条件下的稳定性边界。
这些分析方法和解决方案不仅适用于文中提到的具体案例，也为其他复杂的多环控制系统设计提供了重要的理论指导和实践参考。在实际工程应用中，合理设计各环路的增益分配和频率特性，是确保整个控制系统稳定性和性能的关键。
浙公网安备 33010602011771号