Skip to main content

Why Stable Systems Have Poles Inside the Unit Circle but AR Model Roots Are Outside?


Why Stable Systems Have Poles Inside the Unit Circle but AR Model Roots Are Outside?

Why Stable Systems Have Poles Inside the Unit Circle but AR Model Roots Are Outside?

If you've studied Digital Signal Processing (DSP) or Time Series Analysis, you've probably encountered two statements that seem contradictory:

  • A stable discrete-time system must have its poles inside the unit circle.
  • A stable AR (Auto-Regressive) model requires its roots to lie outside the unit circle.

At first glance, these statements appear to conflict. Fortunately, they are both correct. The apparent contradiction comes from the fact that DSP and statistics define the characteristic polynomial differently.

This article explains why both statements are true and clears up one of the most common sources of confusion in signal processing.


Stability of Continuous-Time Systems

For a continuous-time system, the transfer function is

H(s) = N(s) / D(s)

A system is stable only if all poles satisfy

Re(s) < 0

This means every pole must lie in the left-half of the complex plane.

The reason is straightforward. The natural response contains exponential terms of the form

est = e(σ + jω)t = eσtejωt

When σ is negative, the exponential decays over time, ensuring stability.


Stability of Discrete-Time Systems

For a discrete-time system, the transfer function is

H(z) = B(z) / A(z)

The stability condition changes to

|z| < 1

In other words, all poles must lie inside the unit circle.

Why?

The system response contains terms of the form

zn

Depending on the magnitude of z:

Condition Behavior Result
|z| < 1 zn → 0 Stable
|z| > 1 zn → ∞ Unstable

This is why DSP textbooks always state that stable systems have poles inside the unit circle.


Then Why Do AR Models Require Roots Outside the Unit Circle?

This is where the confusion begins.

An AR model is usually written as

x[n] − a₁x[n−1] − a₂x[n−2] − ... = e[n]

Taking the Z-transform gives

A(z)X(z) = E(z)

where

A(z) = 1 − a₁z−1 − a₂z−2 − ...

The transfer function becomes

H(z) = 1 / A(z)

The poles are obtained by solving

A(z) = 0

These poles must satisfy

|z| < 1

So from the DSP perspective, AR model poles must be inside the unit circle.


Where Does Roots Outside the Unit Circle Come From?

Many statistics and econometrics textbooks define the characteristic polynomial differently:

1 − a₁z − a₂z² − ... − apzp = 0

Notice that this equation uses z instead of z⁻¹.

If we define

w = 1 / z

then every characteristic root is simply the reciprocal of a pole.

Example

Suppose a stable system has a pole at

z = 0.5

The corresponding characteristic root is

w = 1 / 0.5 = 2

Representation Location
DSP Pole Inside Unit Circle (0.5)
Characteristic Root Outside Unit Circle (2)

Both describe exactly the same stable system.


What About MA Models?

An MA (Moving Average) model is written as

x[n] = B(z)e[n]

where

B(z) = 1 + b₁z−1 + b₂z−2 + ...

Unlike AR models, MA models have no non-trivial poles. Therefore, every finite-order MA model is inherently stable.

Then why must MA roots lie outside the unit circle?

The answer is invertibility, not stability.

To recover the input from the output, we use the inverse filter

1 / B(z)

For this inverse filter to remain stable, its poles (which are the zeros of the MA model) must lie inside the unit circle.

Therefore, the zeros of the MA model itself must lie outside the unit circle.


Summary: Stability vs. Invertibility

System Condition Reason
Discrete-Time LTI System Poles inside the unit circle Stability
AR Model (DSP notation) Poles inside the unit circle Stable filter
AR Model (Statistics notation) Characteristic roots outside the unit circle Roots are reciprocals of poles
MA Model Always stable No non-trivial poles
MA Model Zeros outside the unit circle Stable inverse (invertibility)

Summary

  • DSP works directly with the poles of the transfer function, which must lie inside the unit circle.
  • Time-series analysis often discusses the roots of the characteristic polynomial, which are the reciprocals of those poles and therefore lie outside the unit circle.

Why Are MA (Moving Average) Models Always Stable?

One of the key differences between Auto-Regressive (AR) and Moving Average (MA) models is that a finite-order MA model is always stable. This often surprises beginners because AR models require stability conditions, whereas MA models do not.

Understanding the MA Model

A Moving Average model of order q, denoted as MA(q), is defined as

