PR状态方程的MATLAB代码

以下是一个简单的MATLAB代码实现PR状态方程：

% 定义PR状态方程
function x_dot = pr_state_eq(t, x, u)
    % 系统矩阵
    A = [0 1 0;
         0 0 1;
         -0.3 -0.5 -0.2];
    % 控制输入矩阵
    B = [0; 0; 1];
    % 状态向量的导数
    x_dot = A*x + B*u;
end

在这个代码中，`pr_state_eq`函数实现了PR状态方程。它接受三个输入参数：时间 `t`、状态向量 `x` 和控制输入 `u`。它返回状态向量的导数 `x_dot`。

在函数中，我们定义了系统矩阵 `A`、控制输入矩阵 `B`，然后使用这些矩阵计算状态向量的导数。注意，在这个例子中，我们使用了硬编码的系统矩阵和控制输入矩阵。在实际应用中，这些矩阵通常是根据系统模型和控制需求进行计算和设计的。

要使用这个函数，你可以将其保存为一个单独的文件，例如 `pr_state_eq.m`，然后在你的主程序中调用它。例如，如果你想在时间 `t=0` 时计算状态向量的导数，你可以使用以下代码：

% 定义时间和初始状态向量
t = 0;
x = [0; 0; 0];

% 计算控制输入
u = 1;

% 调用PR状态方程计算状态向量的导数
x_dot = pr_state_eq(t, x, u);

在这个例子中，我们传递了时间 `t`、初始状态向量 `x` 和控制输入 `u` 给 `pr_state_eq`函数，并获得状态向量的导数 `x_dot`。