Use of the delta function

  • I
  • Thread starter CAF123
  • Start date
  • #1
CAF123
Gold Member
2,915
88
I've come across the equation $$\int_0^1 dx \frac{dA(x)}{dx} + B = C = \text{finite}$$ in my readings on a certain topic in physics and, in both articles i have read, the following step is taken $$\int_0^1 dx \left( \frac{dA(x)}{dx} + (B-C)\delta(1-x) \right) = \text{finite}$$

For the purposes of my question, I don't think the explicit form of the functions A,B,C are of concern. The reason ##\delta(1-x)## is chosen is associated with the underlying physics which I also think is not important. It seems to me that the replacement $$1 = \int_0^1 \delta(1-x) dx$$ has been made to go from line 1 to line 2 in the above. My questions are:

1) Why is this replacement even correct? The zero of the delta function occurs at the edge of the interval so I don't see why the integral should evaluate to one.
2) I put this integral into wolfram alpha. It returns ##\int_0^1 dx \delta(1-x) = \theta(0)##, where ##\theta(0)## is the step function evaluated at zero. Upon subsequently putting ##\theta(0)## into wolfram alpha it returns 1. So a) how is ##\theta(0)## obtained analytically and b) why is ##\theta(0) = 1##?

Thanks!
 

Answers and Replies

  • #2
Svein
Science Advisor
Insights Author
2,109
677
  1. Integration a delta function is equivalent to an evaluation. Thus [itex]\int_{a}^{b}\delta(x_{0}) dx [/itex] = 1 if a ≤ x0 ≤ b and 0 otherwise.
  2. In the same way [itex]\int_{a}^{b}f(x)\delta(x_{0})dx [/itex] =f(x0) if a ≤ x0 ≤ b and 0 otherwise.
 
  • #3
CAF123
Gold Member
2,915
88
  1. Integration a delta function is equivalent to an evaluation. Thus [itex]\int_{a}^{b}\delta(x_{0}) dx [/itex] = 1 if a ≤ x0 ≤ b and 0 otherwise.
  2. In the same way [itex]\int_{a}^{b}f(x)\delta(x_{0})dx [/itex] =f(x0) if a ≤ x0 ≤ b and 0 otherwise.
Ok, so then what is wrong with the following calculation: $$\int_0^1 \delta(1-x) dx = \int_0^1 \theta(x) \delta(1-x) dx = \theta(1) = 1?$$
(In the first equality I insert a theta function since this equals one for x >0 so I inserted one). The reason I think it is incorrect is because wolfram alpha gives me ##\theta(0)## as the result.
 
  • #4
stevendaryl
Staff Emeritus
Science Advisor
Insights Author
8,401
2,581
  1. Integration a delta function is equivalent to an evaluation. Thus [itex]\int_{a}^{b}\delta(x_{0}) dx [/itex] = 1 if a ≤ x0 ≤ b and 0 otherwise.
  2. In the same way [itex]\int_{a}^{b}f(x)\delta(x_{0})dx [/itex] =f(x0) if a ≤ x0 ≤ b and 0 otherwise.

I think your rules aren't quite right. I think it should be something like (assuming [itex]a < b[/itex]):

[itex]\int_a^b \delta(x-x_0) dx[/itex]
  • ... [itex]=1[/itex] if [itex]a < x_0 < b[/itex]
  • ...[itex]=0[/itex] if [itex]a < b < x_0[/itex] or [itex]x_0 < a < b[/itex]
  • ...[itex]=\frac{1}{2}[/itex] if [itex]x_0 = a[/itex] or [itex]x_0 = b[/itex]
You have to have the last case if you want [itex]\int_a^b \delta(x-x_0) dx + \int_b^c \delta(x-x_0) dx = \int_a^c \delta(x-x_0) dx[/itex]
 
  • #5
stevendaryl
Staff Emeritus
Science Advisor
Insights Author
8,401
2,581
Ok, so then what is wrong with the following calculation: $$\int_0^1 \delta(1-x) dx = \int_0^1 \theta(x) \delta(1-x) dx = \theta(1) = 1?$$
(In the first equality I insert a theta function since this equals one for x >0 so I inserted one). The reason I think it is incorrect is because wolfram alpha gives me ##\theta(0)## as the result.
It's tricky working with discontinuous functions (and the delta function isn't actually even a function). But I would say that:

