Nonhomogeneous Differential Equations — Question 10

PDF ↗

Question 10

On 0≤t≤10\le t\le 1, design a forcing gg for y″+2y′+y=g(t)y''+2y'+y=g(t) so that the trajectory has y(0)=0,y′(0)=0,y(1)=1,y′(1)=0.y(0)=0,\quad y'(0)=0,\qquad y(1)=1,\quad y'(1)=0. Restrict the desired trajectory to polynomials of degree at most three.

Tasks

  1. Find the unique polynomial meeting the four endpoint requirements.

  2. Derive the forcing gg that makes this polynomial an exact solution and verify all four requirements.

  3. Prove that this forcing changes sign exactly once on [0,1][0,1], locating the zero between 3/43/4 and 4/54/5. Can a nonnegative forcing realize this same cubic trajectory?

  4. Keep the designed forcing but perturb the initial data to y(0)=εy(0)=\varepsilon, y′(0)=0y\prime(0)=0. Find the exact position and slope errors at t=1t=1, and determine whether the target state is still reached when ε≠0\varepsilon\ne 0.

Original worksheet page 1: question and worked solution for 3-8-010
Show solutionHide solution

Question 10 – Solution

Strategy. Prescribe a trajectory and compute its forcing directly; then use a homogeneous difference equation to analyze initial-data error.

Step 1: Fit the cubic trajectory. Write y=c0+c1t+c2t2+c3t3y=c_0+c_1t+c_2t^2+c_3t^3. The initial conditions give c0=c1=0c_0=c_1=0. The final conditions give c2+c3=1c_2+c_3=1 and 2c2+3c3=02c_2+3c_3=0, so y*(t)=3t2−2t3.\boxed{y_*(t)=3t^2-2t^3.} The coefficient equations have a unique solution, including among all polynomials of lower degree.

Step 2: Compute the forcing. The derivatives are y*′=6t−6t2y_*'=6t-6t^2 and y*″=6−12ty_*''=6-12t. Thus g=y*″+2y*′+y*=6−9t2−2t3.\boxed{g=y_*''+2y_*'+y_*=6-9t^2-2t^3.} The derivative formulas give zero initial value and slope, and final value 3−2=13-2=1 with slope 6−6=06-6=0. The computed forcing verifies the equation identically.

Step 3: Locate the sign change. We have g(0)=6g(0)=6, g(1)=−5g(1)=-5, and g′=−18t−6t2<0g'=-18t-6t^2<0 for t>0t>0. Therefore there is exactly one zero. More precisely, g(3/4)=3/32>0,g(4/5)=−98/125<0,g(3/4)=3/32>0,\qquad g(4/5)=-98/125<0, so the zero lies strictly between 3/43/4 and 4/54/5. This same cubic fixes gg pointwise, so a nonnegative forcing cannot realize it on the whole interval. This conclusion concerns the specified cubic, not every possible trajectory.

Step 4: Propagate the initial error. The difference h=y−y*h=y-y_* satisfies h″+2h′+h=0h''+2h'+h=0, with h(0)=εh(0)=\varepsilon, h′(0)=0h'(0)=0. Hence h=ε(1+t)e−th=\varepsilon(1+t)e^{-t} and h′=−εte−th'=-\varepsilon te^{-t}. At the endpoint, y(1)−1=2ε/e,y′(1)=−ε/e.\boxed{y(1)-1=2\varepsilon/e,\qquad y'(1)=-\varepsilon/e.} For nonzero ε\varepsilon the target state is missed. The forcing is unchanged; the error evolves entirely through the homogeneous correction.

Original worksheet page 2: question and worked solution for 3-8-010

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