Skip to content

11.2.1 具有退化核的积分方程

如果一个积分方程的核 K(x,y) (kernel K(x,y) of an integral equation) 是两个单变量函数积的有限和,即一个函数仅依赖于 x ,而另一函数仅依赖于 y ,那么它被称为一个退化核(degenerate kernel), 或积核(product kernel).

1. 退化核情形的解

一个具有退化核的第二类弗雷德霍姆积分方程的解导致一个有限维方程组的解. 考虑积分方程

(11.4a)φ(x)=f(x)+λabK(x,y)φ(y)dy,

其中

(11.4b)K(x,y)=α1(x)β1(y)+α2(x)β2(y)++αn(x)βn(y).

诸函数 α1(x),,αn(x)β1(x),,βn(x) 在区间 [a,b] 上被给出,并被假设是连续的. 其次,假设函数 α1(x),,αn(x) 是线性无关的,即,具有常系数 ck 的等式

(11.5)k=1nckαk(x)0

对每个 x[a,b] 成立,仅当 c1=c2==cn=0 时. 否则, K(x,y) 可以被表示成较小数目乘积的和.


①原文将下式中的 f(ξ,y(ξ)) 误为 f(ξ,η(ξ)) . - 译者注


由 (11.4a) 和 (11.4b) 即得

(11.6a)φ(x)=f(x)+λα1(x)abβ1(y)φ(y)dy++λαn(x)abβn(y)φ(y)dy.

式中的这些积分不再是变量 x 的函数,它们是常数. 用 Ak 表示它们:

(11.6b)Ak=abβk(y)φ(y)dy,k=1,,n.

解函数 φ(x) ,如果存在,是扰动函数 f(x) 和诸函数 α1(x),,αn(x) 的一个线性组合之和:

(11.6c)φ(x)=f(x)+λA1α1(x)+λA2α2(x)++λAnαn(x).

2. 解的系数的计算

如下计算诸系数 A1,,An . 用 βk(x) 乘以方程 (11.6c),并在区间 [a,b] 上关于 x 计算其积分:

abβk(x)φ(x)dx=abβk(x)f(x)dx+λA1abβk(x)α1(x)dx+(11.7a)+λAnabβk(x)αn(x)dx.

根据 (11.6b),这个方程的左端等于 Ak . 利用下述记号

(11.7b)bk=abβk(x)f(x)dx,ckj=abβk(x)αj(x)dx,

则对于 k=1,,n ,有

(11.7c)Ak=bk+λck1A1+λck2A2++λcknAn.

有可能不能计算 (11.7b) 中那些积分的精确值. 一旦出现这种情形, 必须用给出在第 1252 页 19.3 中的那些公式之一来计算它们的近似值. 线性方程组 (11.7c) 对于未知值 A1,,An 包含 n 个方程:

(1λc11)A1λc12A2λc1nAn=b1,(11.7d)λc21A1+(1λc22)A2λc2nAn=b2,λcn1A1λcn2A2+(1λcnn)An=bn.

3. 解、本征值、本征函数的分析

从线性方程组的理论知道,(11.7d) 有且仅有一个解 A1,,An ,如果系数矩阵行列式不等于零, 即

(11.8)D(λ)=|(1λc11)λc12λc1nλc21(1λc22)λc2nλcn1λcn2(1λcnn)|0.

显然, D(λ) 不恒为零,因为有 D(0)=1 . 因而,存在一个数 R>0 ,使得当 |λ|<RD(λ)0 . 为了进一步的研究,考虑两种不同的情形.

情形 D(λ)0

积分方程恰好有形如(11.6c)的一个解,并且诸系数 A1,,An 由方程组 (11.7d) 的解给出. 如果 (11.4a) 是一个齐次积分方程,即 f(x)0 ,则 b1=b2= =bn=0 . 此时齐次方程组 (11.7d) 只有平凡解 A1=A2==An=0 . 在这个情形只有函数 φ(x)0 满足积分方程.

情形 D(λ)=0

D(λ) 是一个不高于 n 次的多项式,因而它至多有 n 个根. 对于 λ 的这些值, b1=b2==bn=0 时的齐次方程组 (11.7d) 还有非平凡解,因而除了平凡解 φ(x)0 外,齐次方程组有别的形如

φ(x)=C(A1α1(x)+A2α2(x)++Anαn(x))(C 是一个任意常数 )

