我在极坐标下的二维无规行走中讨论到Rice分布:
(1)
当无规行走发生的“位置”()离原点很远时,值很大,而分布的极值也大至处于处,因此,能看到分布主要特征的范围也在附近。这时,式(1)中的第一类修正贝塞尔函数的值会很大,而自然指数项会很小,但整个式(1)的值仍然是适中的。在MATLAB输入上式计算时,会因为计算过程中涉及到很大的和很小的值,所以尽管最终结果的值是适中的,计算也会溢出。这时可以利用第一类修正贝塞尔函数的渐近展开1。当为定值、很大且时,
(2)
越大,级数会衰减得越快.如果足够大,可近似为。这个近似式并不改变很大的事实。但如果将式(2)代入式(1),
(3)
可见,无论多大,其实只跟与之差有关。式(3)只需要计算。由于Rice分布的范围主要在附近,因此与之差是很小的,因此这一计算不会溢出。
实际计算时,可先判断和的取值是否使式(1)的计算溢出了,若溢出,按式(3)计算时,可根据和的取值大小决定舍去级数的哪些尾项。按照MATLAB的计算极限,能让式(1)溢出的情况,也就必然能使式(3)的所有尾项都舍去了。