Unit 11, lesson 62

Boundary-Value Problems and Solvability

Why a linear boundary-value problem has no solution, one or infinitely many, and how the boundary map, the kernel of the homogeneous problem, resonance and a compatibility condition decide unique solvability.

  • Video 62
  • Solvability and Eigenfunctions

A second-order equation leaves two constants free, and a boundary-value problem supplies two conditions, so it seems the conditions should always pick out one solution. Yet the problem y′′+y=sin⁡xy'' + y = \sin x on [0,π][0, \pi] with y(0)=y(π)=0y(0) = y(\pi) = 0 has no solution at all, and replacing the forcing by cos⁡x\cos x gives infinitely many. Here xx is the independent variable and y′=dy/dxy' = dy/dx. This lesson writes the boundary conditions as a linear system for the free constants, shows how the kernel of the homogeneous problem separates the unique case from the resonant one, and derives the condition the forcing has to meet in the resonant case.

Key terms

Solvability of boundary-value problems

Definition Let pp, qq and ff be continuous on a closed interval [a,b][a, b], and write

L[y]=y′′+p(x) y′+q(x) y.L[y] = y'' + p(x)\,y' + q(x)\,y.

A linear two-point boundary-value problem asks for a twice differentiable function yy on [a,b][a, b] with

L[y]=f  on [a,b],B1[y]=γ1,B2[y]=γ2,\begin{gathered} L[y] = f \ \text{ on } [a, b], \\[0.4em] B_1[y] = \gamma_1, \qquad B_2[y] = \gamma_2, \end{gathered}

where each Bi[y]B_i[y] is a fixed linear combination of the four numbers y(a)y(a), y′(a)y'(a), y(b)y(b) and y′(b)y'(b), and γ1\gamma_1 and γ2\gamma_2 are given numbers. Solvability asks whether a solution exists and, if so, how many there are. Such a problem always has no solution, exactly one or infinitely many. The reason is that every solution of L[y]=fL[y] = f has the form y=yp+c1y1+c2y2y = y_p + c_1y_1 + c_2y_2, where ypy_p is one solution and y1y_1, y2y_2 form a fundamental set of L[y]=0L[y] = 0, so the two boundary conditions become two linear equations for c1c_1 and c2c_2. The classification needs both the equation and the boundary conditions to be linear. A nonlinear problem can have, for example, exactly two solutions.

In plain words The equation leaves two constants free, and the boundary conditions are two linear equations for them. Two linear equations in two unknowns can have no solution, exactly one or infinitely many, and the boundary-value problem inherits exactly these three outcomes.

Example The equation y′′=0y'' = 0 has the general solution y=c1+c2xy = c_1 + c_2x. With y(0)=αy(0) = \alpha and y(1)=βy(1) = \beta the conditions read c1=αc_1 = \alpha and c1+c2=βc_1 + c_2 = \beta, so there is exactly one solution, y=α+(β−α)xy = \alpha + (\beta - \alpha)x. With y′(0)=0y'(0) = 0 and y′(1)=γy'(1) = \gamma both conditions concern c2c_2 alone. They read c2=0c_2 = 0 and c2=γc_2 = \gamma, so for γ≠0\gamma \ne 0 there is no solution, and for γ=0\gamma = 0 every constant function is a solution. Now keep the equation but make one condition nonlinear, y(0)2=1y(0)^2 = 1 and y(1)=0y(1) = 0. The conditions give c1=±1c_1 = \pm 1 and c1+c2=0c_1 + c_2 = 0, so the solutions are exactly y=1−xy = 1 - x and y=x−1y = x - 1. Exactly two solutions can never happen in a linear problem.

Linear boundary map

Definition Let y1y_1 and y2y_2 be a fundamental set of solutions of L[y]=0L[y] = 0. The linear boundary map sends each homogeneous solution y=c1y1+c2y2y = c_1y_1 + c_2y_2 to its pair of boundary values (B1[y],B2[y])(B_1[y], B_2[y]). In terms of the coefficient vector c=(c1,c2)Tc = (c_1, c_2)^{T} it is multiplication by the matrix

