9-4 Eigen Functions (固有函數)

Old Chinese version

If the output of a system has the same type as its input signal, then the input signal is referred to as the eigen function of the system. In this section, we shall demonstrate that the exponential function (including the sinusoidal function as a special case) is an eigen function of the LTI system.

Suppose that the input signal to our LTI system is an exponential function:

x[n]=epn
Then the output can be derived based on convolution:
y[n]=x[n]*h[n]
=Skx[n-k] h[k]
=Skep(n-k) h[k]
=Skepn-pk h[k]
=epnSke-pk h[k]
=epnH(ep)
In other words, the output signal y[n] is equal to the original input signal multiplied by a constant H(ep)(= Ske-pk h[k], which is not a function of time). Hence it can be established that the exponential function is the eigen function of an LTI system.

Based on the above fact, for any given input signal that can be decomposed as a linear combination of exponential functions, we can then apply the superposition principle to find the responses of these exponential functions first, and then linearly combine these responses to obtain the final output. In mathematical notations, we have:

x[n] = aepn + beqn → y[n] = aH(ep) + bH(eq)

Since the exponential function and the sinusoidal function are related by the famous Euler identity:

ejq = cos(q) + j sin(q)
Therefore we can express the sinusoidal function as a linear combination of the exponential function:
cos(q) = (ejq + e-jq)/2 = Re{ejq}
sin(q) = (ejq - e-jq)/(2j) = Im{ejq}
As a result, when the input signal is:
x[n] = cos(wn) = Re{ejwn}
The corresponding response is:
y[n]=Re{ejwnH(ejw)}
=Re{ejwn |H(ejw)| ejq}, q = ∠H(ejw)
=Re{|H(ejw)| ej(wn+q)}
=|H(ejw)| cos(wn + q)
From the above derivation, we can see that the output is still a sinusoidal function, except that the amplitude is multiplied by |H(ejw)|, and the phase is shifted by q = ∠H(ejw).

Some terminologies are explained next.

Similarly, when the input signal is
x[n] = sin(wn) = Im{ejwn}
The output can be derived as follows:
y[n]=Im{ejwnH(ejw)}
=Im{ejwn |H(ejw)| ejq}, q = ∠H(ejw)
=Im{|H(ejw)| ej(wn+q)}
=|H(ejw)| sin(wn + q)
In summary, we have the following important conclusion. For an LTI system, if the input signal is a single sinusoidal function, the output is also a same-frequency sinusoidal function, with the amplitude and phase modified by the system's frequency response H(ejw).

In general, we can use an complex exponential function to subsume the input of sine and cosine functions:

x[n] = ejwn
The corresponding output is
y[n] = ejwnH(ejw)
Alternatively, H(ejw) is also referred to as the the discrete-time Fourier transform of h[n], which can be expressed as the following equation:
H(ejw) = Skh[k]e-jwk
In the real world, we cannot have the complex exponential function as the input to a system. However, the derivation based on complex numbers does hold mathematically and it indeed makes our derivation more concise.
Audio Signal Processing and Recognition (音訊處理與辨識)