Impulse Invariance for Digital Filter Design

Graduate Depth 137 in the knowledge graph I know this Set as goal
filter-design digital-filters discretization impulse-response

Core Idea

Impulse invariance maps analog filter impulse response to digital by sampling: h[n] = T·ha(nT). This preserves the shape of the analog response at sample times but introduces aliasing if the analog filter is not sufficiently bandlimited. Unlike bilinear transform, it does not preserve stability for poles outside the left half-plane.

Explainer

From your study of the Laplace transform, you know how to design analog filters — Butterworth, Chebyshev — specified by poles in the s-plane. From the z-transform, you know how to analyze discrete-time systems. Impulse invariance is one bridge between these two worlds: it constructs a digital filter whose behavior deliberately mimics a known analog prototype, starting from the analog impulse response.

The method begins by expanding the analog transfer function H_a(s) into partial fractions: H_a(s) = Σ A_k / (s − s_k). Each term corresponds to a decaying exponential in the time domain: h_a(t) = Σ A_k e^(s_k t) u(t). Impulse invariance samples this continuous impulse response at intervals T to define the digital impulse response: h[n] = T · h_a(nT) = T · Σ A_k e^(s_k nT) u[n]. Taking the z-transform of each term gives H(z) = T · Σ A_k / (1 − e^(s_k T) z^{-1}). This is the key mapping rule: each analog pole at s = s_k becomes a digital pole at z = e^(s_k T). Poles in the left half s-plane (where Re(s_k) < 0) map to poles inside the unit circle (where |e^(s_k T)| < 1), preserving stability. The frequency axis maps as ω_digital = Ω_analog × T — digital frequencies are a scaled version of analog frequencies, with no warping.

The fatal limitation is aliasing, and your prerequisite on aliasing tells you exactly why it arises. The digital filter's frequency response is a periodic repetition of the analog: H(e^{jω}) = (1/T) Σ_k H_a(j(ω − 2πk)/T). If the analog filter has significant energy above the Nyquist frequency π/T, adjacent copies overlap and corrupt the spectrum. Butterworth and Chebyshev filters have infinite bandwidth — their magnitude falls off polynomially, never reaching exactly zero. Applying impulse invariance to these filters inevitably introduces aliasing. For a well-designed lowpass filter with substantial attenuation before Nyquist (say, −40 dB or better at π/T), the aliasing error may be acceptable. For highpass or wideband filters, the copies wrap around and add in-band, making the result useless.

This is why bilinear transform is often the preferred alternative: it maps the entire left half-plane to the unit disk interior without aliasing by compressing the infinite analog frequency axis onto [−π, π]. The trade-off is frequency warping — a nonlinear compression that distorts the filter shape unless pre-warped during design. Impulse invariance is preferred when you specifically need the digital filter's time-domain samples to match the analog prototype's values — for instance, when digitally simulating an analog physical system where matching the step response at sample instants matters more than matching the stopband at high frequencies. Knowing both methods, and when each is appropriate, is the practical skill.

Practice Questions 5 questions

Prerequisite Chain

