I am interested in characterising the time-dependent response of a system to a sinusoidal input.
The steady state response of the system is also sinusoidal, however I wonder if there is a standard way to represent the information in a manner that is more detailed than the Frequency Response Function to capture the transient information too (as a function of frequency input)?
E.g. The system where input: s(t)=sin(ωt),0≤t<∞ gives output: y(t)=(1+e−t)sin(ωt),0≤t<∞.
The frequency response would be unity, and the impulse response could not be used to capture the transient behaviour?
Answer
As already pointed out in Hilmar's answer, a linear time-invariant (LTI) system is completely described by its impulse response, or - equivalently - by its frequency response, which is the Fourier transform of its impulse response. I think that there is a misunderstanding concerning the Fourier transform, namely that it can only be used to describe steady-state responses to infinitely long sinusoidal input signals, but this is definitely not the case. I think it's illustrative to look at an example:
Let's assume we have a causal real-valued LTI system with frequency response
H(ω)=A(ω)+jB(ω)
where A(ω) and B(ω) are the real and imaginary parts of H(ω), respectively. Now let's define an input signal to this system:
x(t)=u(t)sin(ω0t)
Its Fourier transform is
X(ω)=π2j[δ(ω−ω0)−δ(ω+ω0)]+ω0ω20−ω2
Using A(ω) and B(ω), the imaginary part of the Fourier transform of the output signal y(t) can be written as
ℑ{Y(ω)}=−π2[A(ω0)δ(ω−ω0)−A(−ω0)δ(ω+ω0)]−ω0B(ω)ω2−ω20
Since y(t)=0 for t<0 (because this is true for x(t) and because the system is causal), it is fully described by the imaginary part of its Fourier transform:
y(t)=−2π∫∞0ℑ{Y(ω)}sinωtdω,t>0
From (1) and (2) we get
y(t)=A(ω0)sinω0t+2ω0π∫∞0B(ω)ω2−ω20sinωtdω,t>0
Equation (3) fully describes the response of the LTI system to the causal input signal x(t)=u(t)sinω0t, including transients. This example is supposed to show that the frequency response (or the impulse response) is sufficient for a complete description of the system's response to arbitrary input signals.
EDIT: Derivation of Equation (2):
Let Y(ω)=YR(ω)+jYI(ω) be the Fourier transform of y(t). We need the following basic properties of the Fourier transform (see here):
y(t) real-valued⟹Y(ω)=Y∗(−ω)
If ye(t)=0.5[y(t)+y(−t)] is the even part of y(t) and yo(t)=0.5[y(t)−y(−t)] is the odd part of y(t) then
ye(t)⟺YR(ω)yo(t)⟺jYI(ω)
If y(t) is real-valued the inverse transform can be written as
y(t)=12π∫∞−∞Y(ω)ejωtdω=12π∫∞−∞ℜ{Y(ω)ejωt}dω==12π∫∞−∞[YR(ω)cos(ωt)−YI(ω)sin(ωt)]dω==1π∫∞0[YR(ω)cos(ωt)−YI(ω)sin(ωt)]dω
where the last equality follows from the fact that the integrand is an even function of ω. Consequently we have
ye(t)=1π∫∞0YR(ω)cos(ωt)dωyo(t)=−1π∫∞0YI(ω)sin(ωt)dω
With these preliminaries we can finally state the result using the fact that for causal y(t) the following holds for t>0 (because y(−t)=0 for t>0): y(t)=2ye(t)=2yo(t),t>0
This means
y(t)=2π∫∞0YR(ω)cos(ωt)dω==−2π∫∞0YI(ω)sin(ωt)dω,t>0
The last equality is the same as Equation (2) above.
No comments:
Post a Comment