MHB Numerics on wild oscillating functions

AI Thread Summary
Accurate computation of functions with intense oscillations, such as the integral of sin(x^2), poses challenges due to instability when integrating to infinity. A suggested approach is to split the integral into two parts to manage the infinite limit effectively. One method involves calculating the integral over a finite range and then addressing the improper part using integration by parts to enhance convergence. Additionally, collecting positive and negative contributions separately before summing them may improve accuracy. These strategies aim to achieve better numerical results for oscillating functions.
Theia
Messages
121
Reaction score
1
Hello!

I'd like to ask for a help about how to compute accurately functions which has very intense oscillations. My example is to estimate

$$I = \int_0^{\infty} \sin(x^2) dx= \int_0^{\infty}\frac{\sin(t)}{2\sqrt{t}} dt$$.I tried trapezoid rule over one oscillation at a time, but result is poor. My next though is to collect positive parts and negative parts together and add them later on in the code...

Any comments?

Thank you!
 
Mathematics news on Phys.org
Theia said:
Hello!

I'd like to ask for a help about how to compute accurately functions which has very intense oscillations. My example is to estimate

$$I = \int_0^{\infty} \sin(x^2) dx= \int_0^{\infty}\frac{\sin(t)}{2\sqrt{t}} dt$$.I tried trapezoid rule over one oscillation at a time, but result is poor. My next though is to collect positive parts and negative parts together and add them later on in the code...

Any comments?

Thank you!

I think integrating to infinity is typically unstable, isn't it?
So we should probably split the integral in two parts to eliminate infinity.
That is:
$$I = \int_0^{\infty} \sin(x^2)\,dx = \int_0^{\sqrt{2\pi}} \sin(x^2)\, dx + \int_0^{1/\sqrt{2\pi}} u^{-2}\sin(u^{-2})\, du$$
Perhaps we can integrate the second integral period by period now?
 
Thank you! I computed first

$$\int_0^{2\pi}\frac{\sin t}{2\sqrt{t}}$$,

and then the improper part by using integration by parts to increase $$\sqrt{t}$$ to $$\sqrt{t^3}$$ to obtain faster convergence.
 
Insights auto threads is broken atm, so I'm manually creating these for new Insight articles. In Dirac’s Principles of Quantum Mechanics published in 1930 he introduced a “convenient notation” he referred to as a “delta function” which he treated as a continuum analog to the discrete Kronecker delta. The Kronecker delta is simply the indexed components of the identity operator in matrix algebra Source: https://www.physicsforums.com/insights/what-exactly-is-diracs-delta-function/ by...
Suppose ,instead of the usual x,y coordinate system with an I basis vector along the x -axis and a corresponding j basis vector along the y-axis we instead have a different pair of basis vectors ,call them e and f along their respective axes. I have seen that this is an important subject in maths My question is what physical applications does such a model apply to? I am asking here because I have devoted quite a lot of time in the past to understanding convectors and the dual...
For original Zeta function, ζ(s)=1+1/2^s+1/3^s+1/4^s+... =1+e^(-sln2)+e^(-sln3)+e^(-sln4)+... , Re(s)>1 If we regards it as some function got from Laplace transformation, and let this real function be ζ(x), that means L[ζ(x)]=ζ(s), then: ζ(x)=L^-1[ζ(s)]=δ(x)+δ(x-ln2)+δ(x-ln3)+δ(x-ln4)+... , this represents a series of Dirac delta functions at the points of x=0, ln2, ln3, ln4, ... , It may be still difficult to understand what ζ(x) means, but once it is integrated, the truth is clear...
Back
Top