matlab代码用R_H判据判断系统稳定性

在MATLAB中，Routh-Hurwitz (R_H) 稳定性判据是一种经典的线性系统稳定性分析方法，用于判断连续时间系统的特征根是否都在复平面上的左半平面，即实部为负。这个判据通常应用于数学模型的稳定性检查，特别是对于一阶和二阶系统。

以下是简单的步骤：

1. **设置矩阵**：
设定系统的一阶或二阶传递函数系数，通常形式为 `num` 和 `den`（分子多项式和分母多项式）。

num = [a b]; % a 是常数项，b 是比例系数
den = [1 c d]; % c 和 d 分别是一阶或二阶多项式的系数

2. **构造Routh表**：
使用 `routhTable` 函数创建Routh表，该表列出了系数的递归关系。

RHtable = routhTable(num, den);

3. **应用Routh-Hurwitz规则**：
检查Routh表的第一行和第一列的元素，如果所有元素都是正的，并且最后一列的第一个非零元素也是正的，则系统稳定。

if all(RHtable(1:end-1, 1)) && RHtable(end, 1) > 0
    stability = 'Stable';
else
    stability = 'Unstable';
end

4. **结果输出**：
输出稳定性结论。

disp(['System stability is: ', stability]);