Zulip Chat Archive

Stream: maths

Topic: Non-injective change of variables in integrals


Yury G. Kudryashov (Oct 15 2023 at 20:24):

In #7693 I prove that f(x)dx=nVol(B1)0f(r)rn1dr\int f(\|x\|)\,dx=n\operatorname{Vol}(B_1)\int_0^{\infty}f(r)r^{n-1}\,dr, where xEx\in E, n=dimRE>0n=\dim_{\mathbb R} E>0, f ⁣:RFf\colon \mathbb R\to F, Vol(B1)\operatorname{Vol}(B_1) is the volume of a unit ball in EE.
I wanted to have this formula for any f, not necessarily integrable, so the proof works in two steps: transfer the measure from E{0}E\setminus \{0\} to (0,)×sphere(0,1)(0, \infty)\times\operatorname{sphere}(0, 1), then projection to the first coordinate.

Yury G. Kudryashov (Oct 15 2023 at 20:25):

It would be easier to prove that the push-forward of the Lebesgue measure on E under norm is the Lebesgue measure on real numbers with given density.

Yury G. Kudryashov (Oct 15 2023 at 20:25):

But this way the result will require f to be integrable.

Yury G. Kudryashov (Oct 15 2023 at 20:27):

@Sebastien Gouezel I know that you formalized a very advanced version of the change of variables formula. Is it true in case of a non-injective map (e.g., provided that we know that the push-forward in the Lebesgue measure with some density)?

Yury G. Kudryashov (Oct 15 2023 at 20:28):

If yes, how much different is the proof?

Geoffrey Irving (Oct 15 2023 at 20:49):

(deleted)

Sebastien Gouezel (Oct 16 2023 at 06:07):

If you already know the push-forward, then I guess the formula you want doesn't need anything more. docs#MeasureTheory.integral_map, for instance?

Yury G. Kudryashov (Oct 16 2023 at 16:54):

This lemma assumes that f is AEStronglyMeasurable.

Yury G. Kudryashov (Oct 16 2023 at 16:55):

We have a version without this assumption if φ is a measurable embedding. It is also true for Prod.fst, Prod.snd, and norm.

Yury G. Kudryashov (Oct 16 2023 at 16:57):

I wonder if there is a class of maps that generalizes Prod.fst/Prod.snd/norm. The only difficulty is to prove that AEStronglyMeasurable (f ∘ φ) implies AEStronglyMeasurable f under some assumptions on φ.

Sebastien Gouezel (Oct 16 2023 at 18:10):

Ah, sorry, I had missed the difficulty. Indeed, that's a good question!

Anatole Dedecker (Oct 16 2023 at 20:09):

"Has a measurable section" is the most natural thing that comes to mind and it seems to encompass your three use cases (assuming you mean nnnorm, but I don’t think the result could be true for norm)

Yury G. Kudryashov (Oct 17 2023 at 02:20):

How does a measurable section help?

Yury G. Kudryashov (Oct 17 2023 at 02:20):

Note that we need a statement about NullMeasurableSets, not MeasurableSets.

Yury G. Kudryashov (Oct 17 2023 at 02:21):

If we drop f, then the lemma we need is NullMeasurableSet (φ ⁻¹' s) μ <-> NullMeasurableSet s (μ.map φ).

Yury G. Kudryashov (Oct 17 2023 at 02:21):

(unless I miss something; I didn't formalize the implication)

Felix Weilacher (Oct 17 2023 at 03:59):

Is this different than what was being discussed in this thread?

Felix Weilacher (Oct 17 2023 at 04:00):

For example we showed there that if phi is a measurable map between standard Borel spaces with Borel image, then it is a quotient map for the null-measurable algebras

Felix Weilacher (Oct 17 2023 at 04:10):

And actually I think any measurable map will work.

Yury G. Kudryashov (Oct 17 2023 at 05:41):

Sorry, I forgot that I've already asked this question. I'll read the proof again tomorrow (have a test to grade ASAP).

Xavier Roblot (Oct 17 2023 at 08:33):

Yury G. Kudryashov said:

In #7693 I prove that f(x)dx=nVol(B1)0f(r)rn1dr\int f(\|x\|)\,dx=n\operatorname{Vol}(B_1)\int_0^{\infty}f(r)r^{n-1}\,dr, where xEx\in E, n=dimRE>0n=\dim_{\mathbb R} E>0, f ⁣:RFf\colon \mathbb R\to F, Vol(B1)\operatorname{Vol}(B_1) is the volume of a unit ball in EE.
I wanted to have this formula for any f, not necessarily integrable, so the proof works in two steps: transfer the measure from E{0}E\setminus \{0\} to (0,)×sphere(0,1)(0, \infty)\times\operatorname{sphere}(0, 1), then projection to the first coordinate.

Lately, I have been trying to prove a similar formula (which much less success than you :sweat_smile: ) since it provides a way to compute volume of unit balls in Rn\mathbb{R}^n (and also Cn\mathbb{C}^n) for the norm LpL_p (and also in mixed setting which is the one I really need). Is this something that you are planning to add?

Yury G. Kudryashov (Oct 17 2023 at 16:23):

Could you please give more details?

Yury G. Kudryashov (Oct 17 2023 at 16:24):

I mean, how do you want to compute the volume of the unit ball?

Xavier Roblot (Oct 17 2023 at 19:54):

To compute the volume of the unit ball for the norm p\|\cdot\|_p, you take f(x)=expf(x) = e ^ {-x ^p}, so you get f(xp)dx=(Rexpdx)n\int f(\|x\|_p) dx = (\int_\mathbb{R} e ^ {- |x| ^ p} dx) ^ n where nn is the dimension. Then, you can compute all the integrals (that are equal to values of the Gamma function) and solve for Vol(B1)\mathrm{Vol}(B_1).

I think that works but I might be missing something.

Yury G. Kudryashov (Oct 18 2023 at 03:20):

Where does π\pi come from in this computation?

Yury G. Kudryashov (Oct 18 2023 at 03:20):

From some value of the Gamma function?

Yury G. Kudryashov (Oct 18 2023 at 03:25):

I guess, it comes from Gamma(1/2)

Yury G. Kudryashov (Oct 18 2023 at 04:17):

I meant, in the case p=2.

Xavier Roblot (Oct 18 2023 at 07:10):

Yes, the case p = 2 is done here: https://en.m.wikipedia.org/wiki/Volume_of_an_n-ball#Gaussian_integrals

Xavier Roblot (Oct 18 2023 at 07:11):

I guess I'll formalize these computations since I am the one that needs them.

Xavier Roblot (Oct 30 2023 at 08:35):

#8030


Last updated: Dec 20 2023 at 11:08 UTC