M=(B1[y1]B1[y2]B2[y1]B2[y2]).M = \begin{pmatrix} B_1[y_1] & B_1[y_2] \\ B_2[y_1] & B_2[y_2] \end{pmatrix}.

Because B1B_1 and B2B_2 are linear, the general solution y=yp+c1y1+c2y2y = y_p + c_1y_1 + c_2y_2 of the forced equation meets the boundary conditions exactly when

Mc=d,d=(γ1−B1[yp]γ2−B2[yp]).\begin{gathered} Mc = d, \\[0.4em] d = \begin{pmatrix} \gamma_1 - B_1[y_p] \\ \gamma_2 - B_2[y_p] \end{pmatrix}. \end{gathered}

A different fundamental set multiplies MM on the right by an invertible matrix, so the rank of MM, and in particular whether MM is invertible, does not depend on that choice.

In plain words The boundary map records what the free part of the solution does at the two ends. If its matrix is invertible, any required pair of endpoint values can be reached, and in exactly one way.

Example For y′′+y=fy'' + y = f on [0,π][0, \pi] with y(0)=0y(0) = 0 and y(π)=0y(\pi) = 0, take y1=sin⁡xy_1 = \sin x and y2=cos⁡xy_2 = \cos x. The homogeneous solution Asin⁡x+Bcos⁡xA\sin x + B\cos x has the value BB at 00 and the value −B-B at π\pi, so

M=(010−1).M = \begin{pmatrix} 0 & 1 \\ 0 & -1 \end{pmatrix}.

Both rows constrain BB only, the rank is 11, and MM is singular. On the shorter interval [0,π/2][0, \pi/2] with y(0)=0y(0) = 0 and y(π/2)=0y(\pi/2) = 0, the two values are BB and AA, so

M=(0110),M = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix},

which is invertible, with determinant −1-1.

Kernel of the homogeneous problem

Definition The homogeneous problem keeps the operator and the boundary operators and sets the forcing and the data to zero,

L[y]=0,B1[y]=0,B2[y]=0.\begin{gathered} L[y] = 0, \\[0.4em] B_1[y] = 0, \qquad B_2[y] = 0. \end{gathered}

Its solutions form a vector space, called the kernel of the homogeneous problem. It always contains the zero function. A homogeneous solution c1y1+c2y2c_1y_1 + c_2y_2 lies in it exactly when Mc=0Mc = 0, so the kernel has dimension 2−rank⁡M2 - \operatorname{rank} M, which is 00, 11 or 22. If the forced problem has a solution y∗y^{*}, then its solutions are exactly the functions y∗+uy^{*} + u with uu in the kernel, because the difference of two solutions of the forced problem solves the homogeneous one.

In plain words The kernel collects the shapes that satisfy the equation without forcing and are invisible to the boundary conditions. Any of them can be added to a solution without changing the forcing or the endpoint data, so the kernel measures exactly how uniqueness fails.

Example For y′′+y=0y'' + y = 0 with y(0)=y(π)=0y(0) = y(\pi) = 0, the matrix above forces B=0B = 0 and leaves AA free, so the kernel consists of the functions Asin⁡xA\sin x and has dimension 11. On [0,π/2][0, \pi/2] the invertible matrix forces A=B=0A = B = 0, and the kernel is {0}\{0\}. With the periodic conditions y(0)=y(2π)y(0) = y(2\pi) and y′(0)=y′(2π)y'(0) = y'(2\pi) on [0,2π][0, 2\pi], every function Asin⁡x+Bcos⁡xA\sin x + B\cos x satisfies both conditions, the matrix is the zero matrix and the kernel has dimension 22.

Resonant forcing

Definition A linear boundary-value problem is resonant when its homogeneous problem has a nonzero solution, that is, when MM is singular. Forcing applied to a resonant problem is called resonant forcing. A resonant problem never has exactly one solution. It has none when the vector dd lies outside the column space of MM, and infinitely many when dd lies inside it, since every kernel function can then be added to a solution. Which of the two happens depends on the forcing and the boundary data, not on the operator alone.