的解. 由于 α1(x),,αn(x) 是线性无关的,因此 φ(x) 不恒等于零. D(λ) 的根被称为积分方程的本征值(eigenvalues). 齐次积分方程相应的非零解被称为属于本征值 λ 的本征函数 (eigenfunctions). 几个线性无关的本征函数可以属于同一个本征值. 如果积分方程有一个一般的核, 那么就要考虑齐次积分方程有非平凡解的所有那些本征值 λ . 有些作者把满足 D(λ)=0λ 称为特征数,把 μ=1λ 称为相应于形如 μφ(x)=abK(x,y)φ(y)dy 方程的本征值.

4. 转置积分方程

现在有必要来研究在什么条件下,如果 D(λ)=0 ,齐次积分方程有解? 为此目的, 引进 (11.4a) 的转置积分方程 (transposed integral equation)(或者, 在复情形, 伴随积分方程(adjoint integral equation)):

(11.9a)ψ(x)=g(x)+λabK(y,x)ψ(y)dy.

λ 是一个本征值, φ(x) 是齐次积分方程 (11.4a) 的一个解. 容易证明 λ 也是伴随积分方程的一个本征值. 现在 (11.4a) 的两端乘以齐次伴随积分方程的任一解 ψ(x) , 并在区间 [a,b] 上关于 x 求积分:

(11.9b)abφ(x)ψ(x)dx=abf(x)ψ(x)dx+ab(λabK(x,y)ψ(x)dx)φ(y)dy.

假设 ψ(y)=λabK(x,y)ψ(x)dx ,则有 abf(x)ψ(x)dx=0 .

这就是: 对于某个本征值 λ 齐次积分方程 (11.4a) 有解,当且仅当扰动函数 f(x) 与齐次伴随积分方程的属于同一个 λ 的每个非零解正交(orthogonal). 这个陈述不仅对于具有退化核的积分方程成立, 而且对于具有一般核的积分方程也成立.

A: φ(x)=x+1+1(x2y+xy2xy)φ(y)dy,α1(x)=x2,α2(x)=x,α3(x)= x,β1(y)=y,β2(y)=y2,β3(y)=y.3 个函数 αk(x) 是线性相关的. 因而把原积分方程写为 φ(x)=x+1+1[x2y+x(y2y)]φ(y)dy . 对于这个积分方程,有 α1(x)=x2,α2(x)=x,β1(y)=y,β2(y)=y2y . 如果有解 φ(x) 存在,则它有形式 φ(x)=x+A1x2+A2x .

c11=1+1x3dx=0,c12=1+1x2dx=23,c21=1+1(x4x3)dx=25,c22=1+1(x3x2)dx=23,b1=1+1x2dx=23b2=1+1(x3x2)dx=23.

利用这些值来确定 A1A2 的方程组,则有: A123A2=23,25A1+(1+ 23)A2=23 ,由它即得 A1=1021,A2=27 ,因而 φ(x)=x+1021x227x=1021x2+57x

B: φ(x)=x+λ0πsin(x+y)φ(y)dy ,即 K(x,y)=sin(x+y)=sinxcosy+cosxsiny,φ(x)=x+λsinx0πcosyφ(y)dy+λcosx0πsinyφ(y)dy.

c11=0πsinxcosxdx=0,c12=0πcos2xdx=π2,b1=0πxcosxdx=2,

c21=0πsin2xdx=π2,c22=0πcosxsinxdx=0,b2=0πxsinxdx=π.

利用这些值,方程组 (11.7d) 即为 A1λπ2A2=2,λπ2A1+A2=π . 因为对满足 D(λ)=|1λπ2λπ21|=1λ2π240 的任意 λ,A1A2 的方程组有唯一解. 因而 A1=(λπ222)/(1λ2π24),A2=π(1λ)/(1λ2π24) ,所以此时积分方程的解为 φ(x)=x+λ[(λπ222)sinx+π(1λ)cosx]/(1λ2π24) . 积分方程的本征值为 λ1=2π,λ2=2π .

齐次积分方程 φ(x)=λk0πsin(x+y)φ(y)dy 有形如 φk(x)=λk(A1sinx+ A2cosx)(k=1,2) 的非零解. 对于 λ1=2π ,有 A1=A2 ,因而即得 φ1(x)= A(sinx+cosx) ,其中 A 为任意常数. 类似地,对于 λ2=2π ,有 φ2(x)=B(sinx cosx) ,其中 B 为任意常数.

注 这个先前的解法相当简单, 但它只在退化核的情形有效. 这个方法也可用于在一般核的情形中获得一个好的近似解, 如果可以用退化核足够好地逼近一般核.

version 1.24.0