皮尔西积分

Integrand of Pearcey integral

Pearcey Integral 3D Maple plot

Pearcey Integral Maple contour plot

Pearcey Integral Maple density plot

皮尔西积分(Pearcey Integral)是一种在论述光的传播、光的衍射、分岔理論、突变理论以及关于特殊函数的渐进展开式的研究中常见的多鞍点积分，其定义为

$P(x,y)=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\exp (I*(t^4+x*t^2+y*t))$

=$-(1/2)*\sqrt{(}Pi)*exp(-(1/4*I)*y^2/x)*(limit(exp(I*t^4)*erf(I*x*t/\sqrt{(}-I*x)+(1/2*I)*y/\sqrt{(}-I*x)),t=infinity))/\sqrt{(}-I*x)+(1/2)*\sqrt{(}Pi)*exp(-(1/4*I)*y^2/x)*(limit(exp(I*t^4)*erf(I*x*t/\sqrt{(}-I*x)+(1/2*I)*y/\sqrt{(}-I*x)),t=-infinity))/\sqrt{(}-I*x)+(2*I)*\sqrt{(}Pi)*exp(-(1/4*I)*y^2/x)*(int(exp(I*t^4)*erf(I*x*t/\sqrt{(}-I*x)+(1/2*I)*y/\sqrt{(}-I*x))*t^3,t=-infinity..infinity))/\sqrt{(}-I*x)$

被积分函数的鞍点

对于大的$|x|$,皮尔逊积分的被积分函数左右各有三个鞍点。

在右平面的鞍点是

$rs1:=-(1/3)*\sqrt{(}6)*sinh((1/3)*arcsinh((3/4)*y*\sqrt{(}2)*\sqrt{(}3)/x^{(}3/2)))$

$rs2:=(1/3)*\sqrt{(}6)*sinh((1/3)*arcsinh((3/4)*y*\sqrt{(}2)*\sqrt{(}3)/x^{(}3/2))+(1/3*I)*Pi)$

$s3:=-(1/3*I)*\sqrt{(}6)*cosh((1/3)*arcsinh((3/4)*y*\sqrt{(}2)*\sqrt{(}3)/x^{(}3/2))+(1/6*I)*Pi)$

左平面的鞍点是：

$ls1:=-(1/3)*\sqrt{(}6)*sin((1/3)*arcsin((3/4)*y*\sqrt{(}2)*\sqrt{(}3)/x^{(}3/2))+(1/3)*Pi)$

$ls3:=(1/3)*\sqrt{(}6)*cos((1/3)*arcsin((3/4)*y*\sqrt{(}2)*\sqrt{(}3)/x^{(}3/2))+(1/6)*Pi)$

$ls2:=(1/3)*\sqrt{(}6)*sin((1/3)*arcsin((3/4)*y*\sqrt{(}2)*\sqrt{(}3)/x^{(}3/2)))$

分岔

Bifurcation of Pearcey Integral

皮尔西积分的分岔曲线为

$27*x^2=-8y^3$

斯托克斯曲线

Stokes set of Pearcey Integral

皮尔西积分的斯托克斯曲线为

$y^3=\frac{27}{4}(\sqrt{(}27)-5)x^2$

尖点突变

皮尔西积分的尖点突变区间

在（x,y)平面中，分岔曲线和斯托克斯曲线将平面分化为尖点突变区。

参考文献

• Frank Oliver, NIST Handbook of Mathematical Functions,2010,Cambridge University Press.
• R.B. Paris,D. Kaminski,Hyperasymptotic evaluation of the Pearcey integral via

Hadamard expansions.Journal of Computational and Applied Mathematics 190 (2006) 437–452