In plain words A nonzero kernel function is a free mode that the boundary conditions cannot detect, much like the natural frequency of an undamped oscillator. Forcing can drive that mode in a way the boundaries cannot absorb, and then no solution exists. Resonance alone tells you only that a unique solution is impossible.

Example Take y′′+y=sin⁡xy'' + y = \sin x on [0,π][0, \pi] with y(0)=y(π)=0y(0) = y(\pi) = 0, a problem whose kernel contains sin⁡x\sin x. A particular solution is yp=−12xcos⁡xy_p = -\tfrac{1}{2}x\cos x, because

yp′=−12cos⁡x+12xsin⁡x,yp′′=sin⁡x+12xcos⁡x,\begin{aligned} y_p' &= -\tfrac{1}{2}\cos x + \tfrac{1}{2}x\sin x, \\[0.4em] y_p'' &= \sin x + \tfrac{1}{2}x\cos x, \end{aligned}

and so yp′′+yp=sin⁡xy_p'' + y_p = \sin x. The general solution of the equation is y=−12xcos⁡x+Asin⁡x+Bcos⁡xy = -\tfrac{1}{2}x\cos x + A\sin x + B\cos x. The left condition gives y(0)=B=0y(0) = B = 0. At the right end

y(π)=π2−B=π2,y(\pi) = \frac{\pi}{2} - B = \frac{\pi}{2},

and no choice of AA changes this, because sin⁡π=0\sin\pi = 0. So the boundary-value problem has no solution, although the differential equation alone has a two-parameter family of them. The growing factor xx in ypy_p is the same one that appears when an undamped oscillator is driven at its natural frequency.

Compatibility condition

Definition A compatibility condition is a condition on the forcing and the boundary data that is necessary and sufficient for a resonant problem to have a solution. In terms of the boundary map it says that dd lies in the column space of MM, which amounts to one scalar equation when the kernel has dimension 11 and two when it has dimension 22. For the problem on [0,π][0, \pi]

y′′+y=f(x),y(0)=0,y(π)=0,\begin{gathered} y'' + y = f(x), \\[0.4em] y(0) = 0, \qquad y(\pi) = 0, \end{gathered}

with ff continuous on [0,π][0, \pi], the compatibility condition is

∫0πf(x)sin⁡x dx=0.\int_0^\pi f(x)\sin x\,dx = 0.

The condition is necessary. Multiplying the equation by sin⁡x\sin x and integrating by parts twice gives

