How to use "deval" or any other function non ode solver output

Question

Romio 2023-8-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2004222-how-to-use-deval-or-any-other-function-non-ode-solver-output

编辑： Torsten 2023-8-3

Suppose I have a solution to an ode that I obtained numerically using an accurate method that is close to the exact solution. Call it Y_exact which may be of size 2 x 1000, where 2 is for the two solution components, and 1000 is the number of time points. Suppose I have another solution by a differnt method Y of size 2 x n where n <= 1000. Let tout = 1:dt:n

Now if Y_exact obtained with ode45 or any other solver, I could obtain exact solution at the same time points by doing deval(sol,tout) where sol is the solution structure obtained from ode45.

Now since Y_exact in my case is not obtained with a matlab solver, and is just a vector, is there any way to use deval() or any other way to match Y and Y_exact at the same time points? So that I can compute the error norm(Y - Ytrue,'inf').