Problem in integration function containing gauss hypergeometric function

2 次查看（过去 30 天）
Shweta Rajoria
Shweta Rajoria 2017-12-14
评论： Shweta Rajoria 2018-1-9
I am deriving the coverage probability of heterogeneous network with massive MIMO macro base station. During the process, I am stuck in the in hypergeometric function. The small part of that program is given below.
fun= @(x,y)imag(exp(pi.*1i.*(x).*y.^(-2).*hypergeom(1,[1/2,3/2],-1i.*x./(y.^4))))./x
c= @(y) integral(@(x) fun(x,y), 0.1,1,'ArrayValued',true)
fun1= @(y)((1/2)-((1/pi).*c(y))).*exp(-pi.*y.^2)
c1 = integral(fun1, 0.1,1)
I think iproblem is due to terms(y.^(-2) and and y.^4) in the fun. If any body have any idea please suggest. Any suggestions regarding this will be appreciated.
  7 个评论
显示 5更早的评论
Walter Roberson
Walter Roberson 2018-1-9
编辑：Walter Roberson 2018-1-9
To within the boundaries of double precision, c1 will be nan. The values involved are way way too big to work with in double precision.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Walter Roberson
Walter Roberson 2018-1-8
fun(8/10, 3/10) is about 2.974*10^1912918 .
fun(8/10, 2/10) is about 4.087*10^-24618734963806
fun(8/10, 0.18) is about 2.345*10^207906652598325658
So, fun is very steep, and has a lot of values that overflow floating point representation. You will not be able to integrate that using any normal representation.

类别

在 Help Center 和 File Exchange 中查找有关 Numerical Integration and Differential Equations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by