∫0πf(x)sin⁡x dx=∫0π(y′′+y)sin⁡x dx=[ y′sin⁡x−ycos⁡x ]0π,\begin{aligned} &\int_0^\pi f(x)\sin x\,dx \\[0.4em] &\quad = \int_0^\pi (y'' + y)\sin x\,dx \\[0.4em] &\quad = \Big[\, y'\sin x - y\cos x \,\Big]_0^\pi, \end{aligned}

and the bracket is zero because sin⁡0=sin⁡π=0\sin 0 = \sin\pi = 0 and y(0)=y(π)=0y(0) = y(\pi) = 0. The condition is also sufficient. By variation of parameters, the solution of the equation with y(0)=y′(0)=0y(0) = y'(0) = 0 is

y0(x)=∫0xsin⁡(x−t) f(t) dty_0(x) = \int_0^x \sin(x - t)\,f(t)\,dt

and every solution with y(0)=0y(0) = 0 has the form y0+Asin⁡xy_0 + A\sin x. Since sin⁡(π−t)=sin⁡t\sin(\pi - t) = \sin t, the value y0(π)y_0(\pi) equals the integral of fsin⁡xf\sin x, so when that integral vanishes every function y0+Asin⁡xy_0 + A\sin x solves the problem. The test says that ff is orthogonal to the kernel function sin⁡x\sin x. For self-adjoint problems with zero boundary data the same pattern holds in general, with orthogonality to every kernel function as the test, and this is the content of the Fredholm alternative later in the course.

In plain words Before solving anything, measure how much of the forcing points along the kernel function. A nonzero amount drives the free mode in a way the boundaries cannot absorb, and a zero amount lets solutions through.

Example For f=sin⁡xf = \sin x the test integral is

∫0πsin⁡2x dx=π2≠0,\int_0^\pi \sin^2 x\,dx = \frac{\pi}{2} \ne 0,

which confirms that there is no solution. For f=cos⁡xf = \cos x it is ∫0πsin⁡xcos⁡x dx=0\int_0^\pi \sin x\cos x\,dx = 0, so solutions exist. One of them is u=12xsin⁡xu = \tfrac{1}{2}x\sin x, since

u′=12sin⁡x+12xcos⁡x,u′′=cos⁡x−12xsin⁡x,\begin{aligned} u' &= \tfrac{1}{2}\sin x + \tfrac{1}{2}x\cos x, \\[0.4em] u'' &= \cos x - \tfrac{1}{2}x\sin x, \end{aligned}

so u′′+u=cos⁡xu'' + u = \cos x, and uu vanishes at 00 and at π\pi. Adding the kernel gives all solutions,

y=12xsin⁡x+Asin⁡x,y = \tfrac{1}{2}x\sin x + A\sin x,

one for each real AA. Compatibility restores existence, but it does not restore uniqueness. Nonzero boundary data enter the condition too. With y(0)=αy(0) = \alpha and y(π)=βy(\pi) = \beta the bracket above equals α+β\alpha + \beta, and the condition becomes

∫0πf(x)sin⁡x dx=α+β.\int_0^\pi f(x)\sin x\,dx = \alpha + \beta.

Unique solvability

Definition A linear boundary-value problem is uniquely solvable when it has exactly one solution for every forcing ff continuous on [a,b][a, b] and all boundary values γ1\gamma_1 and γ2\gamma_2. The following statements are equivalent.

  1. The problem is uniquely solvable.
  2. The boundary matrix MM is invertible, that is, det⁡M≠0\det M \ne 0.
  3. The homogeneous problem has only the zero solution.

Indeed, an invertible MM gives exactly one cc for every dd, unique solvability applied to zero forcing and zero data leaves only the zero function, and a zero kernel means Mc=0Mc = 0 only for c=0c = 0. So exactly one of two cases occurs. Either the kernel is {0}\{0\} and every forced problem has exactly one solution, or the kernel is nonzero and every forced problem has none or infinitely many, with both outcomes occurring for suitable data.

In plain words Test the problem with zero forcing and zero endpoint data. If only the zero function survives, every version of the problem has exactly one answer.

Example On [0,π/2][0, \pi/2] the problem y′′+y=1y'' + y = 1 with y(0)=y(π/2)=0y(0) = y(\pi/2) = 0 has the invertible boundary matrix found above, so it has exactly one solution. The constant 11 is a particular solution, and y=1+Asin⁡x+Bcos⁡xy = 1 + A\sin x + B\cos x gives 1+B=01 + B = 0 at 00 and 1+A=01 + A = 0 at π/2\pi/2, so

y=1−sin⁡x−cos⁡x.y = 1 - \sin x - \cos x.

Its second derivative is sin⁡x+cos⁡x\sin x + \cos x, so y′′+y=1y'' + y = 1, and the end values are 1−0−1=01 - 0 - 1 = 0 and 1−1−0=01 - 1 - 0 = 0. In a resonant problem an extra condition can single out one member of a family, but it changes the problem rather than making the original one uniquely solvable. For the family 12xsin⁡x+Asin⁡x\tfrac{1}{2}x\sin x + A\sin x above, imposing the extra requirement ∫0πysin⁡x dx=0\int_0^\pi y\sin x\,dx = 0 and using ∫0πxsin⁡2x dx=π2/4\int_0^\pi x\sin^2 x\,dx = \pi^2/4 gives

π28+π2A=0,A=−π4.\begin{gathered} \frac{\pi^2}{8} + \frac{\pi}{2}A = 0, \\[0.4em] A = -\frac{\pi}{4}. \end{gathered}

Check yourself

  1. Does y′′+y=xy'' + y = x on [0,π][0, \pi] with y(0)=y(π)=0y(0) = y(\pi) = 0 have a solution? Use the compatibility condition, then confirm the answer by solving directly.
Show the answer

This is the resonant problem of this page, so the test is the integral of xsin⁡xx\sin x. Integration by parts gives

∫0πxsin⁡x dx=[−xcos⁡x]0π+∫0πcos⁡x dx=π+0=π,\begin{aligned} \int_0^\pi x\sin x\,dx &= \big[-x\cos x\big]_0^\pi \\[0.4em] &\quad + \int_0^\pi \cos x\,dx \\[0.4em] &= \pi + 0 = \pi, \end{aligned}

which is not zero, so there is no solution. Directly, the function yp=xy_p = x is a particular solution, so every solution of the equation is y=x+Asin⁡x+Bcos⁡xy = x + A\sin x + B\cos x. The left condition gives B=0B = 0, and then y(π)=π−B=π≠0y(\pi) = \pi - B = \pi \ne 0 for every AA.

  1. A classmate says that y′′+y=sin⁡2xy'' + y = \sin 2x with y(0)=y(π)=0y(0) = y(\pi) = 0 has no solution, because the problem is resonant. Is that right? Find all solutions.
Show the answer

No. Resonance rules out exactly one solution, but it allows infinitely many. The compatibility test gives

∫0πsin⁡2xsin⁡x dx=∫0π2sin⁡2xcos⁡x dx=[23sin⁡3x]0π=0,\begin{aligned} &\int_0^\pi \sin 2x\sin x\,dx \\[0.4em] &\quad = \int_0^\pi 2\sin^2 x\cos x\,dx \\[0.4em] &\quad = \Big[\tfrac{2}{3}\sin^3 x\Big]_0^\pi = 0, \end{aligned}

so solutions exist. Trying yp=ksin⁡2xy_p = k\sin 2x gives yp′′+yp=−3ksin⁡2xy_p'' + y_p = -3k\sin 2x, so k=−13k = -\tfrac{1}{3}, and this ypy_p vanishes at 00 and at π\pi. Adding the kernel gives every solution,

y=−13sin⁡2x+Asin⁡x,y = -\tfrac{1}{3}\sin 2x + A\sin x,

for any real AA.

  1. For which b>0b > 0 does the problem y′′+y=fy'' + y = f with y(0)=0y(0) = 0 and y(b)=0y(b) = 0 have exactly one solution for every continuous ff on [0,b][0, b]? For b=2πb = 2\pi, state the compatibility condition.
Show the answer

With the basis sin⁡x\sin x, cos⁡x\cos x, the homogeneous solution Asin⁡x+Bcos⁡xA\sin x + B\cos x has the value BB at 00 and the value Asin⁡b+Bcos⁡bA\sin b + B\cos b at bb, so

M=(01sin⁡bcos⁡b),det⁡M=−sin⁡b.\begin{gathered} M = \begin{pmatrix} 0 & 1 \\ \sin b & \cos b \end{pmatrix}, \\[0.4em] \det M = -\sin b. \end{gathered}

The problem is uniquely solvable exactly when sin⁡b≠0\sin b \ne 0, that is, for every b>0b > 0 that is not an integer multiple of π\pi. For b=2πb = 2\pi the kernel consists of the functions Asin⁡xA\sin x. Every solution with y(0)=0y(0) = 0 is y0+Asin⁡xy_0 + A\sin x with y0y_0 as on this page, and since sin⁡(2π−t)=−sin⁡t\sin(2\pi - t) = -\sin t,

y(2π)=−∫02πf(t)sin⁡t dt.y(2\pi) = -\int_0^{2\pi} f(t)\sin t\,dt.

So the compatibility condition is ∫02πf(x)sin⁡x dx=0\int_0^{2\pi} f(x)\sin x\,dx = 0, and when it holds the solutions are y0+Asin⁡xy_0 + A\sin x for every real AA.

How this lesson connects