周勇攀

(武漢工程大學(xué)計(jì)算機(jī)科學(xué)與工程學(xué)院，湖北 武漢 430074)

?

基于Cauchy主值積分的高振蕩函數(shù)的數(shù)值計(jì)算

周勇攀

(武漢工程大學(xué)計(jì)算機(jī)科學(xué)與工程學(xué)院，湖北 武漢 430074)

摘要：含高振蕩函數(shù)的Cauchy主值積分dx,-1<τ<1，基于解析延拓定理，用最速下降法將其轉(zhuǎn)化成在[0,+∞)上非振蕩且指數(shù)快速下降的積分，再利用Gauss-Laguerre求積法則高效逼近計(jì)算，最后用兩個(gè)數(shù)值實(shí)例來(lái)說(shuō)明該方法的合理性.

關(guān)鍵詞：Cauchy主值積分；最速下降法；Gauss-Laguerre求積法則

0引言

我們討論含高振蕩函數(shù)的Cauchy主值積分

(1.1)

其中，f在包含[-1,1]的一個(gè)充分大的復(fù)區(qū)域內(nèi)解析，如果f在[-1,1]上滿(mǎn)足霍爾德條件[4]，我們知道這個(gè)積分存在[5].對(duì)式(1.1)，用一個(gè)半圓包含奇異點(diǎn)τ的鄰域(圖1)，結(jié)合復(fù)積分方法和最速下降法得到4個(gè)在[0,+∞)上非振蕩且指數(shù)快速下降的積分，再用Gauss-Laguerre求積法則高效逼近計(jì)算[6-7].

圖1　原給異點(diǎn)τ的鄰

圖2　最速下降路徑

1計(jì)算結(jié)果

同理，

由于路徑Γ6是一個(gè)圍住z=τ的半圓，即z-τ=reiθ,0≤θ≤π,則

當(dāng)r→0時(shí)，|z-τ|→0，f(z)在點(diǎn)τ處是連續(xù)的，即|f(z)-f(τ)|→0,

此時(shí)，再由定義的主值積分和式(1.2)，得到

(1.3)

其中，t=ωp.

計(jì)算得到的積分我們就可以用Gauss-Laguerre求積法則來(lái)估計(jì)，含振蕩函數(shù)的Cauchy主值積分估計(jì)得到

(1.4)

其中，xk和wk分別為n階Gauss-Laguerre公式的節(jié)點(diǎn)和權(quán).

其中,ξ1、ξ2、ξ3、ξ4∈C，當(dāng)ω?1，誤差的漸近估計(jì)為O(ω-2n-1)，因此計(jì)算的誤差精度隨著ω的增加而快速提高.

定理1假設(shè)f和g在一個(gè)包含區(qū)間[-1,1]的充分大的復(fù)區(qū)域D內(nèi)解析，并且g的反函數(shù)在D內(nèi)存在，如果下面的條件在D內(nèi)滿(mǎn)足：

?m∈N∶|f(z)|=O(|z|m),?ω0∈R∶|g-1(z)|=O(eω0|z|),|z|→∞,

即對(duì)于x∈[-1,1]，存在一個(gè)函數(shù)F(x)，使得F(x)=∫Γxf(z)eiωg(z)dz，其中Γx是一個(gè)起始于x的一個(gè)路徑，hx(p)是Γx的一個(gè)參數(shù)化表示，p∈[0,∞)，誤差E=F(x)-QF[f,g,hx]的漸近估計(jì)為O(ω-2n-1)，其中QF[f,g,hx]是由n階Gauss-Laguerre求積法則得到的[2].

再用這個(gè)公式，我們能導(dǎo)出誤差的一個(gè)表達(dá)式E=F(x)-QF[f,g,hx];

其中，ξ∈C，誤差的漸近估計(jì)為O(ω-2n-1)得證.

2數(shù)值實(shí)例

表1　Filon方法計(jì)算的絕對(duì)誤差

可以看出，用Filon方法計(jì)算積分，對(duì)于固定的頻率ω，逼近的精度隨插值節(jié)點(diǎn)數(shù)目n的增加而提高，并且頻率越大，精度提高的越快.

表2　利用n點(diǎn)Gauss-Laguerre積分法

計(jì)算的絕對(duì)誤差.

從表2可以看出，對(duì)于含Cauchy核的高振蕩的數(shù)值積分，用數(shù)值最速下降法來(lái)計(jì)算，我們也得到同樣的結(jié)論，即對(duì)于固定的頻率ω，逼近的精度隨Gauss-Laguerre節(jié)點(diǎn)數(shù)目n的增加而提高，并且頻率越大，精度提高的越快；另外與表1對(duì)比，我們只需要取較少的節(jié)點(diǎn)，就可以得到相同的誤差精度.

3參考文獻(xiàn)

[1] Milovanovic G V.Numerical calculation of integrals involving oscillatory and singular Kernels and some applications of quadratures[J].Comput Math Appl,1998,36(8):19-39.

[2] Huybrechs D,Vandewalle S.On the evaluation of highly oscillatory integrals by analytic continuation[J].SIAM J Numer Anal,2006,44(3):1026-1048.

[3] Wang Haiyong ,Xiang Shuhuang.On the evaluation of Cauchy principal value integrals of osc-illatory functions[J].Comput Math Appl,2010,234:95-100.

[4] 路見(jiàn)可.解析函數(shù)邊值問(wèn)題[M].2版.武漢：武漢大學(xué)出版社,2004.

[5] Davis P J ,Rabinowitz P.Methods of Numerical Integration[M].second edition.New York:Acade-mic Press 1984.

[6] Ablowitz M J,Fokas A S.Complex variables:introduction and applications[M].Cambridge,UK:Cambridge University Press,1997.

[7] Wong R.Asymptotic approximation of integrals[M].Philadelphia :SIAM,2001.

[8] Abramowitz M,Sterun I A.Handbook of Mathematical Functions[M].Washington DC:National Bureau of Standard,1964.

[9] Huybrechs D,Olver S.Highly oscillatory quadrature//Engquist B,Fokas T,Hairer E,A.Iserles Highly Oscillatory Problems[M].Cambridge,UK:Cambridge University Press,2009:25-50.

(責(zé)任編輯趙燕)

Value integrals of highly oscillatory functions based on the evaluation of Cauchy principal

ZHOU Yongpan

(School of Computer Science and Engineering,Wuhan Institute of Technology,Wuhan 430074,China)

Abstract:The problem of numerical evaluation of Cauchy principal value integrals of highly oscillatory functions dx,-1<τ<1，had been discussed.Based on analytic continuation and the steepest descent method ,the integrals can be transformed into the problem of integrating on[0,+∞) with the integrand that does not oscillate,and that decays exponentially fast,which can be efficiently computed by using the Gauss-Laguerre rule.The validity of the method has been demonstrated in the provision of two numerical experiments and their results.

Key words:Cauchy principal valueintegrals;steepest descent method;Gauss-Laguerre rule

中圖分類(lèi)號(hào):O241.38；O174.41

文獻(xiàn)標(biāo)志碼:A

DOI:10.3969/j.issn.1000-2375.2016.03.017

文章編號(hào):1000-2375(2016)03-0267-04

作者簡(jiǎn)介：周勇攀(1990-)，男，碩士生

收稿日期:2015-09-15