作者znmkhxrw (QQ)
看板Math
标题[分析] 无限型∫f(x)dx写成无穷级数和之极限
时间Tue Oct 12 00:29:54 2021
想请教一下在Lebesgue积分:
∞ ∞
∫ f(x)dx = lim (1/p)*Σ f(n/p) --(*)
-∞ p→∞ n=-∞
这个式子要成立的话有什麽已知的
充分条件呢?
我自己是有证明一个版本:
-----------------------------------------------------------------
<Theorem>
Let f€L^1(R)
and |f|<=M on R ---(C1)
and f monotonic on [A,∞) and (-∞, -A] for some A>0 ---(C2)
and f is continuous a.e. on R ---(C3)
Then (1) Σ_{n=-∞~∞} f(n/p) exists and finte for all p > 0
(2) (*) holds
-----------------------------------------------------------------
证明的idea是:
(a) 造出阶梯函数f_p(x) := Σ_{n=-∞~∞} f(n/p)*χ_[n/p, (n+1)/p)
造出控制函数g(x) := M , x€[-A,A]
|f(x)| , else
然後就有 |f_p(x)| <= g(x) for all p
P.S.
以上造法仅限於 f在[A,∞)向右递增 & f在(-∞, -A]向左递减
但是根据四种不同的单调排列组合都可以制造出相应的f_p
(b) 因为f在R上a.e.连续, 所以我们有f_p(x)趋近於f(x) a.e. x€R
因为f€L^1(R)所以g€L^1(R)
根据以上两个, LDCT告诉我们下面第一个等号, 第二个等号则by trivial check
∞ ∞ ∞
∫ f(x)dx = lim ∫ f_p(x)dx = lim (1/p)*Σ f(n/p)
-∞ p→∞ -∞ p→∞ n=-∞
==============================================================
总之, (C1)跟(C3)我觉得还OK, 问题就在於
(C2)太严苛了
可是没有单调性我又造不出控制函数QQ
由於我是在推导傅立叶转换与离散时间傅立叶转换的关系时需要(*)式
讯号通常不会有单调性
因此才想问说有没有什麽其他宽松的条件
谢谢!!
(10/13更新)
找到一个符合L^1&(C1)&(C3)的反例:
f(x) = 1 ,for x€Q
0 ,else
则p只要是正整数的话, f(n/p)都是1, 所以级数发散for all p€N
所以(C3)要加强至" 不连续点的数量有限 "才有机会...
(10/13更新2)
全连续, 甚至均匀连续都也有反例, 看来要控制震荡度
(10/14更新)
w大的locally BV也不能...把等腰三角形修改成C^1的话就locally BV了
https://www.desmos.com/calculator/jxkztslgm2
这里我是用sine接起来的, 面积跟等腰三角形一样
所以目前的超强反例是:
f€L^1, C^1, 非负, 趋近於0, 均匀连续, locally BV
而要C^∞我猜也可以, 用test function来接应该就好, 只是面积要再check一下
(10/15) 条件可放宽至以下
<Theorem> (原本)
Let f€L^1(R)
and |f|<=M on R ---(C1)
and f monotonic on [A,∞) and (-∞, -A] for some A>0 ---(C2)
and f is continuous a.e. on R ---(C3)
Then (*)成立
<Theorem> (现在)
Let f€L^1(R)
and f€BV(R)
and f is continuous everywhere
and Σ_{n=-∞~∞}f(n/p) converges
Then (*)成立
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 59.102.225.191 (台湾)
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Math/M.1633969796.A.158.html
1F:推 PPguest : 把(C2)换成non-negative 我猜可能有希望10/12 22:21
谢谢帮想, 不过刚刚写了一下发现(C2)换成non-negative(称作C2')是一样的耶
也就是说
如果(C2)换成(C2')则定理成立, 则把(C2)拿掉也会定理成立
证明如下:
若任何L^1函数满足(C1), (C2'), (C3)则(*)成立
则对於任何f€L^1且满足(C1),(C3), f就能拆成f^+跟f^-且都是L^1且都满足(C2')
如此一来藉由假设能得知f^+跟f^-都满足(*)
接着透过减法就能知道证出f满足(*) 得证
2F:→ willydp : 我猜L^1 + locally BV可以10/13 03:08
确实是震荡的问题, L^1不能控制离散取样的震荡性QQ
3F:→ PPguest : 我觉得(*)等号右边可能要搞定黎曼积分,然後还要解决10/13 10:05
4F:→ PPguest : 从有限区间扩展到实数可能会有的问题10/13 10:06
5F:→ PPguest : (C3)我们要的应该是不连续点的集合0测度,这样反例应10/13 10:07
6F:→ PPguest : 该可以避开10/13 10:08
应该要均匀连续或是有界变量来控制震荡, 全连续也有反例, 刚刚想到一个:
令f在整数点取值为1, 在每个n都造一个底为1/n^2的等腰三角形, 其余都是0, 则f是L^1
但是所有正整数p都会让级数发散QQ
7F:→ PPguest : 不太懂你造的全连续反例,n是什麽?跟(*)里的n有关吗?10/13 18:30
8F:→ PPguest : 等腰三角形的高是?10/13 18:32
9F:→ PPguest : 如果你是说整数点都是等腰三角形的顶点,n是指整数点 10/13 18:36
10F:→ PPguest : (*)等号右边改成黎曼积分的瑕积分会不会就没问题了?10/13 18:39
https://www.desmos.com/calculator/d3gsgfck5u
大致上像这样
11F:→ PPguest : 了解.级数不能爆掉的话f要限制,两边没有趋近於0的应10/13 19:54
12F:→ PPguest : 该是不行10/13 19:54
趋近於0也是有反例, 只要把三角形的高变成1/n, 底变成1/n, 这样f还是L^1
而级数和会是Σ1/n, 发散, 所以还是震荡问题@@
https://www.desmos.com/calculator/3c5gjvfcjj
13F:→ PPguest : 看到z大提到震荡,我想z大应该知道第一个反例,也就是10/13 20:03
14F:→ PPguest : 有理数1,其他0的函数,是处处不连续,不符合(C3) 10/13 20:05
欸对 我在干嘛XDDD 不过现在有处处连续的反例 只能朝震荡发展了QQ
15F:→ PPguest : f两边趋近於0大概是必要条件,直接要求级数收敛会不10/13 20:38
16F:→ PPguest : 会比较快?10/13 20:38
17F:→ PPguest : ^f要满足10/13 20:42
当初有这样想, 但是这条件很难检查(太刻意&太强了), 所以还没往这方向
而且即便收敛, 要说(*)左右相等我也没其他想法
因为最初就是靠控制函数的让 (1) 级数收敛 (2) (*)左右相等
今天如果假设(1)是对的, 还是没其他idea说明(2)
总之我现在是往 "控制震荡後就存在L^1的控制函数"下手
18F:→ PPguest : 目前的想法是考虑f非负,(*)左式用MCT变成有限区间例10/14 10:37
19F:→ PPguest : 如[-k,k]的积分然後对k取极限,接着再把积分变成黎曼10/14 10:37
20F:→ PPguest : 积分,再来想把黎曼积分写成黎曼和的极限,这时会碰到 10/14 10:38
21F:→ PPguest : 两个极限是否能交换,如果可以,大概就变成(*)右式10/14 10:38
22F:→ PPguest : 另外还要处理细部例如黎曼和怎麽取以及边边的部分10/14 10:40
L积分比R积分的好处是容易处理极限交换的问题, 因为R积分极限交换要要求x变数对y变
数是均匀收敛, y变数对x变数逐点收敛, 这更难check...
不过既然locally BV都控制不了震荡性的话, 目前看起来只有单调而已了, 或是你说的直
接让级数收敛当条件看有没有办法逼近积分了
23F:推 Vulpix : 直接BV呢?不要局部看。 10/14 16:08
∞
V大根据你这个idea并且加入f€C^1, 接着利用 Var[a,+∞] = ∫ |f'| dx < ∞
a
目前没什麽进展QQ
就是找不到idea去控制step function的每段的面积...
(10/15 更新)
V大你这个条件OK耶!! 只是要加入P大的条件"级数收敛", 整体叙述我放在文末
24F:→ PPguest : 我想到一个应该是非BV,但级数应该是收敛的例子 10/14 20:49
25F:→ PPguest : 当然这不影响BV很可能有希望 10/14 20:49
26F:→ PPguest : 用1/n^2的整数点,整数点之间造n个震荡 10/14 20:49
27F:→ PPguest : 这样variation应该会像1/n那样跑到无限大 10/14 20:49
28F:→ PPguest : 级数收敛用1/x^2或者p^2/x^2压住(comparison test)10/14 20:50
P大你这个是(1) f€L^1 & 级数收敛for all p
(2) (*)相等
(3) 非BV
(4) f连续且f趋近於0
吗?
(10/15更新) 嗨P大, 级数收敛&V大的全域BV确实就得证了
接着就差好的条件去让级数收敛了
30F:→ PPguest : 之前我想证 若我们有一个跟f有关的级数收敛,则有均10/15 14:12
31F:→ PPguest : 匀收敛,两个极限可以交换(先考虑f非负).初步想想感10/15 14:13
32F:→ PPguest : 觉有希望.後来想说改写成积分形式的,发现基本上就是10/15 14:15
33F:→ PPguest : 在做控制函数,假设是L^1(R),直接用你的证法就解决了 10/15 14:17
34F:→ PPguest : 我考虑的级数换成函数是长这样:10/15 14:19
35F:→ PPguest : Σ_{n=-∞~∞}sup_{[n,n+1]}|f(x)|*χ_[n,n+1) 10/15 14:20
36F:→ PPguest : 如果f不是非负的话,可能要检查一下 10/15 14:20
37F:推 PPguest : 关於我提的例子,如果我没想错的话,(1)是对的 10/15 14:38
38F:→ PPguest : (2)是对的,因为可以找到L^1(R)的控制函数 10/15 14:39
39F:→ PPguest : (3)应该对,虽然以前BV好像都是用在有限区间的函数 10/15 14:41
40F:→ PPguest : (4)是的.0那边用直线连起来,往无限大的部分像1/x^2 10/15 14:42
41F:→ PPguest : (2)我好像讲错,控制函数要和p无关,暂不确定10/15 15:12
42F:→ PPguest : (好多东东差点搞混) ↑我是指用p^2/x^2压住10/15 15:15
43F:→ PPguest : (2)藉由1/x^2平移证明那个控制函数是L^1,所以是对的10/15 15:21
44F:推 PPguest : 直观上我以为加上BV就能解决,毕竟反例应该都挂了10/15 15:28
45F:→ PPguest : f是L^1(R),BV,若f能写成两递增函数相减,那两个递增10/15 15:30
46F:→ PPguest : 函数会L^1(R)吗?10/15 15:31
47F:→ PPguest : 啊,不要理我,应该不行,通常都会爆掉或出问题10/15 15:38
48F:→ PPguest : 直观上觉得有BV的话,就不会有无限多的x让f(x)的值跑10/15 15:47
49F:→ PPguest : 到会让级数发散的地方10/15 15:47
50F:→ PPguest : 我刚想到一个证明,要证若f是L^1(R),BV,连续,则那个10/15 17:55
51F:→ PPguest : 控制函数也会是L^1(R)10/15 17:56
52F:→ PPguest : 先加个f非负好了10/15 17:56
53F:→ PPguest : 假设要证明的东西是错的,那个函数的L积分是∞10/15 17:57
54F:→ PPguest : 目标是要得到variation也会是∞10/15 17:58
55F:→ PPguest : 考虑另一个阶梯函数,也就是把那个控制函数的sup改成 10/15 17:59
56F:→ PPguest : inf,这样子这个函数会是L^1(R)10/15 18:02
57F:→ PPguest : 直观上每个区间sup和inf的差累积起来应该就要爆掉 10/15 18:07
58F:→ PPguest : 细节不确定sup和inf是两端点都有取,会不会造成问题 10/15 18:09
============================================================
帮PPguest贴文把原本的长推文换成写得比较完整一点的文字
1.想到一个证明
本来想证在f€L^1(R),(C1)(C3),非负的条件下,若我们有一个跟f有关的级数收敛,
[kp]
则当k→∞时,(1/p) * Σ f(n/p) 对p均匀收敛。
n=-[kp]
其中k的来源是把整个R上的积分变成[-k,k]区间的积分再取极限,
∞ MCT k
∫ f(x)dx = lim ∫ f(x)dx
-∞ k→∞ -k
p则是来自於把黎曼积分换成黎曼和的极限。
目标是p和k的极限可以交换。
初步想觉得有希望,後来不小心想说把新增的条件改写积分的型式,
结果发现本质上似乎就是做出LDCT的控制函数。
(本来做的控制函数後来发现无法大於等於|f_p(x)|)
-----------------------------------------------------------------
<Theorem>
Let f€L^1(R)
and |f|<=M on R ---(C1)
and f is continuous a.e. on R ---(C3)
and Σ_{n=-∞~∞}sup{|f(x)|, x€[n,n+1)} converges ---(C4)
Then (*)成立
-----------------------------------------------------------------
证明的想法:
大方向和原po最初的证明一样要用LDCT,这里只讲控制函数的部分。
(ⅰ)令 u(x) = Σ_{n=-∞~∞} u_n *χ_[n,n+1), u_n = sup{|f(x)|, x€[n,n+1)}
由(C4)知 u(x)€L^1(R)。
(ⅱ)(造控制函数)
令 g(x) = Σ_{n=-∞~∞} g_n *χ_[n,n+1),
g_n = ︴u_n, if u_(n-1) <= u_n
︴u_(n-1), if u_(n-1) > u_n
(ⅲ)(g(x)€L^1(R))
g(x) = u(x) + Σ_{n=-∞~∞} (g_n-u_n) *χ_[n,n+1)
<= u(x) + Σ_{n=-∞~∞} u_(n-1) *χ_[n,n+1) < +∞
(ⅳ)(|f_p(x)| <= g(x) for p >= 1)
对每个 [n/p, (n+1)/p) 区间,n€Z, n/p 一定落在某个 [N,N+1) 区间中。
若 (n+1)/p <= N+1,则 所有 [n/p, (n+1)/p) 内的 x 都有 |f_p(x)| <= g(x)。
若 (n+1)/p > N+1,则 因为 p >= 1,所以 (n+1)/p <= N+2;
所有 [n/p, N+1) 内的 x 都有 |f_p(x)| <= g(x);
所有 [N+1, (n+1)/p) 内的 x 都有 |f_p(x)| <= u_N <= g_(N+1) = g(x)。
==========================================================================
2.关於原po提的4个问题,
我忘了说 [0, 0.5] 那边用直线连起来,函数值都是1,方便起见只考虑单边。
(1) f 是 L^1(R),因为函数在 [1,∞) 被 1/x^2 给压住,[0,1] 内有界。
级数收敛for all p 也对,用 p^2/x^2 压住 + 级数的comparison test。
(2) (*)会成立。
在 [1,∞) 本来是说用 1/x^2 平移来压住 u(x),於是 u(x)€L^1(R)。
现在发现直接看级数,用1/x^2 压住 (C4)就成立。再用前面证的定理即可。
(3) 非BV没错。
每个正整数点跟前一个整数点之间有 n 个振荡,振幅是 1/n^2,
Σ_{n=1 to ∞} 2 * n * 1/n^2 <= V[f,[0,∞)]
(4) f连续没错,x 往无限大时 f 趋近於 0 没错。
==========================================================================
3.关於BV
直观上觉得有 f€L^1(R)、(C1)、(C3)、BV,应该就会有(*),毕竟反例应该都挂了。
直观的想法是,若有BV,就不会有无限多的 x 让f(x)的值跑到会让级数发散的地方。
-----------------------------------------------------------------
<Theorem>
Let f€L^1(R)
and |f|<=M on R ---(C1)
and f is continuous a.e. on R ---(C3)
and f is BV function.
Then Σ_{n=-∞~∞}sup{|f(x)|, x€[n,n+1)} converges ---(C4)
-----------------------------------------------------------------
证明的想法:
假设(C4)是错的,Σ_{n=-∞~∞}sup{|f(x)|, x€[n,n+1)} = +∞
目标是total variation V[f,R] 要爆掉。
(ⅰ)
令 u(x) = Σ_{n=-∞~∞} u_n *χ_[n,n+1), u_n = sup{|f(x)|, x€[n,n+1)}
l(x) = Σ_{n=-∞~∞} u_n *χ_[n,n+1), l_n = inf{|f(x)|, x€[n,n+1)}
则对每个 x 都有 l(x) <= |f(x)|,於是 l(x)€L^1(R)
由假设(C4)是错的,∫u(x)dx = +∞
(ⅱ)
因为 | |f(x_(i+1))|-|f(x_i)| | <= | f(x_(i+1))- f(x_i) |,
所以 V[|f|,R] <= V[f,R]。
(ⅲ)
因为 ∫u(x)-l(x)dx = +∞,所以Σ_{n=-∞~∞} u_n - l_n = +∞
但 Σ_{n=-∞~∞} u_n - l_n <= V[|f|,R] <= V[f,R],与 f 是BV 矛盾。
==========================================================================
59F:推 PPguest : 另外,也有点好奇z大新版本的定理大概是怎麽证的10/15 21:25
我证单边的模式, 双边类似, 只是f(0)这一项额外提出来给1/p压掉即可
https://www.desmos.com/calculator/3xyysgekfl
60F:→ PPguest : 那串推文可以用"推"分成三部分,完整打可能要晚一点10/16 09:05
61F:→ PPguest : 另外我刚发现第一部分那个控制函数可能无法大於等於10/16 09:07
62F:→ PPguest : |f_p(x)| for all p10/16 09:08
63F:→ PPguest : 等於整个几乎都要打掉(哭10/16 09:11
没关系啦XDD 有空你再寄站内给我我帮你贴上就好, 推长文不太好看
谢谢讨论罗~
64F:→ PPguest : 目前的进展是满足条件下似乎可以构造出LDCT用的控制10/17 20:10
65F:→ PPguest : 函数,但还是要完整写出来才能确定 10/17 20:11
66F:→ PPguest : 另外,我开始好奇如果在有限区间内,条件只有f是L^1, 10/17 20:12
67F:→ PPguest : finite,不连续点的集合0测度,这样是否足够,或者还要10/17 20:14
68F:→ PPguest : 再加什麽条件才能让无上/下界的f能满足(*)10/17 20:16
昨天为了解决" Σ_{n=-∞~∞}f(n/p) converges "这个条件(还是觉得下这个太刻意)
想到" f(x) = O(1/x^q), q>1 " 这个条件可以达成
进而发现这个条件就可以达成我要的控制函数了...
也就是说, 最後我认为最满意的版本是:
-----------------------------------------------
<Theorem>
Let f s.t. |f|<=M on R
and |f|<=g for some monotone g€L^1 (因为g€L^1所以立得g递减到0)
and f is continuous a.e. on R
Then (1) f€L^1
(2) R.H.S of (*) exists for all p > 0
(3) (*) 成立
------------------------------------------------
反而total bounded variation那个版本的证明我发现:
(a) 上界多估了好多东西
(b) 还是无法保证级数的收敛性
好像控制震荡性又不是那麽足够
而最初只想到要由f本身去做假设然後造出控制函数,
就没想到用" 有个单调的控制函数 "当假设即可
※ 编辑: znmkhxrw (59.102.225.191 台湾), 10/17/2021 21:36:05
69F:→ PPguest : 最後这个版本感觉确实比较好用,我也错过了它10/18 16:33
70F:→ PPguest : 好奇你也是用例如平移1单位,并限制p>=1的方式来压住10/18 16:35
71F:→ PPguest : |f_p(x)|吗?10/18 16:35
72F:→ PPguest : 另外你定理写的条件感觉有点怪10/18 16:36
1. 我是正向f(n/p)都往左边造step function然後用g的单调性控制住
即 f(n/p)*χ_((n-1)/p,n/p]
2. 哪边怪怪的呢??
73F:→ PPguest : 定理里面的第二行看起来是只有一个单调函数g 在整个 10/18 20:08
74F:→ PPguest : R上要L^1又要压住|f| 10/18 20:08
没错耶, 概念来自於要让级数存在的话, 我之前有过假设 f(x) = O(1/x^q), q>1
这样就可以让级数存在, 之後我更发现这个 C/x^q 就可以当 f_p(x) 的控制函数
你觉得怪怪的地方我写清楚後应该如下:
Let f s.t. |f|<=M on R
and |f|<=g for some monotone g€L^1
Then |f_p(x)|<=g(x), where f_p(x) = Σ_{n=1~∞} f(n/p)*χ_((n-1)/p,n/p]
75F:→ PPguest : 要不要像(C2)那样分三段,左右两段各有一个单调函数 10/18 22:14
是需要呀, 只是我都写一边而已(正向)XD, 所以左右都要就是
|f|<=g , g€L^1(R), g is monotonic on [A,∞) and (-∞,-A]
※ 编辑: znmkhxrw (59.102.225.191 台湾), 10/18/2021 22:21:06