Q

Matlab Mejier G函數的負參數

2017-03-17 64 views 1 likes

1

是否有可能meijerG函數包含負值（即{-1,0,0}）？我嘗試了Mathematica和Matlab來計算這個meijerG函數，但是他們產生了一個錯誤，這個錯誤代碼是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;

在這裏，我還附上該函數的圖像與文字。

2017-03-17 Baidal Kocham

A

回答

0

從documentation of meijerG：

沒有對參數ai - bj, i = 1, …, n. j = 1, …, m，應該由一個正整數不同[...]。否則，meijerG返回一個錯誤。

複數對任何係數均有效;但在你的情況下，你有a0-b0 = 1這是禁止的。

2017-03-17 08:38:49 Bentoy13

+0

謝謝， Bentoy，爲您解釋。它清除了很多混亂。那麼從你的評論中，我可以推測原始功能是不正確的？上面我也附加了文本中函數的圖像。如果你能對此發表評論並提出一些建議，我會很感激。錯誤？ –

+0

對不起，我對這個特殊功能一無所知。我無法幫助你。也許嘗試尋找meijer G和Bessel K_0函數之間的關係，您可能會發現正確的參數... – Bentoy13

0

我很快看了看。如果一個擴大的log 2（1 + x）的成泰勒級數，替代\γ->的x^2，則積分將是

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

見here的詳細信息。 G是伽馬函數，對於像（k + 1/2）這樣的參數，它通過二項式coeff乘以sqrt（\ pi）來表示，詳情請參閱here。

你的所有，與多項式超過lambdas和b和一些二項式coefs和\ PI等無論是可以概括與否方面的無窮和後 - 我不知道...

2017-04-03 22:02:30

相關問題