Saturday, November 30, 2019

Frequency response equivalent to characterise transient response (to sinusoidal inputs)?


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),0t< gives output: y(t)=(1+et)sin(ωt),0t<.


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-valuedY(ω)=Y(ω)

from which follows YR(ω)=YR(ω) and YI(ω)=YI(ω)


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(ω)

These relation can be found in the link above.


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

digital communications - Understanding the Matched Filter

I have a question about matched filtering. Does the matched filter maximise the SNR at the moment of decision only? As far as I understand, ...