ode15s

12 次查看（过去 30 天）
ahmad
ahmad 2012-1-30
Solving a partial differential equation with ode15s,we know that this ode solver integrates the ode over time direvatives.I need to find the values of these time direvatives,do you have any special command? for example: for i=1:n dudt=uxx(i) end
how can I find these dudt at each time and x(i)? thanks alot

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

采纳的回答

Andrew Newell
Andrew Newell 2012-1-30
If you are solving something like
[t,x] = solver(odefun,tspan,y0)
and you want the time derivative at each point t(i), then that is simply
dxdt = 0*x; % initialize the array
for ii=1:length(t)
dxdt(ii,:) = odefun(t(ii),x(ii,:));
end
(edited to allow for the possibility that x has more than one component at each time t).
  2 个评论
Bård Skaflestad
Bård Skaflestad 2012-1-30
That's more or less what I had in mind. Calling the |odefun| from an |OutputFcn| is just another refinement.
ahmad
ahmad 2012-2-20
So if I want to calculate the integral of square of dudt i don't have any special command in my mind,in fact I used trapz but it keeps saying "the length of x should be the same as the first non-singleton variable" or it says that "x should be a vector".

请先登录，再进行评论。

更多回答（2 个）

Bård Skaflestad
Bård Skaflestad 2012-1-30
If you're going to solve a PDE by the method of lines (i.e., converting it into a system of ODEs), then you will need to provide a spatial discretisation of the PDE along with suitable initial conditions whence ode15s (or others) will likely assist you in the numerical solution.
Please review your favourite text book on discretisation methods.
  4 个评论
显示 2更早的评论
Bård Skaflestad
Bård Skaflestad 2012-1-30
Sorry, I meant |odeset|, not |optimset|.
ahmad
ahmad 2012-1-30
I think that you are good at matlab numerical processes,can you please take a look at my program?I need to call the du/dt values solved by ode solver at every point like du/dt(t,x),is there a special command for it?it is true that its a big amount of data.so how can i do the other calculations using these values?
thanks

请先登录，再进行评论。

Andrew Newell
Andrew Newell 2012-1-30
You could use deval.
EDIT: Here is a modified version of the example at the link provided:
sol = ode45(@vdp1,[0 20],[2 0]);
t = linspace(0,20,100);
[y,dy] = deval(sol,t,1);
plot(t,dy);
  3 个评论
显示 1更早的评论
Andrew Newell
Andrew Newell 2012-1-30
I think it is the time derivative (as in the example I added), but for a polynomial fitting the solution. I provide a more direct solution in a separate answer.
Bård Skaflestad
Bård Skaflestad 2012-1-30
Indeed it is. The second output of |deval| is the derivative (wrt the independent variable) of the continuous output polynomial fitted to the numerical scheme of the particular solver (ODE45 in the above example).
Yet another detail I'd missed while skimming the documentation of |deval|.

请先登录，再进行评论。

类别

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

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by