Мені потрібно чисельно оцінити інтеграл нижче:
де , і . Тут - модифікована функція Бесселя другого роду. У моєму конкретному випадку я , і .x∈R+λ,κ,ν>0Kλ=0,00313κ=0,00825ν=0,33
Я використовую MATLAB, і я спробував вбудовані функції, integral
і quadgk
це дає мені багато помилок (див. Нижче). Я, звичайно, також пробував багато інших речей, таких як інтеграція по частинах та підсумовування інтегралів від до .( k + 1 ) x π
Отже, чи є у вас якісь пропозиції щодо того, який метод я слід спробувати далі?
ОНОВЛЕННЯ (додаткові запитання)
Я читав статтю @Pedro, пов’язану з цим, і не думаю, що це було занадто важко зрозуміти. Однак у мене є кілька питань:
- Чи було б добре використовувати як базисні елементи , в описаному уніваріантному методі Левіна?ψ k
- Чи можу я замість цього просто використовувати метод Філона, оскільки частота коливань фіксована?
Приклад коду
>> integral(@(r) sin(x*r).*sqrt(E(r)),0,Inf)
Warning: Reached the limit on the maximum number of intervals in use. Approximate
bound on error is 1.6e+07. The integral may not exist, or it may be difficult to
approximate numerically to the requested accuracy.
> In funfun\private\integralCalc>iterateScalarValued at 372
In funfun\private\integralCalc>vadapt at 133
In funfun\private\integralCalc at 84
In integral at 89
ans =
3.3197e+06