Systems Modeling - Sinusoidal Inputs

KingDaniel
Messages
44
Reaction score
1

Homework Statement


Finding the complete response (steady-state and transient) is a long and laborious task. My lecturer's notes read (since at our stage of the course, we're mostly interested in the steady-state part of the solution and not so much the transient) :

"The simple method for finding the steady-state part of the response to a sinusoidal input is simply to use the imaginary part of "s", substituting "jω" in place of "s" in the transfer function".

Then he goes on to show how to get the magnitude of the transfer function, G(s) / G(jω), and then on to get the phase.

Please explain what the magnitude of the transfer function has to do with the steady-state part of the solution, yss(t)?

Homework Equations

The Attempt at a Solution


I have no attempt to the solution as I just know that the transfer function, G(s), is y(s)/x(s).
I don't see how to get yss(t) from G(s).
 
Physics news on Phys.org
Consider a sinusoidal input x(t) = x_m \cos(\omega t) u(t) that produces a steady state output y_{ss}(t) and a network that has the transfer function H(s). Then the Laplace transform of the output is: Y(s) = H(s) L \left ( x_m \cos(\omega t) u(t) \right ) = H(s) L \left ( \left ( \frac{x_m}{2} e^{j \omega t} + \frac{x_m}{2} e^{-j \omega t}\right )u(t) \right )
So
Y(s) = \frac{x_m}{2} H(s) \frac{1}{s - j \omega} + \frac{x_m}{2} H(s) \frac{1}{s + j \omega}.

I was taught this in the context of electrical systems. My professor then went on to say that the poles of H(s) influence the transient solution, while the poles of \frac{1}{s \pm j \omega} make up the steady state solution. This feels kind of hand-wavy, but it makes sense when you consider that the transient solution is influenced by the particular system, but the steady state solution (long term) is obviously influenced by the driving force (for instance, consider pushing a spring with a sinusoidally varying force). In this case, if the driving force is a sinusoid, then the laplace transform of the driving force will look like the terms above. So we have Y(s) = \frac{k_1}{s - j \omega} + \frac{k_2}{s + j \omega} + (terms \ involving \ poles \ of \ H(s)) which implies that k_1 = Y(s)(s - j \omega) evaluated at s = j \omega, and by extension, k_2 = k_1^*, or the complex conjugate of k_1. So that means k_1 = \frac{x_m}{2} H(j \omega), k_2 = k_{1}^*.

Finally, you have to remember this little shortcut when taking inverse Laplace transforms: if Y(s) = \frac{k_1}{s + a - i b} + \frac{k_2}{s + a + i b}, then y(t) = 2 |k_1| e^{-a t} \cos(b t + \theta). In this case, a = 0, and \theta is the angle of the complex number k_1.

Putting those things together, when we take the inverse Laplace transform of the equation above, we get y(t) = x_m |H(j \omega)| \cos(\omega t + \theta) + (transient \ terms)
Or, ignoring the transient terms, y_{ss}(t) = x_m |H(j \omega)| \cos(\omega t + \theta) where \theta is the angle of H(j \omega)
 
Thanks @axmls , still seems like a pretty long and laborious task. But basically, the magnitude of the transfer function enables us to get the amplitude of yss(t) and the angle of the transfer function is the phase difference between the input and the output?
 
Hello everyone, I’m considering a point charge q that oscillates harmonically about the origin along the z-axis, e.g. $$z_{q}(t)= A\sin(wt)$$ In a strongly simplified / quasi-instantaneous approximation I ignore retardation and take the electric field at the position ##r=(x,y,z)## simply to be the “Coulomb field at the charge’s instantaneous position”: $$E(r,t)=\frac{q}{4\pi\varepsilon_{0}}\frac{r-r_{q}(t)}{||r-r_{q}(t)||^{3}}$$ with $$r_{q}(t)=(0,0,z_{q}(t))$$ (I’m aware this isn’t...
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
Back
Top