Документация Engee

Подпространственные углы

В этом примере рассматриваются углы между соседними многочленами Лагерра с возмущенной мерой:

где внутреннее произведение определяется .

Для этого свяжем многочлены Лагерра с нормализованными обобщенными многочленами Лагерра, связанными с возмущенной мерой. Из внутреннего произведения коэффициентов связности следует, что:

using FastTransforms, LinearAlgebra

Соседний индекс k и максимальная степень N-1:

k, N = 1, 11
(1, 11)

Параметры связности Лагерра:

α, β = 0.0, 0.125
(0.0, 0.125)

Предварительно вычислим план Лагерра-Лагерра:

P = plan_lag2lag(Float64, N, α, β; norm2=true)
FastTransforms Laguerre--Laguerre plan for 11-element array of Float64

Применим план к тождеству, затем — план сопряженности:

VtV = parent(P*I)
lmul!(P', VtV)
11×11 Matrix{Float64}:
  0.941743    -0.117718    -0.0515016  -0.0321885  -0.0231355  -0.01793    -0.0145681  -0.0122268  -0.0105074  -0.00919399  -0.00815966
 -0.117718     1.07418     -0.125995   -0.0539157  -0.0333201  -0.0237861  -0.0183502  -0.0148608  -0.0124417  -0.0106716   -0.00932328
 -0.0515016   -0.125995     1.14505    -0.131707   -0.0557687  -0.0342413  -0.024336   -0.0187146  -0.0151194  -0.0126345   -0.0108205
 -0.0321885   -0.0539157   -0.131707    1.19444    -0.136106   -0.0572798  -0.0350206  -0.0248132  -0.0190369  -0.0153515   -0.0128093
 -0.0231355   -0.0333201   -0.0557687  -0.136106    1.23272    -0.139703   -0.0585599  -0.0356975  -0.0252354  -0.0193261   -0.015562
 -0.01793     -0.0237861   -0.0342413  -0.0572798  -0.139703    1.26415    -0.142756   -0.059673   -0.0362968  -0.0256145   -0.0195886
 -0.0145681   -0.0183502   -0.024336   -0.0350206  -0.0585599  -0.142756    1.29092    -0.145414   -0.0606594  -0.0368352   -0.0259588
 -0.0122268   -0.0148608   -0.0187146  -0.0248132  -0.0356975  -0.059673   -0.145414    1.31429    -0.147773   -0.0615464   -0.0373246
 -0.0105074   -0.0124417   -0.0151194  -0.0190369  -0.0252354  -0.0362968  -0.0606594  -0.147773    1.33507    -0.149896    -0.062353
 -0.00919399  -0.0106716   -0.0126345  -0.0153515  -0.0193261  -0.0256145  -0.0368352  -0.0615464  -0.149896    1.3538      -0.151829
 -0.00815966  -0.00932328  -0.0108205  -0.0128093  -0.015562   -0.0195886  -0.0259588  -0.0373246  -0.062353   -0.151829     1.37089

Из этой матрицы извлекаются углы:

θ = [acos(VtV[n, n+k]/sqrt(VtV[n, n]*VtV[n+k, n+k])) for n in 1:N-k]
10-element Vector{Float64}:
 1.6881063520094897
 1.6846487354051982
 1.6836556387180577
 1.6831998143073255
 1.6829427202726823
 1.68277949501306
 1.682667534322365
 1.6825864250312665
 1.6825252245704025
 1.682477567030148