在matlab将传递函数离散化的代码

在MATLAB中，将连续传递函数（也称作微分方程模型）离散化通常是为了在计算机上进行仿真和控制系统的数字分析。最常用的一种方法是通过`simulink`工具箱中的State-Space (`ss`) 函数结合`c2d`函数来进行。

假设你已经有了一个状态空间模型 `sys`（例如由 `tf` 或 `ode23s` 函数创建），可以按照以下步骤将其离散化：

1. 确定采样时间 (Ts)：离散化前你需要选择一个合适的采样周期，比如 Ts = 0.1 秒。

2. 将传递函数转换为状态空间形式（如果还没的话）：

   if isa(sys, 'TransferFunction')
       ss_sys = tf2ss(sys);
   end

3. 使用 `c2d` 函数进行离散化：

   discrete_sys = c2d(ss_sys, Ts, 'zoh'); % zoh 表示零阶保持（Zero Order Hold）

这里 'zoh' 是离散化方法的选择，还有其他如 'impulse'、'euler' 等可以选择，取决于你的应用需求。

4. 查看离散化后的系统信息：

   info = bode discrete_sys; % 可视化离散系统的频率响应