提供计算严格适合 SISO 系统的脉冲响应的程序：

输入：

* num：传递函数的分子多项式系数行向量

* den：传递函数的分母多项式系数行向量

* tO：施加单位脉冲输入的时间

* dt：时间步长（小于 1/（最大自然频率））

* tf：脉冲响应计算的最终时间

输出：

* y：脉冲响应

* t：时间点向量