查看“︁冲激不变法”︁的源代码

{{NoteTA
|G1=Signals and Systems
}}
'''冲激不变法'''是利用连续时间滤波器来设计离散时间[[无限冲激响应]]（IIR）滤波器的一种方法，这种方法中对连续时间系统的冲激响应进行采样以产生离散时间系统的冲激响应。离散时间系统的频率响应就会是连续时间系统的频率响应的位移后的拷贝之和；如果连续时间系统的频带大致限制在小于采样的[[奈奎斯特频率]]的范围内，则离散时间系统的频率响应会大致与连续系统的频带相同，低于奈奎斯特频率。

==讨论==

以 <math>T</math> 为采样周期对连续时间系统的冲激响应 <math>h_c(t)</math> 采样得到了离散时间系统的冲激响应 <math>h[n]</math>。

:<math>h[n]=Th_c(nT)\,</math>

因此，该系统的频率响应为

:<math>H(e^{j\omega}) = \sum_{k=-\infty}^\infty{H_c\left(j\frac{\omega}{T} + j\frac{2{\pi}}{T}k\right)}\,</math>

如果连续时间滤波器是大致是带限的（即当 <math>|\Omega| \ge \pi/T</math> 时，<math>H_c(j\Omega) < \delta</math>），则每次采样的频率低于 π（即奈奎斯特频率低于 1/(2''T'') Hz）的离散时间系统的频率响应就会大致为连续时间系统的频率响应：

:<math>H(e^{j\omega}) = H_c(j\omega/T)\,</math> for <math>|\omega| \le \pi\,</math>

===与双线性变换比较===

注意到会出现混叠，包含有奈奎斯特频率以下与超过该频率连续时间滤波器的非零的响应混叠。[[雙線性轉換]]使用不同的映射方法，将连续系统能够达到无限频率的频率响应，映射到在离散时间系统中至多能达到奈奎斯特频率的范围，避免了冲激不变法线性地对频率映射产生的循环混叠，从而可以替代冲激不变法。

===系统函数中极点的效应===

若连续极点位于 <math>s = s_k</math>，系统函数可以用部分分式展开写作

:<math>H_c(s) = \sum_{k=1}^N{\frac{A_k}{s-s_k}}\,</math>

因此，使用拉普拉斯逆变换，得到冲激响应为

:<math>h_c(t) = \begin{cases}
  \sum_{k=1}^N{A_ke^{s_kt}},  & t \ge 0 \\
  0, & \mbox{otherwise}
\end{cases}</math>

对应的离散时间系统的冲激响应定义如下

:<math>h[n] = Th_c(nT)\,</math>

:<math>h[n] = T \sum_{k=1}^N{A_ke^{s_knT}u[n]}\,</math>

对离散时间冲激响应进行Z变换得到下面的离散时间系统函数

:<math>H(z) = T \sum_{k=1}^N{\frac{A_k}{1-e^{s_kT}z^{-1}}}\,</math>

于是连续时间系统函数的极点被搬移到  z = e<sup>s<sub>k</sub>T</sup> 处的极点。零点则不会这样简单地映射过去。{{clarify|date=2013年2月}}

===零极点===

如果系统函数既有零点也有极点，可以用相同方式映射，但结果不再是冲激不变的结果：离散时间冲激响应不是简单地与连续时间冲激响应相等。这种方法名为[[匹配Z变换方法]]，或叫做零极点映射。对于全极点滤波器，两种方法是等价的。

===稳定性和因果性===

由于连续时间系统 ''s'' = ''s<sub>k</sub>'' 处的极点转换到离散时间系统 z = exp(''s<sub>k</sub>T'') 处的极点， ''s'' 左半平面的极点就会映射到 ''z'' 平面的单位圆内部；所以若连续时间滤波器是因果稳定的，则离散时间滤波器也会时因果稳定的。

===修正公式===

当一个因果的连续时间冲激响应在 <math>t=0</math> 处不连续时，上述表达式不一致。<ref>{{Cite journal|title = A correction to impulse invariance|url = http://ieeexplore.ieee.org/xpl/freeabs_all.jsp?tp=&arnumber=870677&isnumber=18861|journal = IEEE Signal Processing Letters|date = 2000-10-01|issn = 1070-9908|pages = 273-275|volume = 7|issue = 10|doi = 10.1109/97.870677|first = L.B.|last = Jackson}}</ref>
这是因为 <math>h_c (0)</math> 应道贡献一半的值给 <math>h[0]</math>。

修正为：

:<math>h[n] = T \left( h_c(nT) - \frac{1}{2}  h_c(0)\delta [n] \right) \,</math>

:<math>h[n] = T \sum_{k=1}^N{A_ke^{s_knT}} \left( u[n] - \frac{1}{2} \delta[n] \right) \,</math>

对离散时间冲激响应进行Z变换得到下面的离散时间系统函数

:<math>H(z) = T \sum_{k=1}^N{\frac{A_k}{1-e^{s_kT}z^{-1}} - \frac{T}{2} \sum_{k=1}^N A_k}.</math>

==参见==

* [[無限脈衝響應濾波器]]
* [[雙線性轉換]]
* [[匹配Z变换方法]]
* 连续时间滤波器：
*:[[切比雪夫滤波器]]
*:[[巴特沃斯滤波器]]
*:[[椭圆函数滤波器]]

==参考文献==

{{Reflist}}

===其他来源===

{{Refbegin}}
* Oppenheim, Alan V. and Schafer, Ronald W. with Buck, John R.  ''Discrete-Time Signal Processing.''    Second Edition.  Upper Saddle River, New Jersey: Prentice-Hall, 1999.
* Sahai, Anant. Course Lecture. Electrical Engineering 123: Digital Signal Processing. University of California, Berkeley. 5 April 2007.

==外部链接==
* [http://www.circuitdesign.info/blog/2009/02/impulse-invariant-transform/ Impulse Invariant Transform at CircuitDesign.info]{{Wayback|url=http://www.circuitdesign.info/blog/2009/02/impulse-invariant-transform/ |date=20160303210334 }} Brief explanation, an example, and application to [[连续信号|Continuous Time]] [[Sigma Delta]] ADC's.
{{Refend}}

<!-- AWB bots, leave this extra space alone! -->

{{DSP}}

{{DEFAULTSORT:冲激不变法}}
[[Category:数字信号处理]]
[[Category:滤波器理论]]