integro-differential equation?

26 visualizzazioni (ultimi 30 giorni)
PWini
PWini il 26 Giu 2016
Commentato: bloodtalon il 16 Giu 2017
I have a problem which I will try to describe in details. Please try to help me, because the exam is coming :) The task is to find the solution of the differential equation as follows:
A*d2v/dt2+B*dv/dt+C*v-P(t)=0
P(t)=integral(fun(t,z)dz)
So I have the integral function implemented into the diff equation, where the integration is not by t-time but some different parameter in this case z-height.
I know, I am not so experienced programmer and MATLAB user but please give me some advice how to write the script. If you need some more info about the nature of that mysterious functions ask me.

Risposte (2)

Roger Stafford
Roger Stafford il 27 Giu 2016
You should solve this using one of the ‘ode’ numerical differential equation functions for a second order equation in the ordinary way, with the exception that the function P(t) must be computed using matlab’s ‘integral’ function, which in turn uses the ‘fun’ function as its integrand. That means it will probably take a considerably longer time for execution because of the repeated full integration times. (I assume you know the limits of integration in P.)
You will of course have to furnish the ‘ode’ function the initial values of v and dv/dt at the initial value of t.
  2 Commenti
PWini
PWini il 27 Giu 2016
Ok, So I calculated P(t) in time domain and the result is a vector of solutions for different time periods. Now how to implement that force into the differential equation formulation? Because for every solution of that diff equation I have a set of results in time domain, but the real problem is I don't know how to connect the each result of integration to that diff equation problem. In my opinion I should obtain the integration with time parameters (in some kind of symbolic function) then solve the diff equetion and substitute the time parameter.
Torsten
Torsten il 27 Giu 2016
If you have a vector for P(t) at different time instances in advance (it seems P(t) does not depend on v in your case), take a look at the example "ODE with Time-Dependent Terms" under
to see how to deal with it.
Best wishes
Torsten.

Accedi per commentare.


Claudio Gelmi
Claudio Gelmi il 30 Gen 2017
Take a look at this MATLAB solver:
"IDSOLVER: A general purpose solver for nth-order integro-differential equations": http://dx.doi.org/10.1016/j.cpc.2013.09.008
MATLAB solver here for download (free): http://cpc.cs.qub.ac.uk/summaries/AEQU_v1_0.html
Best wishes,
Claudio
  1 Commento
bloodtalon
bloodtalon il 16 Giu 2017
Just a link like this isn't helping. Could you actually give some instructions on how to use it and/or some examples, preferably of the one in the question?

Accedi per commentare.

Categorie

Scopri di più su Programming in Help Center e File Exchange

Prodotti

Community Treasure Hunt

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

Start Hunting!

Translated by