利用Simulink中的S函数绘制关于控制系统变量zl和z2的仿真曲线。其满足的 微分方程如下所示： 21=Z2 Z2=Z1Z2 Z1(0)=1 Z2(0)=-1 要求： 绘制出Simulink方框图： 写出S函数中状态变量更新子函数和结果输出子函数的程序； 给出仿真曲线的结果。

为了在Simulink中创建一个S函数来模拟这个系统，我们需要首先定义状态空间模型，然后在S函数中编写状态变量的更新和输出计算。以下是逐步的步骤：

1. **创建Simulink模型**:
- 打开Simulink并在图形界面中添加一个新模型（File > New Model）。
- 创建一个新的S函数图标（ Simscape > S-functions 或者 Library Browser > S-function Block），双击打开编辑器。

2. **定义S函数**:
- 在S函数编辑器中，首先设置基本结构。声明必要的全局变量 `zl` 和 `z2`，它们代表系统的状态。
- 初始化函数（`function z = mySFunction(x, u, p)`）中，处理初始条件 `x`（通常是 `[]` 表示无输入，因为这是仿真开始时的状态）。在这里设置 `zl` 和 `z2` 初始值为给定的值。

   function [zl, z2] = mySFunction(x, ~, p)
       zl = 1;
       z2 = -1;
   end

3. **状态更新函数（State Update Function, dx/dt）**:
- 在 `dx` 函数中，根据系统微分方程 `dzl/dt = z2` 和 `dz2/dt = z1 * z2` 编写状态更新。

   function dz = stateUpdateFcn(~, x, u, ~)
       dz = zeros(2, 1);
       dz(1) = x(2); % dzl/dt = z2
       dz(2) = x(1) * x(2); % dz2/dt = z1 * z2
   end

4. **结果输出函数（Output Function, y）**:
- 在 `output` 函数中，根据需要选择要返回的输出变量。这里可能只需要 `zl` 和 `z2`。

   function y = outputFcn(~, ~, x, ~)
       y = x; % 返回状态变量zl和z2
   end

5. **连接S函数到Simulink**:
- 将刚刚编写的S函数拖放到Simulink模型中，然后将其命名为 "mySFunction"。
- 连接该S函数到适当的位置，例如在状态空间框内，使得 `zl` 和 `z2` 的输出直接传递到其他模块或输出端口。

6. **执行仿真**:
- 设置仿真时间范围（可能是 `tspan = [0, 10];` 或类似的时间步长）。
- 在Simulink模型中运行仿真，观察 `zl` 和 `z2` 变量的变化曲线。

7. **查看结果**:
- 检查 Simulink 仿真窗口中生成的 `zl` 和 `z2` 时间序列图表，确保它们符合预期。