[itex]\int_0^1 \delta(1-x) dx = \int_0^1 \delta(x) = \frac{1}{2}[/itex]

That's the most sensible answer, since
[itex]\int_a^1 \delta(x) dx = 1[/itex] if [itex]a < 0[/itex]
[itex]\int_a^1 \delta(x) dx = 0[/itex] if [itex]a > 0[/itex]

So when [itex]a=0[/itex], we should have:
[itex]\int_0^1 \delta(x) dx = \frac{1}{2}[/itex]

As I said, you need this if you want it to be the case that [itex]\int_{-1}^1 \delta(x) dx =\int_{-1}^0 \delta(x) dx + \int_{0}^1 \delta(x) dx[/itex]
 
  • #6
CAF123
Gold Member
2,915
88
I agree it seems like a plausible definition as you write it but then how to explain

a) the validity of the insertion of the delta from line 1 to 2 in my OP?

b) the fact wolfram alpha gives ##\int_0^1 \delta (1-x) dx = \theta(0) ## which it then, after inputting ##\theta(0)##, says is 1?
 
  • #7
Svein
Science Advisor
Insights Author
2,109
677
o:)Sorry, an error crept into my definitions. The correct definitions should be
  1. Integration a delta function is equivalent to an evaluation. Thus [itex]\int_{a}^{b}\delta(x - x_{0}) dx [/itex] = 1 if a ≤ x0 ≤ b and 0 otherwise.
  2. In the same way [itex]\int_{a}^{b}f(x)\delta(x - x_{0})dx [/itex] =f(x0) if a ≤ x0 ≤ b and 0 otherwise.
 
  • #8
CAF123
Gold Member
2,915
88
Unfortunately I am more confused now than before I posted my question - @Svein and @stevendaryl have given different (as far as I can see) opinions on how to use the delta 'function' when the zero of its argument coincides with one of the endpoints.

Applying what @Svein says to @stevendaryl suggestion of rewriting ##\int_{-1}^{1} \delta(x) dx = \int_{-1}^0 \delta(x) dx + \int_0^1 \delta(x) dx## gives me 1 on the lhs and 2 on the rhs (!) so either: a) we give up the assertion that ##\int_{-1}^{1} \delta(x) dx = \int_{-1}^0 \delta(x) dx + \int_0^1 \delta(x) dx## (perhaps ok if the delta function is not considered to be continuous at 0) or b) the equation ##\int_a^b \delta(x-x_0) dx = 1## if ##a \leq x_0 \leq b##, with ##x_o## strictly being included in the endpoints of the integration interval is wrong. (i.e I would agree that ##\int_a^b \delta(x-x_0) dx = 1## if ##a < x_0 < b## is certainly true).

I understand all these 'fallacies' probably come about from sloppiness in the use of the delta 'function' as a function but hopefully there can be a reconcilation about what has been said.

All my questions in the OP still remain.

Thanks!
 
  • #9
stevendaryl
Staff Emeritus
Science Advisor
Insights Author
8,401
2,581
I agree it seems like a plausible definition as you write it but then how to explain

a) the validity of the insertion of the delta from line 1 to 2 in my OP?

b) the fact wolfram alpha gives ##\int_0^1 \delta (1-x) dx = \theta(0) ## which it then, after inputting ##\theta(0)##, says is 1?
In my opinion, it's more convenient to use
[itex]\theta(x) = 0[/itex] when [itex]x < 0[/itex]
[itex]\theta(x) =\frac{1}{2}[/itex] when [itex]x = 0[/itex]
[itex]\theta(x) = 1[/itex] when [itex]x > 0[/itex]

On the other hand, delta functions and step functions are introduced to make certain calculations easier. The edge cases should be analyzed separately. I think you have to go back to the original problem to see what makes sense.
 
  • Like
Likes jasonRF
  • #10
