チェビシェフ多項式を使用したスペクトル法の難しさ


19

私は論文を理解しようとするのに少し苦労しています。この論文では、スペクトル法を使用して、結合ODEのシステムから得られる固有値を解きます。私の質問の核心に到達するのに十分であるので、私は今1つの方程式だけを書きます。

方程式は

V[r]=e(ν[r]+λ[r])ϵ[r]+p[r][(ϵ[r]+p[r])(eν[r]+λ[r])rW[r]]

派生物を実行して取得する

(Eq1) V=[ϵ+pϵ+p+r(ν+λ)+1]W+rW

論文によると、システムの平衡量)を次の形式のチェビシェフ多項式として展開できるはずです。(ϵ,p,ν,λ

B[r]=Σi=0biTi[y]12b0、ここでは多項式です。Mathematicaで書いたコードを使ってを取得する方法を知っています。また、、の領域はです。b i y = 2 r / R 1 r 0 R Ti[y]biy=2(r/R)1r(0,R)

また、論文では、関数()は 、および一般的にような用語は、F [ r ] = rV,WB[r]F[r]F[r]=(rR)lΣi=0fiTi[y]12f0B[r]F[r]

B[r]F[r]=12(rR)lΣi=0πiTi[y]12π0

ここで 及びのためのと1に等しいため。 Θ K = 0 、K < 0πi=Σj=0[bi+j+Θ(j1)b|i1|]flΘ(k)=0k<0k0

すべてのことを言って、次の平衡関数を作成するとしましょう

ϵ+pϵ+p=B1[r]および、Eq1はr(ν+λ)=B2[r]

(Eq2)((rR)lΣi=0viTi[y]12v0)=[B1[r]+B2[r]+1]((rR)lΣi=0wiTi[y]12w0)+rW

Question1:私は何をします用語?多項式は関数であるため、 X関数[y]のような拡張をどのように行うことができますか?また、方程式の両側でそれらを分割できるように思えますが、その用語の導入のポイントは何でしたか?論文によれば、この用語は、がゼロになるとがゼロになるという境界条件を課すことになっています。[y]B[r]=r(rR)l[y]VWrB[r]=(rR)lV,Wr

* Question2: *私が対処することになっていますどのようにで用語。この論文では、派生語の処理方法について説明していますが、自体についてはどうですか。それを平衡値のように扱い、ような用語にルールを使用するのか、このを用語で表現するのか。それとも、何か他のことをする必要がありますか?r W r B [ r ] F [ r ] r yrrWrB[r]F[r]ry


3
おそらく、参照している論文にリンクできますか?
アロンアーマディア

こんにちはAron、ここにリンク arxiv.org/pdf/gr-qc/0210102.pdfがあります。私が問題を抱えている数値的なものは付録Aに記述されており、上記で調べた式は式(19)です。ありがとう。
tau1777

自体は(したがって)の(線形)関数であることに注意してください。ry rrRr
クリスチャンクラソン

回答:


1

論文を詳細に読まずに質問に答えることが可能かどうかはわかりません。ただし、最初の質問に関しては、です。また、この要素はすべての項を乗算するわけではないため、分割できません。r/R=(y+1)/2

質問2の場合、この方程式はで境界条件を適用するために使用されるため、あなたが言及する用語は消えるべきだと思います。r=0

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.