Evaluating a complex integral

38 ビュー (過去 30 日間)
Michael Devereux
Michael Devereux 2019 年 2 月 14 日
回答済み: Abhishek Hullur 2021 年 8 月 8 日
Hello I'm trying to integrate the following function in MATLAB
but it's returing the wrong answer when I try something like
This is what I have tried so far:
fun = @(t,x,y) exp(1i.*(t.^4+x.*t.^2+y.*t));
P = @(x,y) integral(@(t)fun(t,x,y),-Inf,Inf);
P(1,1)
Any help appreciated and many thanks in advance
  3 件のコメント
Torsten
Torsten 2019 年 2 月 15 日
exp(i*(t^4+x*t^2+y*t)) does not tend to 0 as | t| -> Inf. Thus your integral does not exist (at least in the usual sense).
Michael Devereux
Michael Devereux 2019 年 2 月 15 日
According to WolframAlpha the answer is 1.20759 + 0.601534 i
Keep in mind it's a complex exponential so there is a finite solution. This is know as the Pearcey Integral. I am more concerned that I have entered the formula incorrectly than the actual integral itself. Is this the correct way to approach the problem.

サインインしてコメントする。

採用された回答

Torsten
Torsten 2019 年 2 月 15 日
format long
fun = @(t,x,y) exp(-t.^4 + 1i.*y.*t - x.*t.^2 + 1i*pi*0.125);
P = @(x,y) integral(@(t)fun(t,x*exp(-1i*pi*0.25),y*exp(1i*pi*0.125)),-Inf,Inf);
P(1,1)
Reference:
https://arxiv.org/pdf/1601.03615.pdf
  1 件のコメント
Michael Devereux
Michael Devereux 2019 年 2 月 15 日
Many thanks!

サインインしてコメントする。

その他の回答 (1 件)

Abhishek Hullur
Abhishek Hullur 2021 年 8 月 8 日
. Evaluate around the rectangle with vertices

カテゴリ

Help Center および File ExchangeMathematics についてさらに検索

タグ

製品


リリース

R2018b

Community Treasure Hunt

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

Start Hunting!

Translated by