x[n] = e[n] + b1e[n-1] + b2e[n-2] + ... + bqe[n-q]

where:

  • x[n] is the output signal.
  • e[n] is white noise (input).
  • b1, b2, ..., bq are the MA coefficients.

Taking the Z-transform gives the transfer function

H(z) = 1 + b1z-1 + b2z-2 + ... + bqz-q

Notice the Important Difference

Unlike an AR model, the MA transfer function contains only a numerator polynomial. It can also be written as

H(z) = B(z)

Since there is no denominator polynomial, there are no finite poles that can make the system unstable.

Impulse Response of an MA Model

The impulse response of an MA(q) model is

h[n] = {1, b1, b2, ..., bq, 0, 0, 0, ...}

Notice that the impulse response contains only q + 1 non-zero samples. After that, it becomes exactly zero.

Why Does This Guarantee Stability?

A discrete-time LTI system is BIBO (Bounded-Input Bounded-Output) stable if its impulse response is absolutely summable:

∑ |h[n]| < ∞

Since an MA model has only a finite number of non-zero impulse response samples, the summation becomes

|1| + |b1| + |b2| + ... + |bq|

This sum is always finite for finite coefficient values. Therefore,

∑ |h[n]| < ∞

which satisfies the BIBO stability condition.

What About the Roots of an MA Model?

Although an MA model is always stable, its zeros (also called MA roots) still play an important role. These roots do not affect stability. Instead, they determine whether the model is invertible.

For an MA model to be invertible, all its zeros must lie outside the unit circle. This ensures that the inverse filter

1 / B(z)

has all its poles inside the unit circle and is therefore stable.

Summary
  • An MA model has no denominator polynomial, so it has no finite poles that can cause instability.
  • Its impulse response is of finite duration, making it absolutely summable.
  • Therefore, every finite-order MA model is always BIBO stable.
  • The roots (zeros) of an MA model are related to invertibility, not stability.

Contact Us

Name

Email *

Message *

Popular Posts

Online Simulator for ASK, FSK, and PSK

Interactive Digital Signal Processing (DSP) Tutorial and Simulator for ASK, FSK, and BPSK modulation techniques. Try our new Digital Signal Processing Simulator!   •   Interactive ASK, FSK, and BPSK tools updated for 2025. Start Now Digital Modulation Visualizer: ASK, FSK, & BPSK Simulator Learn and visualize binary modulation techniques (ASK, FSK, BPSK) in real-time with adjustable carrier and sampling parameters. Perfect for DSP students and engineers. 📡 ASK Simulator 📶 FSK Simulator 🎚️ BPSK Simulator 📚 More Topics ASK Modulator FSK Modulator BPSK Modulator More Topics 1. ASK (Amplitude Shift Keying) Simulat...

BER vs SNR for M-ary QAM, M-ary PSK, QPSK, BPSK, ...(MATLAB Code + Simulator)

Bit Error Rate (BER) & SNR Guide Analyze communication system performance with our interactive simulators and MATLAB tools. 📘 Theory 🧮 Simulators 💻 MATLAB Code 📚 Resources BER Definition SNR Formula BER Calculator MATLAB Comparison 📂 Explore M-ary QAM, PSK, and QPSK Topics ▼ 🧮 Constellation Simulator: M-ary QAM 🧮 Constellation Simulator: M-ary PSK 🧮 BER calculation for ASK, FSK, and PSK 🧮 Approaches to BER vs SNR What is Bit Error Rate (BER)? The BER indicates how many corrupted bits are received compared to the total number of bits sent. It is the primary figure of merit f...

MATLAB Code for ASK, FSK, and PSK (with Online Simulator)

MATLAB Code for ASK, FSK, and PSK Comprehensive implementation of digital modulation and demodulation techniques with simulation results. 📘 Theory 📡 ASK Code 📶 FSK Code 🎚️ PSK Code 🕹️ Simulator 📚 Further Reading Amplitude Shift Frequency Shift Phase Shift Live Simulator ASK, FSK & PSK HomePage MATLAB Code MATLAB Code for ASK Modulation and Demodulation COPY % The code is written by SalimWireless.Com clc; clear all; close all; % Parameters Tb = 1; fc = 10; N_bits = 10; Fs = 100 * fc; Ts = 1/Fs; samples_per_bit = Fs * Tb; rng(10); binar...

UGC NET Electronic Science Previous Year Question Papers with Solutions

