Derivative from ode45

Question

jep 2020-4-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/518582-derivative-from-ode45

回答： madhan ravi 2020-4-17

I have a code that solves a second order differential equation denoted in a function called dfq

clear all
close all
y0 = [0;0];
tspan = [0 15];
[t,y] = ode45(@dfq,tspan,y0);
xresp = y(:,1);
plot(t,xresp)
title('Solution Plot')
xlabel('t(s)')
ylabel('x(m)')
grid on

where the response x(t) is delivered and plotted through xresp. I need to plot the derivative of xresp, or x(t), and attempted to use diff to do this but it shortened the size of the vector and could not be plotted. I am wondering if there is a better way to graph the derivative. I also tried to use gradient on xresp, which did not give the correct derivative plot either.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Peter O 2020-4-17

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/518582-derivative-from-ode45#answer_426538

编辑：Peter O 2020-4-17

在 MATLAB Online 中打开

Generally I send the result of ode45 back through my function as a post-processing step.

[t,x] = ode45(@myDerivFcn,tspan,x0)

You just need to be a little cautious of the problem dimensionality. arrayfun can be a neat trick if you've got a one-dimensional system or don't mind manipulating the cell array after the fact. Otherwise, a venerable for loop works:

dx = nan(n_dims_x,length(t))
for ix=1:length(t)
    dx(:,ix) = myDerivFcn(t(ix),x(ix,:).') % Col major fill
end
dx = dx.' % Now it matches the dimensionality of x and has length(t) rows