Stephen Tashi
Science Advisor
7,409
1,381
2) I put this integral into wolfram alpha. It returns ##\int_0^1 dx \delta(1-x) = \theta(0)##, where ##\theta(0)## is the step function evaluated at zero.
Which Wolfram step function did you use?
The Wolfram Language represents the Heaviside generalized function as HeavisideTheta, while using UnitStep to represent the piecewise function Piecewise[[PLAIN]http://mathworld.wolfram.com/images/equations/HeavisideStepFunction/Inline7.gif[PLAIN]http://mathworld.wolfram.com/images/equations/HeavisideStepFunction/Inline8.gif1, [Broken] x >= 0[PLAIN]http://mathworld.wolfram.com/images/equations/HeavisideStepFunction/Inline9.gif[PLAIN]http://mathworld.wolfram.com/images/equations/HeavisideStepFunction/Inline10.gif] [Broken] (which, it should be noted, adopts the convention http://mathworld.wolfram.com/images/equations/HeavisideStepFunction/Inline11.gif instead of the conventional definition [PLAIN]http://mathworld.wolfram.com/images/equations/HeavisideStepFunction/Inline12.gif). [Broken]
 
Last edited by a moderator:
  • #11
CAF123
Gold Member
2,915
88
@Stephen Tashi Oh I see thanks, I didn't know of the difference - I wrote in 'integral from 0 to 1 of delta(1-x)' into wolfram alpha and it gives me theta(0). Next to the output it tells me '##\theta(x)## is the Heaviside step function' but when I click on that to view its properties it takes me to unitstep so I think it is indeed using the Piecewise ##[\left\{ \left\{ 1, x \geq 0\right\} \right\}]## definition, so that ##\theta(0) = 1##. But, if I understand correctly, I could take ##\theta(0) = 1/2## in accordance with what @stevendaryl was saying - so it is just a matter of choice and wolfram just happens to choose one convention?

In any case, how did wolfram derive that ##\int_0^1 dx \delta(1-x) = \theta(0)##?
 
  • #12
Stephen Tashi
Science Advisor
7,409
1,381
In any case, how did wolfram derive that ##\int_0^1 dx \delta(1-x) = \theta(0)##?
I don't use Wolfram alpha, so I don't know. We can also ask why it bothered to express the answer as "##\theta(0)##" instead of just giving a numerical result. Does Wolfram alpha have a feature like "Show me the steps you used"?

A humorous speculation is that it did a change of variable: ##\int_0^1 1 \delta(1-x) dx = \int_1^0 - 1\delta(y) dy = \int_0^1 1 \delta(y)dy ##
 
  • #13
FactChecker
Science Advisor
Gold Member
5,790
2,157
If the delta function is interpreted as a measure at 0, then I think that any integral range which includes 0 should give the full integral value of 1.
 
  • #14
stevendaryl
Staff Emeritus
Science Advisor
Insights Author
8,401
2,581
If the delta function is interpreted as a measure at 0, then I think that any integral range which includes 0 should give the full integral value of 1.
As I said, that would imply that [itex]\int_{-1}^{0} \delta(x) dx + \int_{0}^{1} \delta(x) dx = 2[/itex], while the usual rules for integration would say [itex]\int_{-1}^{0} \delta(x) dx + \int_{0}^{1} \delta(x) dx = \int_{-1}^{+1} \delta(x) dx = 1[/itex]
 
  • Like
Likes FactChecker
  • #15
FactChecker
Science Advisor
Gold Member
5,790
2,157
As I said, that would imply that [itex]\int_{-1}^{0} \delta(x) dx + \int_{0}^{1} \delta(x) dx = 2[/itex], while the usual rules for integration would say [itex]\int_{-1}^{0} \delta(x) dx + \int_{0}^{1} \delta(x) dx = \int_{-1}^{+1} \delta(x) dx = 1[/itex]
Good point. I missed that in your earlier post. But if the measure definition of delta function and Lebesgue integrals are used, including 0 in both integrals should sum to 2.
 
  • #16
CAF123
Gold Member
2,915
88
Hi FactChecker,
Good point. I missed that in your earlier post. But if the measure definition of delta function and Lebesgue integrals are used, including 0 in both integrals should sum to 2.
Could you elaborate on what you mean here?
 
  • #17
FactChecker
Science Advisor
Gold Member
5,790
2,157
Hi FactChecker,


Could you elaborate on what you mean here?
Using general measure theory and Lebesgue integration, you can make a measure that integrates to 1 at the single point x=0. But then if you include 0 in both integrals ∫-10 and ∫0+1, then you have 2 instead of 1. You would need to only include 0 in one integral and not in the other ( like ∫-1-0 + ∫0+1 where -0 means the limit as the integral interval upper value approaches 0 from below.)
 
  • Like
Likes CAF123
  • #18
CAF123
Gold Member
2,915
88
I see, so do you mean if we can make a measure so that ##\int_0^1 dx \delta(x) = 1## then we should define ##\int_{-1}^1 dx \delta(x) = 1 = \int_{-1}^{-0} dx \delta(x) + \int_0^1 dx \delta(x) = 0 + 1?## (I haven't done any measure theory so I just want to see if this is a correct take on things)

Also, going back to the claim in the OP where they go from
$$\int_0^1 dx \frac{dA(x)}{dx} + B = C = \text{finite}$$
to
$$\int_0^1 dx \left( \frac{dA(x)}{dx} + (B-C)\delta(1-x) \right) = \text{finite}$$
I guess all the author has done is defined the measure ##dx## so that ##\int_0^1 dx \delta(1-x)=1##?

Thanks!
 
  • #19
pwsnafu
Science Advisor
1,080
85
I see, so do you mean if we can make a measure so that ##\int_0^1 dx \delta(x) = 1## then we should define ##\int_{-1}^1 dx \delta(x) = 1 = \int_{-1}^{-0} dx \delta(x) + \int_0^1 dx \delta(x) = 0 + 1?## (I haven't done any measure theory so I just want to see if this is a correct take on things)
The Dirac measure is defined so that if a subset ##A \subset \mathbb{R}## contains the point 0 then it has measure 1. Otherwise the measure is 0. I'm going to write this as ##\delta(A)##. Now the Lebesgue integral in our case is (without going into details) ##\int_A f \, d\delta = f(0)\delta(A)##.

Let ##A = [-1,1]##. Lebesgue integration allows us to spit up ##A## but only if the union is disjoint. We can split this into either ##[-1,0] \cup (0,1]## or ##[-1,0) \cup [0,1]## but not ##[-1,0]\cup[0,1]##.

This is an elegant way to define the Dirac, but comes with the downside of the derivatives of the Dirac delta not compadible with "derivatives" as understood in measure theory (known as Radon-Nikodym derivatives).
 
Last edited:
  • Like
Likes FactChecker
  • #20
CAF123
Gold Member
2,915
88
Hi stevendaryl, (I'm still dabbling over this :P)
As I said, that would imply that [itex]\int_{-1}^{0} \delta(x) dx + \int_{0}^{1} \delta(x) dx = 2[/itex], while the usual rules for integration would say [itex]\int_{-1}^{0} \delta(x) dx + \int_{0}^{1} \delta(x) dx = \int_{-1}^{+1} \delta(x) dx = 1[/itex]
Why do we want to have that ##\int_{-1}^1 \delta(x) dx = \int_{-1}^0 \delta(x)dx + \int_0^1 \delta(x) dx## in the first place? The split at zero is where the delta function has its infinite spike so should we expect to be able to write this equation a priori? (it seems like a double counting in some sense)
 
  • #21
stevendaryl
Staff Emeritus
Science Advisor
Insights Author
8,401
2,581
Hi stevendaryl, (I'm still dabbling over this :P)

Why do we want to have that ##\int_{-1}^1 \delta(x) dx = \int_{-1}^0 \delta(x)dx + \int_0^1 \delta(x) dx## in the first place? The split at zero is where the delta function has its infinite spike so should we expect to be able to write this equation a priori? (it seems like a double counting in some sense)
Well, it's just a convenience for working with definite integrals, and is implicit in the use of endpoints, in the first place. If you don't demand that the integral has this property, then you would need to distinguish between an integral over an open interval and an integral over a closed interval. You could certainly do that, but it's a mess.

The real fact is that the delta function is not actually a function, and the various ways of manipulating it are a way of pretending that it is a function without getting into too much trouble.
 
  • Like
Likes FactChecker

Related Threads on Use of the delta function

Replies
1
Views
2K
Replies
3
Views
673
  • Last Post
Replies
4
Views
2K
  • Last Post
Replies
15
Views
2K
Replies
1
Views
1K
Replies
15
Views
6K
Replies
2
Views
8K
Replies
3
Views
957
Replies
9
Views
9K
  • Last Post
Replies
1
Views
1K
Top