3
高度に振動する積分の数値評価
で複素関数論のアプリケーションでは、この上級コースの演習の高い振動不可欠で1ポイントで I(λ)=∫∞−∞cos(λcosx)sinxxdxI(λ)=∫−∞∞cos(λcosx)sinxxdxI(\lambda)=\int_{-\infty}^{\infty} \cos (\lambda \cos x) \frac{\sin x}{x} d x 複素平面で点法を使用して、λλ\lambda大きな値に対して近似する必要があります。 振動性が高いため、この積分は他のほとんどの方法を使用して評価するのは非常に困難です。これらは、異なるスケールでのλ=10λ=10\lambda = 10被積分関数のグラフの2つのフラグメントです。 一次漸近近似は I1(λ)=cos(λ−14π)2πλ−−−√I1(λ)=cos(λ−14π)2πλI_{1}(\lambda) = \cos \left(\lambda-\frac{1}{4} \pi\right) \sqrt{\frac{2 \pi}{\lambda}} そして、さらに(はるかに小さい)改良が用語を追加します I2(λ)=18sin(λ−14π)2πλ3−−−√I2(λ)=18sin(λ−14π)2πλ3I_2(\lambda)=\frac{1}{8} \sin \left(\lambda-\frac{1}{4} \pi\right) \sqrt{\frac{2 \pi}{\lambda^{3}}} λλ\lambda関数としての近似値のグラフは次のようになります。 ここで質問が来ます。近似がどれほど良いかを視覚的に見るために、積分の「実際の値」と比較するか、より正確には独立したアルゴリズムを使用して同じ積分の良い近似と比較したいと思います。サブリーディング補正の小ささにより、これは非常に近いものになると予想されます。 λλ\lambdatanh(sinh)tanh(sinh)\tanh(\sinh) 最後に、実装した重要度サンプルを使用してモンテカルロインテグレーターで運を試しましたが、安定した結果を得ることができませんでした。 λ>1λ>1\lambda > 1