Matlab Mejier G関数の負のパラメータ

meijerG関数に負の値（つまり、{-1,0,0}）が含まれている可能性がありますか？私はこのmeijerG関数を計算するためにMathematicaとMatlabの両方を試しましたが、これはmeijerG is not defined for the given parametersというエラーを生成します。どんな助けも役に立つでしょう。ここでMatlab Mejier G関数の負のパラメータ

D = (0.6); 
lg1 = lg2 = 1; 
G = evalin(symengine, sprintf('meijerG([[0], []], [[-1,0,0], []],%f)',(D/(lg1*lg2)))); 
CD = -((2*D)/(lg1*lg2*(log(4))))*G;

は、私はまた、テキストからの機能の画像を添付しています

は、ここに私のコードです。 documentation of meijerGから

出典

2017-03-17 Baidal Kocham

：

パラメータai - bj, i = 1, …, n. j = 1, …, mのいいえ対は、正の整数だけ異なるべきです[...]。そうでない場合、meijerGはエラーを返します。

複合数は任意の係数に対して有効です。しかし、あなたはa0-b0 = 1が禁止されています。

出典

2017-03-17 08:38:49 Bentoy13

ありがとうございました、 Bentoy、あなたの説明のために。それは多くの混乱を解消した。あなたのコメントから、元の関数が間違っていると推測できますか？私はテキストから関数のイメージを添付しました。私はあなたがそれについてコメントして何かを提案することができれば本当に感謝しますか？エラー？ –

申し訳ありませんが、私はこの特別な機能について何も知らない。私はもっとあなたを助けることはできません。たぶんmeijer GとBessel K_0関数の関係を調べようとすると、正しいパラメータが見つかるかもしれません... – Bentoy13

私はすぐにそれを見ました。 1は、テイラー級数に（1 + x）のLOG2を展開すると、代替\ガンマ>のx^2、その後、積分は

S K0(x) x^m dx = 2^(m-1) G((m+1)/2)^2

だろう詳細はhereを参照してください。 Gはガンマ関数であり、（k + 1/2）のような引数に対しては、二項係数timesrt（\ pi）で表されます。詳細はhereを参照してください。あなたはそれを合計できたか否かlambdasとb上の多項式との用語の無限の合計と、いくつかの二項coefsと\パイなどを持っているすべての後

- 私は知らない...

出典

2017-04-03 22:02:30

Matlab Mejier G関数の負のパラメータ

答えて

関連する問題