Home / Engineering & Other Exams / UGC NET 2026 PYQ ⬇️ Download Papers and Solutions 📋 Exam Pattern 💡 Preparation Tips ❓ FAQs 📊 Exam Highlights: Electronic Science (88) Feature Details Junior Research Fellowship (JRF) ₹37,000 + HRA per month Eligibility M.Sc/M.Tech in Electronics (55%) Validity of Certificate JRF (3 Years) | Lectureship (Lifetime) 📥 Download UGC NET Electronics PDFs Complete collection of previous year question papers, answer keys and explanations for Subject Code 88. Start Downloading 📂 View All Question Papers June 2025 - Question Paper Download PDF June 2025 - Solved Paper + Explanation ...

DFTs-OFDM vs OFDM: Why DFT-Spread OFDM Reduces PAPR Effectively (with MATLAB Code)

Understanding PAPR in DFT-spread OFDM vs. Standard OFDM In modern wireless communications like 4G LTE and 5G NR, managing the Peak-to-Average Power Ratio (PAPR) is critical for hardware efficiency. While OFDM is the gold standard for high-speed data, its high PAPR poses significant challenges for mobile devices. This is where DFTs-OFDM (also known as SC-FDMA) comes in. DFT-spread OFDM (DFTs-OFDM) has lower Peak-to-Average Power Ratio (PAPR) because it "spreads" the data in the frequency domain before applying IFFT, making the time-domain signal behave more like a single-carrier signal rather than a multi-carrier one like OFDM. Deeper Explanation: Aspect OFDM DFTs-OFDM Signal Type Multi-carrier Single-carrier-like Process IFFT of QAM directly QAM → DFT → IFFT PAPR Level High (due to many...

High Level and Low Level Modulation

📘 Overview 📚 Low Level Modulation 📚 High Level Modulation 📚 Simulator for High and Low Level Modulations 📚 MATLAB Code 📚 Further Reading High Level and Low Level Modulation Wireless communication is essential for long-distance data transmission. In wireless systems, modulation is necessary to avoid interference and to reduce antenna size significantly. In the modulation process, we translate the low-frequency baseband signal to a higher frequency by modulating it with a high-frequency carrier signal. For a typical communication system, we generate the carrier signal using a local oscillator. The message signal is modulated with this carrier and then transmitted through the antenna. Low Level Modulation In low-level modulation , the message signal is modulated with the carrier signal at low power levels . After modulation, the signal is passed through a series of linear amplifiers to increase its strength. While this method is easier to imple...

Comparisons among ASK, PSK, and FSK (with MATLAB + Simulator)

Modulation ASK, FSK & PSK Constellation MATLAB Simulink MATLAB Code Comparisons among ASK, PSK, and FSK 📘 Comparisons among ASK, FSK, and PSK 🧮 Online Simulator Bandwidth 🧮 MATLAB Code BER Analysis 📚 Further Reading 📂 View Other Topics on Comparisons among ASK, PSK, and FSK ... 🧮 Comparisons of Noise Sensitivity, Bandwidth, Complexity, etc. 🧮 MATLAB Code for Constellation Diagrams of ASK, FSK, and PSK 🧮 Online Simulator for ASK, FSK, and PSK Generation 🧮 Online Simulator for ASK, FSK, and PSK Constellation 🧮 Some Questions and Answers Comparisons among ASK, PSK, and FSK Comparison among ASK, FSK, and PSK Parameters ASK FSK PSK Variable Characteristics Amplitude ...

MATLAB Code for BER performance of QPSK with BPSK, 4-QAM, 16-QAM, 64-QAM, 256-QAM, etc

📘 Overview 🧮 MATLAB Codes 🧮 Online Simulator for Calculating BER of M-ary PSK and QAM 🧮 QPSK vs BPSK and QAM: A Comparison of Modulation Schemes in Wireless Communication 🧮 Are QPSK and 4-PSK same? 📚 Further Reading   QPSK offers double the data rate of BPSK while maintaining a similar bit error rate at low SNR when Gray coding is used. It shares spectral efficiency with 4-QAM and can outperform 4-QAM or 16-QAM in very noisy channels. QPSK is widely used in practical wireless systems, often alongside QAM in adaptive modulation schemes [Read more...] What is the Gray Code? Gray Code: Gray code is a binary numeral system where two successive values differ in only one bit. This property is called the single-bit difference or unit distance code. It is also known as reflected binary code. Let's convert binary 111 to Gray code: Binary bits: B = 1 1 1 Apply the rule: G[0] = B[0] = 1...