Sturm-Liouville 问题的特征向量和特征值的数值解

计算科学 特征值 数字
2021-12-08 08:48:05

在我的研究中,我必须处理以下问题:

[1DFx]x+fDxcD2F=0

有边界条件

F(0)=0
Fx(L)=0

在哪里f是一个已知常数并且D(x)是一个已知函数x. 我应该澄清一下D(x)在数字上为所有人所知x(我有一个值向量),但它被称为封闭形式的分析函数,因此不容易近似。下标表示导数。

我注意到这是一个 Sturm-Liouville 问题:

[p(x)Fx(x)]x+q(x)F(x)=λr(x)F(x)

p=1/D,q=0, 和r=fDx/D2.

F是特征向量和λ=1/c是特征值。有没有一种简单的方法可以用数字解决它们?谢谢你的帮助。

1个回答

哦!我没有意识到这是多么容易。我在数值方法方面的实践经验以有限差分开始和结束。我发现这是一个特征值问题。我没有意识到,如果我用中心差分算子替换导数,这会简化为一个简单的矩阵特征值问题!

D1D2是居中的差分矩阵,并且与pr定义如上,我可以写:

[D2+diag(px/p)D1]F=λ[diag(r/p)]F

AF=λBF

狄利克雷条件x=0很容易内置到第一行AB和诺伊曼条件x=L也可以用鬼点内置。

作为测试,特征向量和特征值具有解析形式,当D(x)=aexp(2bx),并且我的有限差分方法与我需要的足够精度进行比较。

感谢您的评论。对不起,这太简单了。