范德蒙德矩阵 DG Hestaven

计算科学 插值 正交 不连续-galerkin
2021-12-22 04:02:30

我试图从 Hesthaven 的书中理解节点和模态基础公式(节点间断 Galerkin 方法,Hesthaven,Jan S.,Warburton,Tim)。对于矩阵,它由 Jacobi 多项式计算得到N=2α,β=0

V = Vandermonde1D(2, x)

V =

0.7071   -1.2247    1.5811
0.7071         0   -0.7906
0.7071    1.2247    1.5811

x = (-1,0,1)LGL 点。

在第 48 页的书中,它说对于,表中的 d) 选项是用于具有 Legendre-Gauss-Lobatto 点的正交基础。N=6

现在,如果我从这里运行代码

[x,w,P]=lglnodes(2)

我明白了

x =
-1
 0
 1

这是LGL点和

w =
0.3333
1.3333
0.3333

权重。但后来我得到

P =
1.0000   -1.0000    1.0000
1.0000   -0.0000   -0.5000
1.0000    1.0000    1.0000

在该程序的评论中,它说

计算 Legendre-Gauss-Lobatto 节点、权重和 LGL Vandermonde 矩阵。LGL 节点是(1x2)PN(x)

但这与我从 Hesthaven 的原始程序中得到的 Vandermonde 不同。

现在,我有点明白这个矩阵是什么了。如果你有勒让德多项式P

p0(x)=1,p1(x)=x,p2(x)=12(3x21)

并替换 LGL 点(-1,0,1)我得到这个矩阵。但是如何得到范德蒙德矩阵呢?P

2个回答

您使用的是勒让德多项式的正交化版本,而 Hesthaven 不是。因子归一化,即2n+122(1)+12=1.22472(2)+12=1.5811

您需要使用“归一化”勒让德多项式而不是勒让德多项式。将勒让德多项式除以22n+1.