温州 网站开发/上海百度首页优化
第 4 章 傅里叶级数的一些应用
Fourier series and analogous expansions intervene very
naturally in the general theory of curves and surfaces.
In effect, this theory, conceived from the point of view
of analysis, deals obviously with the study of arbitrary
functions. I was thus led to use Fourier series in sev-
eral questions of geometry, and I have obtained in this
direction a number of results which will be presented
in this work. One notes that my considerations form
only a beginning of a principal series of researches,
which would without doubt give many new results.
(Fourier级数和类似的展开非常自然地介入了曲线和曲面的一般理论。事实上,从分析的角度来看,这个理论显然是处理任意函数的研究。因此,我被引向在几个几何问题中使用Fourier级数,并且我在这个方向上获得了许多成果,这些成果将在本文中展示。有人指出,我的考虑只是一系列主要研究的开始,这些研究无疑会产生许多新的成果。)
----------------------------------------------------------------------------- A. Hurwitz[húəvits], 1902
在前面的章节中,我们介绍了一些关于Fourier分析的基本事实,这是由物理学中出现的问题引起的。弦的运动和热的扩散是自然地导致Fourier级数函数展开的两个例子。我们接下来建议让读者对Fourier分析的更广泛影响有所了解,并说明这些想法如何延伸到数学的其他领域。 尤其要考虑以下三个问题:
I. 在所有位于实数平面 (译注:手写体双线字母)中的长度为𝓁(译注:字母l的数学体细体)的简单封闭曲线中,哪一个封闭曲线围合的面积最大?
II. 已知一个非比数 γ (译注:希腊字母wa-mma小写),关于数例 nγ (n = 1,2,3,...)的小数部分的分布,我们可以说些什么?
III. 存在处处不可微的连续函数吗?
第一个问题本质上显然是几何问题,乍一看似乎与Fourier级数关系不大。第二个问题介于数论和动力系统研究之间,并为我们提供了遍历性(ergodicity)思想的最简单示例。”第三个问题虽然本质上是分析问题,但在最终发现解决方案之前触动了(resisted)许多尝试。 值得注意的是,这三个问题都可以通过使用Fourier级数非常简单直接地得到解决。
在本章的最后一节,我们将回到一个引起我们最初动机的问题。我们考虑圆周上的时间相关的热传导方程。在这里,我们的调查将引导我们找到针对这个圆的重要而神秘的(enigmatic)热内核。然而,在我们应用Poisson求和公式(我们将在下一章中进行应用)之前,将无法完全理解围绕其基本属性的奥秘(mysteries)。
1. 等周不等式(The isoperimetric inequality)
令 Γ(译注:希腊字母wa-mma大写)表示平面上不与其自身相交叉(intersect)的封闭曲线。加外,令 𝓁 表示 Γ 的长度,并且令 𝒜 (译注:字母 A 的数学体细体,LaTex 宏 \mathcal{A})表示 平面中由 Γ 围合的界内区域的面积。现在的问题是,针对已知的曲线长度 𝓁 确定使用面积 𝒜 最大化的曲线 Γ 的表达式(如果这样的曲线存在)。
-----------------------------------------------------图 1. 等周(等周长)问题----------------------------------------
一些实验(experimentation)和深思(reflection)表明这个解应该是一个圆。可以通过以下启发式(heuristic)考虑得出该结论。这条曲线可以被认为是平放在桌面上的一条闭合的弦(string)。如果弦所围合的区域不是凸面(convex)(比如),则可以使部分弦变形(deform)并增加它所围合的面积。此外,通过一些简单的例子,可以说服自己,曲线在某些部分越“扁平(flatter)”,围合区域的效率就越低。因此我们想最大化曲线在每个点的“圆度(roundness)”。
尽管圆是正确的猜想结果,但是,要在数学上精确地表达这种思想却是一件困难的事。
我们给出的等周问题解的关键思想包括针对Fourier级数应用Parseval恒等式。然而,在我们尝试解决这个问题之前,我们必须首先定义简单闭合曲线的概念、它的长度,以及它所围合的区域的面积。
1.1. 曲线、长度及其面积
一条参数化的 γ 曲线是一个映射 。
在平面上,γ 的图像是一个点的集合,我们称其为一条曲线(curve),并用希腊字母 Γ 表示。如果曲线 Γ 其自身不交叉,则曲线很简单,并且,如果其两个端点叠合(coincide),则曲线是围合的(closed)。根据以上参数化表达式,这两个条件转换成:除非满足 ,且有
,否则
。我们可以将 γ 扩展为实数域 ℝ 上周期为 b – a 的周期函数,并将 γ 视为圆周上的函数。此外,我们总是施加某些平滑条件在我们的曲线上,假设 γ 是
类平滑曲线,其一阶导数 γ’ 满足条件 γ’ ≠ 0。这些条件合起来确保了曲线 Γ 在其上每一点具有定义良好的切线,且切线随着曲线上点的变化而连续变化。此外,参数化的 γ 曲线随着 s 从 a 行进到 b 的过程中,会在 Γ 上产生一个可变的方向。
任何一个 类双射映射(bijective mapping)
都会引起 Γ 的另一个参数化表达式,用公式表示为
。 显然,γ 是围合的且简单的这两个条件与所选取的参数无关。此外,如果对于所有的 t ,s ’(t) > 0 ,则我们称这两个参数化表达式 γ 和 η 是等价的;这意味着 η 和 γ 在曲线上产生的方向是相同的。然而,假如 s ’(t) < 0,则 η 产生相反的方向。
如果 γ 的参数化表达式为 , 则曲线 Γ 的长度定义为
。
曲线 Γ 的长度对曲线而言是一个固有的(intrinsic)概念,不取决于其参数表达式。为了理解这个事实,假设 ,则变量变化公式和链式法则意味着
,
这个结果正如预期。
在下面证明这个定理的过程中,我们将使用针对曲线 Γ 的一个特殊参数化类型。假如对于所有的 s,有 |γ (s)| = 1,则我们称 γ 由弧长(arc-length)参数化。这意味着 γ (s) 以固定的速度行进,最终结果便是,Γ 的长度正好等于 b – a 。因此,经过可行的额外转换之后,就得到了一个定义在闭区间[0, 𝓁]上的按弧长参数化的表达式。任意曲线都允许用弧长进行参数化(练习1)。
现在我们转向等周长问题。
尝试给出由简单闭合曲线 Γ 围合区域的面积 𝒜 的准确公式引发了话多棘手的(tricky)问题。在各种简单的面积计算场景中,由下面熟悉的积分公式给出的面积计算公式是显而易见的:
(1) ;
参见,例如,练习3 。因此,在公式化我们的结论的时候,我们将采用容易的权宜之计(expedient),将(1)作为面积的定义。这种方法允许我们快速而简洁地证明等周不等式。在定理证明之后可以找到此简化未解决的问题列表。
1.2. 等周不等式的表述和证明
定理1.1 假如 Γ 是平面 上的一条长度为 𝓁 的闭合曲线,令 𝒜 表示由这条曲线围合的区域的面积。则
,
且当且仅当 Γ 是圆的时候,等号成立。
通过观察,我们首先发现我们可以重新调整(rescale)问题。这意味着,我们可以用一个 δ > 0 的因子按随后的方式更改度量单位。考虑平面 到其自身的映射——将点(x,y)投射到(δx, δy)。观察定义曲线长度的公式可以发现,长度公式表明,如果 Γ 的长为 𝓁 ,则它在这个映射下的图像的长度为 δ𝓁 。因此,这个运算是按比例因子 δ 放大(magnifies)还是缩小(contracts) 曲线长度,这取决于 δ ≥ 1 还是 δ ≤ 1 。类似地,我们可以看到,这个映射按比例因子
放大(或缩小)围合区域的面积。通过取 δ = 2π/𝓁 ,我们可以看到,这足够证明,假如 𝓁 = 2π ,则 𝒜 ≤ π ,且仅当 Γ 是一个圆的时候等式成立。
令 且 γ (s) = (x(s),y(s)) 为用曲线 Γ 的弧长表示参数的参数化表达式,即, 对于所有 s∈[0,2π],有
。这意味着,
(2) 。
因为曲线是闭合的,函数 x(s) 和 y(s) 是 2π 周期的,因此,我们可以考虑其Fourier级数
和
。
则,正如我们在第2章第2节的后面评注的那样,我们有
和
。
应用 Parseval 恒等式到 (2) 得到
(3) 。
现在应用 Parseval 恒等式的双线性形式(第3章,引理1.5)到定义 𝒜 的积分。因为 x(s) 和 y(s) 是实数值函数,我们有 和
, 因此,我们求得
。
接下来,我们观察到
(4) ,
并且,因为 ,我们可以用(3) 求得
,
正如所预期的那样。
当 𝒜 = π 时,我们从以上的论据中可以看出
和
。
因为 且
。从前面我们知道,x(s) 和 y(s) 都是实数值的,因此,
和
。恒等式 (3) 意味着
,因此,我们有(4)中的等式,从而我们一定有
。 我们写成
和
。
这个事实意味着
(译注:由Euler恒乖式推出
), 因此,α – β = kπ ,其中,k 是一个奇数。据此,我们求得
和
, 其中,y(s) 中的符号取决于(k - 1)/2 的奇偶性(parity)。我们发现,在任何情况下,Γ 是圆的时候,等号这个情况总是成立,因此,我们完成了证明。
以上给出的解(归功于1901年的Hurwitz[húəvits]提出的定理,Adolf Hurwitz,1859年3月26日~1919年11月18日,德国数学家)确定非常优雅,但很明显,留下了一些未解决的(answered)重要问题。我们把这些问题列在下面。假设 Γ 是一条简单的闭合曲线。
( ) “由 Γ 围合的区域”是如何定义的?
( ) 这个区域的“面积”的几何定义是什么?这个定义与(1)一致吗?
( ) 这些结论可以推广为最广泛的与这类问题(即这些曲线是“可求长的(rectifiable)”,即,具有有限长度的曲线)相关的简单闭合曲线的通用类型吗?
事实证明,要阐明提出的这些问题,就必然与分析学中其它的重要问题产生联系。我们将在本系列的后续书籍中回到这些问题。
2. Weyl[vail]的等分布定理(Weyl's equidistribution theorem)
我们现在将来自Fourier级数的思想应用到处理非比数(irrational)属性的问题。我们首先简要讨论同余(congruence),这是理解我们的主要定理所需的概念。
2.1. 实数对整数取模(The reals modulo the integers)
假如 x 是实数,令 [x] 为小于等于 x 的最大整数,则称这个量 [x] 为 x 的整数部分(integer part)。而 x 的小数部分(fractional part)定义为 < x > = x - [x] 。特别地,对于任意 x ∈ ℝ , < x > ∈ [0,1) 。例如,2.7 的整数和分数部分分别是 2 和 0.7 ,而 -3.4 的整数和分数部分分别是 -4 和 0.6 。
假如 ,我们可以定义 ℝ 上的一种关系,比如说,定义 x 和 y 的等价关系、或称同余关系(congruence)。则,我们可以写成
或者
。
这意味着,如果两个实数相关一个整数,我们就可以识别它们。观察到,任意实数 x 与 [0,1)中唯一的数(恰好是 < x >,即 x 的小数部分)是同余的。事实上,产生一个实数对整数 ℤ 取模意味着仅关注其小数部分而忽略掉其整数部分。
现在,以实数 γ ≠ 0 开始,并观察序列,γ,2γ ,3γ ,... 。一个引人入胜的(intriguing)问题是,我们要问,如果序列对 ℤ 取模,序列会发生什么,即,如果我们观察其小数部分的序列
<γ>,<2γ> ,<3γ> ,...
会发现什么现象?
下面是一些简单的观察结果:
( ) 假如 γ 是比率数(rational),则仅有有限多个出现在 <nγ> 中的数是不同的。
( ) 假如 γ 是非比数(irrational),则 <nγ> 中所有的数都是不同的。
事实上,对于第 ( i ) 部分,注意到,假如 γ = p/q,在序列中的前面 q 项是
< p/q >,<2 p/q > ,... <(q-1)p/q >,< q p/q> = 0 。
然后,序列开始自重复,因为
<(p + 1) q > = <1 + p/q >= < p/q >,
等等。尽管如此,更精炼的例子,请参看练习6。
此外,对于第( ii )部分,假设并非所有的数都是不同的(译注:即有的数相同)。因此,对于某些 ,我们有
;则
,因此,γ 是比率数,这与假设矛盾。
事实上,可以证明明,假设 γ 是一个非比数,则 <nγ> 在区间[0,1)上是稠密的(dense),这是一个源自Kronecker[króunèkə] 所证明的结论(译注:Kronecker,德国数学家,逻辑学家,1823年12月7日-1891年12月29日)。换句话说,序列 <nγ> 击中[0,1)的每一个子区间(因此,它击中了如此无限的次数)。随着一个更深入的处理序列 <nγ> 的一致分布的定理引理的引出,我们将获得这个事实。
假如对于每一个 (a,b)⊂ [0,1) 的区间,
,
则称区间 [0,1) 中的数列 为等分布的(equidistributed)。其中,#A表示有限集合 A 的基数(cardinality)。这指的是,对于任意大的 N,在 (a,b) 中 n ≤ N 的这部分数
等于区间(a,b)的长度对区间[0,1)的长度的比率。换句话说,序列
均匀地(evenly)扫过整个区间,每个子区间都获得了它的公平份额(fair share)。显然,序列的次序非常重要,正如如下两个例子所示。
例子1: 序列
似乎是均匀分布的,因为它非常均匀地经过区间[0,1)。显然,这不是一个证明,观迎读者给出其证明。对于某种程序上的相似例子,请参看练习 8 中 的情形。
例子2:
令 为区间[0,1)中任意比率数的枚举(enumeration)。则,由
定义序列不是等分布的,因为序列的一半在0处。因此,很明显,序列是稠密的。
现在,我们触及本节的主要定理。
定理 2.1 如果 γ 是非比数(irrational),则由小数部分构成的序列 <γ>,<2γ> ,<3γ> ,... 在区间 [0,1) 上是等分布的。
特别地,<nγ> 在区间 [0,1) 是稠密的,并且,我们将Kronecker定理作为一个推论。在图 2 中,我们演示了当 ,对于 N 的 3 个不同值,点 <γ>,<2γ> ,<3γ> ,...< Nγ> 构成的集合。
----------------------------------------------图 2 当 时的序列<γ>,<2γ> ,<3γ> ,...--------------
固定 (a,b)⊂ [0,1)并令 表示开区间 (a,b) 的特征函数,即,函数在开区间 (a,b)上等于1,在区间 [0,1) - (a,b) 上等于 0 。我们可以将这个函数按周期性(周期为1)扩展到实数域 ℝ 上,且仍然用
表示这个扩展。则,作为这个定义的结果,我们求得
,
并且,定理可以重述为表达式
( 当
时 ) 。
这一步消除了计算小数部分的困难,并将数论简化为分析。
引理2.2 假如 f 是一个周期为 1 的连续函数,且 γ 是一个非比数(irrational),则
当 时 ,
。
引理的证明分 3 步进行。
第 1 步:我们首先检查当 f 是指数函数 1 , , ... ,
, ... 时的这种情况下的极限有效性,假如 f = 1 ,则这个极限显然成立。假如
且 k ≠ 1,则积分等于 0。因为 γ 是非比数,我们有
,因此,
,
当 时,上式趋近于 0 。
第 2 步:很显然,假如 f 和 g 满足引理,则对于任意 A,B ∈ ℂ ,A f + B g 也满足引理。因此,第 1 步意味着对于所有三角多项式,引理为真。
第 3 步:令 ε > 0 。假如 f 是任意周期为 1 的连续周期函数,选取一个三角多项式P 使得 (根据第2章推论5.4,这是可行的)。则,根据第 1 步,对于任意大的 N,我们有
。
因此,
,
因此,引理得证。
现在,我们可以完成定理的证明。选取两个周期为 1 的连续的周期函数 和
,其分别在区间[0,1)上从上限和下限趋近于函数
;
和
都以 1 为界且除了在总长为 2ε 的区间以外与函数
一致(见图3)。
特别地, 且
。
假如
, 则,我们得到
。
------------------------------------------图3 函数 的逼近值--------------------------------------------
因此,
且
。
因为,对于任意 ,
存在则一定等于
,这个表述成立。这就完成了等分布定理的证明。
推论2.3 引理 2.2 的结论对于区间[0,1]上的任意 Riemann 可积函数以及周期为 1 的连续周期函数 f 都成立。
证明:
假设 f 是实数据函数,考虑区间[0,1]的一个划分,比如, 。接下来,定义
( 假如
且
(假如
) 。 则,显然
,并且
。
此外,通过使划分足够好,我们可以确保对于给定的 ,有
。
然而,按照这个定理
,
因为每一个 都区间上特征函数的有限线性组合;类似地,我们有
。
从这两个论断,通过使用前面的逼近论据,我们可以推断出引理的证明。
对于这个引理及其推论,根据简单的动态系统,存在一个有趣的解释。在这个例子中,底层空间是由角度 θ 参数化的圆。我们也考虑一个空间到其自身的映射:在这里,我们按角度 2πγ 选取一个圆的旋转 ρ ,即,变换 。
接下来,我们希望考虑的是,具有底层行为 ρ 的这个空间,如何随着时间演进。换句话说,我们希望考虑 ρ 的迭代,即, 其中,
,
且,在其中,我们认为行为 恰好在 t = n 时刻发生。
对于圆周上的每一个Riemann可积函数 f ,我们也可以将其与旋转 ρ 的对应影响关联起来,并获得一个函数序列
且 。在这个特别的背景下,则这个系统的遍历性(ergodicity)表述为,只要 γ 是非比数,对于每一个 θ,“时间平均(time average)”
存在,并且等于“空间平均(space average)”
。
事实上,一旦我们改变变量 ,则这个论断不过是推论2.3的改述(rephrasing)。
回到等分布序列的问题,我们观察到,定理2.1的证明给出了下面的特征。
Weyl准则(Weyl's criterion) 对于区间[0,1)上的一个实数序列 ,当且仅当对于所有 k ≠ 0 的整数,具有
(当
时 ),
则,这个实数序列是等分布的。这个定理的一个方向在上面得到了有效证明,相反的方向可以在练习 7 中找到。特别地,我们发现,为了理解序列 的等分布属性,估算对应的“指数和(exponential sum)”
就已足够。例如, 使用Weyl准则可以证明,只要 γ 是非比数,序列
是等分布的。此例及其它例子可以在练习8和9,以及问题2和3中找到。
作为最后的评述,我们提到序列 <nγ> 的分布属性的一个很好的几何解释。假设正方形的边是反射镜,并且光线从正方形内部的一点射出。光线会走出怎样的路径?
----------------------------------------------图 4. 一束光在正方形中的反射--------------------------------------
为了求解这个问题,主要思想是考虑通过连续反射初始正方形在其边上形成的平面的网格。选取一个合适的坐标轴,在正方形中由光追踪的路径对应平面上的直线 。其结果,读者可以观察到,这个路径要么将是闭合且是周期的,要么将在正方形中是稠密的。这些场景中的第一个发生的条件是,当且仅当光线(相对于正方形的一侧确定)的初始方向的斜率 γ 是比率数(rational),而这个密度(density)可由Kronecker定理推导出。人们从等分布定理中能得出什么更有说服力的结论呢?
3. 一个连续但无处可微的函数(A continuous but nowhere differentiable function)
存在很多在某一点连续但不可微的典型函数,比如说,。要构造一个在任何给定的有限点集上不可微,甚至在包含可数个点的恰当集合上也是不可微的这样一个连续函数,几乎同样容易。一个更“微妙(subtle)”的问题是,是否存在一个无处可微的连续函数。在 1861年,Riemann 猜出了这个定义为
(5)
的函数是无处可微的。他之所以考虑这个函数,是因为它与将在第 5 章介绍的 ζ 函数有密切联系。Riemann 从未给出证明,但在他的一次演讲中提到了这个例子。这引发了 Weierstrass 的兴趣,他在试图寻找证明时遇到了连续但无处可微函数的第一个例子。比如说,0 < b < 1 ,且 a 是一个大于 1 的整数。在 1872年,他证明了,假如 ab > 1 + 3π/2 , 则函数
是无处可微的。
但是,如果不对Riemann原函数(original function)作最后说明(word),这个故是不完整的。在1916年,Hardy[há:di]( 1877 年 2 月 7 日 – 1947 年 12 月 1 日,英国数学家)证明了,R 在所有 π 的非比数倍数点处不可微。 然而,直到很久后的1969年,Gerver(美国数学家)才完全解决了这个问题;他首先证明,事实上,函数 R 在所有形如πp/q (p 和 q 是奇数整数)的 π 的比率数倍数点处是可微的,然后证明,R 在所有余下的情况下是不可微的。
在这一节,我们证明下面的定理。
定理3.1 假如 0 < α < 1 ,则, 函数
是连续但无处可微的。
此函数的连续性是显而易见的,因为级数绝对收敛。我们所需要的 f 的关键属性在于,它具有很多消失的Fourier系数。一个跳过很多项的Fourier级数,就像以上给出的一样,或者像 W(x)一样,这样的级数被称为缺项(lacunary[ ləˈkjuːnərɪ])Fourier级数。
定理的证明实际上是求Fourier级数和的3个方法的故事。第一种就去,根据部分和公式 , 存在常规的收敛性。第二种方法,存在Cesàro可和公式
且
是 Féjer 核。第三种方法,很显然与第二种有关,涉及定义为
的延迟均值(delayed means)。
因此, 。这些方法用可视化最好理解,如图5所示。
----------------------------------------图 5. 3种Fourier级数的求和方法--------------------------------------
假如 。则:
假如
,
以项
乘以 1 的形式出现;假如
,
以项
乘以 0 的形式出现。
假如
,
以
乘以
的形式出现;假如
,
以
乘以 0 的形式出现。
假如
,
以项
乘以 1 的形式出现;假如
,
以项
乘以
的形式出现;假如
,
以项
乘以 0 的形式出现。
例如,注意
。
其它诊断的证明与此相似。
延迟均值有两个重要的特征。在一方面,它们的属性与Cesàro均值的(好)特征紧密相关。在另一方面,对于那些像 f 一样具有缺项属性的级数而言,延迟均值在本质上等于部分和。特别地,注意,对于我们的函数
(6) ,
其中, 是形如
且
的最大整数,通过考察图5 及 f 的定义,可以看出这是显而易见的。
我们回到严格意义上的定理证明,并采用反证法论证;即,我们假设 f 在某个 点存在导数
。
引理 3.2 令 g 为在 点可微的任意连续函数。则 Cesàro均值满足
, 因此,
。(译注:第2章提到过,这种 “O” 记法表示左侧以右侧的常数倍为界。)
证明:
首先,我们有
,
其中, 是 Fejér 核。因为
是周期的,我们有
且意味着
。
从假设 g 在 点可微我们得到
。
现在,观穿到 满足两个估算
和
。
对于第一个不等式,我们记得, 是一个次数为 N 的三角多项式,其系数以 1 为界。因此,
是一个次数为 N 的三角多项式,其系数不大于 N 。因此
。
对于第二个不等式,我们记得
。
微分这个表达式,我们得到两项:
。
然后,假如我们使用两个事实, 和
(对于| t | ≤ π) ,我们得到对
的预期估算。使用这些估计,我们求得
。
一旦我们调用 的定义,就完成了对这个引理的证明。
引理 3.3 假如 , 则
。
这个引理从前面的观察(6)可以推出,因为 和
。
现在,按第一个引理,我们有
,
且第二个引理也意味着
。
这与预期相悖,因为 比
增长得更快。
关于我们的这个函数 ,按顺序另加几条评述。
相对于上面的函数 R 和 W,这个函数是复数值的,因此, 无处可微并不意味着其实数部分和虚数部也同样无处可微。然而,通过少许修改我们的证明便可证明,事实上,
的实数部分
,
以及其虚数部分,两者都是无处可微的。为了理解这一点,首先通过第一个证明观察到,引理3.2具有下列泛化性:假如 g 是一个在 点可微的连续函数,则
(只要
) 。
则,我们继续以 推进,如上,注意到,
;其结果就是,假设 F 在
点可微,则当
且
时,我们得到
。
为了得出矛盾的逻辑,我们仅需选取 h 使得其满足 即可;通过设置 δ 使其等于从
到形如
的最近的数的距离(因此,
) ,并且取
,这样即可完成证明。
显然,当 α > 1 时,函数 是连续可微的,因为级数可以逐项微分。最后,我们已经证明了当 α < 1 时函数无处可微,事实上这可以推广到 α = 1的情况,只需通过论据的合适改进即可(见第5章问题8)。事实上,使用这些更详尽的方法,读者也可以证明,对于 Weierstrass 函数,假如 ab ≥ 1 ,则其是无处可微的。
4. 圆周上的热传导方程(The heat equation on the circle)
作为最后的例释(illustration),我们回到最初由Fourier考虑的热传导(diffusion)问题。
假如给我们一个环上位于 t = 0 时刻的初始温度分布,并要求我们描述环上 t > 0 的某时刻的温度。
这个环是一个单位圆的模型。圆上的一点用其角度 θ = 2πx 描述,其中变量 x 位于0 到 1 之间。假如用 u(x,t ) 表示用角度 θ 描述的一点在 t 时刻的温度,则类似第 1 章中给出来的考量结果,u 满足微分方程
(7) 。
常量 c 是一个正的物理常量,其取决于制造环的材料(参见第1章第2.1节)。在重新调节时间变量 t 之后,我们可以假设 c = 1 ,假如 f 是我们的初始数据,我们施加一个已知条件
。
为了解这个问题,我们分离变量并寻找形如
的特解。
然后,则 u 的表达式插入热传导方程,我们得到
。
因此,等式两侧都是常量,比如说,设其等于 λ 。因为 A 一定是周期为 1 的周期函数,我们发现,唯一的可能性是 , 其中,n ∈ ℤ 。则 A 是指数函数
和
的线性组合,而 B ( t )是函数
的倍数。通过叠加这些解,导出方程
(8) ,
其中,设置 t = 0 , 我们发现 是 f 的 Fourier系数。
注意,当 f 是 Riemann可积的时候,系数 是有界的。因为因子
以极快的速度趋近于 0 , 其定义的 u 这个级数收敛。事实上,在这种情况下,u 是二次可微的,且作为方程 (7) 的解。
针对边界条件,很自然地出现的问题如下:当 t ⟶ 0 的时候,u(x,t) ⟶ f ( x ) ?如果是,又是在什么意义上?简单地应用Parseval恒等式证明,在均值平方的意义上,这个极限成立(练习11)。为了更好地理解方程解(8)的属性,我们将其写为
,
其中, 圆的热传导核(heat kernel for the circle),由
(8)
给出,其中,周期为 1 的函数的卷积定义为
。
热传导核和(第2章中给出的)Poisson核之间的类比在练习12中给出。然而,与Poisson核的情况不同,热传导核不存在基本公式。因此,证明它是一个好核(在第2章好核的意义上)。证明并不明显,需要使用著名的Poisson求和公式,这将在第 5 章中介绍。作为一个推论,我们也将会发现, 处处为正,从其定义表达式(9)来看,这也是一个不明显的事实。假如我们以处处小于 0 初始温度分布 f 开始。则对于所有的 t ,预期 u(x,t ) ≤ 0 在物理上是合理的,因为热量由热向冷传导(travels)。现在,
。
假如 对于某个
点为负值,则我们可以选取由附近的
所支持的 f ≤ 0 ,这将意味着
, 这就导致了出现矛盾。
内容来源:
<<Fourier Analysis: An Introduction>> E.M. Stein & R. Shakarchi
相关文章:
Fourier分析导论——第4章——Fourier级数的一些应用(E.M. Stein R. Shakarchi)
第 4 章 傅里叶级数的一些应用 Fourier series and analogous expansions intervene very naturally in the general theory of curves and surfaces. In effect, this theory, conceived from the point of view of analysis, deals obviously with the study of arbitra…...

c语言使用fdk_aac库对aac音频解码为pcm
//示例为adts的aac流数据(adts数据可以每一包都可以独立解析不需要拼凑) //解码数据的采样率同解码前的采样率,如果不满足需求,需要对数据进行重采样 #include <aacdecoder_lib.h>int m_fd -1; int m_fd2 -1;void aac2pc…...

zustand管理工具--React
npm i zustand 1.函数参数必须返回一个对象 对象内部编写状态数据和方法 2.set是用来修改数据的专门方法必须调用它来修改数据 import { useEffect } from "react"; import { create } from "zustand";// 1. 创建store const goodsStore create((set) …...

Elasticsearch内存分析
文章目录 Elasticsearch JVM内存由哪些部分组成Indexing BufferNode Query CacheShard Request CacheField Data CacheSegments Cache查询 非堆内存内存压力mat分析es的jvm缓存监控 Elasticsearch JVM内存由哪些部分组成 官方建议Elasticsearch设置堆内存为32G,因为…...

Alert警告提示(antd-design组件库)简单使用
1.Alert警告提示 警告提示,展现需要关注的信息。 2.何时使用 当某个页面需要向用户显示警告的信息时。 非浮层的静态展现形式,始终展现,不会自动消失,用户可以点击关闭。 组件代码来自: 警告提示 Alert - Ant Design 3…...

Linux提权方法总结
1、内核漏洞提权 利用内核漏洞提取一般三个环节:首先对目标系统进行信息收集,获取系统内核信息及版本信息 第二步,根据内核版本获取对应的漏洞以及exp 第三步,使用exp对目标进行攻击,完成提权 注:此处可…...

力扣第300题 最长递增子序列 c++ 动态规划题 附Java代码
题目 300. 最长递增子序列 中等 相关标签 数组 二分查找 动态规划 给你一个整数数组 nums ,找到其中最长严格递增子序列的长度。 子序列 是由数组派生而来的序列,删除(或不删除)数组中的元素而不改变其余元素的顺序。例…...

Si3262 集成低功耗SOC 三合一智能门锁应用芯片
Si3262 是一款G度集成的低功耗 SOC 芯片,其集成了基于 RISC-V 核的低功耗MCU 和工作在 13.56MHz 的非接触式读写器模块。 读写器模块支持 ISO/IEC 14443 A/B/MIFARE 协议,支持自动载波侦测功能(ACD)。无需外W其他电路,…...

linux rsyslog介绍
Rsyslog网址:https://www.rsyslog.com/ Rsyslog is the rocket-fast system for log processing. It offers high-performance, great security features and a modular design. While it started as a regular syslogd, rsyslog has evolved into a kind of swis…...

项目部署之安装和配置Canal
1.Canal介绍 Canal是阿里巴巴的一个开源项目,基于java实现,整体已经在很多大型的互联网项目生产环境中使用,包括阿里、美团等都有广泛的应用,是一个非常成熟的数据库同步方案,基础的使用只需要进行简单的配置即可。 …...

基于Skywalking的全链路跟踪实现
在前文“分布式应用全链路跟踪实现”中介绍了分布式应用全链路跟踪的几种实现方法,本文将重点介绍基于Skywalking的全链路实现,包括Skywalking的整体架构和基本概念原理、Skywalking环境部署、SpringBoot和Python集成Skywalking监控实现等。 1、Skywalki…...

Spark Core
Spark Core 本文来自 B站 黑马程序员 - Spark教程 :原地址 第一章 RDD详解 1.1 为什么需要RDD 分布式计算需要 分区控制shuffle控制数据存储、序列化、发送数据计算API等一系列功能 这些功能,不能简单的通过Python内置的本地集合对象(如…...

[算法日志]图论: 广度优先搜索(BFS)
[算法日志]图论: 广度优先搜索(BFS) 广度优先概论 广度优先遍历也是一种常用的遍历图的算法策略,其思想是将本节点相关联的节点都遍历一遍后才切换到相关联节点重复本操作。这种遍历方式类似于对二叉树节点的层序遍历,即先遍历完子节点后…...

Xilinx FPGA SPIx4 配置速度50M约束语句(Vivado开发环境)
qspi_50m.xdc文件: set_property BITSTREAM.GENERAL.COMPRESS TRUE [current_design] set_property BITSTREAM.CONFIG.SPI_BUSWIDTH 4 [current_design] set_property BITSTREAM.CONFIG.CONFIGRATE 50 [current_design] set_property CONFIG_VOLTAGE 3.3 [curren…...

Linux Shell和权限
目录 Shell命令及运行原理 权限 1.文件基本属性 2.文件权限值的表示方法 3.文件访问权限的相关设置方法 3.(1)chmod 组名修改 3.(2)chmod 二进制修改 3.(3)chown 3.(4)chgrp 3.(5)umask 4.目录权限 Shell命令及运行原理 Linux的操作系统,狭义上是…...

Git同时配置Gitee和GitHub
Git同时配置Gitee和GitHub 一、删除原先ssh密钥二、生成密钥 这里的同时配置是针对于之前配置过单个gitee或者github而言的,如果需要看git从安装开始的配置,则可以看这一篇文章 git安装配置教程 一、删除原先ssh密钥 在C盘下用户/用户名/.ssh文件下找到…...

IGP高级特性简要介绍(OSPF-上篇)
OSPF高级特性 一、OSPF_提升故障收敛及网络恢复速度 1.FRR与BFD快速恢复故障 1.1 FRR 在传统转发模式下,当到达同一个目的网络存在多条路由时,路由器总是选择最优路由使用,并且下发到FIB表指导数据转发。 当最优路由故障时,需…...

Oracle-Ogg集成模式降级为经典模式步骤
前言: Ogg集成模式降级为经典模式的场景比较少,因为降级为经典模式会导致无法支持压缩表同步,XA事务,多线程模式,PDB模式同步等功能,除非遇到集成模式暂时无法解决的bug或者环境不支持集成模式,比如DG备库环…...

链表面试OJ题(1)
今天讲解两道链表OJ题目。 1.链表的中间节点 给你单链表的头结点 head ,请你找出并返回链表的中间结点。 如果有两个中间结点,则返回第二个中间结点。 示例 输入:head [1,2,3,4,5] 输出:[3,4,5] 解释:链表只有一个…...

[极客大挑战 2019]Upload 1
题目环境: 根据题目和环境可知此题目是一道文件上传漏洞 编写一句话木马脚本<?php eval($_POST[shell]);?>将脚本文件更改为jpg图片文件我这里是flag.jpg上传文件并burpsuite抓包Repeater重放 报错一句话木马里面有<?字符 换一种一句话木马继续编写木马…...

OpenFeign讲解+面试题
一:OpenFeign是什么? 是一个声明式的web客户端,只需要创建一个接口,添加注解即可完成微服务之间的调用 二:调用微服务的方式? ribbon restTemplate方式调用openFeign通过接口注解的方式调用 三:如何使用OpenFeign&…...

嬴图 | LLM+Graph:大语言模型与图数据库技术的协同
前言 2022年11月以来,大语言模型席卷全球,在自然语言任务中表现卓越。尽管存在一系列伦理、安全等方面的担心,但各界对该技术的热情和关注并未减弱。 本文不谈智能伦理方面的问题,仅集中于Ulitpa嬴图在应用中的一些探索与实践&a…...

微信小程序下载文件和转发文件给好友总结
这段时间公司让我负责小程序的一些功能开发,回想上次开发小程序还是在上一次,这次开发小程序主要实现的功能就是转发文件给好友和下载文件,总结一下这次遇到的各种问题和解决方法。 下载文件 首先正常下载 wx.downloadFile({url: https://img.haihaina.cn/月度支出表.xls,…...

一文掌握 Apache SkyWalking
Apache SkyWalking SkyWalking是一个开源可观测平台,用于收集、分析、聚合和可视化来自服务和云原生基础设施的数据。SkyWalking 提供了一种简单的方法来保持分布式系统的清晰视图,甚至跨云。它是一种现代APM,专为云原生、基于容器的分布式系…...

外贸网站优化常用流程和一些常识
外贸网站google排名,总以为是单个网页标签的优化过程。 显然,这些观点都是错误的,九凌网络是做谷歌优化服务,九凌网络跟大家分享外贸网站Google优化常用流程和一些常识需要做以下几个步骤: 第一步:网站诊断࿰…...

Hive的时间操作函数
目录 前言函数使用介绍实际使用判断该天是星期几判断该天对应的周(包含一周开始和结束) 前言 hive 里面的时间函数有很多,今天单讲dayofweek函数,背景:有时候不仅要出日报,还要出周报,需要很多…...

【Web安全】CORS跨域资源共享漏洞
文章目录 前言一、漏洞概述二、漏洞原理三、CORS响应头类型四、漏洞挖掘五、修复建议前言 本篇文章主要介绍CORS跨域漏洞产生的原理,漏洞复现过程,挖掘手段以及如何进行修复,文章难免会有失误,烦请留下宝贵建议,谢谢! 一、漏洞概述 跨域资源共享(CORS)是一种浏览器机制…...

IntelliJ IDEA 如何修改默认Maven仓库地址
在使用idea过程中,每次新建项目或者打开项目时,maven仓库地址都会变为默认地址。如何修改默认地址,让其保持不变,如下这种方式可以简单快捷的设置。 1.打开idea,取消项目自动加载 2.点击 Customize,然后再点击 All se…...

Vue3 <script setup>是什么?作用?
结论先行: <script setup> 是 Vue3 的语法糖,简化了组合式 API 的写法,实现了 “顶层的绑定”。例如: ① 声明的属性和方法无需 return,就可以直接在模板使用; ② 引入组件的时候,会自…...

2.9 CSS 响应式布局
1.媒体:media 媒体类型: all:检测所有设备。screen:检测电子屏幕,包括:电脑屏幕、平板屏幕、手机屏幕等。print:检测打印机 媒体特性: width:检测视口宽度。max-width:…...