如何用MATLAB编写代码求解控制系统控制器的K与T参数？

前提说明 #

由于原问题未给出被控对象的具体传递函数，以下以典型二阶被控对象 G(s) = 1/(s(s+1)) 为例，控制器采用超前校正形式 Gc(s) = K(Ts+1) 展开实现。若你的实际系统结构不同，只需替换对应传递函数即可。

核心思路 #

闭环系统的特征方程需与目标极点 s = -2±j2 对应的特征方程匹配：

• 目标特征方程：(s+2+j2)(s+2-j2) = s² + 4s + 8 = 0
• 由开环传递函数推导闭环特征方程，通过系数相等建立方程组求解K和T。

方法一：符号计算法（依赖Symbolic Math Toolbox） #

直接通过符号变量建立方程，求解参数：

% 定义符号变量
syms K T

% 目标特征方程的系数
target_coeff = [1 4 8];

% 根据系统结构推导闭环特征方程系数
% 开环传递函数：G_open(s) = K(Ts+1)/(s(s+1))
% 闭环特征方程：s² + (1+KT)s + K = 0
sys_coeff = [1 1+K*T K];

% 建立系数相等的等式
eq1 = sys_coeff(2) == target_coeff(2);
eq2 = sys_coeff(3) == target_coeff(3);

% 求解方程组
sol = solve([eq1, eq2], [K, T]);

% 输出结果
disp('求解得到的控制器参数：');
disp(['K = ', num2str(sol.K)]);
disp(['T = ', num2str(sol.T)]);

方法二：数值求解法（依赖Optimization Toolbox） #

针对复杂系统，无法直接写出符号方程时，使用fsolve进行数值求解：

% 构造目标极点对应的特征多项式
target_poly = conv([1 2+2j], [1 2-2j]);

% 定义误差函数：输入[K,T]，输出特征方程在目标极点处的误差
fun = @(x) polyval([1 1+x(1)*x(2) x(1)], roots(target_poly)) - 0;

% 设置参数初始猜测值（可根据系统特性调整）
x0 = [1, 1];

% 求解方程组
sol = fsolve(fun, x0);

% 输出结果
disp('求解得到的控制器参数：');
disp(['K = ', num2str(sol(1))]);
disp(['T = ', num2str(sol(2))]);

结果验证 #

求出参数后，可验证闭环极点是否符合要求：

% 代入求解得到的参数（根据所用方法选择赋值方式）
K = sol.K; % 符号法用此行，数值法替换为 K = sol(1);
T = sol.T; % 数值法替换为 T = sol(2);

% 构建控制器与被控对象传递函数
Gc = tf(K*[T 1], [1]);
G = tf(1, [1 1 0]);
G_open = Gc * G;

% 构建闭环系统
G_close = feedback(G_open, 1);

% 获取并打印闭环极点
poles = pole(G_close);
disp('验证得到的闭环极点：');
disp(poles);

注意事项 #

1. 若你的被控对象传递函数不同，需同步修改代码中特征方程系数或传递函数定义
2. fsolve的求解结果受初始猜测值影响，建议根据系统物理意义设置合理初始值
3. 确保已安装对应MATLAB工具箱（符号工具箱/优化工具箱）

内容的提问来源于stack exchange，提问作者NITIN KUMAR