how to interpolate the data over uniform time grid

3 次查看（过去 30 天）
PetronasAMG
PetronasAMG 2021-9-26
回答： Nipun 2024-5-30
Since numerical solution obtained from ode45 is a non uniform in time, I need to obtain the frequency power spectrim of the solution, but first, i need to interpolate the date over unifrom time grid. How would i able to do so? I am completely stuck after inputting vlause from ODE45! Please explain the steps!
function dxdt = frs(t,x)
dxdt = zeros(3,1);
dx1dt = -((4/5)*x(2)+x(1));
dx2dt = x(3) + ((4/9)*x(2));
dx3dt = (1/8)+((x(1)-2.5)*x(2));
end
tspan = [0:10:100];
xi = [1.0;4.0;2.0];
[t,x]= ode45(@frs,tspan,xi);
from here I do not know how to process with the question please help!
  2 个评论
Walter Roberson
Walter Roberson 2021-9-26
Since numerical solution obtained from ode45 is a non uniform in time
That is only the case if you pass in a two-element tspan. When you pass in a tspan that is uniform in time, the way you do, then the output from ode45 will only be at those times, and so the output would be uniform in time.
Mathieu NOE
Mathieu NOE 2021-9-27
you can also interpolate the output data of your ode solution over a new time vector (linear)
use interp1
maybe this will be faster than forcing ode to work on fixed , refined time vector

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Nipun
Nipun 2024-5-30
Hi [First Name]
I understand that you want to interpolate data over a uniform time grid. Here’s how you can do it using interp1 in MATLAB:
% Given data
time = [0, 1, 4, 5, 6]; % Non-uniform time points
data = [10, 20, 15, 25, 30]; % Corresponding data points
% Uniform time grid
uniformTime = linspace(min(time), max(time), 100);
% Interpolate data over uniform time grid
uniformData = interp1(time, data, uniformTime, 'linear');
% Display results
plot(time, data, 'o', uniformTime, uniformData, '-');
xlabel('Time');
ylabel('Data');
legend('Original Data', 'Interpolated Data');
For more details on interpolation, refer to the MathWorks documentation on interp1: https://www.mathworks.com/help/matlab/ref/interp1.html
Hope this helps.
Regards,
Nipun

类别

在 Help Center 和 File Exchange 中查找有关 MATLAB 的更多信息

标签

产品


版本

R2021b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by