Second Order Unit Step Response

Second order Unit Step Response 1. Unit Step Response We will use the example of an undamped harmonic oscillator with in­ put f (t) modeled by .. mx + kx = f (t). The unit step response is the solution to this equation with input u(t) and rest initial conditions x (t) = 0 for t < 0. That is, it is the solution to the initial value problem (IVP) .. mx + kx = u(t), . x (0− ) = 0, x (0− ) = 0. This could be an undamped spring-mass system with mass m and spring constant k. The mass is at rest at equilibrium until time t = 0 when a steady force starts to act on it. Force represents a change in momentum over time. A finite force F (t) can only cause an ininitesimal change in momentum (i.e. F (t) dt) at a time. Therefore, the mass does not change position abruptly, nor does it change velocity instantaneously. Because of this we should expect a solution which is continuous with continuous derivative. Only the acceleration experi­ ences a discontinuity. For t < 0 we already know that x (t) = 0. For t > 0 the DE is .. mx + kx = 1. This has a constant particular solution x (t) = 1/k, and a general homoge­ neous solution √ xh (t) = c1 cos(ωn t) + c2 sin(ωn t), where ωn = k/m. Putting the two together gives the general solution x (t) = 1/k + c1 cos(ωn t) + c2 sin(ωn t) for t > 0. . . . The continuity of x and x implies x (0) = x (0− ) = 0 and x (0) = x (0− ) = 0. This allows us to find c1 and c2 . 0 = x (0) = 1/k + c1 . 0 = x (0) = c2 ω n ⇒ c1 = −1/k ⇒ c2 = 0. Second order Unit Step Response OCW 18.03SC The unit step response for this system is (in both cases and u-format) � 0 for t < 0 1 x (t) = 1 = (1 − cos(ωn t))u(t). k for t > 0. k (1 − cos( ωn t )) As in the first order case, we will sometimes denote this v(t). The claim that we get a continuous response is true, but may feel a bit unjustified. Let’s redo the above example very carefully without making this assumption. It will take more work, but we will get the same answer. In cases format the equation for the IVP is � .. . 0 for t < 0 mx + kx = x (0− ) = 0, x (0− ) = 0. 1 for t > 0, Solving the two pieces we get � c1 cos(ωn t) + c2 sin(ωn t) x (t) = 1/k + c3 cos(ωn t) + c4 sin(ωn t) (1) for t < 0 for t > 0. . The pre-initial conditions x (0− ) = x (0− ) = 0 easily imply c1 = c2 = 0. So our solution looks like � 0 for t < 0 x (t) = 1/k + c3 cos(ωn t) + c4 sin(ωn t) for t > 0. To find c3 and c4 we substitute x (t) into equation (1). . To measure the jumps we compute x (0+ ) = 1/k + c3 and x (0+ ) = c4 ωn . We use this as we compute derivatives of x. � 0 for t < 0 . x (t) = (1/k + c3 )δ(t) + −c3 ωn sin(�ωn t) + c4 ωn cos(ωn t) for t > 0. .. 0 −c3 ωn2 cos(ωn t) − c4 ωn2 sin(ωn t) x (t) = (1/k + c3 )δ� (t) + c4 ωn δ(t) + for t < 0 for t > 0. Since the right-hand side of equation (1) does not have any delta functions .. or any δ� (t) the coefficients in front of these terms in the formula for x must be 0: 1/k + c3 = 0 c4 ω n = 0 ⇒ c3 = −1/k ⇒ c4 = 0. 2 Second order Unit Step Response OCW 18.03SC In the end, we have exactly the same solution as above for the unit step response. To summarize: the continuity assumptions follow because any jumps . in x (t) or x (t) would result in delta functions when x is substituted into equation (1). The generalized derivative δ� (t) is not something we’ve seen before. It is often called a doublet. There is an entire theory of these and other gener­ alized functions, but we will only use δ(t) in this course. Figure 1 shows the graph of the unit step response (with k = 1 and m = 0.5. 2/k t .. Figure 1. The unit step response for the system mx + kx = u(t). If we added some damping the homogeneous part of the solution would go to 0 and the unit step response would go asymptotically to 1/k. The Meaning of the Phrase ’Unit Step Response’ As we noted in the first order case, the unit step response is the response of the system to a unit step input. For example, if our system is .. . . mx + bx + kx = by and we consider y to be the input, then the unit step response is the solution to .. . . .. . mx + bx + kx = bu(t) equivalently mx + bx + kx = bδ(t). 3