Dirac Delta Function — Question 7

PDF ↗

Question 7

Use causal one-sided Laplace transforms. Write Ha(t)=0H_a(t)=0 for t<at<a and 11 for t≥at\ge a. The unit impulse satisfies ∫ϕ(t)δ(t−a)dt=ϕ(a)\int\phi(t)\delta(t-a)\,dt=\phi(a) for continuous ϕ\phi near aa. Interpret equations between impulses and through their jump conditions; use right-hand values at jumps. Write [v]a=v(a+)−v(a−)[v]_a=v(a^+)-v(a^-) for a jump.

An unstable second-order system is given one corrective impulse: y″−y=Jδ(t−a),y(0)=1,y′(0)=0,a>0.y''-y=J\delta(t-a),\qquad y(0)=1,\quad y'(0)=0,\qquad a>0. Choose JJ to make the response decay to zero after the impulse. A realized impulse may differ by a constant error Δ\Delta.

Tasks

  1. Derive Y(s)Y(s) and y(t)y(t) for arbitrary JJ, including the initial response.

  2. Find the unique decaying choice J*J_*, its explicit tail and the global maximum of the selected response.

  3. Derive the response error for strength J*+ΔJ_*+\Delta and a sharp bound on |Δ||\Delta| ensuring error at most ϵ>0\epsilon>0 throughout a≤t≤a+La\le t\le a+L, where L>0L>0.

  4. Compare exact real transform domains for the selected and perturbed strengths. Explain the effect of cancellation on the system’s stability.

Original worksheet page 1: question and worked solution for 4-8-007
Show solutionHide solution

Question 7 – Solution

Strategy. Cancel the growing mode’s coefficient, then inspect the remaining mode and the error produced by imperfect cancellation.

Step 1: Transform and invert. For sufficiently large ss, (s2−1)Y−s=Je−as,Y=s+Je−ass2−1,y=cosh⁡t+JHa(t)sinh⁡(t−a).(s^2-1)Y-s=Je^{-as},\qquad Y=\frac{s+Je^{-as}}{s^2-1},\qquad \boxed{y=\cosh t+JH_a(t)\sinh(t-a).} The added term starts at zero with derivative JJ, so it preserves displacement and produces the correct velocity jump. Away from aa both terms solve the homogeneous equation.

Step 2: Select and interpret the decaying mode. After aa, the ete^t coefficient is (1+Je−a)/2(1+Je^{-a})/2. Decay requires J*=−ea,y*(t)={cosh⁡t,0≤t<a,cosh⁡ae−(t−a),t≥a.\boxed{J_*=-e^a,\qquad y_*(t)=\begin{cases}\cosh t,&0\le t<a,\\ \cosh a\,e^{-(t-a)},&t\ge a. \end{cases}} The pre-impact velocity is sinh⁡a\sinh a and the post-impact velocity is −cosh⁡a-\cosh a; their difference is −ea-e^a, as required. The response increases strictly before aa and decreases strictly after it. Its unique global maximum is y*(a)=cosh⁡a\boxed{y_*(a)=\cosh a}.

Step 3: Quantify the sensitivity. Linearity or direct subtraction gives y(t)−y*(t)=ΔHa(t)sinh⁡(t−a).\boxed{y(t)-y_*(t)=\Delta H_a(t)\sinh(t-a).} Because sinh⁡τ\sinh\tau increases strictly on [0,L][0,L], the desired error bound is equivalent to |Δ|≤ϵsinh⁡L.\boxed{|\Delta|\le\frac{\epsilon}{\sinh L}.} This is necessary as well as sufficient, and equality attains the bound at a+La+L. The selected solution still approaches zero only asymptotically; the impulse does not put it at the zero state.

Step 4: Check convergence and stability. The selected tail is a positive multiple of e−te^{-t}, so its exact real transform domain is s>−1\boxed{s>-1}. The apparent singularity at s=1s=1 in its transform is removable. For every Δ≠0\Delta\ne 0, the error contains a nonzero growing mode and the actual domain becomes s>1\boxed{s>1}. The differential equation’s unstable homogeneous mode is unchanged. The plot uses a=1a=1, with a dashed curve for Δ=0.05\Delta=0.05.

See the diagram in the original worksheet below.

Original worksheet page 2: question and worked solution for 4-8-007

Original worksheet layout. Use Enlarge or open the PDF for a closer view.