How can I overcome errors in RK4 due to linearly interpolated data?

Liferider
Messages
41
Reaction score
0
I am simulating an ODE where the differential function is a function of sampled data points and are having trouble obtaining data in between sampled points (for the step computations of RK4). The first thing that I tried was to linearly interpolate the data, but that introduced large accumulating errors in the simulation. Is there a way to circumvent this?

One possibility of course, is to use Euler integration with the appropriate step length, but it would be nice to use a higher order scheme if possible.
 
Physics news on Phys.org
Don't bother, just very sensitive regarding step length.
 
There is the following linear Volterra equation of the second kind $$ y(x)+\int_{0}^{x} K(x-s) y(s)\,{\rm d}s = 1 $$ with kernel $$ K(x-s) = 1 - 4 \sum_{n=1}^{\infty} \dfrac{1}{\lambda_n^2} e^{-\beta \lambda_n^2 (x-s)} $$ where $y(0)=1$, $\beta>0$ and $\lambda_n$ is the $n$-th positive root of the equation $J_0(x)=0$ (here $n$ is a natural number that numbers these positive roots in the order of increasing their values), $J_0(x)$ is the Bessel function of the first kind of zero order. I...
Are there any good visualization tutorials, written or video, that show graphically how separation of variables works? I particularly have the time-independent Schrodinger Equation in mind. There are hundreds of demonstrations out there which essentially distill to copies of one another. However I am trying to visualize in my mind how this process looks graphically - for example plotting t on one axis and x on the other for f(x,t). I have seen other good visual representations of...
Back
Top