离散距离过滤器

打开直播脚本

该示例显示了几种用于离散化陷波滤波器的技术的比较。虽然控制系统组件通常在连续时间内设计，但通常必须离散化以实现数字计算机和嵌入式处理器。

连续时间陷波滤波器

陷波滤波器被设计成在特定频率下抑制该频率的增益来拒绝特定频率的信号内容。对于此示例，我们考虑以下陷波滤波器：

$H （ S. ） = \frac{{S.}^{2} + 0. 。 5. S. + 1 0. 0.}{{S.}^{2} + 5. S. + 1 0. 0.}$

您可以绘制此过滤器的频率响应b命令：

H = TF（[1 0.5 100]，[1 5 100]）;BODE（H），网格

图包含2个轴。轴1包含类型线的对象。该对象表示H.轴2包含类型线的对象。该对象代表H.

该凹口滤波器在频率w = 10 rad / s处提供20dB衰减。

选择离散化方法

您可以将连续时间系统离散C2D.命令。控制系统工具箱™支持几种离散化算法，包括：金宝app

零阶持有
一阶持有
冲动不变
Tustin（双线性近似）
带频率预警的Tustin
匹配的杆和零

哪种方法选择取决于应用程序和要求。

零级和一阶保持方法和脉冲不变方法非常适合于时域中的离散近似。例如，ZOH离散化的阶跃响应与每个时间步骤（独立于采样率）匹配连续时间步长响应：

ts = 0.1;HDZ = C2D（H，TS，'ZOH'）;步骤（h，'B'，HDZ，'r'）， 传奇（'连续'那'在10 Hz离散化）

图包含轴。轴包含2个类型的型号。这些物体代表连续，在10 Hz离散化。

同样，脉冲不变离散化与原始系统相同的脉冲响应：

g = tf（[1 -3]，[1 2 10]）;gd = c2d（g，ts，'IMP'）;脉冲（g，'B'，gd，'r'） 传奇（'连续'那'在10 Hz离散化）

图包含轴。轴包含2个类型的型号。这些物体代表连续，在10 Hz离散化。

相比之下，突然突变和匹配的方法在频域中倾向于更好地在奈奎斯特频率附近引入较少的增益和相位失真。例如，使用ZOH，Tustin和匹配的算法比较连续时间陷波滤波器的Bode响应及其离散化：

HDT = C2D（H，TS，'tustin'）;HDM = C2D（H，TS，'比赛'）;BODE（H，'B'，HDZ，'r'，HDT，'M'，HDM，'G'，{1 100}），网格传奇（'连续'那'ZOH'那'tustin'那'匹配'）

图包含2个轴。轴1包含4个类型的线。这些物体代表连续，ZOH，汀汀，匹配。轴2包含4个类型线的物体。这些物体代表连续，ZOH，汀汀，匹配。

该比较表明匹配方法提供了陷波滤波器的最精确的频域近似。然而，您可以通过指定等于陷阱频率的预调高频率来进一步提高Tustin算法的准确性。这确保了靠近W = 10 rad / s的准确匹配：

HDP = C2D（H，TS，'prewarp'，10）;BODE（H，'B'，HDT，'M'，HDP，'G'，{1 100}），网格传奇（'连续'那'tustin'那'突然与预警'）

图包含2个轴。轴1包含3个类型的线。这些物体代表连续，汀汀，汀汀与预警。轴2包含3个类型的线。这些物体代表连续，汀汀，汀汀与预警。

选择采样率

采样率越高，连续和离散响应之间的匹配越近。但采样率或等同地，采样间隔的程度多么小？作为经验的规则，如果您希望连续和离散的模型将紧密地匹配一些频率浪，确保奈奎斯特频率（采样率PI）至少是两次浪。对于陷波滤波器，您需要保持10 RAD / S附近的形状，因此奈奎斯特频率应超出20 rad / s，这在最多pi / 20 = 0.16 s的采样周期。

要确认此选择，请将匹配的离散化与采样期0.1,0.15和0.3进行比较：

HD1 = C2D（H，0.1，'M'）;HD2 = C2D（H，0.15，'M'）;HD3 = C2D（H，0.3，'M'）;BODE（H，'B'，HD1，'r'，HD2，'M'，HD3，'G'，{1 100}），网格传奇（'连续'那'ts = 0.1'那'ts = 0.15'那'ts = 0.3'）