从频率空间来看光滑性
我们现在学习一个从频率的观点看函数的重要的例子:
给定 T 上定义可积函数 f. 如下两个叙述是等价的:
1) | f 是光滑函数, 即 f∈C∞(T); |
2) | f(k) 的具有任意可能的多项式衰减, 即对任意的 N⩾0, 总存在常数 CN, 使得对任意的 k∈Z, 我们有∣∣f(k)∣∣⩽(1+∣k∣)NCN. |
函数 Fourier 系数的衰减速度决定了函数的光滑性, 反之亦然. 换句话说, 频率空间的衰减等价于物理空间的光滑性. 这是学习 Fourier 分析必须要记住的原则之一.
证明. 我们先证明 1)⇒2). 不妨假设 k=0, 证明的关键就是利用光滑性, 我们可以连续做多次分部积分并且分部积分得到的在端点 0 和 2π 处的值会相互抵消 (利用周期性) : f(k)=2π1∫02πf(x)e−ikxdx=2π1(−ik)1∫02πf(x)(e−ikx)′dx=2π1ik1∫02πf(x)′e−ikxdx ⋯ ⋯=2π1(ik)N1∫02πf(N)(x)e−ikxdx.所以, ∣f(k)∣⩽2π1∣k∣N1∫02π∣f(N)(x)∣dx⩽∣k∣N∥f(N)∥L∞.
证明 2)⇒1) 的方法是使用 Lebesgue 控制收敛定理或者上学期证明的关于求导和求和可交换的命题. 首先, 根据条件 2) (选 N=2 即可) , 我们知道 k∈Z∑∣f(k)∣ 是绝对收敛的, 所以, 根据刚刚证明的结论: f(x)=k∈Z∑f(k)eikx是连续函数. 另外, 利用 N=3, 我们知道函数项级数k∈Z∑(f(k)eikx)′=k∈Z∑ikf(k)eikx是绝对收敛的, 因为k∈Z∑∣ikf(k)∣⩽k∈Z∑(1+∣k∣)3C3∣k∣<∞.
根据 Lebesgue 控制收敛的第二个推论, 我们有
f(x)′=(k∈Z∑f(k)eikx)′=绝对收敛,所以是连续函数k∈Z∑ikf(k)eikx.所以,
f∈C1(T). 我们可以继续这个过程, 从而证明对任意的
ℓ∈Z⩾1,
f∈Cℓ(T).
从证明中可见, 对函数光滑 f 求导数等价于对每个频率 k 的 Fourier 系数乘以 −ik. 把这个现象和线性代数联系起来是很有启发性的: 我们考虑线性空间 L2(T), 那么, dxd:L2(T)→L2(T)是线性算子 (当然, 你可能会关心 dxd 是否是良好定义的, 让我们假设这个成立 (其实不成立) ) . 我们现在考虑测度空间 (Z,P(Z),μ), 其中, μ 就是数元素个数的测度, 我们令 ℓ2(Z)=L2(Z,P(Z),μ). 那么, 根据关于 L2 空间上的 Fourier 级数的基本定理, 我们得到了 (连续的) 映射 (同构) : F:L2(T)→ℓ2(Z), f↦(f(k))k∈Z.我们把这个同构看做是在 L2(T) 选取了一个基. 那么, 上面求导可以用下面的交换图标来表达:
换句话说, 换了这组基后, 微分算子 dxd 被对角化了 (变成了成在每个分量上乘一个数) , 这把求微分的操作变成了代数操作.
仿照线性代数中对一个正定的矩阵可以定义它的开方的思路, 我们对一个光滑函数的 Fourier 系数分别乘以 ∣k∣21, 这可以定义求二分之一次导数. 这实际上是所谓的拟微分算子的基本想法 (如果下学期时间足够的话, 我们会讲到) .
我们上学期已经证明如下的 Riemann-Lebesgue 引理, 这个引理是 Fourier 分析中的核心定理:
对任意的 f∈L1(T), 我们都有∣k∣→∞limf(k)=0.
证明. 上面的定理表明, 如果
f∈C∞(T), 那么命题成立, 再利用光滑函数在
L1(T) 中稠密即可. 我们不再重复证明的细节.
实际上, 我们可以证明, 对任意的序列 {ak}k∈Z, 其中 ak=0, 其中k→∞limak=∞,总存在某个可积函数 f∈L1(T), 使得k→∞limak∣f(k)∣=∞.换句话说, 一个 L1 函数的 Fourier 系数的衰减可以任意得慢. 这个证明需要用到泛函分析中的 Banach-Steinhaus 定理, 我们这里仅仅陈述结论而不给证明. (在本次的作业中, 我们有一个较弱的版本, 证明上面的极限对一个子序列是成立的)
证明中分部积分的技巧实际上还证明了如下的结论: 如果 f∈Cm(T), 其中 m⩾1 是正整数, 那么存在常数 C, 使得对任意的 k∈Z, 我们有∣∣f(k)∣∣⩽(1+∣k∣)mCN.也就是说, 每多一阶可微性, Fourier 系数的衰减也快一阶.
Dirichlet 核和 Féjer 核
我们上周就提过: 可以对 T 上的两个 L1 函数定义卷积: f∗g:T→C, x↦(f∗g)(x)=∫02πf(x−y)g(y)2πdy,与 Rn 上的证明完全一致, 利用 Fubini 定理, 我们有∥f∗g∥L1⩽∥f∥L1∥g∥L1.这表明, 我们有如下的乘法结构L1(T)×L1(T)⟶∗L1(T).仍然根据 Fubini 定理, 我们可以说明卷积满足交换律和结合律:
1) | 对于几乎处处的 x∈T, 我们有 (f∗g)(x)=(g∗f)(x). |
2) | 假设 f,g,h∈L1(Rn), 那么 ((f∗g)∗h)(x)=(f∗(g∗h))(x) 几乎处处成立. |
类似地, 我们可以证明: 对于 p=1,2 或者 ∞, 有L1(T)×Lp(T)⟶∗Lp(T),即对任意的 f∈L1(T), g∈Lp(T), 函数(f∗g)(x)=∫Tf(x−y)g(y)dy在 Lp(Rn) 中是良好定义的 (几乎处处) . 特别地, 我们还有如下的不等式: ∥f∗g∥Lp(T)⩽∥f∥L1(T)∥g∥Lp(T).
卷积的定义不仅仅是概念上的推广, 它有着很重要的实际意义. 我们来看下面的例子:
假设 f∈L1(T), 考虑 Fourier 分析中最基本的函数ek(x)=eikx.我们来计算 f∗ek: (f∗ek)(x)=∫02πf(y)eik(x−y)2πdy=eikx∫02πf(y)e−iky2πdy=f(k)eikx.所以, 与 eikx 做卷积就给出了 f 向 eikx 这个方向上的投影.
特别地, 我们在 Fourier 技术中感兴趣的部分和就有如下的表达式:
SN(f)(x)=−N⩽k⩽N∑f(k)eikx=f∗(−N⩽k⩽N∑ek).我们定义 Dirichlet 核 DN(x) 如下, 其中, N⩾0 是整数: DN(x)=−N⩽k⩽N∑ek=eix−1ei(N+1)x−e−iNx=ei2x−e−i2xei(N+21)x−e−i(N+21)x=sin(21x)sin((N+21)x).所以, 对于 x∈R, 我们定义 Dirichlet 核 DN(x): DN(x)=⎩⎨⎧sin(21x)sin((N+21)x), 2N+1, x=2kπ;x=2kπ.这也是定义在 T 上的函数.
我们就证明了下面的引理:
对任意的 f∈L1(T), 它的 Fourier 级数的部分和 SN(f) 可以表达为SN(f)(x)=−N⩽k⩽N∑f(k)eikx=(f∗DN)(x).
在传统的 Fourier 分析中, 我们还考虑一种被称作 Fejér 和的部分和. 按照定义, 它是部分和序列 {SN}N⩾0 的 Cesàro 求和 (前 n 项的平均值) , 即σN(f)(x)=NS0(x)+S1(x)+⋯+SN−1(x)我们将会看到, Cesàro 求和会让 Fourier 级数的收敛速度变得很快!
我们现在来计算 Féjer 核的公式. 按照定义, 我们有σN(f)(x)=(f∗N1k=0∑N−1Dk)(x)=(f∗FN)(x).其中, Féjer 核 FN 被定义成FN(x)=N1k=0∑N−1Dk(x).我们可以进一步计算 FN(x) 的解析表达式: FN(x)=N1k=0∑N−1sin(21x)sin((k+21)x)=N1sin(21x)1k=0∑N−1Im(e2ixeikx)=N1sin(21x)1Im(e2ixeix−1eiNx−1)=N1sin(2x)1Im(ei2Nxei2x−e−i2xei2Nx−e−i2Nx)=N1(sin(2x)sin(2Nx))2所以, 对于 x∈R, 我们定义如下的 Féjer 核: FN(x)=⎩⎨⎧N1(sin(2x)sin(2Nx))2, N, x=2kπ;x=2kπ.
证明, FN(x) 的积分是 2π: ∫02πFN(x)2πdx=1.
这是一个有趣的联系, 最好的证明是观察到
∫−ππDN(x)2πdx=1, 这是因为 (不要用解析表达式)
DN=−N⩽k⩽N∑ek.
我们之前卷积的核函数是 χε, 与之类似, 我们的 Féjer 核 FN 满足下面的性质:
1) | FN 是正的, 即对任意的 x∈R, 我们有 FN(x)⩾0; |
2) | FN 是 T 上的概率密度, 即∫−ππFN(x)2πdx=1; |
3) | 把 FN 视作是 [−π,π] 上的函数, 那么, 当 N→∞ 时, FN 集中在 x=0 附近, 即对任意的 δ>0, 当 N→∞ 时, 我们有∫δ⩽∣x∣⩽πFN(x)2πdx⟶0, N→∞. 证明是直截了当的: 当 ∣x∣∈[δ,π] 时, 由于 ∣∣sin(2x)∣∣⩽2∣x∣, 我们有∣FN(x)∣⩽N1(δ2)2(sin(2Nx))2⩽Nδ21.所以, 当 N→0 时, 我们有∫δ⩽∣x∣⩽πFN(x)2πdx⩽∫δ⩽∣x∣⩽πNδ242πdx→0. |
传统的 Fourier 分析把 FN 的性质提炼成所谓的好的积分核的性质:
给定 R 上的一族连续函数 {KN(x)}N⩾0, 我们可以认为 x∈T. 假设它们满足下面的性质:
1) | KN 在 L1(T) 中有界: 存在常数 M>0, 使得对任意的 N⩾0, 都有∫−ππ∣KN∣2πdx⩽M; |
2) | KN 是 T 上的概率密度, 即 ∫−ππKN(x)2πdx=1; |
3) | 当 N→∞ 时, {KN}⩾0 集中在 x=0 附近, 即对任意的 δ>0, 当 N→∞ 时, 我们有∫δ⩽∣x∣⩽π∣KN(x)∣2πdx→0. |
我们就称 {KN(x)} 是 T 上一族好的积分核.
按照刚才的构造, Féjer 核 {FN}N⩾0 是一族好的积分核.
假设 {KN(x)}N⩾1 是一族好的积分核, 那么对任意的 f∈C0(T), 我们都有一致收敛: N→∞lim∥f∗KN−f∥L∞=0,其中f∗KN(x)=∫−ππf(x−y)KN(y)2πdy.
证明. 利用
KN 是概率密度函数, 我们有
∣(f∗KN)(x)−f(x)∣=∣∣∫−ππ(f(x−y)−f(x))KN(y)2πdy∣∣⩽A∣∣∫∣y∣⩽δ(f(x−y)−f(x))KN(y)2πdy∣∣ +B∣∣∫δ⩽∣y∣⩽π(f(x−y)−f(x))KN(y)2πdy∣∣.由于
f 在
[0,2π] 上面连续, 所以一致连续性. 从而, 对任意的
ε>0, 存在
δ>0, 当
∣y∣<δ 时, 对任意的
x, 我们都有
∣f(x−y)−f(x)∣<ε. 所以, 根据好的积分核的性质 1), 我们有
A⩽ε∫∣y∣⩽δ∣KN(y)∣2πdy⩽εM.对于
B, 我们用好的积分核的性质 3):
B⩽2∥f∥L∞∫δ⩽∣y∣⩽π∣KN(y)∣2πdy.所以, 存在
N0 依赖于
δ 从而依赖于
ε, 使得当
N⩾N0 时, 对任意的
x∈T, 我们有
B⩽ε.综合上述, 对任意的
ε>0, 存在
N0, 使得当
N⩾N0 时, 对任意的
x∈T, 我们有
∣(f∗KN)(x)−f(x)∣⩽(M+1)ε.命题得到证明.
对任意的 f∈C0(T), 我们都有一致收敛N→0lim∥f∗FN−f∥L∞=0.
特别地, 我们注意到
FN(x)=N1∣k∣⩽N∑(N−∣k∣)eikx,所以,
σN(f)(x)=FN∗f(x)=∣k∣≤N∑(1−N∣k∣)f(k)eikx.这是个有限的三角级数, 它将一致收敛到
f, 所以, 我们实际上再次证明了周期连续函数情形下的 Stone-Weierstrass 定理.
Féjer 核的一个应用: L1 函数的原函数与 Fatou 的反例
根据 Riemann-Lebesgue 引理, f∈L1(T), 那么 f(k)=0. 那么, 是不是每个极限是 0 序列 {ak}k∈Z, 它必为某个可积函数的 Fourier 系数呢, 即是否存在 f∈L1(T), 使得对任意的 k, 我们都有 f=ak?
给定 f∈L1(T) 是周期函数, 我们假设 f(0)=0. (f(0)=0 等价于 f 的积分为零) 我们定义 R 上的函数: F:R→C, F(x)=∫0xf(t)dt.根据 f 是可积的, F(x) 是连续函数 (比如, 可以用 Lebesgue 控制收敛来证明) . 由于 f 在一个周期上的积分消失, 所以 F(x) 也是以 2π 为周期的. 我们可以计算 F(x) 的 Fourier 系数:
给定 f∈L1(T), 我们假设 f(0)=0. 那么, 以 2π 为周期的连续函数 F(x)=∫0xf(t)dt 的 Fourier 系数可以用 f 的信息表示如下: F(k)={−∫02πxf(x)2πdx, ik1f(k), k=0;k=0.
如果 f 是连续函数, 那么 F 是 C1 的, 所以, F′=f. 此时, 由于求导在频率空间来看是乘以 −ik, 所以当 k=0 时, 命题时显然的 (可以直接分部积分) .
证明. 当
k=0 时, 利用 Fubini 定理, 我们有
F(0)=∫02π(∫0xf(t)dt)2πdx=∫[0,2π]2f(t)1t⩽x(t,x)2πdxdt=2π1∫02π(∫t2πf(t)dx)dt=2π1∫02π(2π−t)f(t)dt.由于
f=0, 所以
F(0)=−∫02πxf(x)2πdx.当
k=0 时, 我们有
F(k)=∫02π(∫0xf(t)dt)e−ikx2πdx=2π1∫02π(∫t2πe−ikxf(t)dx)dt=2π1∫02π−ik1(1−e−ikt)f(t)dt=ik1f(k).命题成立.
在同样的假设下, 由于 F 是连续的, 根据 Féjer 核的理论, 当 N→∞ 时, FN∗F 一致收敛到 F. 另外, 由于FN∗F(x)=∣k∣⩽N∑(1−N∣k∣)F(k)eikx,所以, 在 x=0 处, 我们有FN∗F(0)=∣k∣⩽N∑(1−N∣k∣)F(k)=−2π1∫02πxf(x)dx+1⩽∣k∣⩽N∑(1−N∣k∣)ik1f(k)=−2π1∫02πxf(x)dx−i1⩽∣k∣⩽N∑kf(k)+ANNi1⩽k⩽N∑f(k)−BNNi−N⩽k⩽−1∑f(k).根据 Riemann-Lebesgue 引理, 我们知道 k→∞limf(k)=0, 所以, N→∞limAN=0 (这是第一学期关于数列极限的标准习题: 如果一个数列的极限是 0, 那么它所对应的 Cesàro 和 (前 n 项的平均) 的极限也是 0) ; 类似地, N→∞limBN=0. 由于 F(0)=0, 所以, 当 N→∞ 时, 上面的等式给出了N→∞lim⎝⎛1⩽∣k∣⩽N∑kf(k)⎠⎞=2πi∫02πxf(x)dx.特别地, 这表明对于一个 L1 函数的 Fourier 系数, 数列 {kf(k)}k∈Z 具有一定的 “可求和性”. Fatou 根据这个性质构造了如下的反例:
对于如下定义的数列 {ak}k⩾1, 其中ak=⎩⎨⎧0, 2ilogk1, 2ilog∣k∣−1, ∣k∣⩽1;k⩾2;k⩽−2.不存在 f∈L1(T), 使得对任意的 k∈Z, f(k)=ak.
证明. 我们用反证法: 如果不然, 那么,
f 的 Fourier 级数的部分和为
SN(f)(x)=∣k∣⩽N∑akeikx(=k⩾2∑logksin(kx)).上面的计算表明,
⎩⎨⎧1⩽∣k∣⩽N∑kak⎭⎬⎫N⩾1 是有极限的, 然而, 根据
ak 的定义, 我们有
1⩽∣k∣⩽N∑kak=i12⩽k⩽N∑klogk1.我们可以用上学期最后学习的面积方法来用
(log(x))−1 的积分估计上面求和的大小, 这表明
2⩽k⩽N∑klogk1∼loglogN→∞.然而, 上面的积分应该收敛到
2πi∫02πxf(x)dx, 矛盾.
尽管 Féjer 核的理论更简洁, 我们却更想知道是否 SN(f) 能够足够好地逼近 f, 因为这是最自然的部分和. 这使得我们要研究 Dirichlet 核函数 DN 的性质. 实际上, 尽管 DN(x) 的积分是 1, 但是, 当 N 很大的时候, 它的振荡很厉害, 所有有很多正的部分对积分的贡献可以被负的部分对积分的贡献消掉, 而 DN(x) 并不满足好的积分核的定义中的第一条. 实际上, 我们有
当 N→∞ 时, 我们有∫−ππ∣DN(x)∣2πdx=π24log(N)+O(1).
证明. 我们重新写下 DN(x) 的表达式DN(x)=⎩⎨⎧sin(21x)sin((N+21)x), 2N+1, x=2kπ;x=2kπ.证明的基本想法是把分母中的 sin 函数替换成更容易控制的函数. 我们首先说明: 对任意的 x∈[0,∞), 我们有x−61x3⩽sin(x)⩽x.后面一个不等式是熟知的; 为了说明前一个, 我们知道1−21x2⩽cos(x), x⩾0.对 x 积分即可.
所以, 对 π⩾x⩾0, 我们有x2⩽sin(2x)1⩽x21−241x21⩽x2(1+2×241x2)⩽x2+6x.即x2⩽sin(2x)1⩽x2+6x.
从而, 通过把
[−π,π] 上的积分变成
[0,π] 上的积分, 我们有
0⩽2∫0π∣DN(x)∣2πdx−2∫0π2x∣sin((N+21)x)∣∣2πdxIN⩽6π1∫0πx∣sin((N+21)x)∣∣dx⩽6π1∫0πxdx=12π=O(1).所以, 我们只需要对
IN=π2∫0πx∣sin((N+21)x)∣∣dx的增长进行估计即可. 此时, 我们有
IN=π2∫0(N+21)πx∣sinx∣dx=π2∫0Nπx∣sinx∣dx+O(N−1)π2∫Nπ(N+21)πx∣sinx∣dx=π2k=0∑N−1∫kπ(k+1)πx∣sinx∣dx+O(1)=π2k=1∑N−1∫0πx+kπ∣sinx∣dx+O(1).我们进一步放缩分母, 我们有
(k+1)π1⩽x+kπ1⩽kπ1.从而
π2k=1∑N−1∫0π(k+1)π∣sinx∣dx⩽IN−O(1)⩽π2k=1∑N−1∫0πkπ∣sinx∣dx.然而, 左右两边我们可以直接计算, 从而
π24k=2∑Nk1dx⩽IN−O(1)⩽π24k=1∑N−1k1.最终利用
1∑N−1k1=logN+γ+O(N1),其中
γ 是 Euler 常数. 命题得证.
Dirichlet 核函数的研究将是我们之后学习的重点. 我们先证明 Fourier 级数的部分和或者说 Dirichlet 核的局部化的性质. 这是一个非常有启发意义的命题, 它的证明是 Riemann-Lebesgue 的一个漂亮的应用:
对任意的 x0∈R 任意的小正数 δ>0, 如果函数 f,g∈L1(T) 满足f∣∣(x0−δ,x0+δ)=g∣∣(x0−δ,x0+δ),那么, 我们有N→∞lim(SN(f)(x0)−SN(g)(x0))=0.
据此, 为了研究 f 在某点 x0 处的 Fourier 级数是否收敛, 只要将注意力集中在这个点的邻域就好! 这是绝对是不平凡的结论: 按照定义, Fourier 系数是依赖于 f 在整个 T 的积分的, 然而 Fourier 级数的收敛性却是局部的!
证明. 首先, 由于
DN(x) 是偶函数, 通过变量替换, 我们有
SN(f)(x0)−SN(g)(x0)=∫−ππDN(y)(f(x0−y)−g(x0−y))2πdy=∫0πDN(y)(f(x0−y)+f(x0+y)−g(x0−y)−g(x0+y))2πdy=∫δπDN(y)(f(x0−y)+f(x0+y)−g(x0−y)−g(x0+y))2πdy.根据
DN(x) 的解析表达式, 我们有
SN(f)(x0)−SN(g)(x0)=∫δπsin((N+21)y)F(y)sin(2x)(f(x0−y)+f(x0+y)−g(x0−y)−g(x0+y))2πdy=∫−ππsin((N+21)y)⋅F(y)1δ⩽∣y∣⩽π(y)dy.由于在
x⩾δ 时,
sin(2x)⩾sin(2δ) 有正的下界, 所以
F(x)1δ⩽∣x∣⩽π(x) 实际上在
L1(T) 中. 根据 Riemann-Lebesgue 引理, 当
N→∞ 时, 我们有
∣SN(f)(x0)−SN(g)(x0)∣→0.命题得证.
对上面的证明稍加修改, 我们就可以证明: 对任意的 h∈L1(T) 和任意的小正数 δ>0, 当 N→∞ 时, 我们有∫δ⩽∣x∣⩽πDN(x)h(x)2πdx∫δ⩽∣x∣⩽πFN(x)h(x)2πdx→0,=O(N1).我们观察到, 尽管上面的的极限都消失, 但是后者有固定的衰减速度而前者并没有 (由 Riemann-Lebesgue 引理给出的衰减是没有衰减速率的) .