Counting to 10Counting to 20Understanding ZeroThe Number ZeroCounting to FiveOne-to-One CorrespondenceCombining Small Groups Within 5Addition Within 10Addition Within 20Two-Digit Addition Without RegroupingTwo-Digit Addition with RegroupingAddition Within 100Repeated Addition as MultiplicationMultiplication Facts Within 100Division as Equal SharingDivision as Grouping (Measurement Division)Division: Grouping (Repeated Subtraction) ModelDivision: Fair Sharing ModelDivision as Equal SharingDivision as GroupingBasic Division FactsDivision Facts Within 100Two-Digit by One-Digit DivisionDivision with RemaindersRemainders and Quotients in DivisionDivision Word ProblemsIntroduction to Long DivisionFactors and MultiplesPrime and Composite NumbersEquivalent FractionsRelating Fractions and DecimalsDecimal Place ValueReading and Writing DecimalsComparing and Ordering DecimalsAdding and Subtracting DecimalsMultiplying DecimalsDividing DecimalsDividing FractionsMixed Number ArithmeticOrder of OperationsInteger Order of OperationsVariable ExpressionsCombining Like TermsOne-Step EquationsTwo-Step EquationsSolving Multi-Step EquationsEquations with Variables on Both SidesAngle Pairs: Complementary, Supplementary, and VerticalParallel Lines and TransversalsCorresponding AnglesAlternate Interior AnglesTriangle Angle Sum TheoremExterior Angle TheoremTriangle Inequality TheoremSimilar Triangles: AA SimilaritySimilar Triangles: SSS and SAS SimilarityProportions in Similar TrianglesRight Triangle Trigonometry IntroductionTrigonometric Ratios ReviewRadian MeasureConverting Between Degrees and RadiansThe Unit CircleGraphing Sine and CosineGraphing Tangent and Reciprocal Trigonometric FunctionsDerivatives of Trigonometric FunctionsAntiderivativesIterated Integrals and Fubini's TheoremDouble Integrals in Cartesian CoordinatesDouble Integrals over Rectangular RegionsDouble Integrals in Polar CoordinatesDouble Integrals: Definition and SetupIterated Integrals and Fubini's TheoremDouble Integrals over Rectangular RegionsDouble Integrals over General RegionsApplications of Double Integrals: Area, Mass, and MomentsTriple Integrals in Cartesian CoordinatesTriple Integrals in Cylindrical and Spherical CoordinatesChange of Variables and the Jacobian DeterminantApplications of Triple Integrals: Volume and MassVector Fields and Their RepresentationsLine Integrals of Vector FieldsGreen's TheoremSurface Integrals and Flux of Vector FieldsSurface Integrals and Flux of Vector FieldsDivergence Theorem: Flux and OutflowDivergence TheoremElectric FluxGauss's LawConductors in Electrostatic EquilibriumCapacitance and CapacitorsDielectricsDielectric Constant and Relative PermittivityElectric Field Inside Dielectric MaterialsDielectric Materials and PolarizationDielectric Susceptibility and PermittivityEnergy Density in Electric FieldsElectric Current and Current DensityElectrical Resistance and ResistivityOhm's Law and Circuit ElementsElectromotive Force (EMF) and BatteriesKirchhoff's Circuit Laws: Voltage and CurrentDC Circuit Network Analysis MethodsTransient Response in RC CircuitsRC CircuitsLC and RLC CircuitsAC Circuits: FundamentalsImpedance and ReactanceAC Power and ResonanceElectromagnetic WavesFrequency-Dependent Permittivity and DispersionElectromagnetic Waves in Anisotropic MediaBirefringence and DichroismWave Plates: Quarter-Wave and Half-Wave PlatesCircular and Elliptical Polarization ProductionPolarization States: Linear, Circular, and EllipticalLinear Superposition of WavesSuperposition Principle in ElectrostaticsElectric Field Lines and VisualizationElectric Potential and Potential EnergyElectric Potential and VoltageIdeal Voltage and Current SourcesSeries, Parallel, and Combined Resistor NetworksVoltage Divider Principle and ApplicationsKirchhoff's Voltage and Current LawsNodal Analysis MethodLinearity, Superposition, and ScalingAC Steady-State Circuit AnalysisAC Circuit Analysis Using PhasorsAC Power AnalysisResonance in RLC CircuitsFrequency Response and Bode PlotsNotch Filters and Resonator DesignAnti-Aliasing Filters and Pre-Sampling DesignDecimation, Anti-Aliasing, and DownsamplingAliasing, Anti-Aliasing Filters, and Signal ReconstructionDigital Signal Processing FundamentalsFIR Filter Design and RealizationImpulse Invariance for Digital Filter Design

Longest path: 138 steps · 791 total prerequisite topics

Prerequisites (6)

Leads To (0)

No topics depend on this one yet.