Solving an ODE when the function is given as discrete values -matlab-

by

You can use interp1 to create an interpolated lookup table function:

fx = [0.5644 0.6473 0.7258 0.7999 0.8697 0.9353 0.9967 1.0540 1.1072 1.1564 ...
      1.2016 1.2429 1.2803 1.3138 1.3435 1.3695 1.3917 1.4102 1.4250 1.4362 ...
      1.4438 1.4477 1.4482 1.4450 1.4384 1.4283 1.4147 1.3977 1.3773 1.3535 ...
      1.3263 1.2957 1.2618 1.2246 1.1841 1.1403 1.0932 1.0429 0.9893 0.9325 0.8725];
x = 0:0.25:10
f = @(xq)interp1(x,fx,xq);

Then you should be able to use ode45 as normal:

tspan = [0 1];
x0 = 2;
xout = ode45(@(t,x)f(x),tspan,x0);

Note that you did not specify what values of of x your function (fx here) is evaluated over so I chose zero to ten. You’ll also not want to use the copy-and-pasted values from the command window of course because they only have four decimal places of accuracy. Also, note that because ode45 required the inputs t and then x, I created a separate anonymous function using f, but f can created with an unused t input if desired.

More Related Contents:

  1. Matlab ode solvers: changing state and specified time
  2. MATLAB: Is it possible to have two event values whilst using ode45?
  3. Absolute error of ODE45 and Runge-Kutta methods compared with analytical solution
  4. Plotting logarithmic values in separate graphs [closed]
  5. How to draw a Histogram in Matlab
  6. For loop to split matrix to equal sized sub-matrices
  7. How to normalize a histogram in MATLAB?
  8. Difference between accessing cell elements using curly braces and parentheses
  9. Understanding Matlab FFT example
  10. Nearest-neighbor interpolation algorithm in MATLAB
  11. How to assign values to a MATLAB matrix on the diagonal?
  12. Find location of current m-file in MATLAB
  13. MATLAB – extracting rows of a matrix
  14. Is there a foreach in MATLAB? If so, how does it behave if the underlying data changes?
  15. MATLAB is running out of memory but it should not be
  16. Suppressing Output MATLAB
  17. Matrix “Zigzag” Reordering
  18. Fast Algorithms for Finding Pairwise Euclidean Distance (Distance Matrix)
  19. How do I get real integer overflows in MATLAB/Octave?
  20. How do you handle resources in MATLAB in an exception safe manner? (like “try … finally”)
  21. How can I close files that are left open after an error?
  22. quiver3 arrow color corresponding to magnitude
  23. MatLab accumarray unexpectedly changing ordering
  24. How can I plot a 3D-plane in Matlab?
  25. Plotting Implicit Algebraic equations in MATLAB
  26. imresize – trying to understand the bicubic interpolation
  27. How to use Hash Tables (dictionaries) in MATLAB? [duplicate]
  28. How to draw a line on an image in matlab?
  29. Matlab: arrayfun, cellfun, spfun and structfun vs. simple for-loop
  30. Automatically maximize a figure

Leave a Comment